Layer- and Direction-Specific Material Properties, Extreme Extensibility and Ultimate Material Strength of Human Abdominal Aorta and Aneurysm: A Uniaxial Extension Study

What does the slope of stress-stretch curves tell us about vascular tissue response?

We examined a group of 50 uniaxial stress-stretch curves obtained from human ascending aortic aneurysm tissues. The curves were believed to be associated with elastic response because the stress is monotonically increasing in all curves, and so is the slope. However, 26 curves exhibit exponential-like slope while the remaining 24 curves have sigmoid slopes. We hypothesized that the slope patterns stemmed from collage waviness distribution. A structural constitutive model was introduced to describe the responses. The model employed a unimodal density function to describe the waviness distribution, from which a two-phase response ensued. In the first phase the slope is quasi-exponential, and in the second phase the slope is sigmoid. The model fitted all 50 curves perfectly well. An exponential model was also introduced for a comparison. The model fitted the curves of quasi-exponential slope generally well, but performed worse over the curves of sigmoid slope. The work suggests that the slope may encode significant information about collagen waviness, and underscores a limitation of exponential-based models.

In Vivo Quantification of Ascending Thoracic Aortic Aneurysm Wall Stretch Using MRI: Relationship to Repair Threshold Diameter and Ex Vivo Wall Failure Behavior

Ascending thoracic aortic aneurysms (aTAAs) can lead to life-threatening dissection and rupture. Recent studies have highlighted aTAA mechanical properties as relevant factors associated with progression. The aim of this study was to quantify in vivo aortic wall stretch in healthy participants and aTAA patients using displacement encoding with stimulated echoes (DENSE) magnetic resonance imaging. Moreover, aTAA wall stretch between surgical and nonsurgical patients was investigated. Finally, DENSE measurements were compared to reference-standard mechanical testing on aTAA specimens from surgical repairs. In total, 18 subjects were recruited, six healthy participants and 12 aTAA patients, for this prospective study. Electrocardiogram-gated DENSE imaging was performed to measure systole-diastole wall stretch, as well as the ratio of aTAA stretch to unaffected descending thoracic aorta stretch. Free-breathing and breath-hold DENSE protocols were used. Uniaxial tensile testing-measured indices were correlated to DENSE measurements in five harvested specimens. in vivo aortic wall stretch was significantly lower in aTAA compared to healthy subjects (1.75±1.44% versus 5.28±1.92%, respectively, P = 0.0004). There was no correlation between stretch and maximum aTAA diameter (P = 0.56). The ratio of aTAA to unaffected thoracic aorta wall stretch was significantly lower in surgical candidates compared to nonsurgical candidates (0.993±0.011 versus 1.017±0.016, respectively, P = 0.0442). Finally, in vivo aTAA wall stretch correlated to wall failure stress and peak modulus of the intima (P = 0.017 and P = 0.034, respectively), while the stretch ratio correlated to whole-wall thickness failure stretch and stress (P = 0.013 and P = 0.040, respectively). Aortic DENSE has the potential to assess differences in aTAA mechanical properties and progressions.

Comparison of Biomechanical and Microstructural Properties of Aortic Graft Materials in Aortic Repair Surgeries

Mechanical mismatch between native aortas and aortic grafts can induce graft failure. This study aims to compare the mechanical and microstructural properties of different graft materials used in aortic repair surgeries with those of normal and dissected human ascending aortas. Five types of materials including normal aorta (n = 10), dissected aorta (n = 6), human pericardium (n = 8), bovine pericardium (n = 8) and Dacron graft (n = 5) were collected to perform uniaxial tensile testing to determine their material stiffness, and ultimate strength/stretch. The elastin and collagen contents in four tissue groups except for Dacron were quantified by histological examinations, while the material ultrastructure of five material groups was visualized by scanning electron microscope. Statistical results showed that three graft materials including Dacron, human pericardium and bovine pericardium had significantly higher ultimate strength and stiffness than both normal and dissected aortas. Human and bovine pericardia had significantly lower ultimate stretch than native aortas. Histological examinations revealed that normal and diseased aortic tissues had a significantly higher content of elastic fiber than two pericardial tissues, but less collagen fiber content. All four tissue groups exhibited lamellar fiber ultrastructure, with aortic tissues possessing thinner lamella. Dacron was composed of densely coalesced polyethylene terephthalate fibers in thick bundles. Aortic graft materials with denser fiber ultrastructure and/or higher content of collagen fiber than native aortic tissues, exhibited higher ultimate strength and stiffness. This information provides a basis to understand the mechanical failure of aortic grafts, and inspire the design of biomimetic aortic grafts.

Exploiting light-based 3D-printing for the fabrication of mechanically enhanced, patient-specific aortic grafts

