A mechanistic model on the role of "radially-running" collagen fibers on dissection properties of human ascending thoracic aorta

Effects of Age, Elastin Density, and Glycosaminoglycan Accumulation on the Delamination Strength of Human Thoracic and Abdominal Aortas

Aortic dissection is a life-threatening condition caused by layer separation. Despite extensive research, the relationship between the aortic wall's structural integrity and dissection risk remains unclear. Glycosaminoglycan (GAG) accumulation and elastin loss are suspected to play significant roles. We investigated how age-related changes in aortic structure affect dissection susceptibility. Peeling tests were performed on longitudinal and circumferential thoracic (TA) and abdominal aortic (AA) strips from 35 donors aged 13-76 years (mean 38±15 years, 34% female). GAG, elastin, collagen, and smooth muscle cell (SMC) contents were assessed using bidirectional histology. Young TAs resisted longitudinal peeling better than circumferential, with delamination strengths of 65.4 mN/mm and 44.2 mN/mm, respectively. Delamination strength decreased with age in both directions, more rapidly longitudinally, equalizing at ∼20-25 mN/mm in older TAs. Delamination strength in AAs was 22% higher than in TAs. No sex differences were observed. GAG density increased, while elastin density decreased by 2.5% and 4% per decade, respectively. Collagen density did not change with age, while SMC density decreased circumferentially. GAGs partially mediated the reduction in longitudinal delamination strength due to aging, while circumferential strength reduction was not mediated by changes in either GAG or elastin densities. This study explains why aortic dissections are more common in TAs, especially in older individuals, and why they typically propagate spirally. TAs exhibit lower delamination strength compared to AAs and experience strength reduction with age, a phenomenon linked to increased GAG accumulation and elastin loss. These findings enhance our understanding of the pathophysiological mechanisms behind aortic dissection. STATEMENT OF SIGNIFICANCE: This work explores the age-dependent relationships between delamination strength in human aortas and wall structural content. We investigated 35 human aortas from donors aged 13 to 76 years, providing new insights into the biomechanical and histological factors that influence aortic dissection risk. Our findings elucidate how variations in elastin, glycosaminoglycan, collagen, and smooth muscle cell densities impact the structural integrity of the aorta, contributing significantly to the understanding of aortic dissection mechanisms.

Traction-separation law parameters for the description of age-related changes in the delamination strength of the human descending thoracic aorta

Aortic dissection is a life-threatening disease that consists in the development of a tear in the wall of the aorta. The initial tear propagates as a discontinuity leading to separation within the aortic wall, which can result in the creation of a so-called false lumen. A fatal threat occurs if the rupture extends through the whole thickness of the aortic wall, as blood may then leak. It is generally accepted that the dissection, which can sometime extend along the entire length of the aorta, propagates via a delamination mechanism. The aim of the present paper is to provide experimentally validated parameters of a mathematical model for the description of the wall's cohesion. A model of the peeling experiment was built in Abaqus. The delamination interface was described by a piecewise linear traction-separation law. The bulk behavior of the aorta was assumed to be nonlinearly elastic, anisotropic, and incompressible. Our simulations resulted in estimates of the material parameters for the traction-separation law of the human descending thoracic aorta, which were obtained by minimizing the differences between the FEM predictions and the delamination force given by the regression of the peeling experiments. The results show that the stress at damage initiation, Tc, should be understood as an age-dependent quantity, and under the assumptions of our model this dependence can be expressed by linear regression as Tc =  - 13.03·10-4·Age + 0.2485 if the crack front advances in the axial direction, and Tc =  - 7.58·10-4·Age + 0.1897 if the crack front advances in the direction of the aortic circumference (Tc [MPa], Age [years]). Other model parameters were the stiffness K and the separation at failure, δf-δc (K = 0.5 MPa/mm, δf-δc = 0.1 mm). The material parameters provided by our study can be used in numerical simulations of the biomechanics of dissection propagation through the aorta especially when age-associated phenomena are studied.

Atherosclerotic Calcifications Have a Local Effect on the Peel Behavior of Human Aortic Media

Aortic dissections, characterized by the propagation of a tear through the layers of the vessel wall, are critical, life-threatening events. Aortic calcifications are a common comorbidity in both acute and chronic dissections, yet their impact on dissection mechanics remains unclear. Using micro-computed tomography (CT) imaging, peel testing, and finite element modeling, this study examines the interplay between atherosclerotic calcifications and dissection mechanics. Samples cut from cadaveric human thoracic aortas were micro-CT imaged and subsequently peel-tested to map peel tension curves to the location of aortic calcifications. Empirical mode decomposition separated peel tension curves into high and low-frequency components, with high-frequency effects corresponding to interlamellar bonding mechanics and low-frequency effects to peel tension fluctuations. Finally, we used an idealized finite element model to examine how stiff calcifications affect aortic failure mechanics. Results showed that atherosclerosis influences dissection behavior on multiple length scales. Experimentally, atherosclerotic samples exhibited higher peel tensions and greater variance in the axial direction. The variation was driven by increased amplitudes of low-frequency tension fluctuations in diseased samples, indicating that more catastrophic propagations occur near calcifications. The simulations corroborated this finding, suggesting that the low-frequency changes resulted from the presence of a stiff calcification in the vessel wall. There were also modifications to the high-frequency peel mechanics, a response likely attributable to alterations in the microstructure and interlamellar bonding within the media. Considered collectively, these findings demonstrate that dissection mechanics are modified in aortic media nearby and adjacent to aortic calcifications.

