Experimental characterization of rupture in human aortic aneurysms using a full-field measurement technique

Experimental Protocols to Test Aortic Soft Tissues: A Systematic Review

Experimental protocols are fundamental for quantifying the mechanical behaviour of soft tissue. These data are crucial for advancing the understanding of soft tissue mechanics, developing and calibrating constitutive models, and informing the development of more accurate and predictive computational simulations and artificial intelligence tools. This paper offers a comprehensive review of experimental tests conducted on soft aortic tissues, employing the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) methodology, based on the Scopus, Web of Science, IEEE, Google Scholar and PubMed databases. This study includes a detailed overview of the test method protocols, providing insights into practical methodologies, specimen preparation and full-field measurements. The review also briefly discusses the post-processing methods applied to extract material parameters from experimental data. In particular, the results are analysed and discussed providing representative domains of stress-strain curves for both uniaxial and biaxial tests on human aortic tissue.

On the Importance of Including Cohesive Zone Models in Modelling Mixed-Mode Aneurysm Rupture

INTRODUCTION

The precise mechanism of rupture in abdominal aortic aneurysms (AAAs) has not yet been uncovered. The phenomenological failure criterion of the coefficient of proportionality between von Mises stress and tissue strength does not account for any mechanistic foundation of tissue fracture. Experimental studies have shown that arterial failure is a stepwise process of fibrous delamination (mode II) and kinking (mode I) between layers. Such a mechanism has not previously been considered for AAA rupture.

METHODS

In the current study we consider both von Mises stress in the wall, in addition to interlayer tractions and delamination using cohesive zone models. Firstly, we present a parametric investigation of the influence of a range of AAA anatomical features on the likelihood of elevated interlayer traction and delamination.

RESULTS

We observe in several cases that the location of peak von Mises stress and tangential traction coincide. Our simulations also reveal however, that peak von Mises and intramural tractions are not coincident for aneurysms with Length/Radius less than 2 (short high-curvature aneurysms) and for aneurysms with symmetric intraluminal thrombus (ILT). For an aneurysm with (L/R = 2.0), the peak σ vm moves slightly towards the origin while the peak T t is near the peak bulge with a separation distance of ~ 17 mm. Additionally, we present three patient-specific AAA models derived directly from CT scans, which also illustrate that the location of von Mises stress does not correlate with the point of interlayer delamination.

CONCLUSION

This study suggests that incorporating cohesive zone models into clinical based FE analyses may capture a greater proportion of ruptures in-silico.

Mechanics of ascending aortic aneurysms based on a modulus of elasticity dependent on aneurysm diameter and pressure

The mechanical stresses and strains are examined, in ascending thoracic aortic aneurysm (aTAA) models, in a patient-specific aTAA as well as in healthy thoracic aortic models, via Finite Element Analysis. The aneurysms are assumed spherical, 1.5 mm thick, with diameters between 47 mm and 80 mm, eccentrically positioned. The geometry and wall thickness distribution of the aorta along its length are based on open literature data for an average patient age of 66.25 years, accounting for the Body Surface Area (BSA) parameter. The vessel wall material is assumed isotropic and incompressible, with its Young's modulus varying with the aneurysm diameter and the applied intraluminal pressure (120 mmHg to 240 mmHg). In the aTAAs, peak stresses were found to increase nonlinearly with aneurysm diameter (for a given pressure) tending to reach a plateau, appearing at the proximal area of the aneurysm, whereas lower stresses were found at its distal part and even smaller at the aneurysm maximum diameter. Regarding the patient-specific aTAA model, the peak stresses appeared at the distal part of the aneurysm where a tear of the intima layer was detected during surgical intervention. Peak strains exhibited for each pressure a maximum at a certain aneurysm diameter beyond which they dropped so that essentially the vessel wall's distensibility was thus reduced. Examining more than 100 geometry cases and employing a failure stress criterion, the rupture diameter thresholds were estimated to be 65, 52.5, 50 and 47.5 mm for a pressure of 120, 160, 200 and 240 mmHg respectively.

Intervertebral Disc Elastography to Relate Shear Modulus and Relaxometry in Compression and Bending

Intervertebral disc degeneration is the most recognized cause of low back pain, characterized by the decline of tissue structure and mechanics. Image-based mechanical parameters (e.g., strain, stiffness) may provide an ideal assessment of disc function that is lost with degeneration but unfortunately remains underdeveloped. Moreover, it is unknown whether strain or stiffness of the disc may be predicted by MRI relaxometry (e.g. T1 or T2), an increasingly accepted quantitative measure of disc structure. In this study, we quantified T1 and T2 relaxation times and in-plane strains using displacement-encoded MRI within the disc under physiological levels of compression and bending. We then estimated shear modulus in orthogonal image planes and compared these values to relaxation times and strains within regions of the disc. Intratissue strain depended on the loading mode, and shear modulus in the nucleus pulposus was typically an order of magnitude lower than the annulus fibrosis, except in bending, where the apparent stiffness depended on the loading. Relative shear moduli estimated from strain data derived under compression generally did not correspond with those from bending experiments, with no correlations in the sagittal plane and only 4 of 15 regions correlated in the coronal plane, suggesting that future inverse models should incorporate multiple loading conditions. Strain imaging and strain-based estimation of material properties may serve as imaging biomarkers to distinguish healthy and diseased discs. Additionally, image-based elastography and relaxometry may be viewed as complementary measures of disc structure and function to assess degeneration in longitudinal studies.