Despite polyester vascular grafts being routinely used in life-saving aortic aneurysm surgeries, they are less compliant than the healthy, native human aorta. This mismatch in mechanical behaviour has been associated with disruption of haemodynamics contributing to several long-term cardiovascular complications. Moreover, current fabrication approaches mean that opportunities to personalise grafts to the individual anatomical features are limited. Various modifications to graft design have been investigated to overcome such limitations; yet optimal graft functionality remains to be achieved. This study reports on the development and characterisation of an alternative vascular graft material. An alginate:PEGDA (AL:PE) interpenetrating polymer network (IPN) hydrogel has been produced with uniaxial tensile tests revealing similar strength and stiffness (0.39 ± 0.05 MPa and 1.61 ± 0.19 MPa, respectively) to the human aorta. Moreover, AL:PE tubular conduits of similar geometrical dimensions to segments of the aorta were produced, either via conventional moulding methods or stereolithography (SLA) 3D-printing. While both fabrication methods successfully demonstrated AL:PE hydrogel production, SLA 3D-printing was more easily adaptable to the fabrication of complex structures without the need of specific moulds or further post-processing. Additionally, most 3D-printed AL:PE hydrogel tubular conduits sustained, without failure, compression up to 50% their outer diameter and returned to their original shape upon load removal, thereby exhibiting promising behaviour that could withstand pulsatile pressure in vivo. Overall, these results suggest that this AL:PE IPN hydrogel formulation in combination with 3D-printing, has great potential for accelerating progress towards personalised and mechanically-matched aortic grafts.

Material strengths of shear-induced platelet aggregation clots and coagulation clots

Arterial occlusion by thrombosis is the immediate cause of some strokes, heart attacks, and peripheral artery disease. Most prior studies assume that coagulation creates the thrombus. However, a contradiction arises as whole blood (WB) clots from coagulation are too weak to stop arterial blood pressures (> 150 mmHg). We measure the material mechanical properties of elasticity and ultimate strength for Shear-Induced Platelet Aggregation (SIPA) type clots, that form under stenotic arterial hemodynamics in comparison with coagulation clots. The ultimate strength of SIPA clots averaged 4.6 ± 1.3 kPa, while WB coagulation clots had a strength of 0.63 ± 0.3 kPa (p < 0.05). The elastic modulus of SIPA clots was 3.8 ± 1.5 kPa at 1 Hz and 0.5 mm displacement, or 2.8 times higher than WB coagulation clots (1.3 ± 1.2 kPa, p < 0.0001). This study shows that the SIPA thrombi, formed quickly under high shear hemodynamics, is seven-fold stronger and three-fold stiffer compared to WB coagulation clots. A force balance calculation shows a SIPA clot has the strength to resist arterial pressure with a short length of less than 2 mm, consistent with coronary pathology.

Tri-layered constitutive modelling unveils functional differences between the pig ascending and lower thoracic aorta

The arterial wall's tri-layered macroscopic and layer-specific microscopic structure determine its mechanical properties, which vary at different arterial locations. Combining layer-specific mechanical data and tri-layered modelling, this study aimed to characterise functional differences between the pig ascending (AA) and lower thoracic aorta (LTA). AA and LTA segments were obtained for n=9 pigs. For each location, circumferentially and axially oriented intact wall and isolated layer strips were tested uniaxially and the layer-specific mechanical response modelled using a hyperelastic strain energy function. Then, layer-specific constitutive relations and intact wall mechanical data were combined to develop a tri-layered model of an AA and LTA cylindrical vessel, accounting for the layer-specific residual stresses. AA and LTA behaviours were then characterised for in vivo pressure ranges while stretched axially to in vivo length. The media dominated the AA response, bearing>2/3 of the circumferential load both at physiological (100 mmHg) and hypertensive pressures (160 mmHg). The LTA media bore most of the circumferential load at physiological pressure only (57±7% at 100 mmHg), while adventitia and media load bearings were comparable at 160 mmHg. Furthermore, increased axial elongation affected the media/adventitia load-bearing only at the LTA. The pig AA and LTA presented strong functional differences, likely reflecting their different roles in the circulation. The media-dominated compliant and anisotropic AA stores large amounts of elastic energy in response to both circumferential and axial deformations, which maximises diastolic recoiling function. This function is reduced at the LTA, where the adventitia shields the artery against supra-physiological circumferential and axial loads.

Mechanical and histological characteristics of aortic dissection tissues

AIMS

This study investigated the association between the macroscopic mechanical response of aortic dissection (AoD) flap, its fibre features, and patient physiological features and clinical presentations.

METHODS

Uniaxial test was performed with tissue strips in both circumferential and longitudinal directions from 35 patients with (AoD:CC) and without (AoD:w/oCC) cerebral/coronary complications, and 19 patients with rheumatic or valve-related heart diseases (RH). A Bayesian inference framework was used to estimate the expectation of material constants (C₁, D₁, and D₂) of the modified Mooney-Rivlin strain energy density function. Histological examination was used to visualise the elastin and collagen in the tissue strips and image processing was performed to quantify their area percentages, fibre misalignment and waviness.

RESULTS

The elastin area percentage was negatively associated with age (p = 0.008), while collagen increased about 6% from age 40 to 70 (p = 0.03). Elastin fibre was less dispersed and wavier in old patients and no significant association was found between patient age and collagen fibre dispersion or waviness. Features of fibrous microstructures, either elastin or collagen, were comparable between AoD:CC and AoD:w/oCC group. Elastin and collagen area percentages were positively correlated with C₁ and D₂, respectively, while the elastin and collagen waviness were negatively correlated with C₁ and D₂, respectively. Elastin dispersion was negatively correlated to D₂. Multivariate analysis showed that D₂ was an effective parameter which could differentiate patient groups with different age and clinical presentations, as well as the direction of tissue strip.

CONCLUSION

Fibre dispersion and waviness in the aortic dissection flap changed with patient age and clinical presentations, and these can be captured by the material constants in the strain energy density function.

