Over length quantification of the multiaxial mechanical properties of the ascending, descending and abdominal aorta using Digital Image Correlation

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.

Location specific multi-scale characterization and constitutive modeling of pig aorta

The mechanical and structural behavior of the aorta depend on physiological functions and vary from proximal to distal. Understanding the relation between regionally varying mechanical and multi-scale structural response of aorta can be helpful to assess the disease outcomes. Therefore, this study investigated the variation in mechanical and multi-scale structural properties among the major segments of aorta such as ascending aorta (AA), descending aorta (DA) and abdominal aorta (ABA), and established a relation between mechanical and multi-structural parameters. The obtained results showed significant increase in anisotropy and nonlinearity from proximal to distal aorta. The change in periphery length and radii between load and stress free configuration was also found increasing far from the heart. Opening angle was significantly large for ABA than AA and DA (AA/DA vs ABA; p = 0.001). Mean circumferential residual stretch (ratio of mean periphery length at load and stress free configurations) was found decreasing between AA and DA, and then increasing between DA to ABA and its value was significantly more for ABA (AA vs DA; p = 0.041, AA vs ABA; p = 0.001, DA vs ABA; p = 0.001). The waviness of collagen fibers, collagen fiber content, collagen fibril diameter and total protein content were found significantly increasing from proximal to distal. Pearson correlation test showed a significant linear correlation between variation in mechanical and multi-scale structural parameters over the aortic length. Residual stretch was found positively correlated with collagen fiber content (r = 0.82) whereas opening angel was found positively correlated with total protein content (TPC) (r = 0.76).

Image-Based Finite Element Modeling Approach for Characterizing In Vivo Mechanical Properties of Human Arteries

Mechanical properties of the arterial walls could provide meaningful information for the diagnosis, management and treatment of cardiovascular diseases. Classically, various experimental approaches were conducted on dissected arterial tissues to obtain their stress-stretch relationship, which has limited value clinically. Therefore, there is a pressing need to obtain biomechanical behaviors of these vascular tissues in vivo for personalized treatment. This paper reviews the methods to quantify arterial mechanical properties in vivo. Among these methods, we emphasize a novel approach using image-based finite element models to iteratively determine the material properties of the arterial tissues. This approach has been successfully applied to arterial walls in various vascular beds. The mechanical properties obtained from the in vivo approach were compared to those from ex vivo experimental studies to investigate whether any discrepancy in material properties exists for both approaches. Arterial tissue stiffness values from in vivo studies generally were in the same magnitude as those from ex vivo studies, but with lower average values. Some methodological issues, including solution uniqueness and robustness; method validation; and model assumptions and limitations were discussed. Clinical applications of this approach were also addressed to highlight their potential in translation from research tools to cardiovascular disease management.

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

An ultrastructural 3D reconstruction method for observing the arrangement of collagen fibrils and proteoglycans in the human aortic wall under mechanical load

An insight into changes of soft biological tissue ultrastructures under loading conditions is essential to understand their response to mechanical stimuli. Therefore, this study offers an approach to investigate the arrangement of collagen fibrils and proteoglycans (PGs), which are located within the mechanically loaded aortic wall. The human aortic samples were either fixed directly with glutaraldehyde in the load-free state or subjected to a planar biaxial extension test prior to fixation. The aortic ultrastructure was recorded using electron tomography. Collagen fibrils and PGs were segmented using convolutional neural networks, particularly the ESPNet model. The 3D ultrastructural reconstructions revealed a complex organization of collagen fibrils and PGs. In particular, we observed that not all PGs are attached to the collagen fibrils, but some fill the spaces between the fibrils with a clear distance to the collagen. The complex organization cannot be fully captured or can be severely misinterpreted in 2D. The approach developed opens up practical possibilities, including the quantification of the spatial relationship between collagen fibrils and PGs as a function of the mechanical load. Such quantification can also be used to compare tissues under different conditions, e.g., healthy and diseased, to improve or develop new material models. STATEMENT OF SIGNIFICANCE: The developed approach enables the 3D reconstruction of collagen fibrils and proteoglycans as they are embedded in the loaded human aortic wall. This methodological pipeline comprises the knowledge of arterial mechanics, imaging with transmission electron microscopy and electron tomography, segmentation of 3D image data sets with convolutional neural networks and finally offers a unique insight into the ultrastructural changes in the aortic tissue caused by mechanical stimuli.