Estimation of anisotropic properties of CMR patient-specific left ventricle using the virtual field method

Left ventricle (LV) myocardial dysfunction has been recently investigated using the estimation of isotropic myocardial stiffness from magnetic resonance imaging (MRI). However, Myocardium is known to have a 3D complex geometry with anisotropic stiffness. The assessment of the anisotropy properties characterizes structural changes in myocardium as a consequence of heart failure (HF). From image data, the virtual field method (VFM) can determine material stiffness in a non-invasive manner. In the present work, the objective is to compare two inverse identification methods, given the isotropic and anisotropic models in the characterization of properties of myocardium in acute lymphoblastic leukemia (ALL) survivors using VFM and MRI. Two types of VFM approach are presented. Using the first, the virtual displacements (VFs) allow whole-field LV to be imposed into VFM formulation and caused to directly estimate two independent parameters from isotropic constitutive relation. With the second, anisotropic parameters are estimated using piece-wise (Finite element-based) VFM. The resulting values showed significant differences between the subjects in comparative study of leukemia survivors, and variance in estimated parameters by two different VFM approach. This approach would be an efficient tool to characterize early cardiac dysfunction. This work elucidates the benefits and shortcomings of using VFM to determine anisotropic parameters of LV myocardium in linear elastic and of using the FEM application to generate meshes of patient-specific LVs from MRI images.

Experimental Study of the Propagation Process of Dissection Using an Aortic Silicone Phantom

BACKGROUND

The mortality of acute aortic dissection (AD) can reach 65~70%. However, it is challenging to follow the progress of AD formation. The purpose of this work was to observe the process of dissection development using a novel tear-embedded silicone phantom.

METHODS

Silicone phantoms were fabricated by embedding a torn area and primary tear feature on the inner layer. CT scanning and laser lightening were conducted to observe the variations in thickness and volume of the true lumen (TL) and false lumen (FL) during development.

RESULTS

The model with a larger interlayer adhesion damage required a lower pressure to trigger the development of dissection. At the initiation stage of dissection, the volume of TL increased by 25.5%, accompanied by a 19.5% enlargement of tear size. The force analysis based on the change of tear size verified the deduction of the process of interlaminar separation from the earlier studies.

CONCLUSIONS

The primary tear and the weakening adhesion of the vessel layers are key factors in AD development, suggesting that some forms of primary damage to the arterial wall, in particular, the lumen morphology of vessels with straight inner lumen, should be considered as early risk predictors of AD.

Material Identification on Thin Shells Using the Virtual Fields Method, Demonstrated on the Human Eardrum

Characterization of material parameters from experimental data remains challenging, especially on biological structures. One of such techniques allowing for the inverse determination of material parameters from measurement data is the virtual fields method (VFM). However, application of the VFM on general structures of complicated shape has not yet been extensively investigated. In this paper, we extend the framework of the VFM method to thin curved solids in three-dimensional, commonly denoted shells. Our method is then used to estimate the Young's modulus and hysteretic damping of the human eardrum. By utilizing Kirchhoff plate theory, we assume that the behavior of the shell varies linearly through the thickness. The total strain of the shell can then be separated in a bending and membrane strain. This in turn allowed for an application of the VFM based only on data of the outer surface of the shell. We validated our method on simulated and experimental data of a human eardrum made to vibrate at certain frequencies. It was shown that the identified material properties were accurately determined based only on data from the outer surface and are in agreement with literature. Additionally, we observed that neither the bending nor the membrane strain in an human eardrum can be neglected and both contribute significantly to the total strain found experimentally.

Multi-view digital image correlation systems for in vitro testing of arteries from mice to humans

BACKGROUND

Digital image correlation (DIC) methods are increasingly used for non-contact optical assessment of geometry and deformation in soft tissue biomechanics, thus providing the full-field strain estimates needed for robust inverse material characterization. Despite the well-known flexibility and ease of use of DIC, issues related to spatial resolution and depth-of-field remain challenging in studies of quasi-cylindrical biological samples such as arteries.

OBJECTIVE

After demonstrating that standard surrounding multi-view DIC systems are inappropriate for such usage, we submit that both the optical setup and the data analysis need to be specifically designed with respect to the size of the arterial sample of interest. Accordingly, we propose novel and optimized DIC systems for two distinct ranges of arterial diameters: less than 2.5 mm (murine arteries) and greater than 10 mm (human arteries).

METHODS

We designed, set up, and validated a four-camera panoramic-DIC system for testing murine arteries and a multi-biprism DIC system for testing human arteries. Both systems enable dynamic 360-deg measurements with refraction correction over the entire surface of submerged samples in their native geometries.

RESULTS

Illustrative results for 3D shape and full-surface deformation fields were obtained for a mouse infrarenal aorta and a latex cylinder of size similar to the human infrarenal aorta.

CONCLUSION

Results demonstrated the feasibility and accuracy of both proposed methods in providing quantitative information on the regional behavior of arterial samples tested in vitro under physiologically relevant loading.

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.

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

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

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.