STATEMENT OF SIGNIFICANCE

Aortic dissection (AoD) is a severe cardiovascular disease. Understanding the mechanical property of intimal flap is essential for its risk evaluation. In this study, mechanical testing and histology examination were combined to quantify the relationship between mechanical presentations and microstructure features. A Bayesian inference framework was employed to estimate the expectation of the material constants in the modified Mooney-Rivlin constitutive equation. It was found that fibre dispersion and waviness in the AoD flap changed with patient age and clinical presentations, and these could be captured by the material constants. This study firstly demonstrated that the expectation of material constants can be used to characterise tissue microstructures and differentiate patients with different clinical presentations.

Cushioned on the way up, controlled on the way down during resuscitative endovascular balloon occlusion of the aorta (REBOA): investigating a novel compliant balloon design for optimizing safe overinflation combined with partial REBOA ability

Background

There are a variety of devices capable of performing resuscitative endovascular balloon occlusion of the aorta (REBOA), with most containing compliant balloon material. While compliant material is ideal for balloon inflation due to its “cushioning” effect, it can be problematic to “control” during deflation. The COBRA-OS (Control Of Bleeding, Resuscitation, Arterial Occlusion System) was designed to optimize inflation and deflation of its compliant balloon and was tested in vitro and in vivo with respect to its overinflation and partial REBOA abilities.

Methods

For overinflation, the COBRA-OS was inflated in three differently sized inner diameter (ID) vinyl tubes until balloon rupture. It was then overinflated in six harvested swine aortas and in all three REBOA zones of three anesthetized swine. For partial REBOA, the COBRA-OS underwent incremental deflation in a pulsatile benchtop aortic model and in zone 1 of three anesthetized swine.

Results

For overinflation, compared with the known aortic rupture threshold of 4 atm, the COBRA-OS exceeded this value in only the smallest of the vinyl tubes: 8 mm ID tube, 6.5 atm; 9.5 mm ID tube, 3.5 atm; 13 mm ID tube, 1.5 atm. It also demonstrated greater than 500% overinflation ability without aortic damage in vitro and caused no aortic damage when inflated to maximum inflation volume in vivo. For partial REBOA, the COBRA-OS was able to provide a titration window of between 3 mL and 4 mL in both the pulsatile vascular model (3.4±0.12 mL) and anesthetized swine (3.8±0.35 mL).

Discussion

The COBRA-OS demonstrated the ability to have a cushioning effect during inflation combined with titration control on deflation in vitro and in vivo. This study suggests that despite its balloon compliance, both safe overinflation and partial REBOA can be successfully achieved with the COBRA-OS.

Level of evidence

Basic science.

Micromechanical and Ultrastructural Properties of Abdominal Aortic Aneurysms

Abdominal aortic aneurysms are a common condition of uncertain pathogenesis that can rupture if left untreated. Current recommended thresholds for planned repair are empirical and based entirely on diameter. It has been observed that some aneurysms rupture before reaching the threshold for repair whilst other larger aneurysms do not rupture. It is likely that geometry is not the only factor influencing rupture risk. Biomechanical indices aiming to improve and personalise rupture risk prediction require, amongst other things, knowledge of the material properties of the tissue and realistic constitutive models. These depend on the composition and organisation of the vessel wall which has been shown to undergo drastic changes with aneurysmal degeneration, with loss of elastin, smooth muscle cells, and an accumulation of isotropically arranged collagen. Most aneurysms are lined with intraluminal thrombus, which has an uncertain effect on the underlying vessel wall, with some authors demonstrating a reduction in wall stress and others a reduction in wall strength. The majority of studies investigating biomechanical properties of ex vivo abdominal aortic aneurysm tissues have used low-resolution techniques, such as tensile testing, able to measure the global material properties at the macroscale. High-resolution engineering techniques such as nanoindentation and atomic force microscopy have been modified for use in soft biological tissues and applied to vascular tissues with promising results. These techniques have the potential to advance the understanding and improve the management of abdominal aortic aneurysmal disease.

Mechanical properties of tissue-mimicking composites formed by material jetting additive manufacturing

Capitalizing on features including high resolution, smooth surface finish, large build volume, and simultaneous multi-color/multi-material printing, material jetting additive manufacturing enables the fabrication of full-scale anatomic models. The ability to print materials that resemble relevant, compliant tissues has especially motivated implementation of material jetting for patient-specific surgical planning or training models. In an effort to broaden the material selection for the material jetting process, and to provide materials that more closely mimic the functional needs for a wider variety of tissues, a composite material system is explored that uses non-curing fluid dispersed into a photo-curable medium. The material properties of the composites are examined through both thermal and mechanical analysis (dynamic mechanical analysis, Shore hardness testing, puncture testing, and tensile testing). Higher contributions of non-curing fluid generally reduce part strength and stiffness, and exponential and second-order polynomial expressions are appropriate fits for many of the mechanical properties as functions of non-curing fluid concentration. Through the fundamental exploration of the impact of an added diluent on material properties, the study advances knowledge on the process-property relationship for multi-material jetting. Additionally, better understanding of the mechanical property space offered by these materials will expand the capabilities of material jetting in the context of biomedical applications. The collection of mechanical properties serve as reference data sets to facilitate quicker screening for tissue-mimicking, medical models.

Vessel structural stress mediates aortic media degeneration in bicuspid aortopathy: New insights based on patient-specific fluid-structure interaction analysis