Correlation between age, location, orientation, loading velocity and delamination strength in the human aorta

Aortic dissection is a biomechanical phenomenon associated with a failure of internal cohesion, which manifests itself through the delamination of the aortic wall. The goal of this study is to deepen our knowledge of the delamination strength of the aorta. To achieve this, 661 peeling experiments were carried out with strips of the human aorta collected from 46 cadavers. The samples were ordered into groups with respect to (1) anatomical location, (2) orientation of the sample, and (3) extension rate used within the experiment. The obtained results are in accordance with the hypothesis that delamination resistance is not sensitive to the extension rates 0.1, 1, 10, and 50 mms^-1. We arrived at this conclusion for all positions along the aorta investigated in our study. These were the thoracic ascending (AAs), thoracic descending (ADs), and the abdominal aorta (AAb), simultaneously considering both the longitudinal (L) as well as the circumferential (C) orientations of the samples. On the other hand, our results showed that the delamination strength differs significantly with respect to the anatomical position and orientation of the sample. The medians of the delamination strength were as follows, 4.1 in AAs-L, 3.2 in AAs-C, 3.1 in ADs-L, 2.4 in ADs-C, AAb-L in 3.6, and 2.7 in AAb-C case (all values are in 0.01·Nmm^-1). This suggests that resistance to crack propagation should be an anisotropic property and that the aorta is inhomogeneous along its length from the point of view of delamination resistance. Finally, correlation analysis proved that the delamination strength of the human aorta significantly decreases with age.

Regional delamination strength in the human aorta underlies the anatomical localization of the dissection channel

Aortic dissection is a life-threatening event, during which a primary tear propagates along the aorta causing catastrophic delamination of the inner (intima with most of the media) from the outer layers (leftover media with adventitia). Our understanding of mode-I fracture resistance at different aortic regions is incomplete, although the anatomical localization of the dissection channel may be assigned to this factor. To determine whether the susceptibility to dissection propagation varied with aortic region, the average and standard deviation of peel tension (indices of adhesive strength between layers when pulled apart and its fluctuation) were measured in 24 cadaveric subjects. Measurements were made in the inner and outer quadrants of 9 consecutive regions. Strong regional heterogeneity was established that was age-related based on the following evidence: (1) the average and standard deviation of peel tension peaked in the ascending aorta, decreasing to almost constant values in the descending thoracic aorta, but increasing across the abdominal aorta; (2) axial differences were more pronounced in the inner quadrant, with differences among quadrants reaching significance proximally; (3) the average peel tension was greatly impaired from <40 to 40-60 but much less to >60-year-old subjects at most regions/quadrants, leading to non-uniform axial variations in all age groups; (4) gender affected little the data. This comprehensive series of delamination tests explains the clinical observation of most dissections initiating in the ascending aorta to extend distally and of few dissections initiating in the descending thoracic aorta to extend proximally, while supporting the increased vulnerability in aged subjects.

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.

Association between single nucleotide polymorphisms of tropoelastin gene and aortic dissection

OBJECTIVES

To evaluate the relation between single nucleotide polymorphisms (SNPs) of tropoelastin gene and aortic dissection (AD) via identifying SNPs in the tropoelastin gene, and to detect the level of tropoelastin mRNA, elastin and elastic fibers.

METHODS

The specimens of the AD group (n=96) and the control group (n=95), including their blood and aortic wall tissues, were collected. DNA was extracted from the blood samples in the 2 groups, and the SNPs in the tropoelastin gene were examined by the MassARRAY genotyping technique, and their haplotypes were constructed by PHASE software. The expression of tropoelastin mRNA and elastin in the aortic tunica media was respectively detected by real-time PCR or Western blotting. Elastin Van Gieson (EVG) staining was used to observe the shape of aortic tunica media and clarify the distribution of elastic fibers. The frequency of genotypes and haplotypes of SNP loci in the tropoelastin gene was analyzed and compared between the 2 groups, and the expression of tropoelastin mRNA, elastin and elastic fibers were also compared.

RESULTS

Seven SNP loci of the tropoelastin gene were detected in these samples. Among them, 5 SNP loci were polymorphic. The frequency of 3 SNP loci[rs2071307 (G/A), rs34945509 (C/T) and rs17855988 (G/C)] was significantly different between the AD group and the control group (all P<0.05). There were significantly different in the haplotypes frequency of rs2071307 (G/A), rs34945509 (C/T) and rs17855988 (G/C) between the 2 groups (all P<0.01). Real-time PCR and Western blotting showed that the relative expression of tropoelastin mRNA and elastin in the aortic tunica media in the AD group was significantly lower than that in the control group (P<0.05). EVG staining showed that the aortic tunica media was torn, the morphology and structure of elastic fibers were broken, cracked, and disordered in the AD group, while the aortic tunica media was in complete structure and well arrangement.The elastic fibers were presented closely and orderly in the control group.

CONCLUSIONS

The polymorphisms of rs2071307 (G/A), rs34945509 (C/T), and rs17855988(G/C) in the tropoelastin gene may eventually affect the synthesis of elastic fibers and they may play an important role in the occurrence of AD.

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.

Review of Current Advances in the Mechanical Description and Quantification of Aortic Dissection Mechanisms