Investigation of Pathophysiological Aspects of Aortic Growth, Remodeling, and Failure Using a Discrete-Fiber Microstructural Model

Aortic aneurysms are inherently unpredictable. One can never be sure whether any given aneurysm may rupture or dissect. Clinically, the criteria for surgical intervention are based on size and growth rate, but it remains difficult to identify a high-risk aneurysm, which may require intervention before the cutoff criteria, versus an aneurysm than can be treated safely by more conservative measures. In this work, we created a computational microstructural model of a medial lamellar unit (MLU) incorporating (1) growth and remodeling laws applied directly to discrete, individual fibers, (2) separate but interacting fiber networks for collagen, elastin, and smooth muscle, (3) active and passive smooth-muscle cell mechanics, and (4) failure mechanics for all three fiber types. The MLU model was then used to study different pathologies and microstructural anomalies that may play a role in vascular growth and failure. Our model recapitulated many aspects of arterial remodeling under hypertension with no underlying genetic syndrome including remodeling dynamics, tissue mechanics, and failure. Syndromic effects (smooth muscle cell (SMC) dysfunction or elastin fragmentation) drastically changed the simulated remodeling process, tissue behavior, and tissue strength. Different underlying pathologies were able to produce similarly dilatated vessels with different failure properties, providing a partial explanation for the imperfect nature of aneurysm size as a predictor of outcome.

Null strain analysis of submerged aneurysm analogues using a novel 3D stereomicroscopy device

To measure the inhomogeneous 3D-strain fields present during inflation-extension testing of physiologically submerged micro-aneurysms, a Stereo Digital Image Correlation (StereoDIC) microscopy system is developed that revolves 15° stereo-angle cameras around a centrally-mounted target. Calibration is performed using submerged dot patterns and system accuracy verified using strain and deformation analyses for rigid body motions of speckle-patterned, micro-aneurysmal surrogates. In terms of the Green-Lagrange strain tensor and the 3D displacement fields, the results are stable even after 120 minutes, with maxima in both strain bias and strain standard deviation less than 2E-03 for all components, and micron-level displacement standard deviation.

Ascending aorta mechanics in bicuspid aortopathy: controversy or fact?

Bicuspid aortic valve is the most common congenital cardiovascular defect, often associated with proximal aortic dilatation, and the ideal management strategy is debated. The inconsistency in previous and present guideline recommendations emphasizes the insufficiency of the maximal diameter as the sole criterion for prophylactic repair. Our ability to guide clinical decisions may improve through an understanding of the mechanical properties of ascending thoracic aortic aneurysms in bicuspid compared to tricuspid aortic valve patients and non-aneurysmal aortas, because dissection and rupture are aortic wall mechanical failures. Such an understanding of the mechanical properties has been attempted by several authors, and this article addresses whether there is a controversy in the accumulated knowledge. The available mechanical studies are briefly reviewed, discussing factors such as age, sex, and the region of mechanical examination that may be responsible for the lack of unanimity in the reported findings. The rationale for acquiring layer-specific properties is presented along with the main results from our recent study. No mechanical vulnerability of ascending thoracic aortic aneurysms was evidenced in bicuspid aortic valve patients, corroborating present conservative guidelines concerning the management of bicuspid aortopathy. Weakening and additional vulnerability was evidenced in aged patients and those with coexisting valve pathology, aortic root dilatation, hypertension, and hyperlipidemia. Discussion of these results from age- and sex-matched subjects, accounting for the region- and layer-specific aortic heterogeneity, in relation to intact wall results and histologic confirmation, helps to reconcile previous findings and affords a universal interpretation of ascending aorta mechanics in bicuspid aortopathy.

Thermo-Mechanical Behaviour of Human Nasal Cartilage

The aim of this study was to undergo a comprehensive analysis of the thermo-mechanical properties of nasal cartilages for the future design of a composite polymeric material to be used in human nose reconstruction surgery. A thermal and dynamic mechanical analysis (DMA) in tension and compression modes within the ranges 1 to 20 Hz and 30 °C to 250 °C was performed on human nasal cartilage. Differential scanning calorimetry (DSC), as well as characterization of the nasal septum (NS), upper lateral cartilages (ULC), and lower lateral cartilages (LLC) reveals the different nature of the binding water inside the studied specimens. Three peaks at 60–80 °C, 100–130 °C, and 200 °C were attributed to melting of the crystalline region of collagen matrix, water evaporation, and the strongly bound non-interstitial water in the cartilage and composite specimens, respectively. Thermogravimetric analysis (TGA) showed that the degradation of cartilage, composite, and subcutaneous tissue of the NS, ULC, and LLC take place in three thermal events (~37 °C, ~189 °C, and ~290 °C) showing that cartilage releases more water and more rapidly than the subcutaneous tissue. The water content of nasal cartilage was estimated to be 42 wt %. The results of the DMA analyses demonstrated that tensile mode is ruled by flow-independent behaviour produced by the time-dependent deformability of the solid cartilage matrix that is strongly frequency-dependent, showing an unstable crystalline region between 80–180 °C, an amorphous region at around 120 °C, and a clear glass transition point at 200 °C (780 kJ/mol). Instead, the unconfined compressive mode is clearly ruled by a flow-dependent process caused by the frictional force of the interstitial fluid that flows within the cartilage matrix resulting in higher stiffness (from 12 MPa at 1 Hz to 16 MPa at 20 Hz in storage modulus). The outcomes of this study will support the development of an artificial material to mimic the thermo-mechanical behaviour of the natural cartilage of the human nose.

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.

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.