This study aimed to assess the relationship between local mechanical stimuli and regional aortic tissue degeneration using fluid-structure interaction (FSI) analysis in patients with bicuspid aortic valve (BAV) disease. Nine patients underwent ascending aortic replacement were recruited. Tissues were collected to evaluate the pathology features in four regions, greater curvature (GC-region), posterior (P-region), anterior (A-region), and lesser curvature (LC-region). FSI analysis was performed to quantify vessel structural stress (VSS) and flow-induced parameters, including wall shear stress (WSS), oscillatory shear index (OSI), and particle relative residence time (RRT). The correlation between these biomechanical metrics and tissue degeneration was analyzed. Elastin in the medial layer and media thickness were thinnest and the gap between fibers was biggest in the GC-region, followed by the P-region and A-region, while the elastin and media thickness were thickest and the gap smallest in the LC-region. The collagen deposition followed a pattern with the biggest in the GC-region and least in the LC-region. There is a strong negative correlation between mean or peak VSS and elastin thickness in the arterial wall in the GC-region (r = -0.917; p = 0.001 and r = -0.899; p = 0.001), A-region (r = -0.748; p = 0.020 and r = -0.700; p = 0.036) and P-region (r = -0.773; p = 0.014 and r = -0.769; p = 0.015), and between mean VSS and fiber distance in the A-region (r = -0.702, p = 0.035). Moreover, strong negative correlation between mean or peak VSS and media thickness was also observed. No correlation was found between WSS, OSI, and RRT and aortic tissue degeneration in these four regions. These findings indicate that increased VSS correlated with local elastin degradation and aortic media degeneration, implying that it could be a potential biomechanical parameter for a refined risk stratification for patients with BAV.

Synthetic hydrogels as blood clot mimicking wound healing materials

Excessive bleeding-or hemorrhage-causes millions of civilian and non-civilian casualties every year. Additionally, wound sequelae, such as infections, are a significant source of chronic morbidity, even if the initial bleeding is successfully stopped. To treat acute and chronic wounds, numerous wound healing materials have been identified, tested, and adopted. Among them are topical dressings, such as gauzes, as well as natural and biomimetic materials. However, none of these materials successfully mimic the complex and dynamic properties of the body's own wound healing material: the blood clot. Specifically, blood clots exhibit complex mechanical and biochemical properties that vary across spatial and temporal scales to guide the wound healing response, which make them the ideal wound healing material. In this manuscript, we review blood clots' complex mechanical and biochemical properties, review current wound healing materials, and identify opportunities where new materials can provide additional functionality, with a specific focus on hydrogels. We highlight recent developments in synthetic hydrogels that make them capable of mimicking a larger subset of blood clot features: as plugs and as stimuli for tissue repair. We conclude that future hydrogel materials designed to mimic blood clot biochemistry, mechanics, and architecture can be combined with exciting platelet-like particles to serve as hemostats that also promote the biological wound healing response. Thus, we believe synthetic hydrogels are ideal candidates to address the clear need for better wound healing materials.

Computation of a probabilistic and anisotropic failure metric on the aortic wall using a machine learning-based surrogate model

Scalar-valued failure metrics are commonly used to assess the risk of aortic aneurysm rupture and dissection, which occurs under hypertensive blood pressures brought on by extreme emotional or physical stress. To compute failure metrics under an elevated blood pressure, a classical patient-specific computer model consists of multiple computation steps involving inverse and forward analyses. These classical procedures may be impractical for time-sensitive clinical applications that require prompt feedback to clinicians. In this study, we developed a machine learning-based surrogate model to directly predict a probabilistic and anisotropic failure metric, namely failure probability (FP), on the aortic wall using aorta geometries at the systolic and diastolic phases. Ascending thoracic aortic aneurysm (ATAA) geometries of 60 patients were obtained from their CT scans, and biaxial mechanical testing data of ATAA tissues from 79 patients were collected. Finite element simulations were used to generate datasets for training, validation, and testing of the ML-surrogate model. The testing results demonstrated that the ML-surrogate can compute the maximum FP failure metric, with 0.42% normalized mean absolute error, in 1 s. To compare the performance of the ML-predicted probabilistic FP metric with other isotropic or deterministic metrics, a numerical case study was performed using synthetic "baseline" data. Our results showed that the probabilistic FP metric had more discriminative power than the deterministic Tsai-Hill metric, isotropic maximum principal stress, and aortic diameter criterion.

Development of an FEA framework for analysis of subject-specific aortic compliance based on 4D flow MRI