Aortic dissection is a life-threatening event associated with a very poor outcome. A number of complex phenomena are involved in the initiation and propagation of the disease. Advances in the comprehension of the mechanisms leading to dissection have been made these last decades, thanks to improvements in imaging and experimental techniques. However, the micro-mechanics involved in triggering such rupture events remains poorly described and understood. It constitutes the primary focus of the present review. Towards the goal of detailing the dissection phenomenon, different experimental and modeling methods were used to investigate aortic dissection, and to understand the underlying phenomena involved. In the last ten years, research has tended to focus on the influence of microstructure on initiation and propagation of the dissection, leading to a number of multiscale models being developed. This review brings together all these materials in an attempt to identify main advances and remaining questions.

Mechanical and structural contributions of elastin and collagen fibers to interlamellar bonding in the arterial wall

The artery relies on interlamellar structural components, mainly elastin and collagen fibers, for maintaining its integrity and resisting dissection propagation. In this study, the contribution of arterial elastin and collagen fibers to interlamellar bonding was studied through mechanical testing, multiphoton imaging and finite element modeling. Steady-state peeling experiments were performed on porcine aortic media and the purified elastin network in the circumferential (Circ) and longitudinal (Long) directions. The peeling force and energy release rate associated with mode-I failure are much higher for aortic media than for the elastin network. Also, longitudinal peeling exhibits a higher energy release rate and strength than circumferential peeling for both the aortic media and elastin. Multiphoton imaging shows the recruitment of both elastin and collagen fibers within the interlamellar space and points to in-plane anisotropy of fiber distributions as a potential mechanism for the direction-dependent phenomena of peeling tests. Three-dimensional finite element models based on cohesive zone model (CZM) of fracture were created to simulate the peeling tests with the interlamellar energy release rate and separation distance at damage initiation obtained directly from peeling test. Our experimental results show that the separation distance at damage initiation is 80 μm for aortic media and 40 μm for elastin. The damage initiation stress was estimated from the model for aortic media (Circ: 60 kPa; Long: 95 kPa) and elastin (Circ: 9 kPa; Long: 14 kPa). The interlamellar separation distance at complete failure was estimated to be 3 - 4 mm for both media and elastin. Furthermore, elastin and collagen fibers both play an important role in bonding of the arterial wall, while collagen has a higher contribution than elastin to interlamellar stiffness, strength and toughness. These results on microstructural interlamellar failure shed light on the pathological development and progression of aortic dissection.

Effect of localized tendon remodeling on supraspinatus tear propagation

Rotator cuff tear propagation is multifactorial and may be due to localized changes in mechanical properties from tendon remodeling based on the inhomogeneous stresses experienced by a tendon with a tear. The objective of this study was to investigate the effect of localized tendon remodeling on tear propagation for simulated supraspinatus tendon tears. A validated computational model of a supraspinatus tendon using subject-specific geometry and material properties with a 1 cm wide anterior tear was used. The medial edge of the supraspinatus tendon was displaced 5 mm to induce tear propagation and cohesive elements were used to model tear propagation. Four remodeling scenarios were investigated: (1) Baseline (no remodeling), (2) Positive remodeling (increased fiber stiffness) and (3) Negative remodeling (decreased fiber stiffness) at tear tips, and (4) Negative remodeling along the medial-lateral tear edge. Output parameters included the amount of tear propagation, critical load to propagate the tear, and maximum principal stress at the tear tips. Positive remodeling at the tear tips resulted in the largest amount of tear propagation (18.4 mm), highest peak maximum principal stress (25.2 MPa), and lowest critical load to propagate the tear (249N). Conversely, negative remodeling at the tear tips resulted in the least amount of tear propagation (16 mm), lowest peak maximum principal stress (17.6 MPa) and highest critical load to propagate the tear (278N). Overall, remodeling at the tear tips has the greatest effect on tear propagation. Therefore, a better method for clinicians to measure tendon stiffness at the tear tips would be helpful to improve outcome of patients.

Avalanches and power law behavior in aortic dissection propagation

Aortic dissection is a devastating cardiovascular disease known for its rapid propagation and high morbidity and mortality. The mechanisms underlying the propagation of aortic dissection are not well understood. Our study reports the discovery of avalanche-like failure of the aorta during dissection propagation that results from the local buildup of strain energy followed by a cascade failure of inhomogeneously distributed interlamellar collagen fibers. An innovative computational model was developed that successfully describes the failure mechanics of dissection propagation. Our study provides the first quantitative agreement between experiment and model prediction of the dissection propagation within the complex extracellular matrix (ECM). Our results may lead to the possibility of predicting such catastrophic events based on microscopic features of the ECM.

Segmental Variations in the Peel Characteristics of the Porcine Thoracic Aorta

Aortic dissection occurs predominantly in the thoracic aorta and the mechanisms for the initiation and propagation of the tear in aortic dissection are not well understood. We study the tearing characteristics of the porcine thoracic aorta using a peeling test and we estimate the peeling energy per unit area in the ascending and the descending segments. The stretch and the peel force per unit width undergone by the peeled halves of a rectangular specimen are measured. We find that there can be significant variation in the stretch within the specimen and the stretch between the markers in the specimen varies with the dynamics of peeling. We found that in our experiment the stretch achieved in the peeled halves was such that it was in the range of the stretch at which the stress-stretch curve for the uniaxial experiment starts deviating from linearity. Higher peeling energy per unit area is required in the ascending aorta compared to the descending aorta. Longitudinal specimens required higher peeling energy per unit area when compared to the circumferential specimens.