Does the Knowledge of the Local Thickness of Human Ascending Thoracic Aneurysm Walls Improve Their Mechanical Analysis?

Ascending thoracic aortic aneurysm (ATAA) ruptures are life threatening phenomena which occur in local weaker regions of the diseased aortic wall. As ATAAs are evolving pathologies, their growth represents a significant local remodeling and degradation of the microstructural architecture and thus their mechanical properties. To address the need for deeper study of ATAAs and their failure, it is required to analyze the mechanical behavior at the sub-millimeter scale by making use of accurate geometrical and kinematical measurements during their deformation. For this purpose, we propose a novel methodology that combined an accurate tool for thickness distribution measurement of the arterial wall, digital image correlation to assess local strain fields and bulge inflation to characterize the physiological and failure response of flat unruptured human ATAA specimens. The analysis of the heterogeneity of the local thickness and local physiological stress and strain was carried out for each investigated subject. At the subject level, our results state the presence of a non-consistent relationship between the local wall thickness and the local physiological strain field and high heterogeneity of the variables. At the inter-subject level, thicknesses were studied in relation to physiological strain and stress and load at rupture. The rupture pressure was correlated with neither the average thickness nor the lowest thickness of the specimens. Our results confirm that intrinsic material strength (hence structure) differs a lot from a subject to another and even within the same subject.

Inflation and rupture of vaginal tissue

Around 80% of women experience vaginal tears during labour when the diameter of the vagina must increase to allow the passage of a full-term baby. Current techniques for evaluating vaginal tears are qualitative and often lead to an incorrect diagnosis and inadequate treatment, severely compromising the quality of life of women. In order to characterize the failure properties of the vaginal tissue, whole vaginal tracts from rats (n = 18) were subjected to free-extension inflation tests until rupture using a custom-built experimental set-up. The resulting deformations were measured using the digital image correlation technique. Overall, the strain and changes in curvature in the hoop direction were significantly larger relative to the axial direction. At a failure pressure of 110 ± 23 kPa (mean ± s.d.), the hoop and axial stresses were computed to be 970 ± 340 kPa and 490 ± 170 kPa, respectively. Moreover, at such pressure, the hoop and axial strains were found to be 12.8 ± 4.4 % and 6.4 ± 3.7 % , respectively. Rupture of the vaginal specimens always occurred in the hoop direction by tearing along the axial direction. This knowledge about the rupture properties of the vaginal tissue will be crucial for the development of clinical approaches for preventing and mitigating vaginal tearing and the associated short- and long-term traumatic conditions.

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.

Material characterization of cardiovascular biomaterials using an inverse finite-element method and an explicit solver

The ability to accurately model soft tissue behavior, such as that of heart valve tissue, is essential for developing reliable numerical simulations and determining patient-specific care options. Although several material models can predict soft tissue behavior, complications may arise when these models are implemented into finite element (FE) programs, due to the addition of an arbitrary penalty parameter for numerically enforcing material incompressibility. Herein, an inverse methodology was developed in MATLAB to use previously published stress-strain data from experimental planar equibiaxial testing of five biomaterials used in heart valve cusp replacements, in conjunction with commercial explicit FE solver LS-DYNA, to optimize the material parameters and the penalty parameter for an anisotropic hyperelastic strain energy function. A two-parameter optimization involving the scaling constant of the strain energy function and the penalty parameter proved sufficient to produce acceptable material responses when compared with experimental behaviors under the same testing conditions, as long as analytically derived material constants were available for the other non-optimized parameters and the actual tissue thickness was not much less than 1 mm. Variations in the penalty parameter had a direct effect on the accuracy of the simulated responses, with a practical range determined to be 5×10⁸-9×10⁸ times the scaling constant of the strain energy function.

A Uniaxial Testing Approach for Consistent Failure in Vascular Tissues

Although uniaxial tensile testing is commonly used to evaluate failure properties of vascular tissue, there is no established protocol for specimen shape or gripping method. Large percentages of specimens are reported to fail near the clamp and can potentially confound the studies, or, if discarded will result in sample waste. The objective of this study is to identify sample geometry and clamping conditions that can achieve consistent failure in the midregion of small arterial specimens, even for vessels from older individuals. Failure location was assessed in 17 dogbone specimens from human cerebral and sheep carotid arteries using soft inserts. For comparison with commonly used protocols, an additional 22 rectangular samples were tested using either sandpaper or foam tape inserts. Midsample failure was achieved in 94% of the dogbone specimens, while only 14% of the rectangular samples failed in the midregion, the other 86% failing close to the clamps. Additionally, we found midregion failure was more likely to be abrupt, caused by cracking or necking. In contrast, clamp failure was more likely to be gradual and included a delamination mode not seen in midregion failure. Hence, this work provides an approach that can be used to obtain consistent midspecimen failure, avoiding confounding clamp-related artifacts. Furthermore, with consistent midregion failure, studies can be designed to image the failure process in small vascular samples providing valuable quantitative information about changes to collagen and elastin structure during the failure process.