Biomechanical characterization and constitutive modeling of the layer-dissected residual strains and mechanical properties of abdominal porcine aorta

We analyze the residual stresses and mechanical properties of layer-dissected infrarenal abdominal aorta (IAA). We measured the axial pre-stretch and opening angle and performed uniaxial tests to study and compare the mechanical behavior of both intact and layer-dissected porcine IAA samples under physiological loads. Finally, some of the most popular anisotropic hyperelastic constitutive models (GOH and microfiber models) were proposed to estimate the mechanical properties of the abdominal aorta by least-square fitting of the recorded in-vitro uniaxial test results. The results show that the residual stresses are layer dependent. In all cases, we found that the OA in the media layer is lower than in the whole artery, the intima and the adventitia. For the axial pre-stretch, we found that the adventitia and the media were slightly stretched in the environment of the intact arterial strip, whereas the intima appears to be compressed. Regarding the mechanical properties, the media seems to be the softest layer over the whole deformation domain showing high anisotropy, while the intima and adventitia exhibit considerable stiffness and a lower anisotropy response. Finally, all the hyperelastic anisotropic models considered in this study provided a reasonable approximation of the experimental data. The GOH model showed the best fitting.

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.

Influence of Neighborhood Size and Cross-Correlation Peak-Fitting Method on Location Accuracy

A known technique to obtain subpixel resolution by using object tracking through cross-correlation consists of interpolating the obtained correlation function and then refining peak location. Although the technique provides accurate results, peak location is usually biased toward the closest integer coordinate. This effect is known as the peak-locking error and it strongly limits this calculation technique’s experimental accuracy. This error may differ depending on the scene and algorithm used to fit and interpolate the correlation peak, but in general, it may be attributed to a sampling problem and the presence of aliasing. Many studies in the literature analyze this effect in the Fourier domain. Here, we propose an alternative analysis on the spatial domain. According to our interpretation, the peak-locking error may be produced by a non-symmetrical sample distribution, thus provoking a bias in the result. According to this, the peak interpolant function, the size of the local domain and low-pass filters play a relevant role in diminishing the error. Our study explores these effects on different samples taken from the DIC Challenge database, and the results show that, in general, peak fitting with a Gaussian function on a relatively large domain provides the most accurate results.

Unraveling the multilayer mechanical response of aorta using layer-specific residual stresses and experimental properties

To test the capability of the multilayer model, we used previously published layer-specific experimental data relating to the axial pre-stretch, the opening angle, the fiber distribution obtained by polarized light microscopy measurements, and the uniaxial and biaxial response of the porcine descending and abdominal aorta. We fitted the mechanical behavior of each arterial layer using Gasser, Holzapfel and Ogden strain energy function using the dispersion parameter κ as phenomenological parameter obtained during the fitting procedure or computed from the experimental fiber distribution. A multilayer finite element model of the whole aorta with the dimensions of the circumferential and longitudinal strips were then built using layer-specific material parameters previously fitted. This model was used to capture the whole aorta response under uniaxial and biaxial stress states and to reproduce the response of the whole aorta to internal pressure. Our results show that a model based on a multilayer structure without residual stresses is unable to render the uniaxial and biaxial mechanical response of the aorta (R²=0.6954 and R²=0.8582 for descending thoracic aorta (DTA) and infrarenal abdominal aorta (IAA), respectively). Only an appropriate multilayer model that includes layer-specific residual stresses can reproduce the response of the whole aorta (R²=0.9787 and R²=0.9636 for DTA and IAA respectively). In addition, a multilayer model without residual stresses produces the same stress distribution as a monolayer model without residual stresses where the maximal value of circumferential and longitudinal stresses appears at the inner radius of the intima. Finally, if layer-specific residual stresses are not available, there is less error the stress distribution using a monolayer model with residual stresses that a multilayer model without residual stresses.

Intramural Glycosaminoglycans Distribution vs. Residual Stress in Porcine Ascending Aorta: a Computational Study. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020;2020:2816-2819. [PMID: 33018592 DOI: 10.1109/embc44109.2020.9176381]