This paper presents a subject-specific in-silico framework in which we uncover the relationship between the spatially varying constituents of the aorta and the non-linear compliance of the vessel during the cardiac cycle uncovered through our MRI investigations. A microstructurally motivated constitutive model is developed, and simulations reveal that internal vessel contractility, due to pre-stretched elastin and actively generated smooth muscle cell stress, must be incorporated, along with collagen strain stiffening, in order to accurately predict the non-linear pressure-area relationship observed in-vivo. Modelling of elastin and smooth muscle cell contractility allows for the identification of the reference vessel configuration at zero-lumen pressure, in addition to accurately predicting high- and low-compliance regimes under a physiological range of pressures. This modelling approach is also shown to capture the key features of elastin digestion and SMC activation experiments. The volume fractions of the constituent components of the aortic material model were computed so that the in-silico pressure-area curves accurately predict the corresponding MRI data at each location. Simulations reveal that collagen and smooth muscle volume fractions increase distally, while elastin volume fraction decreases distally, consistent with reported histological data. Furthermore, the strain at which collagen transitions from low to high stiffness is lower in the abdominal aorta, again supporting the histological finding that collagen waviness is lower distally. The analyses presented in this paper provide new insights into the heterogeneous structure-function relationship that underlies aortic biomechanics. Furthermore, this subject-specific MRI/FEA methodology provides a foundation for personalised in-silico clinical analysis and tailored aortic device development. STATEMENT OF SIGNIFICANCE: This study provides a significant advance in in-silico medicine by capturing the structure/function relationship of the subject-specific human aorta presented in our previous MRI analyses. A physiologically based aortic constitutive model is developed, and simulations reveal that internal vessel contractility must be incorporated, along with collagen strain stiffening, to accurately predict the in-vivo non-linear pressure-area relationship. Furthermore, this is the first subject-specific model to predict spatial variation in the volume fractions of aortic wall constituents. Previous studies perform phenomenological hyperelastic curve fits to medical imaging data and ignore the prestress contribution of elastin, collagen, and SMCs and the associated zero-pressure reference state of the vessel. This novel MRI/FEA framework can be used as an in-silico diagnostic tool for the early stage detection of aortic pathologies.

The influence of sample geometry and size on porcine aortic material properties from uniaxial tensile tests using custom-designed tissue cutters, clamps and molds

The aim of this study was to identify the influence of specimen geometry and size on the results of aortic uniaxial tensile tests using custom-designed tissue cutters, clamps and molds. Six descending thoracic aortas from pigs were used for rectangular sample tests, in which the circumferential and axial specimens had widths of 6 mm, 8 mm and 10 mm. The other six aortas were used for the dog-bone-shaped sample tests and were punched into circumferential, axial and oblique specimens with widths of 2 mm, 4 mm and 6 mm. We performed uniaxial tensile tests on the specimens and compared the test results. The results showed that mid-sample failure occurred in 85.2% of the dog-bone-shaped specimens and in 11.1% of the rectangular samples, which could be caused by Saint-Venant’s principle. Therefore, rectangular specimens were not suitable for aortic uniaxial tensile testing performed until rupture. The results also showed that the size effect of the aorta conformed to Weibull theory, and dog-bone-shaped specimens with a width of 4 mm were the optimal choice for aortic uniaxial tensile testing performed until rupture.

A whole blood thrombus mimic: Constitutive behavior under simple shear

Deep vein thrombosis and pulmonary embolism affect 300,000-600,000 patients each year in the US. To better understand the highly mechanical pathophysiology of pulmonary embolism, we set out to develop an in-vitro thrombus mimic and to test this mimic under large deformation simple shear. In addition to reporting on the mechanics of our mimics under simple shear, we explore the sensitivity of their mechanics to coagulation conditions and blood storage time, and compare three hyperelastic material models for their ability to fit our data. We found that thrombus mimics made from whole blood demonstrate strain-stiffening, a negative Poynting effect, and hysteresis when tested quasi-statically to 50% strain under simple shear. Additionally, we found that the stiffness of these mimics does not significantly vary with coagulation conditions or blood storage times. Of the three hyperelastic constitutive models that we tested, the Ogden model provided the best fits to both shear stress and normal stress. In conclusion, we developed a robust protocol to generate regularly-shaped, homogeneous thrombus mimics that lend themselves to simple shear testing under large deformation. Future studies will extend our model to include the effect of maturation and explore its fracture properties toward a better understanding of embolization.

A Novel Anisotropic Failure Criterion With Dispersed Fiber Orientations for Aortic Tissues

Accurate failure criteria play a fundamental role in biomechanical analyses of aortic wall rupture and dissection. Experimental investigations have demonstrated a significant difference of aortic wall strengths in the circumferential and axial directions. Therefore, the isotropic von Mises stress and maximum principal stress, commonly used in computational analysis of the aortic wall, are inadequate for modeling of anisotropic failure properties. In this study, we propose a novel stress-based anisotropic failure criterion with dispersed fiber orientations. In the new failure criterion, the overall failure metric is computed by using angular integration (AI) of failure metrics in all directions. Affine rotations of fiber orientations due to finite deformation are taken into account in an anisotropic hyperelastic constitutive model. To examine fitting capability of the failure criterion, a set of off-axis uniaxial tension tests were performed on aortic tissues of four porcine individuals and 18 human ascending thoracic aortic aneurysm (ATAA) patients. The dispersed fiber failure criterion demonstrates a good fitting capability with the off-axis testing data. Under simulated biaxial stress conditions, the dispersed fiber failure criterion predicts a smaller failure envelope comparing to those predicted by the traditional anisotropic criteria without fiber dispersion, which highlights the potentially important role of fiber dispersion in the failure of the aortic wall. Our results suggest that the deformation-dependent fiber orientations need to be considered when wall strength determined from uniaxial tests are used for in vivo biomechanical analysis. More investigations are needed to determine biaxial failure properties of the aortic wall.

A Dual-VENC Four-Dimensional Flow MRI Framework for Analysis of Subject-Specific Heterogeneous Nonlinear Vessel Deformation