A review of techniques for visualising soft tissue microstructure deformation and quantifying strain Ex Vivo

Many biological tissues have a complex hierarchical structure allowing them to function under demanding physiological loading conditions. Structural changes caused by ageing or disease can lead to loss of mechanical function. Therefore, it is necessary to characterise tissue structure to understand normal tissue function and the progression of disease. Ideally intact native tissues should be imaged in 3D and under physiological loading conditions. The current published in situ imaging methodologies demonstrate a compromise between imaging limitations and maintaining the samples native mechanical function. This review gives an overview of in situ imaging techniques used to visualise microstructural deformation of soft tissue, including three case studies of different tissues (tendon, intervertebral disc and artery). Some of the imaging techniques restricted analysis to observational mechanics or discrete strain measurement from invasive markers. Full-field local surface strain measurement has been achieved using digital image correlation. Volumetric strain fields have successfully been quantified from in situ X-ray microtomography (micro-CT) studies of bone using digital volume correlation but not in soft tissue due to low X-ray transmission contrast. With the latest developments in micro-CT showing in-line phase contrast capability to resolve native soft tissue microstructure, there is potential for future soft tissue mechanics research where 3D local strain can be quantified. These methods will provide information on the local 3D micromechanical environment experienced by cells in healthy, aged and diseased tissues. It is hoped that future applications of in situ imaging techniques will impact positively on the design and testing of potential tissue replacements or regenerative therapies. LAY DESCRIPTION: The soft tissues in our bodies, such as tendons, intervertebral discs and arteries, have evolved to have complicated structures which deform and bear load during normal function. Small changes in these structures can occur with age and disease which then leads to loss of function. Therefore, it is important to image tissue microstructure in 3D and under functional conditions. This paper gives an overview of imaging techniques used to record the deformation of soft tissue microstructures. Commonly there are compromises between obtaining the best imaging result and retaining the samples native structure and function. For example, invasive markers and dissecting samples damages the tissues natural structure, and staining or clearing (making the tissue more transparent) can distort tissue structure. Structural deformation has been quantified from 2D imaging techniques (digital image correlation) to create surface strain maps which help identify local tissue mechanics. When extended to 3D (digital volume correlation), deformation measurement has been limited to bone samples using X-ray micro-CT. Recently it has been possible to image the 3D structure of soft tissue using X-ray micro-CT meaning that there is potential for internal soft tissue mechanics to be mapped in 3D. Future application of micro-CT and digital volume correlation will be important for soft tissue mechanics studies particularly to understand normal function, progression of disease and in the design of tissue replacements.

Numerical aspects of anisotropic failure in soft biological tissues favor energy-based criteria: A rate-dependent anisotropic crack phase-field model

A deeper understanding to predict fracture in soft biological tissues is of crucial importance to better guide and improve medical monitoring, planning of surgical interventions and risk assessment of diseases such as aortic dissection, aneurysms, atherosclerosis and tears in tendons and ligaments. In our previous contribution (Gültekin et al., 2016) we have addressed the rupture of aortic tissue by applying a holistic geometrical approach to fracture, namely the crack phase-field approach emanating from variational fracture mechanics and gradient damage theories. In the present study, the crack phase-field model is extended to capture anisotropic fracture using an anisotropic volume-specific crack surface function. In addition, the model is equipped with a rate-dependent formulation of the phase-field evolution. The continuum framework captures anisotropy, is thermodynamically consistent and based on finite strains. The resulting Euler-Lagrange equations are solved by an operator-splitting algorithm on the temporal side which is ensued by a Galerkin-type weak formulation on the spatial side. On the constitutive level, an invariant-based anisotropic material model accommodates the nonlinear elastic response of both the ground matrix and the collagenous components. Subsequently, the basis of extant anisotropic failure criteria are presented with an emphasis on energy-based, Tsai-Wu, Hill, and principal stress criteria. The predictions of the various failure criteria on the crack initiation, and the related crack propagation are studied using representative numerical examples, i.e. a homogeneous problem subjected to uniaxial and planar biaxial deformations is established to demonstrate the corresponding failure surfaces whereas uniaxial extension and peel tests of an anisotropic (hypothetical) tissue deal with the crack propagation with reference to the mentioned failure criteria. Results favor the energy-based criterion as a better candidate to reflect a stable and physically meaningful crack growth, particularly in complex three-dimensional geometries with a highly anisotropic texture at finite strains.

Computational predictions of damage propagation preceding dissection of ascending thoracic aortic aneurysms