Ex Vivo Mechanical Tests and Multiscale Computational Modeling Highlight the Importance of Intramural Shear Stress in Ascending Thoracic Aortic Aneurysms

Ascending thoracic aortic aneurysms (ATAAs) are anatomically complex in terms of architecture and geometry, and both complexities contribute to unpredictability of ATAA dissection and rupture in vivo. The goal of this work was to examine the mechanism of ATAA failure using a combination of detailed mechanical tests on human tissue and a multiscale computational model. We used (1) multiple, geometrically diverse, mechanical tests to characterize tissue properties; (2) a multiscale computational model to translate those results into a broadly usable form; and (3) a model-based computer simulation of the response of an ATAA to the stresses generated by the blood pressure. Mechanical tests were performed in uniaxial extension, biaxial extension, shear lap, and peel geometries. ATAA tissue was strongest in circumferential extension and weakest in shear, presumably because of the collagen and elastin in the arterial lamellae. A multiscale, fiber-based model using different fiber properties for collagen, elastin, and interlamellar connections was specified to match all of the experimental data with one parameter set. Finally, this model was used to simulate ATAA inflation using a realistic geometry. The predicted tissue failure occurred in regions of high stress, as expected; initial failure events involved almost entirely interlamellar connections, consistent with arterial dissection-the elastic lamellae remain intact, but the connections between them fail. The failure of the interlamellar connections, paired with the weakness of the tissue under shear loading, is suggestive that shear stress within the tissue may contribute to ATAA dissection.

Particle-based computational modelling of arterial disease

Accumulated glycosaminoglycans (GAGs) can sequester water and induce swelling within the intra-lamellar spaces of the medial layer of an artery. It is increasingly believed that stress concentrations caused by focal swelling can trigger the damage and delamination that is often seen in thoracic aortic disease. Here, we present computational simulations using an extended smoothed particle hydrodynamics approach to examine potential roles of pooled GAGs in initiating and propagating intra-lamellar delaminations. Using baseline models of the murine descending thoracic aorta, we first calculate stress distributions in a healthy vessel. Next, we examine increases in mechanical stress in regions surrounding GAG pools. The simulations show that smooth muscle activation can partially protect the wall from swelling-associated damage, consistent with experimental observations, but the wall can yet delaminate particularly in cases of smooth muscle dysfunction or absence. Moreover, pools of GAGs located at different but nearby locations can extend and coalesce, thus propagating a delamination. These findings, combined with a sensitivity study on the input parameters of the model, suggest that localized swelling can alter aortic mechanics in ways that eventually can cause catastrophic damage within the wall. There is, therefore, an increased need to consider roles of GAGs in aortic pathology.

Biomechanics of aortic wall failure with a focus on dissection and aneurysm: A review

Aortic dissections and aortic aneurysms are fatal events characterized by structural changes to the aortic wall. The maximum diameter criterion, typically used for aneurysm rupture risk estimations, has been challenged by more sophisticated biomechanically motivated models in the past. Although these models are very helpful for the clinicians in decision-making, they do not attempt to capture material failure. Following a short overview of the microstructure of the aorta, we analyze the failure mechanisms involved in the dissection and rupture by considering also traumatic rupture. We continue with a literature review of experimental studies relevant to quantify tissue strength. More specifically, we summarize more extensively uniaxial tensile, bulge inflation and peeling tests, and we also specify trouser, direct tension and in-plane shear tests. Finally we analyze biomechanically motivated models to predict rupture risk. Based on the findings of the reviewed studies and the rather large variations in tissue strength, we propose that an appropriate material failure criterion for aortic tissues should also reflect the microstructure in order to be effective. STATEMENT OF SIGNIFICANCE: Aortic dissections and aortic aneurysms are fatal events characterized by structural changes to the aortic wall. Despite the advances in medical, biomedical and biomechanical research, the mortality rates of aneurysms and dissections remain high. The present review article summarizes experimental studies that quantify the aortic wall strength and it discusses biomechanically motivated models to predict rupture risk. We identified contradictory observations and a large variation within and between data sets, which may be due to biological variations, different sample sizes, differences in experimental protocols, etc. Based on the findings of the reviewed literature and the rather large variations in tissue strength, it is proposed that an appropriate criterion for aortic failure should also reflect the microstructure.

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.

A combined experimental-numerical lamellar-scale approach of tensile rupture in arterial medial tissue using X-ray tomography

Aortic dissection represents a serious cardio-vascular disease and life-threatening event. Dissection is a sudden delamination event of the wall, possibly leading to rupture within a few hours. Current knowledge and practical criteria to understand and predict this phenomenon lack reliable models and experimental observations of rupture at the lamellar scale. In an attempt to quantify rupture-related parameters, the present study proposes an analytical model that reproduces a uniaxial test on medial arterial samples observed under X-ray tomography. This model is composed of several layers that represent the media of the aortic wall, each having proper elastic and damage properties. Finite element models were created to validate the analytical model using user-defined parameters. Once the model was validated, an inverse analysis was used to fit the model parameters to experimental curves of uniaxial tests from a published study. Because this analytical model did not consider delamination strength between layers, a finite element model that included this phenomenon was also developed to investigate the influence of the delamination on the stress-strain curve through a sensitivity analysis. It was shown that shear delamination strength between layers, i.e. mode II separation, is essential in the rupture process observed experimentally.

Isotropic Failure Criteria Are Not Appropriate for Anisotropic Fibrous Biological Tissues