Advancement of subject-specific in silico medicine requires new imaging protocols tailored to specific anatomical features, paired with new constitutive model development based on structure/function relationships. In this study, we develop a new dual-velocity encoding coefficient (VENC) 4D flow MRI protocol that provides unprecedented spatial and temporal resolution of in vivo aortic deformation. All previous dual-VENC 4D flow MRI studies in the literature focus on an isolated segment of the aorta, which fail to capture the full spectrum of aortic heterogeneity that exists along the vessel length. The imaging protocol developed provides high sensitivity to all blood flow velocities throughout the entire cardiac cycle, overcoming the challenge of accurately measuring the highly unsteady nonuniform flow field in the aorta. Cross-sectional area change, volumetric flow rate, and compliance are observed to decrease with distance from the heart, while pulse wave velocity (PWV) is observed to increase. A nonlinear aortic lumen pressure-area relationship is observed throughout the aorta such that a high vessel compliance occurs during diastole, and a low vessel compliance occurs during systole. This suggests that a single value of compliance may not accurately represent vessel behavior during a cardiac cycle in vivo. This high-resolution MRI data provide key information on the spatial variation in nonlinear aortic compliance, which can significantly advance the state-of-the-art of in-silico diagnostic techniques for the human aorta.

Investigating the Mechanical Behavior of Clot Analogues Through Experimental and Computational Analysis

With mechanical thrombectomy emerging as the new standard of care for stroke treatment, clot analogues provide an extremely useful tool in the testing and design of these treatment devices. The aim of this study is to characterise the mechanical behavior of thrombus analogues as a function of composition. Platelet-contracted clot analogues were prepared from blood mixtures of various hematocrits. Mechanical testing was performed whereby clots were subjected to unconfined compression between two rigid plates. Two loading protocols were imposed: cyclic compression for 10 cycles at a constant strain-rate magnitude; stress-relaxation at a constant applied compressive strain. A hyper-viscoelastic constitutive law was identified and calibrated based on the experimental mechanical test data. Scanning electron microscopy (SEM) investigated the clot microstructure at various time-points. Clot analogue composition was found to strongly affect the observed mechanical behavior. The SEM found that the microstructure of the clot analogues was affected by the storage solution and age of the clot. The proposed hyper-viscoelastic constitutive model was found to successfully capture the material test data. The results presented in this study are of key importance to the evaluation and future development mechanical thrombectomy devices and procedures.

Structure, mechanical properties, and modeling of cyclically compressed pulmonary emboli

Pulmonary embolism occurs when blood flow to a part of the lungs is blocked by a venous thrombus that has traveled from the lower limbs. Little is known about the mechanical behavior of emboli under compressive forces from the surrounding musculature and blood pressure. We measured the stress-strain responses of human pulmonary emboli under cyclic compression, and showed that emboli exhibit a hysteretic stress-strain curve. The fibrin fibers and red blood cells (RBCs) are damaged during the compression process, causing irreversible changes in the structure of the emboli. We showed using electron and confocal microscopy that bundling of fibrin fibers occurs due to compression, and damage is accumulated as more cycles are applied. The stress-strain curves depend on embolus structure, such that variations in composition give quantitatively different responses. Emboli with a high fibrin component demonstrate higher normal stress compared to emboli that have a high RBC component. We compared the compression response of emboli to that of whole blood clots containing various volume fractions of RBCs, and found that RBCs rupture at a certain critical stress. We describe the hysteretic response characteristic of foams, using a model of phase transitions in which the compressed foam is segregated into coexisting rarefied and densified phases whose fractions change during compression. Our model takes account of the rupture of RBCs in the compressed emboli and stresses due to fluid flow through their small pores. Our results can help in classifying emboli as rich in fibrin or rich in red blood cells, and can help in understanding what responses to expect when stresses are applied to thrombi in vivo.

In vivo measurement of blood clot mechanics from computational fluid dynamics based on intravital microscopy images

Ischemia which leads to heart attacks and strokes is one of the major causes of death in the world. Whether an occlusion occurs or not depends on the ability of a growing thrombus to resist flow forces exerted on its structure. This manuscript provides the first known in vivo measurement of how much stress a clot can withstand, before yielding to the surrounding blood flow. Namely, Lattice-Boltzmann Method flow simulations are performed based on 3D clot geometries, which are estimated from intravital microscopy images of laser-induced injuries in cremaster microvasculature of live mice. In addition to reporting the blood clot yield stresses, we also show that the thrombus "core" does not experience significant deformation, while its "shell" does. This indicates that the shell is more prone to embolization. Therefore, drugs should be designed to target the shell selectively, while leaving the core intact to minimize excessive bleeding. Finally, we laid down a foundation for a nondimensionalization procedure which unraveled a relationship between clot mechanics and biology. Hence, the proposed framework could ultimately lead to a unified theory of thrombogenesis, capable of explaining all clotting events. Thus, the findings presented herein will be beneficial to the understanding and treatment of heart attacks, strokes and hemophilia.

Bayesian Inference-Based Estimation of Normal Aortic, Aneurysmal and Atherosclerotic Tissue Mechanical Properties: From Material Testing, Modeling and Histology

OBJECTIVE

Mechanical properties of healthy, aneurysmal, and atherosclerotic arterial tissues are essential for assessing the risk of lesion development and rupture. Strain energy density function (SEDF) has been widely used to describe these properties, where material constants of the SEDF are traditionally determined using the ordinary least square (OLS) method. However, the material constants derived using OLS are usually dependent on initial guesses.