In this computational modelling work, we explored the mechanical roles that various glycosaminoglycans (GAGs) distributions may play in the porcine ascending aortic wall, by studying both the transmural residual stress as well as the opening angle in aortic ring samples. A finite element (FE) model was first constructed and validated against published data generated from rodent aortic rings. The FE model was then used to simulate the response of porcine ascending aortic rings with different GAG distributions prescribed through the wall of the aorta. The results indicated that a uniform GAG distribution within the aortic wall did not induce residual stresses, allowing the aortic ring to remain closed when subjected to a radial cut. By contrast, a heterogeneous GAG distribution led to the development of residual stresses which could be released by a radial cut, causing the ring to open. The residual stresses and opening angle were shown to be modulated by the GAG content, gradient, and the nature of the transmural distribution.

Finite strain parametric HFGMC micromechanics of soft tissues

A micromechanical analysis is offered for the prediction of the global behavior of biological tissues. The analysis is based on the isotropic-hyperelastic behavior of the individual constituents (Collagen and Elastin), their volume fractions, and takes into account their detailed interactions. The present analysis predicts the instantaneous tensors from which the effective current first tangent tensor is established, thus providing the overall anisotropic constitutive behavior of the composite and the resulting field distribution in the composite. This is in contradistinction with the macroanalysis in which the composite internal energy, which involves unknown functions that depend on several strain invariants, must be proposed. The offered micromechanical analysis forms a generalization to the finite strain high-fidelity generalized method of cells (HFGMC) based on the homogenization technique for periodic composites to the parametric finite strain. This involves an arbitrary discretization of the repeating unit-cell of the periodic composites. Results are given for the response of the human abdominal aorta, which consists of three layered tissues: intima, media, and adventitia, all of which are composed out of the Collagen and Elastin. The isotropic-hyperelastic constituents (Mooney-Rivlin and Yeoh) of the composites are calibrated by utilizing available experimental data which describe the response of the tissue. Validation of the results is performed by comparison of the predicted Cauchy stress and stretches with the experimental measurements. In addition, results are given in the form of Cauchy stress and deformation gradient field distributions in the constituents of several tissues.

Variation of Axial Residual Strains Along the Course and Circumference of Human Aorta Considering Age and Gender

Our understanding of aortic biomechanics is customarily limited by lack of information on the axial residual stretches of the vessel in both humans and experimental animals that would facilitate the identification of its actual zero-stress state. The aim of this study was thus to acquire hitherto unreported quantitative knowledge of axial opening angle and residual stretches in different segments and quadrants of the human aorta according to age and gender. Twenty-three aortas were harvested during autopsy from the aortic root to the iliac bifurcation and were divided into ≥12 segments and 4 quadrants. Morphometric measurements were taken in the excised/curled configuration of rectangular strips considered to be under zero-stress using image-analysis software to study the axial/circumferential variation of axial opening angle, internal/external residual stretch, and thickness of the aortic wall. The measured data demonstrated: (1) an axial opening angle peak at the arch branches, decreasing toward the ascending and to a near-constant value in the descending thoracic aorta, and increasing in the abdominal aorta; (2) the variation of residual stretches resembled that of opening angle, but axial differences in external residual stretch were more prominent; (3) wall thickness showed a progressive diminution along the vessel; (4) the highest opening angle/residual stretches were found in the inner quadrant and the lowest in the outer quadrant; (5) the anterior was the thinnest quadrant throughout the aorta; (6) age caused thickening but greatly reduced axial opening angle/residual stretches, without differences between males and females.

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.

Parametric Evaluation of Errors Using Isolated Dots for Movement Measurement by Image Cross-Correlation

Digital Image Correlation (DIC) is a common tool for assessing the movement of objects in a scene. Among others, one of the most popular techniques consists of tracking a dotted texture imitating speckle patterns. In this work, we analyzed the individual dots that form this pattern in order to propose an optimum size, shape, and dynamic range that allows minimizing the tracking error. Tracking was accomplished by using normalized cross-correlation with peak interpolation in order to obtain subpixel accuracy. For the models here used, we show that dot radii of 30–40 px with 150 gray levels are enough to obtain an accurate subpixel tracking resolution. Also, we show that 0.002 px is the performance limit of this technique, being this limit in accordance with the experimentally achievable subpixel limit.