The von Mises (VM) stress is a common stress measure for finite element models of tissue mechanics. The VM failure criterion, however, is inherently isotropic, and therefore may yield incorrect results for anisotropic tissues, and the relevance of the VM stress to anisotropic materials is not clear. We explored the application of a well-studied anisotropic failure criterion, the Tsai–Hill (TH) theory, to the mechanically anisotropic porcine aorta. Uniaxial dogbones were cut at different angles and stretched to failure. The tissue was anisotropic, with the circumferential failure stress nearly twice the axial (2.67 ± 0.67 MPa compared to 1.46 ± 0.59 MPa). The VM failure criterion did not capture the anisotropic tissue response, but the TH criterion fit the data well (R2 = 0.986). Shear lap samples were also tested to study the efficacy of each criterion in predicting tissue failure. Two-dimensional failure propagation simulations showed that the VM failure criterion did not capture the failure type, location, or propagation direction nearly as well as the TH criterion. Over the range of loading conditions and tissue geometries studied, we found that problematic results that arise when applying the VM failure criterion to an anisotropic tissue. In contrast, the TH failure criterion, though simplistic and clearly unable to capture all aspects of tissue failure, performed much better. Ultimately, isotropic failure criteria are not appropriate for anisotropic tissues, and the use of the VM stress as a metric of mechanical state should be reconsidered when dealing with anisotropic tissues.

Effects of Tendon Degeneration on Predictions of Supraspinatus Tear Propagation

Rotator cuff tendons undergo degeneration with age, which could have an impact on tear propagation. The objective of this study was to predict tear propagation for different levels of tissue degeneration using an experimentally validated finite element model of a supraspinatus tendon. It was hypothesized that greater amounts of degeneration will result in tear propagation at lower loads than tendons with less degeneration. Using a previously-validated computational model of supraspinatus tendon, 1-cm tears were introduced in the anterior, middle, and posterior thirds of the tendon. Cohesive elements were assigned subject-specific failure properties to model tear propagation, and tendon degeneration ranging from "minimal" to "severe" was modeled by modifying its mechanical properties. Tears in tendons with severe degeneration required the smallest loads to propagate (122-207 N). Posterior tears required greater loads compared to middle and anterior tears at all levels of degeneration. Stress and strain required for tear propagation decreased substantially with degeneration, ranging from 8.5 MPa and 32.6% strain for minimal degeneration and 0.6 MPa and 4.5% strain for severe degeneration. Overall, this work indicates that greater amounts of tendon degeneration lead to greater risk of tear propagation, supporting the need for early detection and treatment of rotator cuff tears.

Biomechanical relevance of the microstructure in artery walls with a focus on passive and active components

The microstructure of arteries, consisting, in particular, of collagen, elastin, and vascular smooth muscle cells, plays a very significant role in their biomechanical response during a cardiac cycle. In this article, we highlight the microstructure and the contributions of each of its components to the overall mechanical behavior. We also describe the changes of the microstructure that occur as a result of abdominal aortic aneurysms and disease, such as atherosclerosis. We also focus on how the passive and active constituents are incorporated into a mathematical model without going into detail of the mathematical formulation. We conclude by mentioning open problems toward a better characterization of the biomechanical aspects of arteries that will be beneficial for a better understanding of cardiovascular pathophysiology.

Effects of tear size and location on predictions of supraspinatus tear propagation

Rotator cuff tears remain a significant clinical problem with a high incidence rate and severe clinical burden. Previous computational models developed to study rotator cuff tears have not modeled tissue damage and tear propagation. The objective of this study was to predict tear propagation for various combinations of tear size and location using an experimentally validated finite element model of supraspinatus tendon. It was hypothesized that larger rotator cuff tears propagate at lower loads than smaller tears, and that posterior tears require higher loads to propagate than anterior tears. Using a previously validated computational model of supraspinatus tendon, tears of size 0.5-1.5 cm were introduced to the tendon geometry in the anterior, middle, and posterior tendon thirds. Cohesive elements were assigned subject-specific failure properties and used to model tissue damage and tear propagation. A displacement of 5 mm was applied to the medial tendon edge to induce tear propagation. Model outputs included critical load required to propagate the tear, and principal stress and maximum principal strain at the anterior and posterior tear tips. For all tear sizes, posterior tears required the highest loads to propagate (247-567 N). Anterior tears generally required the least load to propagate (171-280 N). Stress and strain were larger on the articular side (maximum 33.9% articular strain vs 27.8% bursal strain). Overall, larger tears located in the anterior supraspinatus tendon that interrupt the rotator cable are most at risk for tear propagation, and should be carefully followed by clinicians when considering treatment options.

Structural modeling reveals microstructure-strength relationship for human ascending thoracic aorta

High lethality of aortic dissection necessitates accurate predictive metrics for dissection risk assessment. The not infrequent incidence of dissection at aortic diameters <5.5 cm, the current threshold guideline for surgical intervention (Nishimura et al., 2014), indicates an unmet need for improved evidence-based risk stratification metrics. Meeting this need requires a fundamental understanding of the structural mechanisms responsible for dissection evolution within the vessel wall. We present a structural model of the repeating lamellar structure of the aortic media comprised of elastic lamellae and collagen fiber networks, the primary load-bearing components of the vessel wall. This model was used to assess the role of these structural features in determining in-plane tissue strength, which governs dissection initiation from an intimal tear. Ascending aortic tissue specimens from three clinically-relevant patient populations were considered: non-aneurysmal aorta from patients with morphologically normal tricuspid aortic valve (CTRL), aneurysmal aorta from patients with tricuspid aortic valve (TAV), and aneurysmal aorta from patients with bicuspid aortic valve (BAV). Multiphoton imaging derived collagen fiber organization for each patient cohort was explicitly incorporated in our model. Model parameters were calibrated using experimentally-measured uniaxial tensile strength data in the circumferential direction for each cohort, while the model was validated by contrasting simulated tissue strength against experimentally-measured strength in the longitudinal direction. Orientation distribution, controlling the fraction of loaded collagen fibers at a given stretch, was identified as a key feature governing anisotropic tissue strength for all patient cohorts.