METHODS

To avoid such dependencies, Bayesian inference-based estimation was used to fit experimental stress-stretch curves of 312 tissue strips from 8 normal aortas, 19 aortic aneurysms, and 21 carotid atherosclerotic plaques to determine the constants, C₁, D₁, and D₂ of the modified Mooney-Rivlin SEDF.

RESULTS

Compared with OLS, material constants varied much less with prior in the Bayesian inference-based estimation. Moreover, fitted material constants differed amongst distinct tissue types. Atherosclerotic tissues associated with the biggest D₂, an indicator of the rate of increase in stress during stretching, followed by aneurysmal tissues and those from normal aortas. Histological analyses showed that C₁ and D₂ were associated with elastin content and details of the collagen configuration, specifically, waviness and dispersion, in the structure.

CONCLUSION

Bayesian inference-based estimation robustly determines material constants in the modified Mooney-Rivlin SEDF and these constants can reflect the inherent physiological and pathological features of the tissue structure.

SIGNIFICANCE

This study suggested a robust procedure to determine the material constants in SEDF and demonstrated that the obtained constants can be used to characterize tissues from different types of lesions, while associating with their inherent microstructures.

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.

Resuscitative endovascular balloon occlusion of the aorta: rupture risk and implications for blind inflation

Background

Resuscitative endovascular balloon occlusion of the aorta (REBOA) is a torso hemorrhage control technique. To expedite deployment, inflation is frequently performed as a blind technique with minimal imaging, which carries a theoretical risk of aortic injury. The objective of this study was to examine the relationship between balloon inflation, deformation and the risk of aortic rupture.

Methods

Compliant balloon catheters were incrementally inflated in segments of cadaveric swine aorta. Serial longitudinal and circumferential measurements were recorded, along with the incidence of aortic rupture.

Results

Fourteen cadaveric swine aorta segments were tested with mean (±SD) baseline aortic diameter (mm) of 14.2±3.4. Rupture occurred in three aortas. The mean baseline diameters (mm) of the aortic segments that were ruptured were significantly smaller than those that did not rupture (8.9±1.2 vs 15.6±1.9; P<0.001). The maximal circumferential stretch ratios were significantly higher in the aorta segments that ruptured compared with those that did not (1.9±0.1 vs 1.5±0.1; P<0.001). The maximal amount of balloon longitudinal deformation was 80 mm (116% longer than the intended working length).

Conclusions

Inflation of aortic balloon catheters carries an inherent risk of aortic injury, which may be minimized through an understanding of the intrinsic characteristics of the aorta and compliant balloons. Smaller diameter aortic segments undergoing overinflation, particularly beyond a circumferential stretch ratio of 1.8, are at risk of aortic rupture.

Level of evidence

Level II.

Microstructure and mechanics of healthy and aneurysmatic abdominal aortas: experimental analysis and modelling

Soft biological tissues such as aortic walls can be viewed as fibrous composites assembled by a ground matrix and embedded families of collagen fibres. Changes in the structural components of aortic walls such as the ground matrix and the embedded families of collagen fibres have been shown to play a significant role in the pathogenesis of aortic degeneration. Hence, there is a need to develop a deeper understanding of the microstructure and the related mechanics of aortic walls. In this study, tissue samples from 17 human abdominal aortas (AA) and from 11 abdominal aortic aneurysms (AAA) are systematically analysed and compared with respect to their structural and mechanical differences. The collagen microstructure is examined by analysing data from second-harmonic generation imaging after optical clearing. Samples from the intact AA wall, their individual layers and the AAA wall are mechanically investigated using biaxial stretching tests. A bivariate von Mises distribution was used to represent the continuous fibre dispersion throughout the entire thickness, and to provide two independent dispersion parameters to be used in a recently proposed material model. Remarkable differences were found between healthy and diseased tissues. The out-of-plane dispersion was significantly higher in AAA when compared with AA tissues, and with the exception of one AAA sample, the characteristic wall structure, as visible in healthy AAs with three distinct layers, could not be identified in AAA samples. The collagen fibres in the abluminal layer of AAAs lost their waviness and exhibited rather straight and thick struts of collagen. A novel set of three structural and three material parameters is provided. With the structural parameters fixed, the material model was fitted to the mechanical experimental data, giving a very satisfying fit although there are only three material parameters involved. The results highlight the need to incorporate the structural differences into finite-element simulations as otherwise simulations of AAA tissues might not be good predictors for the actual in vivo stress state.

STRESS ANALYSIS OF INTERNAL CAROTID ARTERY WITH LOW STENOSIS LEVEL: THE EFFECT OF MATERIAL MODEL AND PLAQUE GEOMETRY

Stress concentration in carotid stenosis has been proven to assist plaque morphology in disease diagnosis and vulnerability. This work focuses on numerical analysis of stress and strain distribution in the cross-section of internal carotid artery using a 2D structure-only method. The influence of four different idealized plaque geometries (circle, ellipse, oval and wedge) is investigated. Numerical simulations are implemented utilizing linear elastic model along with four hyperelastic constitutive laws named neo-Hookean, Ogden, Yeoh and Mooney–Rivlin. Each case is compared to the real geometry. Results show significant strength of oval and wedged geometries in predicting stress and strain values. Our results emphasize that Yeoh and Ogden hyperelastic materials are more reliable in stress prediction with errors less than 3%. The same concept is observed in locating critical stresses where oval and wedged plaque geometries are the most accurate models. Similar results are observed in predicting maximum principal elastic strain with errors less than 1%. However, the strain distribution in idealized plaque models showed a considerable difference in comparison with real geometry.