Dissections of ascending thoracic aortic aneurysms (ATAAs) cause significant morbidity and mortality worldwide. They occur when a tear in the intima-media of the aorta permits the penetration of the blood and the subsequent delamination and separation of the wall in 2 layers, forming a false channel. To predict computationally the risk of tear formation, stress analyses should be performed layer-specifically and they should consider internal or residual stresses that exist in the tissue. In the present paper, we propose a novel layer-specific damage model based on the constrained mixture theory, which intrinsically takes into account these internal stresses and can predict appropriately the tear formation. The model is implemented in finite-element commercial software Abaqus coupled with user material subroutine. Its capability is tested by applying it to the simulation of different exemplary situations, going from in vitro bulge inflation experiments on aortic samples to in vivo overpressurizing of patient-specific ATAAs. The simulations reveal that damage correctly starts from the intimal layer (luminal side) and propagates across the media as a tear but never hits the adventitia. This scenario is typically the first stage of development of an acute dissection, which is predicted for pressures of about 2.5 times the diastolic pressure by the model after calibrating the parameters against experimental data performed on collected ATAA samples. Further validations on a larger cohort of patients should hopefully confirm the potential of the model in predicting patient-specific damage evolution and possible risk of dissection during aneurysm growth for clinical applications.

Comparison of in vivo vs. ex situ obtained material properties of sheep common carotid artery

Patient-specific biomechanical modelling can improve preoperative surgical planning. This requires patient-specific geometry as well as patient-specific material properties as input. The latter are, however, still quite challenging to estimate in vivo. This study focuses on the estimation of the mechanical properties of the arterial wall. Firstly, in vivo pressure, diameter and thickness of the arterial wall were acquired for sheep common carotid arteries. Next, the animals were sacrificed and the tissue was stored for mechanical testing. Planar biaxial tests were performed to obtain experimental stress-stretch curves. Finally, parameters for the hyperelastic Mooney-Rivlin and Gasser-Ogden-Holzapfel (GOH) material model were estimated based on the in vivo obtained pressure-diameter data as well as on the ex situ experimental stress-stretch curves. Both material models were able to capture the in vivo behaviour of the tissue. However, in the ex situ case only the GOH model provided satisfactory results. When comparing different fitting approaches, in vivo vs. ex situ, each of them showed its own advantages and disadvantages. The in vivo approach estimates the properties of the tissue in its physiological state while the ex situ approach allows to apply different loadings to properly capture the anisotropy of the tissue. Both of them could be further enhanced by improving the estimation of the stress-free state, i.e. by adding residual circumferential stresses in vivo and by accounting for the flattening effect of the tested samples ex vivo. • Competing interests: none declared • Word count: 4716.

Nondestructive mechanical characterization of developing biological tissues using inflation testing

One of the hallmarks of biological soft tissues is their capacity to grow and remodel in response to changes in their environment. Although it is well-accepted that these processes occur at least partly to maintain a mechanical homeostasis, it remains unclear which mechanical constituent(s) determine(s) mechanical homeostasis. In the current study a nondestructive mechanical test and a two-step inverse analysis method were developed and validated to nondestructively estimate the mechanical properties of biological tissue during tissue culture. Nondestructive mechanical testing was achieved by performing an inflation test on tissues that were cultured inside a bioreactor, while the tissue displacement and thickness were nondestructively measured using ultrasound. The material parameters were estimated by an inverse finite element scheme, which was preceded by an analytical estimation step to rapidly obtain an initial estimate that already approximated the final solution. The efficiency and accuracy of the two-step inverse method was demonstrated on virtual experiments of several material types with known parameters. PDMS samples were used to demonstrate the method's feasibility, where it was shown that the proposed method yielded similar results to tensile testing. Finally, the method was applied to estimate the material properties of tissue-engineered constructs. Via this method, the evolution of mechanical properties during tissue growth and remodeling can now be monitored in a well-controlled system. The outcomes can be used to determine various mechanical constituents and to assess their contribution to mechanical homeostasis.

A Methodology for Verifying Abdominal Aortic Aneurysm Wall Stress

An abdominal aortic aneurysm (AAA) is a permanent focal dilatation of the abdominal aorta of at least 1.5 times its normal diameter. Although the criterion of maximum diameter is still used in clinical practice to decide on a timely intervention, numerical studies have demonstrated the importance of other geometric factors. However, the major drawback of numerical studies is that they must be validated experimentally before clinical implementation. This work presents a new methodology to verify wall stress predicted from the numerical studies against the experimental testing. To this end, four AAA phantoms were manufactured using vacuum casting. The geometry of each phantom was subject to microcomputed tomography (μCT) scanning at zero and three other intraluminal pressures: 80, 100, and 120 mm Hg. A zero-pressure geometry algorithm was used to calculate the wall stress in the phantom, while the numerical wall stress was calculated with a finite-element analysis (FEA) solver based on the actual zero-pressure geometry subjected to 80, 100, and 120 mm Hg intraluminal pressure loading. Results demonstrate the moderate accuracy of this methodology with small relative differences in the average wall stress (1.14%). Additionally, the contribution of geometric factors to the wall stress distribution was statistically analyzed for the four phantoms. The results showed a significant correlation between wall thickness and mean curvature (MC) with wall stress.

Investigation of inhomogeneous and anisotropic material behavior of porcine thoracic aorta using nano-indentation tests