Failure of the Porcine Ascending Aorta: Multidirectional Experiments and a Unifying Microstructural Model

The ascending thoracic aorta is poorly understood mechanically, especially its risk of dissection. To make better predictions of dissection risk, more information about the multidimensional failure behavior of the tissue is needed, and this information must be incorporated into an appropriate theoretical/computational model. Toward the creation of such a model, uniaxial, equibiaxial, peel, and shear lap tests were performed on healthy porcine ascending aorta samples. Uniaxial and equibiaxial tests showed anisotropy with greater stiffness and strength in the circumferential direction. Shear lap tests showed catastrophic failure at shear stresses (150-200 kPa) much lower than uniaxial tests (750-2500 kPa), consistent with the low peel tension (∼60 mN/mm). A novel multiscale computational model, including both prefailure and failure mechanics of the aorta, was developed. The microstructural part of the model included contributions from a collagen-reinforced elastin sheet and interlamellar connections representing fibrillin and smooth muscle. Components were represented as nonlinear fibers that failed at a critical stretch. Multiscale simulations of the different experiments were performed, and the model, appropriately specified, agreed well with all experimental data, representing a uniquely complete structure-based description of aorta mechanics. In addition, our experiments and model demonstrate the very low strength of the aorta in radial shear, suggesting an important possible mechanism for aortic dissection.

Layer-dependent role of collagen recruitment during loading of the rat bladder wall

In this work, we re-evaluated long-standing conjectures as to the source of the exceptionally large compliance of the bladder wall. Whereas these conjectures were based on indirect measures of loading mechanisms, in this work we take advantage of advances in bioimaging to directly assess collagen fibers and wall architecture during biaxial loading. A custom biaxial mechanical testing system compatible with multiphoton microscopy was used to directly measure the layer-dependent collagen fiber recruitment in bladder tissue from 9 male Fischer rats (4 adult and 5 aged). As for other soft tissues, the bladder loading curve was exponential in shape and could be divided into toe, transition and high stress regimes. The relationship between collagen recruitment and loading curves was evaluated in the context of the inner (lamina propria) and outer (detrusor smooth muscle) layers. The large extensibility of the bladder was found to be possible due to folds in the wall (rugae) that provide a mechanism for low resistance flattening without any discernible recruitment of collagen fibers throughout the toe regime. For more extensible bladders, as the loading extended into the transition regime, a gradual coordinated recruitment of collagen fibers between the lamina propria layer and detrusor smooth muscle layer was found. A second important finding was that wall extensibility could be lost by premature recruitment of collagen in the outer wall that cut short the toe region. This change was correlated with age. This work provides, for the first time, a mechanistic understanding of the role of collagen recruitment in determining bladder extensibility and capacitance.

Multiphoton microscopy observations of 3D elastin and collagen fiber microstructure changes during pressurization in aortic media

Elastin and collagen fibers play important roles in the mechanical properties of aortic media. Because knowledge of local fiber structures is required for detailed analysis of blood vessel wall mechanics, we investigated 3D microstructures of elastin and collagen fibers in thoracic aortas and monitored changes during pressurization. Using multiphoton microscopy, autofluorescence images from elastin and second harmonic generation signals from collagen were acquired in media from rabbit thoracic aortas that were stretched biaxially to restore physiological dimensions. Both elastin and collagen fibers were observed in all longitudinal-circumferential plane images, whereas alternate bright and dark layers were observed along the radial direction and were recognized as elastic laminas (ELs) and smooth muscle-rich layers (SMLs), respectively. Elastin and collagen fibers are mainly oriented in the circumferential direction, and waviness of collagen fibers was significantly higher than that of elastin fibers. Collagen fibers were more undulated in longitudinal than in radial direction, whereas undulation of elastin fibers was equibiaxial. Changes in waviness of collagen fibers during pressurization were then evaluated using 2-dimensional fast Fourier transform in mouse aortas, and indices of waviness of collagen fibers decreased with increases in intraluminal pressure. These indices also showed that collagen fibers in SMLs became straight at lower intraluminal pressures than those in EL, indicating that SMLs stretched more than ELs. These results indicate that deformation of the aorta due to pressurization is complicated because of the heterogeneity of tissue layers and differences in elastic properties of ELs, SMLs, and surrounding collagen and elastin.

A structural finite element model for lamellar unit of aortic media indicates heterogeneous stress field after collagen recruitment