Influence of overlapping pattern of multiple overlapping uncovered stents on the local mechanical environment: A patient-specific parameter study

BACKGROUND

Multiple overlapping uncovered stents (MOUS) system has shown potentials in managing complex aortic aneurysms with side branches involvement. It promotes the development of thrombus by modulating local flow pattern that reduces the wall tension, while maintaining patency of side branches. However the modulation of local hemodynamic parameters depends on various factors that have not been assessed comprehensively.

METHODS

Aneurysm 3D geometry was reconstructed based on CT images. One-way fluid-structure interaction analysis was performed to quantify structural stress concentration in the wall, and changes of blood velocity, wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT) and pressure in the sac due to the stent deployment.

RESULTS

High structural stress concentration due to stent deployment was found in the landing zone and it increased linearly with the number of stents deployed. The wall tension in the sac was unaffected by the stent deployment. Stress within the wall was insensitive to the different overlapping pattern. After one stent was deployed, the mean flow velocity in the sac reduced by 36.4%. The deployment of the 2nd stent further reduced the mean sac velocity by 10%. WSS decreased while both OSI and RRT increased after stent deployment, however pressure in the sac remained nearly unchanged. Except for the cases with complete stents struts alignment, different overlapping pattern had little effect on flow parameters.

CONCLUSIONS

Mechanical parameters modulated by the MOUS are insensitive to different overlapping pattern suggesting that endovascular procedure can be performed with less attention to the overlapping pattern.

Non-contrast 3D black blood MRI for abdominal aortic aneurysm surveillance: comparison with CT angiography

OBJECTIVES

Management of abdominal aortic aneurysms (AAAs) is based on diameter. CT angiography (CTA) is commonly used, but requires radiation and iodinated contrast. Non-contrast MRI is an appealing alternative that may allow better characterization of intraluminal thrombus (ILT). This study aims to 1) validate non-contrast MRI for measuring AAA diameter, and 2) to assess ILT with CTA and MRI.

METHOD

28 patients with AAAs (diameter 50.7 ± 12.3 mm) underwent CTA and non-contrast MRI. MRI was acquired at 3 T using 1) a conventional 3D gradient echo (GRE) sequence and 2) a 3D T1-weighted black blood fast-spin-echo sequence. Two radiologists independently measured the AAA diameter. The ratio of signal of ILT and adjacent psoas muscle (ILTr = signalILT/signalMuscle) was quantified.

RESULTS

Strong agreement between CTA and non-contrast MRI was shown for AAA diameter (intra-class coefficient > 0.99). Both approaches had excellent inter-observer reproducibility (ICC > 0.99). ILT appeared homogenous on CTA, whereas MRI revealed compositional variations. Patients with AAAs ≥5.5 cm and <5.5 cm had a variety of distributions of old/fresh ILT types.

CONCLUSIONS

Non-contrast 3D black blood MRI provides accurate and reproducible AAA diameter measurements as validated by CTA. It also provides unique information about ILT composition, which may be linked with elevated risk for disease progression.

KEY POINTS

• Non-contrast MRI is an appealing alternative to CTA for AAA management. • Non-contrast MRI can accurately measure AAA diameters compared to CTA. • MRI affords unique characterization of intraluminal thrombus composition.

The influence of constitutive law choice used to characterise atherosclerotic tissue material properties on computing stress values in human carotid plaques

Calculating high stress concentration within carotid atherosclerotic plaques has been shown to be complementary to anatomical features in assessing vulnerability. Reliability of stress calculation may depend on the constitutive laws/strain energy density functions (SEDFs) used to characterize tissue material properties. Different SEDFs, including neo-Hookean, one-/two-term Ogden, Yeoh, 5-parameter Mooney–Rivlin, Demiray and modified Mooney–Rivlin, have been used to describe atherosclerotic tissue behavior. However, the capacity of SEDFs to fit experimental data and the difference in the stress calculation remains unexplored. In this study, seven SEDFs were used to fit the stress–stretch data points of media, fibrous cap, lipid and intraplaque hemorrhage/thrombus obtained from 21 human carotid plaques. Semi-analytic solution, 2D structure-only and 3D fully coupled fluid-structure interaction (FSI) analyses were used to quantify stress using different SEDFs and the related material stability examined. Results show that, except for neo-Hookean, all other six SEDFs fitted the experimental points well, with vessel stress distribution in the circumferential and radial directions being similar. 2D structural-only analysis was successful for all seven SEDFs, but 3D FSI were only possible with neo-Hookean, Demiray and modified Mooney–Rivlin models. Stresses calculated using Demiray and modified Mooney–Rivlin models were nearly identical. Further analyses indicated that the energy contours of one-/two-term Ogden and 5-parameter Mooney–Rivlin models were not strictly convex and the material stability indictors under homogeneous deformations were not always positive. In conclusion, considering the capacity in characterizing material properties and stabilities, Demiray and modified Mooney–Rivlin SEDF appear practical choices for mechanical analyses to predict the critical mechanical conditions within carotid atherosclerotic plaques.