This study investigates the inhomogeneity and anisotropy of porcine descending thoracic aorta in three dimensions using a custom-made nano-indentation technique and a quasi-linear viscoelastic modeling approach. The indentation tests were conducted in axial, circumferential, and radial orientations with about 100 μm spatial resolution. The ratio of the elastic moduli obtained in different orientations was used to quantify the tissue local anisotropy. The distal sections were generally stiffer than the proximal ones in both axial and circumferential indentations. Four distinct layers were identified across the thickness with significantly different mechanical properties. The stiffness of the medial quadrant was significantly lower than all other quadrants in axial indentation. The anisotropic behavior of the tissue was more pronounced in the lateral quadrant of the distal sections. The results of this study can be used to better understand the mechanisms of aorta deformation and improve the spatial accuracy of computational models of aorta.

Verification of a virtual fields method to extract the mechanical properties of human optic nerve head tissues in vivo

We aimed to verify a custom virtual fields method (VFM) to estimate the patient-specific biomechanical properties of human optic nerve head (ONH) tissues, given their full-field deformations induced by intraocular pressure (IOP). To verify the accuracy of VFM, we first generated 'artificial' ONH displacements from predetermined (known) ONH tissue biomechanical properties using finite element analysis. Using such deformations, if we are able to match back the known biomechanical properties, it would indicate that our VFM technique is accurate. The peripapillary sclera was assumed anisotropic hyperelastic, while all other ONH tissues were considered isotropic. The simulated ONH displacements were fed into the VFM algorithm to extract back the biomechanical properties. The robustness of VFM was also tested against rigid body motions and noise added to the simulated displacements. Then, the computational speed of VFM was compared to that of a gold-standard stiffness measurement method (inverse finite element method or IFEM). Finally, as proof of principle, VFM was applied to IOP-induced ONH deformation data (obtained from one subject's eye imaged with OCT), and the biomechanical properties of the prelamina and lamina cribrosa (LC) were extracted. From given ONH displacements, VFM successfully matched back the biomechanical properties of ONH tissues with high accuracy and efficiency. For all parameters, the percentage errors were less than 0.05%. Our method was insensitive to rigid body motions and was also able to recover the material parameters in the presence of noise. VFM was also found 125 times faster than the gold-standard IFEM. Finally, the estimated shear modulus for the prelamina and the LC of the studied subject's eye were 33.7 and 63.5 kPa, respectively. VFM may be capable of measuring the biomechanical properties of ONH tissues with high speed and accuracy. It has potential in identifying patient-specific ONH biomechanical properties in the clinic if combined with optical coherence tomography.

A phase-field approach to model fracture of arterial walls: Theory and finite element analysis

This study uses a recently developed phase-field approach to model fracture of arterial walls with an emphasis on aortic tissues. We start by deriving the regularized crack surface to overcome complexities inherent in sharp crack discontinuities, thereby relaxing the acute crack surface topology into a diffusive one. In fact, the regularized crack surface possesses the property of Gamma-Convergence, i.e. the sharp crack topology is restored with a vanishing length-scale parameter. Next, we deal with the continuous formulation of the variational principle for the multi-field problem manifested through the deformation map and the crack phase-field at finite strains which leads to the Euler-Lagrange equations of the coupled problem. In particular, the coupled balance equations derived render the evolution of the crack phase-field and the balance of linear momentum. As an important aspect of the continuum formulation we consider an invariant-based anisotropic constitutive model which is additively decomposed into an isotropic part for the ground matrix and an exponential anisotropic part for the two families of collagen fibers embedded in the ground matrix. In addition we propose a novel energy-based anisotropic failure criterion which regulates the evolution of the crack phase-field. The coupled problem is solved using a one-pass operator-splitting algorithm composed of a mechanical predictor step (solved for the frozen crack phase-field parameter) and a crack evolution step (solved for the frozen deformation map); a history field governed by the failure criterion is successively updated. Subsequently, a conventional Galerkin procedure leads to the weak forms of the governing differential equations for the physical problem. Accordingly, we provide the discrete residual vectors and a corresponding linearization yields the element matrices for the two sub-problems. Finally, we demonstrate the numerical performance of the crack phase-field model by simulating uniaxial extension and simple shear fracture tests performed on specimens obtained from a human aneurysmatic thoracic aorta. Model parameters are obtained by fitting the set of novel experimental data to the predicted model response; the finite element results agree favorably with the experimental findings.

Biaxial rupture properties of ascending thoracic aortic aneurysms

Although hundreds of samples obtained from ascending thoracic aortic aneurysms (ATAA) of patients undergoing elective surgical repair have already been characterized biomechanically, their rupture properties were always derived from uniaxial tensile tests. Due to their bulge shape, ATAAs are stretched biaxially in vivo. In order to understand the biaxial rupture of ATAAs, our group developed a novel methodology based on bulge inflation and full-field optical measurements. The objective of the current paper is threefold. Firstly, we will review the failure properties (maximum stress, maximum stretch) obtained by bulge inflation testing on a cohort of 31 patients and compare them with failure properties obtained by uniaxial tension in a previously published study. Secondly, we will investigate the relationship between the failure properties and the age of patients, showing that patients below 55years of age display significantly higher strength. Thirdly, we will define a rupture risk based on the extensibility of the tissue and we will show that this rupture risk is strongly correlated with the physiological elastic modulus of the tissue independently of the age, ATAA diameter or the aortic valve phenotype of the patient.

STATEMENT OF SIGNIFICANCE

Despite their medical importance, rupture properties of ascending thoracic aortic aneurysms (ATAA) subjected to biaxial tension were inexistent in the literature. In order to address this lack, our group developed a novel methodology based on bulge inflation and full-field optical measurements. Here we report rupture properties obtained with this methodology on 31 patients. It is shown for the first time that rupture occurs when the stretch applied to ATAAs reaches the maximum extensibility of the tissue and that this maximum extensibility correlates strongly with the elastic properties. The outcome is a better detection of at-risk individuals for elective surgical repair.

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

Mechanical analysis has the potential to provide complementary information to aneurysm morphology in assessing its vulnerability. Reliable calculations require accurate material properties of individual aneurysmal components. Quantification of extreme extensibility and ultimate material strength of the tissue are important if rupture is to be modelled. Tissue pieces from 11 abdomen aortic aneurysm (AAA) from patients scheduled for elective surgery and from 8 normal aortic artery (NAA) from patients who scheduled for kidney/liver transplant were collected at surgery and banked in liquid nitrogen with the use of Cryoprotectant solution to minimize frozen damage. Prior to testing, specimen were thawed and longitudinal and circumferential tissue strips were cut from each piece and adventitia, media and thrombus if presented were isolated for the material test. The incremental Young’s modulus of adventitia of NAA was direction-dependent at low stretch levels, but not the media. Both adventitia and media had a similar extreme extensibility in the circumferential direction, but the adventitia was much stronger. For aneurysmal tissues, no significant differences were found when the incremental moduli of adventitia, media or thrombus in both directions were compared. Adventitia and media from AAA had similar extreme extensibility and ultimate strength in both directions and thrombus was the weakest material. Adventitia and media from AAA were less extensible compared with those of NAA, but the ultimate strength remained similar. The material properties, including extreme extensibility and ultimate strength, of both healthy aortic and aneurysmal tissues were layer-dependent, but not direction-dependent.

Biomechanical evaluation of ascending aortic aneurysms

The biomechanical properties of ascending aortic aneurysms were investigated only in the last decade in a limited number of studies. Indeed, in recent years, there has been a growing interest in this field in order to identify new predictive parameters of risk of dissection, which may have clinical relevance. The researches performed so far have been conducted according to the methods used in the study of abdominal aortic aneurysms. In most cases, uniaxial or biaxial tensile tests were used, while in a smaller number of studies other methods, such as opening angle, bulge inflation, and inflation-extension tests, were used. However, parameters and protocols of these tests are at present very heterogeneous in the studies reported in the literature, and, therefore, the results are not comparable and are sometimes conflicting. The purpose of this review then thence to provide a comprehensive analysis of the experimental methodology for determination of biomechanical properties in the specific field of aneurysms of the ascending aorta to allow for better comparison and understanding of the results.

Nontraumatic acute aortic emergencies: Part 2, Pre- and postsurgical complications related to aortic aneurysm in the emergency clinical setting

OBJECTIVE

The aim of this article is to illustrate the imaging findings and spectrum of disease entities affecting the aorta. The clinical presentation and assessment of acute aortic pathology can be elusive or deceptive, making the diagnosis challenging. The widespread availability of advanced cross-sectional imaging technology in the emergency setting puts the radiologist at the forefront of accurate and timely diagnosis.

CONCLUSION

Cross-sectional imaging plays a pivotal role in the diagnosis and delineation of aortic pathology. Awareness of the imaging findings and complications can help in swift and accurate diagnosis.

Biomechanics of porcine renal arteries and role of axial stretch

It is known that arteries experience significant axial stretches in vivo. Several authors have shown that the axial force needed to maintain an artery at its in vivo axial stretch does not change with transient cyclical pressurization over normal ranges. However, the axial force phenomenon of arteries has never been explained with microstructural considerations. In this paper we propose a simple biomechanical model to relate the specific axial force phenomenon of arteries to the predicted load-dependent average collagen fiber orientation. It is shown that (a) the model correctly predicts the authors' experimentally measured biaxial behavior of pig renal arteries and (b) the model predictions are in agreement with additional experimental results reported in the literature. Finally, we discuss the implications of the model for collagen fiber orientation and deposition in arteries.

Compared prediction of the experimental failure of a thin fibrous tissue by two macroscopic damage models

Several models for fibrous biological tissues have been proposed in the past, taking into account the fibrous microstructure through different homogenization methods. The aim of this paper is to compare theoretically and experimentally two existing homogenization methods - the Angular Integration method and the Generalized Structure Tensor method - by adapting them to a damage model for a planar fibrous tissue made of linear elastic and brittle fibers. The theoretical implementation of the homogenization methods reveals some differences once damage starts in the fibrous tissue; in particular, the anisotropy of the tissue evolves differently. The experimental aspect of this work consists in identifying the parameters of the damage model, with both homogenization methods, using inflation tests until rupture on a biological membrane. The numerical identification method is based on the simulation of the tests with the real geometry of the samples and the real boundary conditions computed by Stereo Digital Image Correlation. The identification method is applied to human liver capsule. The collagen fibers Young's modulus (19±6MPa) as well as their ultimate longitudinal strain (33±4%) are determined; no significant difference was observed between the two methods. However, by using the experimental boundary conditions, we could observe that the damage progression is faster for the Angular Integration version of the model.