Incorporation of collagen structural information into the study of biomechanical behavior of ascending thoracic aortic (ATA) wall tissue should provide better insight into the pathophysiology of ATA. Structurally motivated constitutive models that include fiber dispersion and recruitment can successfully capture overall mechanical response of the arterial wall tissue. However, these models cannot examine local microarchitectural features of the collagen network, such as the effect of fiber disruptions and interaction between fibrous and non-fibrous components, which may influence emergent biomechanical properties of the tissue. Motivated by this need, we developed a finite element based three-dimensional structural model of the lamellar units of the ATA media that directly incorporates the collagen fiber microarchitecture. The fiber architecture was computer generated utilizing network features, namely fiber orientation distribution, intersection density and areal concentration, obtained from image analysis of multiphoton microscopy images taken from human aneurysmal ascending thoracic aortic media specimens with bicuspid aortic valve (BAV) phenotype. Our model reproduces the typical J-shaped constitutive response of the aortic wall tissue. We found that the stress state in the non-fibrous matrix was homogeneous until the collagen fibers were recruited, but became highly heterogeneous after that event. The degree of heterogeneity was dependent upon local network architecture with high stresses observed near disrupted fibers. The magnitude of non-fibrous matrix stress at higher stretch levels was negatively correlated with local fiber density. The localized stress concentrations, elucidated by this model, may be a factor in the degenerative changes in aneurysmal ATA tissue.

Mechanical assessment of arterial dissection in health and disease: Advancements and challenges

Arterial dissection involves a complex series of coupled biomechanical events. The past two decades have witnessed great advances in the understanding of the intrinsic mechanism for dissection initiation, and hence in the development of novel therapeutic strategies for surgical repair. This is due in part to the profound advancements in characterizing emerging behaviors of dissection using state-of-the-art tools in experimental and computational biomechanics. In addition, researchers have identified the important role of the microstructure in determining the tissue׳s fracture modality during dissection propagation. In this review article, we highlight a variety of approaches in terms of biomechanical measurements, computational modeling and histological/microstructural analysis used to characterize a dissection that propagates in healthy and diseased arteries. Notable findings with quantitative mechanical data are reviewed. We conclude by discussing some unsolved problems that are of interest for future research.

Constitutive modeling of ascending thoracic aortic aneurysms using microstructural parameters

Ascending thoracic aortic aneurysm (ATAA) has been associated with diminished biomechanical strength and disruption in the collagen fiber microarchitecture. Additionally, the congenital bicuspid aortic valve (BAV) leads to a distinct extracellular matrix structure that may be related to ATAA development at an earlier age than degenerative aneurysms arising in patients with the morphological normal tricuspid aortic valve (TAV). The purpose of this study was to model the fiber-reinforced mechanical response of ATAA specimens from patients with either BAV or TAV. This was achieved by combining image-analysis derived parameters of collagen fiber dispersion and alignment with tensile testing data. Then, numerical simulations were performed to assess the role of anisotropic constitutive formulation on the wall stress distribution of aneurysmal aorta. Results indicate that both BAV ATAA and TAV ATAA have altered collagen fiber architecture in the medial plane of experimentally-dissected aortic tissues when compared to normal ascending aortic specimens. The study findings highlight that differences in the collagen fiber distribution mostly influences the resulting wall stress distribution rather than the peak stress. We conclude that fiber-reinforced constitutive modeling that takes into account the collagen fiber defect inherent to the aneurysmal ascending aorta is paramount for accurate finite element predictions and ultimately for biomechanical-based indicators to reliably distinguish the more from the less 'malignant' ATAAs.

Crack Propagation and Its Shear Mechanisms in the Bovine Descending Aorta

Aortic dissection and rupture may involve circumferential shear stress in the circumferential-longitudinal plane. Inflation of bovine descending aortic ring specimens provides evidence of such shear from the non-uniform circumferential distortion of radial lines drawn on the circumferential-radial ring face. Delamination without tensile peeling induces cracks that propagate nearly circumferentially in the circumferential-longitudinal plane from the root of a radial cut representing rupture initiation in a ring. Translational shear deformation tests of small rectangular aortic wall blocks in the circumferential and longitudinal direction measure the consequences of such shear on substructures in the aortic wall, in particular the collagen fibers. The two directions of shear deformation produce no statistical difference in the shear stress response of the wall. Possibly, the interfiber connections between collagen fibers are put into tension by either translational shear deformation so that the stress measured reflects the tensile response of these connections. Wall rupture may involve failure of these connections; such failure is supported by the voids parallel to the collagen fibers observed in a histological study after translational shear. Further, interstitial fluid is redistributed by shear as evidenced by the measured weight loss of a set of specimens during the translational shear of blocks. Because the mass changes, mathematical modeling of aortic tissue in vitro as incompressible is an approximation. These observations suggest that no simple modification of classical rupture theories, whether based on energy functions, stress or strain, suffices to predict the rupture of hydrated soft biological tissue that has complex substructures.

The Effect of Size and Location of Tears in the Supraspinatus Tendon on Potential Tear Propagation

Rotator cuff tears are a common problem in patients over the age of 50 yr. Tear propagation is a potential contributing factor to the failure of physical therapy for treating rotator cuff tears, thus requiring surgical intervention. However, the evolution of tears within the rotator cuff is not well understood yet. The objective of this study is to establish a computational model to quantify initiation of tear propagation in the supraspinatus tendon and examine the effect of tear size and location. A 3D finite element (FE) model of the supraspinatus tendon was constructed from images of a healthy cadaveric tendon. A tear of varying length was placed at six different locations within the tendon. A fiber-reinforced Mooney–Rivlin material model with spatial variation in material properties along the anterior–posterior (AP) axis was utilized to obtain the stress state of the computational model under uniaxial stretch. Material parameters were calibrated by comparing computational and experimental stress–strain response and used to validate the computational model. The stress state of the computational model was contrasted against the spatially varying material strength to predict the critical applied stretch at which a tear starts propagating further. It was found that maximum principal stress (as well as the strain) was localized at the tips of the tear. The computed critical stretch was significantly lower for the posterior tip of the tear than for the anterior tip suggesting a propensity to propagate posteriorly. Onset of tear propagation was strongly correlated with local material strength and stiffness in the vicinity of the tear tip. Further, presence of a stress-shielded zone along the edges of the tear was observed. This study illustrates the complex interplay between geometry and material properties of tendon up to the initiation of tear propagation. Future work will examine the evolution of tears during the propagation process as well as under more complex loading scenarios.

Structure-based constitutive model can accurately predict planar biaxial properties of aortic wall tissue

Structure-based constitutive models might help in exploring mechanisms by which arterial wall histology is linked to wall mechanics. This study aims to validate a recently proposed structure-based constitutive model. Specifically, the model's ability to predict mechanical biaxial response of porcine aortic tissue with predefined collagen structure was tested. Histological slices from porcine thoracic aorta wall (n=9) were automatically processed to quantify the collagen fiber organization, and mechanical testing identified the non-linear properties of the wall samples (n=18) over a wide range of biaxial stretches. Histological and mechanical experimental data were used to identify the model parameters of a recently proposed multi-scale constitutive description for arterial layers. The model predictive capability was tested with respect to interpolation and extrapolation. Collagen in the media was predominantly aligned in circumferential direction (planar von Mises distribution with concentration parameter bM=1.03 ± 0.23), and its coherence decreased gradually from the luminal to the abluminal tissue layers (inner media, b=1.54 ± 0.40; outer media, b=0.72 ± 0.20). In contrast, the collagen in the adventitia was aligned almost isotropically (bA=0.27 ± 0.11), and no features, such as families of coherent fibers, were identified. The applied constitutive model captured the aorta biaxial properties accurately (coefficient of determination R(2)=0.95 ± 0.03) over the entire range of biaxial deformations and with physically meaningful model parameters. Good predictive properties, well outside the parameter identification space, were observed (R(2)=0.92 ± 0.04). Multi-scale constitutive models equipped with realistic micro-histological data can predict macroscopic non-linear aorta wall properties. Collagen largely defines already low strain properties of media, which explains the origin of wall anisotropy seen at this strain level. The structure and mechanical properties of adventitia are well designed to protect the media from axial and circumferential overloads.

Effect of aneurysm on biomechanical properties of "radially-oriented" collagen fibers in human ascending thoracic aortic media

We recently reported a mechanistic model to link micro-architectural information to the delamination strength (Sd) of human ascending thoracic aorta (ATA). That analysis demonstrated that the number density (N) and failure energy (Uf) of the radially-oriented collagen fibers contribute to the Sd of both aneurysmal (ATAA) and non-aneurysmal (CTRL-ATA) aortic tissue. Among the set of ATAA samples, we studied specimens from patients displaying bicuspid (BAV) and tricuspid aortic valve (TAV) morphologic phenotypes. Results from our prior work were based on the assumption that the Uf was independent of dissection direction. In the current study, we excluded that assumption and hypothesized that Uf correlates with the Sd of ATAA. To test the hypothesis, we used previously-reported experimentally-determined Sd measurements and N of radially-oriented collagen fibers as input in our validated mechanistic model to calculate Uf for BAV-ATAA, TAV-ATAA and CTRL-ATA tissue specimens. The results of our analysis revealed that Uf is significantly lower for both BAV-ATAA and TAV-ATAA compared to CTRL-ATA cases, and does not differ between BAV-ATAA and TAV-ATAA. Furthermore, we found that Uf is consistent between circumferential-radial and longitudinal-radial planes in either of BAV-ATAA, TAV-ATAA or CTRL-ATA specimens. These findings employ a novel mechanistic model to increase our understanding of the putative interrelationship between biomechanical properties, extracellular matrix biology, and failure energy of aortic dissection.

Varying degrees of nonlinear mechanical behavior arising from geometric differences of urogynecological meshes

Synthetic polypropylene meshes were designed to restore pelvic organ support for women suffering from pelvic organ prolapse; however, the FDA released two notifications regarding potential complications associated with mesh implantation. Our aim was to characterize the structural properties of Restorelle and UltraPro subjected to uniaxial tension along perpendicular directions, and then model the tensile behavior of these meshes utilizing a co-rotational finite element model, with an imbedded linear or fiber-recruitment local stress-strain relationship. Both meshes exhibited a highly nonlinear stress-strain behavior; Restorelle had no significant differences between the two perpendicular directions, while UltraPro had a 93% difference in the low (initial) stiffness (p=0.009) between loading directions. Our model predicted that early alignment of the mesh segments in the loading direction and subsequent stretching could explain the observed nonlinear tensile behavior. However, a nonlinear stress-strain response in the stretching regime, that may be inherent to the mesh segment, was required to better capture experimental results. Utilizing a nonlinear fiber recruitment model with two parameters A and B, we observed improved agreement between the simulations and the experimental results. An inverse analysis found A=120 MPa and B=1.75 for Restorelle (RMSE=0.36). This approach yielded A=30 MPa and B=3.5 for UltraPro along one direction (RMSE=0.652), while the perpendicular orientation resulted in A=130 MPa and B=4.75 (RMSE=4.36). From the uniaxial protocol, Restorelle was found to have little variance in structural properties along these two perpendicular directions; however, UltraPro was found to behave anisotropically.