Biaxial mechanical properties of intact and layer-dissected human carotid arteries at physiological and supraphysiological loadings

Nonlinear Elasticity of Blood Vessels and Vascular Grafts

The transplantation of vascular grafts has emerged as a prevailing approach to address vascular disorders. However, the development of small-diameter vascular grafts is still in progress, as they serve in a more complicated mechanical environment than their counterparts with larger diameters. The biocompatibility and functional characteristics of small-diameter vascular grafts have been well developed; however, mismatch in mechanical properties between the vascular grafts and native arteries has not been accomplished, which might facilitate the long-term patency of small-diameter vascular grafts. From a point of view in mechanics, mimicking the nonlinear elastic mechanical behavior exhibited by natural blood vessels might be the state-of-the-art in designing vascular grafts. This review centers on elucidating the nonlinear elastic behavior of natural blood vessels and vascular grafts. The biological functionality and limitations associated with as-reported vascular grafts are meticulously reviewed and the future trajectory for fabricating biomimetic small-diameter grafts is discussed. This review might provide a different insight from the traditional design and fabrication of artificial vascular grafts.

Significance of Dynamic Axial Stretching on Estimating Biomechanical Behavior and Properties of the Human Ascending Aorta

Contrary to most vessels, the ascending thoracic aorta (ATA) not only distends but also elongates in the axial direction. The purpose of this study is to investigate the biomechanical behavior of the ascending thoracic aorta (ATA) in response to dynamic axial stretching during the cardiac cycle. In addition, the implications of neglecting this dynamic axial stretching when estimating the constitutive model parameters of the ATA are investigated. The investigations were performed through in silico simulations by assuming a Gasser-Ogden-Holzapfel (GOH) constitutive model representative of ATA tissue material. The GOH model parameters were obtained from biaxial tests performed on four human ATA tissues in a previous study. Pressure-diameter curves were simulated as synthetic data to assess the effect of neglecting dynamic axial stretching on estimating constitutive model parameters. Our findings reveal a significant increase in axial stress (~ 16%) and stored strain energy (~ 18%) in the vessel when dynamic axial stretching is considered, as opposed to assuming a fixed axial stretch. All but one artery showed increased volume compliance while considering a dynamic axial stretching condition. Furthermore, we observe a notable difference in the estimated constitutive model parameters when dynamic axial stretching of the ATA is neglected, compared to the ground truth model parameters. These results underscore the critical importance of accounting for axial deformations when conducting in vivo biomechanical characterization of the ascending thoracic aorta.

Interstitial flow, pressure and residual stress in the aging carotid artery model in FEBio

Vascular smooth muscle cells (VSMCs) are subject to interstitial flow-induced shear stress, which is a critical parameter in cardiovascular disease progression. Transmural pressure loading and residual stresses alter the hydraulic conductivity of the arterial layers and modulate the interstitial fluid flux through the arterial wall. In this paper, a biphasic multilayer model of a common carotid artery (CCA) with anisotropic fiber-reinforced soft tissue and strain-dependent permeability is developed in FEBio software. After the verification of the numerical predictions, age-related arterial thickening and stiffening effects on arterial deformation and interstitial flow are computed under physiological geometry and physical parameters. We found that circumferential residual stress shifts outward in each layer and its gradient increases up to 6 times with aging. Internally pressurized CCA displays nonlinear deformation. In the aged artery, the circumferential stress becomes greater on the media layer (82-158 kPa) and lower on the intima and adventitia (19-23 kPa and 25-28 kPa, respectively). The radial compression of the intima reduces the total hydraulic conductivity by 48% in the young and 16% in the aged arterial walls. Consequently, the average radial interstitial flux increases with pressure by 14% in the young and 91% in the aged arteries. Accordingly, the flow shear stress experienced by the VSMCs becomes more significant for aged arteries, which may accelerate cardiovascular disease progression compared to young arteries.

Quantitative and large-format histochemistry to characterize peripheral artery compositional gradients

The femoropopliteal artery (FPA) is a long, flexible vessel that travels down the anteromedial compartment of the thigh as the femoral artery and then behind the kneecap as the popliteal artery. This artery undergoes various degrees of flexion, extension, and torsion during normal walking movements. The FPA is also the most susceptible peripheral artery to atherosclerosis and is where peripheral artery disease manifests in 80% of cases. The connection between peripheral artery location, its mechanical flexion, and its physiological or pathological biochemistry has been investigated for decades; however, histochemical methods remain poorly leveraged in their ability to spatially correlate normal or abnormal extracellular matrix and cells with regions of mechanical flexion. This study generates new histological image processing pipelines to quantitate tissue composition across high-resolution FPA regions-of-interest or low-resolution whole-section cross-sections in relation to their anatomical locations and flexions during normal movement. Comparing healthy ovine femoral, popliteal, and cranial-tibial artery sections as a pilot, substantial arterial contortion was observed in the distal popliteal and cranial tibial regions of the FPA which correlated with increased vascular smooth muscle cells and decreased elastin content. These methods aim to aid in the quantitative characterization of the spatial distribution of extracellular matrix and cells in large heterogeneous tissue sections such as the FPA. RESEARCH HIGHLIGHTS: Large-format histology preserves artery architecture. Elastin and smooth muscle content is correlated with distance from heart and contortion during flexion. Cell and protein analyses are sensitive to sectioning plane and image magnification.

Genetic deletion of MMP12 ameliorates cardiometabolic disease by improving insulin sensitivity, systemic inflammation, and atherosclerotic features in mice

BACKGROUND

Matrix metalloproteinase 12 (MMP12) is a macrophage-secreted protein that is massively upregulated as a pro-inflammatory factor in metabolic and vascular tissues of mice and humans suffering from cardiometabolic diseases (CMDs). However, the molecular mechanisms explaining the contributions of MMP12 to CMDs are still unclear.

METHODS

We investigated the impact of MMP12 deficiency on CMDs in a mouse model that mimics human disease by simultaneously developing adipose tissue inflammation, insulin resistance, and atherosclerosis. To this end, we generated and characterized low-density lipoprotein receptor (Ldlr)/Mmp12-double knockout (DKO) mice fed a high-fat sucrose- and cholesterol-enriched diet for 16-20 weeks.

RESULTS

DKO mice showed lower cholesterol and plasma glucose concentrations and improved insulin sensitivity compared with LdlrKO mice. Untargeted proteomic analyses of epididymal white adipose tissue revealed that inflammation- and fibrosis-related pathways were downregulated in DKO mice. In addition, genetic deletion of MMP12 led to alterations in immune cell composition and a reduction in plasma monocyte chemoattractant protein-1 in peripheral blood which indicated decreased low-grade systemic inflammation. Aortic en face analyses and staining of aortic valve sections demonstrated reduced atherosclerotic plaque size and collagen content, which was paralleled by an improved relaxation pattern and endothelial function of the aortic rings and more elastic aortic sections in DKO compared to LdlrKO mice. Shotgun proteomics revealed upregulation of anti-inflammatory and atheroprotective markers in the aortas of DKO mice, further supporting our data. In humans, MMP12 serum concentrations were only weakly associated with clinical and laboratory indicators of CMDs.

CONCLUSION

We conclude that the genetic deletion of MMP12 ameliorates obesity-induced low-grade inflammation, white adipose tissue dysfunction, biomechanical properties of the aorta, and the development of atherosclerosis. Therefore, therapeutic strategies targeting MMP12 may represent a promising approach to combat CMDs.

External airway splint placement for severe pediatric tracheobronchomalacia

OBJECTIVE

To present external airway splinting with bioabsorbable airway supportive devices (ASD) for severe, life-threatening cases of pediatric tracheomalacia (TM) or tracheobronchomalacia (TBM).

METHODS

A retrospective cohort was performed for 5 pediatric patients with severe TM or TBM who underwent ASD placement. Devices were designed and 3D-printed from a bioabsorbable material, polycaprolactone (PCL). Pre-operative planning included 3-dimensional airway modeling of tracheal collapse and tracheal suture placement using nonlinear finite element (FE) methods. Pre-operative modeling revealed that triads along the ASD open edges and center were the most effective suture locations for optimizing airway patency. Pediatric cardiothoracic surgery and otolaryngology applied the ASDs by suspending the trachea to the ASD with synchronous bronchoscopy. Respiratory needs were trended for all cases. Data from pediatric patients with tracheostomy and diagnosis of TM or TBM, but without ASD, were included for discussion.

RESULTS

Five patients (2 Females, 3 Males, ages 2-9 months at time of ASD) were included. Three patients were unable to wean from respiratory support after vascular ring division; all three weaned to room air post-ASD. Two patients received tracheostomies prior to ASD placement, but continued to experience apparent life-threatening events (ALTE) and required ventilation with supraphysiologic ventilator settings. One patient weaned respiratory support successfully after ASD placement. The last patient died post-ASD due to significant respiratory co-morbidity.

CONCLUSION

ASD can significantly benefit patients with severe, unrelenting tracheomalacia or tracheobronchomalacia. Proper multidisciplinary case deliberation and selection are key to success with ASD. Pre-operative airway modeling allows proper suture placement to optimally address the underlying airway collapse.

Modeling planar response of vascular tissues using quadratic functions of effective strain

Simulation-based studies of the cardiovascular structure such as aorta have become increasingly popular for many biomedical problems such as predictions of aneurysm rupture. A critical step in these simulations is the development of constitutive models that accurately describe the tissue's mechanical behavior. In this work, we present a new constitutive model, which explicitly accounts for the gradual recruitment of collagen fibers. The recruitment is considered using an effective stretch, which is a continuum-scale kinematic variable measuring the uncrimped stretch of the tissue in an average sense. The strain energy of a fiber bundle is described by a quadratic function of the effective strain. Constitutive models formulated in this manner are applied to describe the responses of ascending thoracic aortic aneurysm and porcine thoracic aorta tissues. The heterogeneous properties of the ATAA tissue are extracted from bulge inflation test data, and then used in finite element analysis to simulate the inflation test. The descriptive and predictive capabilities are further assessed using planar testing data of porcine thoracic aortic tissues. It is found that the constitutive model can accurately describe the stress-strain relations. In particular, the finite element simulation replicates the displacement, strain, and stress distributions with excellent fidelity.

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).

Inverse material parameter estimation of patient-specific finite element models at the carotid bifurcation: The impact of excluding the zero-pressure configuration and residual stress

The carotid bifurcation experiences a complex loading environment due to its anatomical structure. Previous in-vivo material parameter estimation methods often use simplified model geometries, isotropic hyperelastic constitutive equations or neglect key aspects of the vessel, such as the zero-pressure configuration or residual stress, all of which have independently been shown to alter the stress environment of the vessel wall. Characterizing the location of high stress in the vessel wall has often been proposed as a potential indicator of structural weakness. However, excluding the afore-mentioned zero-pressure configuration, residual stress and patient-specific material parameters can lead to an incorrect estimation of the true stress values observed, meaning that stress alone as a risk indicator of rupture is insufficient. In this study, we investigate how the estimated material parameters and overall stress distributions in geometries of carotid bifurcations, extracted from in-vivo MR images, alter with the inclusion of the zero-pressure configuration and residual stress. This approach consists of the following steps: (1) geometry segmentation and hexahedral meshing from in-vivo magnetic resonance images (MRI) at two known phases; (2) computation of the zero-pressure configuration and the associated residual stresses; (3) minimization of an objective function built on the difference between the stress states of an "almost true" stress field at two known phases and a "deformed" stress field by altering the input material parameters to determine patient-specific material properties; and (4) comparison of the stress distributions throughout these carotid bifurcations for all cases with estimated material parameters. This numerical approach provides insights into the need for estimation of both the zero-pressure configuration and residual stress for accurate material property estimation and stress analysis for the carotid bifurcation, establishing the reliability of stress as a rupture risk metric.

Tunable elastomer materials with vascular tissue-like rupture mechanics behavior

Purpose:Laboratory models of human arterial tissues are advantageous to examine the mechanical response of blood vessels in a simplified and controllable manner. In the present study, we investigated three silicone-based materials for replicating the mechanical properties of human arteries documented in the literature.Methods:We performed uniaxial tensile tests up to rupture on Sylgard184, Sylgard170 and DowsilEE-3200 under different curing conditions and obtained their True (Cauchy) stress-strain behavior and Poisson's ratios by means of digital image correlation (DIC). For each formulation, we derived the constitutive parameters of the 3-term Ogden model and designed numerical simulations of tubular models under a radial pressure of 250mmHg.Results:Each material exhibits evident non-linear hyperelasticity and dependence on the curing condition. Sylgard184 is the stiffest formulation, with the highest shear moduli and ultimate stresses at relative low strains (µ₁₈₄=0.52-0.88MPa, σ₁₈₄=15.90-16.54MPa, ε₁₈₄=0.72-0.96). Conversely, Sylgard170 and DowsilEE-3200 present significantly lower shear moduli and ultimate stresses that are closer to data reported for arterial tissues (µ₁₇₀=0.33-0.7MPa σ₁₇₀=2.61-3.67MPa, ε₁₇₀=0.69-0.81; µ_dow=0.02-0.09MPa σ_dow=0.83-2.05MPa, ε_dow=0.91-1.05). Under radial pressure, all formulations except DowsilEE-3200 at 1:1 curing ratio undergo circumferential stresses that remain in the elastic region with values ranging from 0.1 to 0.18MPa. Conclusion: Sylgard170 and DowsilEE-3200 appear to better reproduce the rupture behavior of vascular tissues within their typical ultimate stress and strain range. Numerical models demonstrate that all three materials achieve circumferential stresses similar to human common carotid arteries (Sommer et al. 2010), making these formulations suited for cylindrical laboratory models under physiological and supraphysiological loading.

Impact of physiological loads of arterial wall on nucleus deformation in endothelial cells: A computational study

INTRODUCTION

Computational modeling can enhance the understanding of cell mechanics. To achieve this, finite element models of endothelial cells were proposed with shapes mimicking their natural state inside the endothelium within the cardiovascular system. Implementing the recently proposed bendo-tensegrity concept, these models consider flexural (buckling) as well as tensional/compressional behavior of microtubules and also incorporate the waviness of intermediate filaments.

MATERIALS AND METHODS

Four different models were created (flat and domed hexagons, both regular and elongated in the direction of blood flow) and loaded by biaxial deformation, blood pressure, and shear load from blood flow - natural physiological conditions of the arterial endothelium - aiming to investigate the "in situ" mechanical response of the cell.

RESULTS

The impact of individual components of loads on the nucleus deformation (more specifically on the first principal strain) potentially influencing mechanotransduction was investigated and the role of the cytoskeleton and its constituents in the mechanical response of the endothelial cell was assessed. The results show (i) the impact of pulsating blood pressure on cyclic deformations of the nucleus, which increase substantially with decreasing axial pre-stretch of the cell, (ii) the importance of relatively low shear stresses in the cell response and nucleus deformation.

CONCLUSION

Not only the pulsatile blood pressure but also the wall shear stress may induce significant deformation of the nucleus and thus trigger remodelation processes in endothelial cells.

The effects of viscoelasticity on residual strain in aortic soft tissues

Residual stress is thought to play a critical role in modulating stress distributions in soft biological tissues and in maintaining the mechanobiological stress environment of cells. Residual stresses in arteries and other tissues are classically assessed through opening angle experiments, which demonstrate the continuous release of residual stresses over hours. These results are then assessed through nonlinear biomechanical models to provide estimates of the residual stresses in the intact state. Although well studied, these analyses typically focus on hyperelastic material models despite significant evidence of viscoelastic phenomena over both short and long timescales. In this work, we extended the state-of-the-art structural tensor model for arterial tissues from Holzapfel and Ogden for fractional viscoelasticity. Models were tuned to capture consistent levels of hysteresis observed in biaxial experiments, while also minimizing the fractional viscoelastic weighting and opening angles to correctly capture opening angle dynamics. Results suggest that a substantial portion of the human abdominal aorta is viscoelastic, but exhibits a low fractional order (i.e. more elastically). Additionally, a significantly larger opening angle in the fully unloaded state is necessary to produce comparable hysteresis in biaxial testing. As a consequence, conventional estimates of residual stress using hyperelastic approaches over-estimate their viscoelastic counterparts by a factor of 2. Thus, a viscoelastic approach, such as the one illustrated in this study, in combination with an additional source of rate-controlled viscoelastic data is necessary to accurately analyze the residual stress distribution in soft biological tissues. STATEMENT OF SIGNIFICANCE: Residual stress plays a crucial role in achieving a homeostatic stress environment in soft biological tissues. However, the analysis of residual stress typically focuses on hyperelastic material models despite evidence of viscoelastic behavior. This work is the first attempt at analyzing the effects of viscoelasticity on residual stress. The application of viscoelasticity was crucial for producing realistic opening dynamics in arteries. The overall residual stresses were estimated to be 50% less than those from using hyperelastic material models, while the opening angles were projected to be 25% more than those measured after 16 hours, suggesting underestimated residual strain. This study highlights the importance viscoelasticity in the analysis of residual stress even in weakly dissipative materials like the human aorta.

Physiological basis for longitudinal motion of the arterial wall

As opposed to arterial distension in the radial plane, longitudinal wall motion (LWM) is a multiphasic and bidirectional displacement of the arterial wall in the anterograde (i.e., in the direction of blood flow) and retrograde (i.e., opposing direction of blood flow) directions. While initially disregarded as imaging artifact, LWM has been consistently reported in ultrasound investigations in the last decade and is reproducible beat-to-beat, albeit with large inter-individual variability across healthy and diseased populations. Emerging literature has sought to examine the mechanistic control of LWM to explain the shape and variability of the motion pattern but lacks considerations for key foundational vascular principles at the level of the arterial wall ultrastructure. The purpose of this review is to summarize the potential factors that underpin the causes and control of arterial LWM, spanning considerations from the arterial extracellular matrix to systems-level integrative theories. First, an overview of LWM and relevant aspects wall composition will be discussed, including major features of the multiphasic pattern, arterial wall extracellular components, tunica fiber orientations, and arterial longitudinal pre-stretch. Second, current theories on the systems-level physiological mechanisms driving LWM will be discussed in the context of available evidence including experimental human research, porcine studies, and mathematical models. Throughout, we discuss implications of these observations with suggestions for future priority research areas.

Intramural Distributions of GAGs and Collagen vs. Opening Angle of the Intact Porcine Aortic Wall

The heterogeneity and contribution of collagen and elastin to residual stresses have been thoroughly studied, but more recently, glycosaminoglycans (GAGs) also emerged as potential regulators. In this study, the opening angle of aortic rings (an indicator of circumferential residual stresses) and the mural distributions of sulfated GAGs (sGAG), collagen, and elastin were quantified in the ascending, aortic arch and descending thoracic regions of 5- to 6-month-old pigs. The opening angle correlated positively with the aortic ring's mean radius and thickness, with good and moderate correlations respectively. The correlations between the sGAG, collagen, elastin, and collagen:sGAG ratio and the opening angle were evaluated to identify aortic compositional factors that could play roles in regulating circumferential residual stresses. The total collagen:sGAG ratio displayed the strongest correlation with the opening angle (r = - 0.715, p < 0.001), followed by the total sGAG content which demonstrated a good correlation (r = 0.623, p < 0.001). Additionally, the intramural gradients of collagen, sGAG and collagen:sGAG correlated moderately with the opening angle. We propose that, in addition to the individual role sGAG play through their content and intramural gradient, the interaction between collagen and sGAG should be considered when evaluating circumferential residual stresses in the aorta.

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 a Collagen Fibre Remodelling Rupture Risk Metric for Potentially Vulnerable Carotid Artery Atherosclerotic Plaques

Atherosclerotic plaque rupture in carotid arteries can lead to stroke which is one of the leading causes of death or disability worldwide. The accumulation of atherosclerotic plaque in an artery changes the mechanical properties of the vessel. Whilst healthy arteries can continuously adapt to mechanical loads by remodelling their internal structure, particularly the load-bearing collagen fibres, diseased vessels may have limited remodelling capabilities. In this study, a local stress modulated remodelling algorithm is proposed to explore the mechanical response of arterial tissue to the remodelling of collagen fibres. This stress driven remodelling algorithm is used to predict the optimum distribution of fibres in healthy and diseased human carotid bifurcations obtained using Magnetic Resonance Imaging (MRI). In the models, healthy geometries were segmented into two layers: media and adventitia and diseased into four components: adventitia, media, plaque atheroma and lipid pool (when present in the MRI images). A novel meshing technique for hexahedral meshing of these geometries is also demonstrated. Using the remodelling algorithm, the optimum fibre patterns in various patient specific plaques are identified and the role that deviations from these fibre configurations in plaque vulnerability is shown. This study provides critical insights into the collagen fibre patterns required in carotid artery and plaque tissue to maintain plaque stability.

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.

Stretch and stress distributions in the human artery based on two-layer model considering residual stresses

The objective is to know the stress distributions in the arterial walls under residual stresses based on two-layer model. Human common carotid arteries were analysed to show stress distributions at physiological and supraphysiological intraluminal pressures. The analyses for the loaded states were performed with stretch ratios with reference to a Riemannian stress-free configuration which is a 3D non-Euclidean manifold due to the nonzero Riemann curvature tensor. The experimental data obtained by other literature were used for the common carotid arteries to analyse the stretch and stress distributions in the arterial wall although kinematics is different from the literature. The stretches and stresses were calculated for the unloaded state, i.e. the residual stretches and stresses. And those at the axial stretch ratio 1.1 with reference to the unloaded state were calculated at the intraluminal pressures 16, 50, and 100 kPa. The stresses increased from the inner surface to the outer surface at all pressures analysed. These results suggest that in the human arteries the mechanical loads are mainly supported with the adventitia even though the media and intima play an important role to control of physiological functions.

Stress distribution in the walls of major arteries: implications for atherogenesis

Background

There is a correlation between the sites of atheroma development and stress points in the arterial system. Generally, pulse pressure results in stresses acting on the vascular vessel, including longitudinal stress, radial or normal stress, tangential stress or hoop stress and shear stress. This paper explores the relationship between arterial wall shear stress and pulsatile blood pressure with the aim of furthering the understanding of atherogenesis and plaque progression.

Methods

We computed the magnitude of the shear stresses within the carotid bifurcation geometry of a patient and calculated the increase in shear stress levels that would occur when the blood pressure and pulse pressures rise during exertion. We also determined in which layer of the artery wall the maximum shear stress is located, and computed the shear stress at different levels within the media. We used the theory of laminate analysis, (Classical Laminate Plate Theory), to analyse the stress distribution on the carotid artery wall. Computational Fluid Dynamics (CFD) analysis was used on anatomy based on a CT angiogram of the carotid bifurcation of a patient with a 90% stenosis on the right side and 10% on the left. The pulsatile non-Newtonian blood flow with a resting blood pressure of 120/80 mmHg and an exertion pressure of 200/100 mmHg was simulated and the resultant forces were transferred to an ANSYS Composite PrepPost (ACP) model for wall shear stress analysis. A multilayer elastic, anisotropic, and inhomogeneous arterial wall (intima, internal elastic lamina, media, external elastic lamina, and adventitial layers) was modelled and the shear stress magnitudes and change over time between the layers was calculated.

Results

Shear stress in the individual composite layers is far greater than that acting on the endothelium (less than 5 Pa). At rest, the maximum variation of shear stress in the arterial wall occurs in the intima (138 Pa) and adventitia (135 Pa). The medial layer has the lowest variation of shear stress. Under severe exertion, the maximum shear stress magnitude in the intimal layer and the adjacent medial layer is near the ultimate stress level. The maximum/minimum shear stress ratios during the cardiac cycle vary most widely in the innermost part of the media, adjacent to the intima, with a four-fold ratio increase. This compares with a less than two-fold increase in all the other layers including the intima and adventitia, making the inner media the most vulnerable layer to mechanical injury.

Conclusions

This study showed that the magnitude of exertion-induced shear stress approaches the ultimate stress limit in the intima and the immediate adjacent medial layer. The variation in stress is maximal in the inner layer of the media. These findings correlate the site of atheroma development with the most vulnerable site for injury in the media and emphasise the impact of pulse pressure. Further biological studies are required to ascertain whether this leads to injury that initiates atheroma that then precipitates an injury/healing cycle.

A comparison of passive and active wall mechanics between elastic and muscular arteries of juvenile and adult rats

The elastic abdominal aorta and muscular femoral artery are susceptible to aneurysm and atherosclerosis, respectively. The vessel wall mechanics should be an important element for the difference. The objective of the study is to demonstrate a comparison of vessel wall mechanics between elastic and muscular arteries of juvenile and adult rats to show the changes of mechanical properties relevant to aging. The passive and active mechanical tests, theoretical analysis, and histological evaluation were carried out to investigate mechanical properties of vessel walls in the abdominal aorta and carotid and femoral arteries of young and adult rats. There are stiffening femoral artery, unchanged carotid artery, and distensible abdominal aorta in adult rats as compared with the young. The opening angle has values of 54 ± 13°, 82 ± 13°, and 94 ± 13° in the abdominal aorta and carotid and femoral arteries of adult rats, respectively, as well as 80 ± 22°, 93 ± 19°, and 100 ± 23° in the young. The findings are explained by the significantly reduced width of collagen fibers in the abdominal aorta, relatively unchanged width in the carotid artery, and significantly increased width in the femoral artery of adult rats as compared with the young. In conjunction with available literatures, we concluded that inconsistency for nonlinear age-related changes of artery wall mechanics occurs between arteries of different types, which may be a risk factor for the occurrence of abdominal aorta aneurysm and femoral artery atherosclerosis.

Effects of the Haemodynamic Stimulus on the Location of Carotid Plaques Based on a Patient-Specific Mechanobiological Plaque Atheroma Formation Model

In this work, we propose a mechanobiological atheroma growth model modulated by a new haemodynamic stimulus. To test this model, we analyse the development of atheroma plaques in patient-specific bifurcations of carotid arteries for a total time of 30 years. In particular, eight geometries (left or right carotid arteries) were segmented from clinical images and compared with the solutions obtained computationally to validate the model. The influence of some haemodynamical stimuli on the location and size of plaques is also studied. Plaques predicted by the mechanobiological models using the time average wall shear stress (TAWSS), the oscillatory shear index (OSI) and a new index proposed in this work are compared. The new index predicts the shape index of the endothelial cells as a combination of TAWSS and OSI values and was fitted using data from the literature. The mechanobiological model represents an evolution of the one previously proposed by the authors. This model uses Navier-Stokes equations to simulate blood flow along the lumen in the transient mode. It also employs Darcy's law and Kedem-Katchalsky equations for plasma and substance flow across the endothelium using the three-pore model. The mass balances of all the substances that have been considered in the model are implemented by convection-diffusion-reaction equations, and finally the growth of the plaques has been computed. The results show that by using the new mechanical stimulus proposed in this study, prediction of plaques is, in most cases, better than only using TAWSS or OSI with a minimal and maximal errors on stenosis ratio of 2.77 and 32.89 %, respectively. However, there are a few geometries in which haemodynamics cannot predict the location of plaques, and other biological or genetic factors would be more relevant than haemodynamics. In particular, the model predicts correctly eleven of the fourteen plaques presented in all the geometries considered. Additionally, a healthy geometry has been computed to check that plaque is not developed with the model in this case.

Age-related carotid extra-media thickening is associated with increased blood pressure and arterial stiffness

Ageing results in higher blood pressure and arterial stiffening leading to increased cardiovascular disease (CVD) risk. The extra-media thickness (EMT) is a composite measure of the arterial adventitia, perivascular adipose tissue (PVAT) and the jugular vein, but the association among EMT and ageing, blood pressure and arterial stiffness is largely unknown. We hypothesized that EMT is associated with ageing, blood pressure and arterial stiffness. Fifty (18M/32F, age range 20-79 years.) individuals underwent measures of EMT (media-adventitia border to jugular lumen interface) via ultrasonography, blood pressure (brachial; carotid), arterial stiffness (carotid beta-stiffness, distensibility and Young's modulus) and body composition (dual X-ray absorptiometry). Independent two-tailed t-tests compared characteristics between young and middle-aged/older adults (MA/O). Bivariate correlations assessed the relation between EMT, ageing, blood pressure and arterial stiffness endpoints. Partial correlations were used to adjust for age, sex, body mass index (BMI) and body fat percentage. Compared to young, MA/O adults had higher EMT, blood pressure, arterial stiffness, BMI and BF% (all, p ≤ 0.05). Carotid EMT was positively correlated with age (r = 0.46), brachial SBP (r = 0.32), carotid SBP (r = 0.42), PP (r = 0.42), beta-stiffness (r = 0.37) and Young's elastic modulus (r = 0.43) (all, p ≤ 0.05), and negatively correlated with carotid distensibility (r = -0.34, p ≤ 0.05). All endpoints correlated with EMT remained after adjusting for BMI, BF% and sex (p ≤ 0.05). These data suggest EMT is a clinically relevant target that may be associated with age-related CVD risk in humans, yet further investigation is required to elucidate the role of EMT in age-related increases in blood pressure and arterial stiffness.

Potential damage in pulmonary arterial hypertension: An experimental study of pressure-induced damage of pulmonary artery

Pulmonary arterial hypertension (PAH) is associated with elevated pulmonary arterial pressure. PAH prognosis remains poor with a 15% mortality rate within 1 year, even with modern clinical management. Previous clinical studies proposed wall shear stress (WSS) to be an important hemodynamic factor affecting cell mechanotransduction, growth and remodeling, and disease progress in PAH. However, WSS in vivo is typically at most 2.5 Pa and a doubt has been cast whether WSS alone can drive disease progress. Furthermore, our current understanding of PAH pathology largely comes from small animals' studies in which caliber enlargement, a hallmark of PAH in humans, is rarely reported. Therefore, a large-animal experiment on pulmonary arteries (PAs) is needed to validate whether increased pressure can induce enlargement of PAs caliber. In this study, we use an inflation testing device to characterize the mechanical behavior, both nonlinear elastic behavior and irreversible damage of porcine arteries. The parameters of elastic behavior are estimated from the inflation test at a low-pressure range before and after over-pressurization. Then, histological images are qualitatively examined for medial and adventitial layers. This study sheds light on the relevance of pressure-induced damage mechanism in human PAH.

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

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

Comparison of morphometric, structural, mechanical, and physiologic characteristics of human superficial femoral and popliteal arteries

Peripheral arterial disease differentially affects the superficial femoral (SFA) and the popliteal (PA) arteries, but their morphometric, structural, mechanical, and physiologic differences are poorly understood. SFAs and PAs from 125 human subjects (age 13-92, average 52±17 years) were compared in terms of radii, wall thickness, and opening angles. Structure and vascular disease were quantified using histology, mechanical properties were determined with planar biaxial extension, and constitutive modeling was used to calculate the physiologic stress-stretch state, elastic energy, and the circumferential physiologic stiffness. SFAs had larger radii than PAs, and both segments widened with age. Young SFAs were 5% thicker, but in old subjects the PAs were thicker. Circumferential (SFA: 96→193°, PA: 105→139°) and longitudinal (SFA: 139→306°, PA: 133→320°) opening angles increased with age in both segments. PAs were more diseased than SFAs and had 11% thicker intima. With age, intimal thickness increased 8.5-fold, but medial thickness remained unchanged (620μm) in both arteries. SFAs had 30% more elastin than the PAs, and its density decreased ~50% with age. SFAs were more compliant than PAs circumferentially, but there was no difference longitudinally. Physiologic circumferential stress and stiffness were 21% and 11% higher in the SFA than in the PA across all ages. The stored elastic energy decreased with age (SFA: 1.4→0.4kPa, PA: 2.5→0.3kPa). While the SFA and PA demonstrate appreciable differences, most of them are due to vascular disease. When pathology is the same, so are the mechanical properties, but not the physiologic characteristics that remain distinct due to geometrical differences.

Constitutive modeling using structural information on collagen fiber direction and dispersion in human superficial femoral artery specimens of different ages

Arterial mechanics plays an important role in vascular pathophysiology and repair, and advanced imaging can inform constitutive models of vascular behavior. We have measured the mechanical properties of 14 human superficial femoral arteries (SFAs) (age 12-70, mean 48±19 years) using planar biaxial extension, and determined the preferred collagen fiber direction and dispersion using multiphoton microscopy. The collagen fiber direction and dispersion were evaluated using second-harmonic generation imaging and modeled using bivariate von Mises distributions. The microstructures of elastin and collagen were assessed using two-photon fluorescence imaging and conventional bidirectional histology. The mechanical and structural data were used to describe the SFA mechanical behavior using two- and four-fiber family invariant-based constitutive models. Older SFAs were stiffer and mechanically more nonlinear than younger specimens. In the adventitia, collagen fibers were undulated and diagonally-oriented, while in the media, they were straight and circumferentially-oriented. The media was rich in collagen that surrounded the circumferentially-oriented smooth muscle cells, and the elastin was present primarily in the internal and external elastic laminae. Older SFAs had a more circumferential collagen fiber alignment, a decreased circumferential-radial fiber dispersion, but the same circumferential-longitudinal fiber dispersion as younger specimens. Both the two- and the four-fiber family constitutive models were able to capture the experimental data, and the fits were better for the four-fiber family formulation. Our data provide additional details on the SFA intramural structure and inform structurally-based constitutive models.

Vascular Grafts with Tailored Stiffness and a Ligand Environment via Multiarmed Polymer Sheath for Expeditious Regeneration

The bypass graft is the mainstream of surgical intervention to treat vascular diseases. Ideal bypass materials, yet to be developed, require mechanical properties, availability, clinically feasible manufacturing logistics, and bioactivities with precise physicochemical cues defined to guide cell activities for arterial regeneration. Such needs instigated our fabrication of vascular grafts, which consist of coaxial, nanostructured fibers exhibiting a polycaprolactone (PCL) core and a photoclickable, 4-arm thiolated polyethylene glycol-norbornene (PEG-NB) sheath. The graft strength and bioactivity were modulated by the PCL concentration and the peptides (RGD, transforming growth factor β-1 or TGF-β1) conjugated to thiol-ene of PEG-NB, respectively. Structural, physical, and mechanical characterizations demonstrated that the fibrous grafts mimicked the key features of the native extracellular matrix, including a crosslinked fiber network for structural stability, viscoelasticity emulating arteries, hydration property, and high porosity for cell infiltration. Meanwhile, these grafts displayed strength and toughness exceeding or meeting surgical criteria. Furthermore, the grafts with higher PCL concentration (3 vs 1.8%) showed thicker fibers, lower porosity and pore size, and increased elastic and storage moduli. Graft bioactivity was determined by the mesenchymal stem cell (MSC) behaviors on the grafts and arterial regeneration in vivo using interposition grafting. Results showed that the cell adhesion and proliferation increased with the RGD density (25 vs 5 mM). After 1 week implantation, all peptide-functionalized PCL/PEG-NB grafts with or without MSC preseeding, as opposed to PCL grafts, showed expeditious endothelial lining, abundant vascular cell infiltration, and matrix production. Compared to RGD grafts, RGD/TGF-β1 grafts enhanced MSC differentiation into smooth muscle cells in vitro and developed thicker smooth muscle cell layers in vivo. Overall, the versatile porous vascular grafts offer superior properties and tunability for future translation.

Mechanical, structural, and physiologic differences in human elastic and muscular arteries of different ages: Comparison of the descending thoracic aorta to the superficial femoral artery

Elastic and muscular arteries differ in structure, function, and mechanical properties, and may adapt differently to aging. We compared the descending thoracic aortas (TA) and the superficial femoral arteries (SFA) of 27 tissue donors (average 41±18 years, range 13-73 years) using planar biaxial testing, constitutive modeling, and bidirectional histology. Both TAs and SFAs increased in size with age, with the outer radius increasing more than the inner radius, but the TAs thickened 6-fold and widened 3-fold faster than the SFAs. The circumferential opening angle did not change in the TA, but increased 2.4-fold in the SFA. Young TAs were relatively isotropic, but the anisotropy increased with age due to longitudinal stiffening. SFAs were 51% more compliant longitudinally irrespective of age. Older TAs and SFAs were stiffer, but the SFA stiffened 5.6-fold faster circumferentially than the TA. Physiologic stresses decreased with age in both arteries, with greater changes occurring longitudinally. TAs had larger circumferential, but smaller longitudinal stresses than the SFAs, larger cardiac cycle stretch, 36% lower circumferential stiffness, and 8-fold more elastic energy available for pulsation. TAs contained elastin sheets separated by smooth muscle cells (SMCs), collagen, and glycosaminoglycans, while the SFAs had SMCs, collagen, and longitudinal elastic fibers. With age, densities of elastin and SMCs decreased, collagen remained constant due to medial thickening, and the glycosaminoglycans increased. Elastic and muscular arteries demonstrate different morphological, mechanical, physiologic, and structural characteristics and adapt differently to aging. While the aortas remodel to preserve the Windkessel function, the SFAs maintain higher longitudinal compliance.

Fatigue of textiles used in vascular surgery: Application to carotid endarterectomy

Fatigue is a material-based phenomenon playing a significant role in the mechanical behavior of components and structures. Although fatigue has been well studied for traditional materials, such as metals, its underlying mechanisms are not thoroughly understood in novel applications such as the case of textiles used as patches to close the arteriotomy in carotid endarterectomy. The latter is a type of vascular surgery for the treatment of carotid artery disease in which after an arteriotomy and removal of atherosclerotic plaque closure is made with a patch sutured on the artery. Completion of the operation signals the initiation of complex mechanical and hemodynamic phenomena. Fatigue performance of the patch eventually determines the successful outcome of carotid endarterectomy. In this study, we evaluate with a two-fold approach the mechanics of patch angioplasty in carotid endarterectomy. First, an analytical model for the fatigue behavior of textiles is developed, considering the microstructure and geometry of the fabric. Then, the surgical procedure is simulated and a finite element analysis of the endarterectomized and patched carotid artery is employed. Stress fields are calculated, while deformation at the site of patch angioplasty indicates a potential cause for the formation of aneurismal degeneration after the surgery. Such analysis can provide a better understanding in the establishment of follow-up protocols.

Consideration of stiffness of wall layers is decisive for patient-specific analysis of carotid artery with atheroma

The paper deals with the impact of chosen geometric and material factors on maximal stresses in carotid atherosclerotic plaque calculated using patient-specific finite element models. These stresses are believed to be decisive for the plaque vulnerability but all applied models suffer from inaccuracy of input data, especially when obtained in vivo only. One hundred computational models based on ex vivo MRI are used to investigate the impact of wall thickness, MRI slice thickness, lipid core and fibrous tissue stiffness, and media anisotropy on the calculated peak plaque and peak cap stresses. The investigated factors are taken as continuous in the range based on published experimental results, only the impact of anisotropy is evaluated by comparison with a corresponding isotropic model. Design of Experiment concept is applied to assess the statistical significance of these investigated factors representing uncertainties in the input data of the model. The results show that consideration of realistic properties of arterial wall in the model is decisive for the stress evaluation; assignment of properties of fibrous tissue even to media and adventitia layers as done in some studies may induce up to eightfold overestimation of peak stress. The impact of MRI slice thickness may play a key role when local thin fibrous cap is present. Anisotropy of media layer is insignificant, and the stiffness of fibrous tissue and lipid core may become significant in some combinations.

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.

Assessment of mechanical indicators of carotid plaque vulnerability: Geometrical curvature metric, plaque stresses and damage in tissue fibres

Stroke is a major cause of death worldwide. The rupture of atherosclerotic carotid plaques is the leading single cause of stroke. Currently there is no accepted clinical measure to quantitatively assess the risk of carotid plaque rupture. Structural analyses of vulnerable plaques, using finite element (FE) analysis, have retrospectively found that regions of high stress tend to be the site of plaque rupture. The current study proposes a new clinical measure, based on plaque geometry, to assess the risk of carotid plaque rupture. This measure, named the weighted curvature difference, is based on the curvature of both the lumen and intima-media boundary, and the local plaque thickness. A series of idealized and realistic, 2-D and 3-D geometries are used to systematically assess this novel geometrical metric. The areas predicted to be at high risk of rupture using this geometrical metric are compared with areas of high stress, obtained from both isotropic and anisotropic material models. These results are also compared with areas in diseased carotid arteries that are predicted to have high damage accumulation in collagen fibres using a continuum damage model. Results show the new geometrical metric consistently predicts the locations of high stress in all of the vessel geometries examined. The drawbacks of using lumen curvature only as a risk measure are highlighted; particularly in the case of outward remodelled vessels. Weighted curvature difference shows great potential to be used as a metric to efficiently distinguish the rupture prone areas in a diseased vessels in a way that is independent of material properties.

A novel technique for the assessment of mechanical properties of vascular tissue

Accurate estimation of mechanical properties of the different atherosclerotic plaque constituents is important in assessing plaque rupture risk. The aim of this study was to develop an experimental set-up to assess material properties of vascular tissue, while applying physiological loading and being able to capture heterogeneity. To do so, a ring-inflation experimental set-up was developed in which a transverse slice of an artery was loaded in the radial direction, while the displacement was estimated from images recorded by a high-speed video camera. The performance of the set-up was evaluated using seven rubber samples and validated with uniaxial tensile tests. For four healthy porcine carotid arteries, material properties were estimated using ultrasound strain imaging in whole-vessel-inflation experiments and compared to the properties estimated with the ring-inflation experiment. A 1D axisymmetric finite element model was used to estimate the material parameters from the measured pressures and diameters, using a neo-Hookean and Holzapfel–Gasser–Ogden material model for the rubber and porcine samples, respectively. Reproducible results were obtained with the ring-inflation experiment for both rubber and porcine samples. Similar mean stiffness values were found in the ring-inflation and tensile tests for the rubber samples as 202 kPa and 206 kPa, respectively. Comparable results were obtained in vessel-inflation experiments using ultrasound and the proposed ring-inflation experiment. This inflation set-up is suitable for the assessment of material properties of healthy vascular tissue in vitro. It could also be used as part of a method for the assessment of heterogeneous material properties, such as in atherosclerotic plaques.

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.

A comparative study of methods used to generate the arterial fiber structure in a clinically relevant numerical analysis

The advanced constitutive material models of artery wall require the definition of the mean collagen fiber directions in the material configuration. There are several proposed methods; however, it is unclear how much does the fiber structures obtained by these methods differ one from the other and how much this difference may affect the results of the structural analysis of a clinically relevant scenario. Therefore, in this paper, we address this issue by presenting the results of the comparative study of our developed and currently state-of-the-art fiber definition methods. In addition, we present the verification of our developed numerical model that incorporates the extended Holzapfel-Gasser-Ogden (HGO) constitutive material model and the generalized prestressing algorithm (GPA). In the case of the patient-specific internal carotid artery (ICA), the percentage error of the mean fiber directions defined by different methods does not exceed 17.73% (at least 0.05%, at most 81.82%) and has negligible effect on the stress levels, as the percentage error of the mean circumferential Cauchy stress does not exceed 0.1%. Both fiber definition methods produce comparable fiber structure, but our proposed method has an advantage, as it does not depend on method and software used to model the arterial wall mechanics.

Vulnerability analysis on the interaction between Asymmetric stent and arterial layer

The utilization of Asymmetric stent for recovering atherosclerotic diseases, particularly non-symmetric obstruction, is a quite challenging breakthrough treatment. In terms of eccentric plaque, the non-uniform stiffness of arterial layer causes the increasingly complex issues of vulnerability. This study investigated the vulnerability of the interaction between the Asymmetric stent and the surrounding arterial layer using structural transient dynamic analysis in ANSYS. Four combinations of stent deployment, i.e. the Sinusoidal stent expanded by the offset balloon, the Sinusoidal stent expanded by the ordinary cylindrical balloon, the Asymmetric stent expanded by the offset balloon, and the Asymmetric stent expanded by the ordinary cylindrical balloon, are generated for this comparative study. Multilayer material properties from recent in vitro experiments are adopted for the surrounding arterial layer, such as a fibrous cap, lipid core, diseased-healthy intima, and diseased-healthy media. In order to address plaque vulnerability, the Cauchy stresses and Hencky strains are used for stress measure because of convenience in comparison with the uniaxial/biaxial tension test data. The location-specific threshold value from the diseased human carotid artery is adopted for rupture criteria. The simulation indicated that as regards the eccentric plaque, the plaque vulnerability is caused by the plaque shape and components rather than caused by the geometrical structure of the stent or balloon expansion method. Nevertheless, the non-symmetric inflation of balloon, which leads against the plaque, contributed to an increase in the vulnerability of fibrous cap of fibroatheroma plaque.

Mechanical stresses associated with flattening of human femoropopliteal artery specimens during planar biaxial testing and their effects on the calculated physiologic stress-stretch state

Planar biaxial testing is commonly used to characterize the mechanical properties of arteries, but stresses associated with specimen flattening during this test are unknown. We quantified flattening effects in human femoropopliteal arteries (FPAs) of different ages and determined how they affect the calculated arterial physiologic stress-stretch state. Human FPAs from 472 tissue donors (age 12-82 years, mean 53 ± 16 years) were tested using planar biaxial extension, and morphometric and mechanical characteristics were used to assess the flattening effects. Constitutive parameters for the invariant-based model were adjusted to account for specimen flattening and used to calculate the physiologic stresses, stretches, axial force, circumferential stiffness, and stored energy for the FPAs in seven age groups. Flattened specimens were overall 12 ± 4% stiffer longitudinally and 19 ± 11% stiffer circumferentially when biaxially tested. Differences between the stress-stretch curves adjusted and non-adjusted for the effects of flattening were relatively constant across all age groups longitudinally, but increased with age circumferentially. In all age groups, these differences were smaller than the intersubject variability. Physiologic stresses, stretches, axial force, circumferential stiffness, and stored energy were all qualitatively and quantitatively similar when calculated with and without the flattening effects. Stresses, stretches, axial force, and stored energy reduced with age, but circumferential stiffness remained relatively constant between 25 and 65 years of age suggesting a homeostatic target of 0.75 ± 0.02 MPa. Flattening effects associated with planar biaxial testing are smaller than the intersubject variability and have little influence on the calculated physiologic stress-stretch state of human FPAs.

Biaxial Mechanical Characterizations of Atrioventricular Heart Valves

Extensive biaxial mechanical testing of the atrioventricular heart valve leaflets can be utilized to derive optimal parameters used in constitutive models, which provide a mathematical representation of the mechanical function of those structures. This presented biaxial mechanical testing protocol involves (i) tissue acquisition, (ii) the preparation of tissue specimens, (iii) biaxial mechanical testing, and (iv) postprocessing of the acquired data. First, tissue acquisition requires obtaining porcine or ovine hearts from a local Food and Drug Administration-approved abattoir for later dissection to retrieve the valve leaflets. Second, tissue preparation requires using tissue specimen cutters on the leaflet tissue to extract a clear zone for testing. Third, biaxial mechanical testing of the leaflet specimen requires the use of a commercial biaxial mechanical tester, which consists of force-controlled, displacement-controlled, and stress-relaxation testing protocols to characterize the leaflet tissue's mechanical properties. Finally, post-processing requires the use of data image correlation techniques and force and displacement readings to summarize the tissue's mechanical behaviors in response to external loading. In general, results from biaxial testing demonstrate that the leaflet tissues yield a nonlinear, anisotropic mechanical response. The presented biaxial testing procedure is advantageous to other methods since the method presented here allows for a more comprehensive characterization of the valve leaflet tissue under one unified testing scheme, as opposed to separate testing protocols on different tissue specimens. The proposed testing method has its limitations in that shear stress is potentially present in the tissue sample. However, any potential shear is presumed negligible.

Anisotropic residual stresses in arteries

The paper provides a deepened insight into the role of anisotropy in the analysis of residual stresses in arteries. Residual deformations are modelled following Holzapfel and Ogden (Holzapfel and Ogden 2010, J. R. Soc. Interface 7, 787-799. ( doi:10.1098/rsif.2009.0357 )), which is based on extensive experimental data on human abdominal aortas (Holzapfel et al. 2007, Ann. Biomed. Eng. 35, 530-545. ( doi:10.1007/s10439-006-9252-z )) and accounts for both circumferential and axial residual deformations of the individual layers of arteries-intima, media and adventitia. Each layer exhibits distinctive nonlinear and anisotropic mechanical behaviour originating from its unique microstructure; therefore, we use the most general form of strain-energy function (Holzapfel et al. 2015, J. R. Soc. Interface 12, 20150188. ( doi:10.1098/rsif.2015.0188 )) to derive residual stresses for each layer individually. Finally, the systematic experimental data (Niestrawska et al. 2016, J. R. Soc. Interface 13, 20160620. ( doi:10.1098/rsif.2016.0620 )) on both mechanical and structural properties of the different layers of the human abdominal aorta facilitate our discussion on (i) the importance of anisotropy in modelling residual stresses; (ii) the variability of residual stresses within the same class of tissue, the abdominal aorta; (iii) the limitations of conventional opening angle method to account for complex residual deformations; and (iv) the effect of residual stresses on the loaded configuration of the aorta mimicking in vivo conditions.

An investigation into the role of different constituents in damage accumulation in arterial tissue and constitutive model development

Carotid atherosclerotic plaque rupture is one of the leading causes of stroke. Treatments for atherosclerosis can induce tissue damage during the deployment of an intravascular device or through external tissue clamping during surgery. In this paper, a constituent specific study was performed to investigate the role of the ground matrix and collagen fibres of arterial tissue in response to supra-physiological loads. Cyclic mechanical tests were conducted on intact and collagenase-digested strips of porcine common carotid arteries. Using these tests, four passive damage-relevant phenomena were studied, namely (i) Mullins effect, (ii) hysteresis, (iii) permanent set and (iv) matrix failure and fibre rupture. A constitutive model was also developed to capture all of these damage-relevant phenomena using a continuum damage mechanics approach. The implemented constitutive model was fit to experimental results for both intact and digested samples. The results of this work demonstrate the important role of the ground matrix in the permanent deformation of the arterial tissue under high loads. Supra-physiological load-induced tissue damage may play a key role in vascular remodelling in arteries with atherosclerosis or following interventional procedures.

Biomechanical properties and microstructure of neonatal porcine ventricles

Neonatal heart disorders represent a major clinical challenge, with congenital heart disease alone affecting 36,000 new-borns annually within the European Union. Surgical intervention to restore normal function includes the implantation of synthetic and biological materials; however, a lack of experimental data describing the mechanical behaviour of neonatal cardiac tissue is likely to contribute to the relatively poor short- and long-term outcome of these implants. This study focused on characterising the mechanical behaviour of neonatal cardiac tissue using a porcine model, to enhance the understanding of how this differs to the equivalent mature tissue. The biomechanical properties of neonatal porcine cardiac tissue were characterised by uniaxial tensile, biaxial tensile, and simple shear loading modes, using samples collected from the anterior and posterior walls of the right and left ventricles. Histological images were prepared using Masson's trichrome staining, to enable assessment of the microstructure and correlation with tissue behaviour. The mechanical tests demonstrated that the neonatal cardiac tissue is non-linear, anisotropic, viscoelastic and heterogeneous. Our data provide a baseline describing the biomechanical behaviour of immature porcine cardiac tissue. Comparison with published data also indicated that the neonatal porcine cardiac tissue exhibits one-half the stiffness of mature porcine tissue in uniaxial extension testing, one-third in biaxial extension testing, and one-fourth stiffness in simple shear testing; hence, it provides an indication as to the relative change in characteristics associated with tissue maturation. These data may prove valuable to researchers investigating neonatal cardiac mechanics.

Mechanical response of human subclavian and iliac arteries to extension, inflation and torsion

Peripheral vascular trauma due to injuries of the upper and lower limbs are life-threatening, and their treatment require rapid diagnosis and highly-qualified surgical procedures. Experienced surgeons have recognized that subclavian arteries, affected by injuries of the upper limbs, require a more careful handling due to fragility than common iliac arteries, which are may be affected by injures of the lower limbs. We investigated these two artery types with comparable diameter to evaluate the differences in the biomechanical properties between subclavian and iliac arteries. Human subclavian and common iliac arteries of 14 donors either from the right or the left side (age: 63 yrs, SD: 19,9 female and 5 male) were investigated. Extension-inflation-torsion experiments at different axial strains (0-20%), transmural pressures (0-200 mmHg) and torsion (±25°) on preconditioned arterial tubes were performed. Residual stresses in both circumferential and axial direction were determined. Additionally, the microstructure of the tissues was determined via second-harmonic generation imaging and by histological investigations. At physiological conditions (p_i=13.3 kPa, λ_z=1.1) common iliac arteries revealed higher Cauchy stresses in circumferential and axial directions but a more compliant response in the circumferential direction than subclavian arteries. Both arteries showed distinct stiffer behavior in circumferential than in axial direction. Circumferential stiffness of common iliac arteries at physiological conditions increased significantly with aging (r=-0.67,p=0.02). The median inversion stretches, where the axial force is basically independent of the transmural pressure, were determined to be 1.05 for subclavian arteries and 1.11 for common iliac arteries. Both arteries exhibited increased torsional stiffness, when either axial prestretch or inflation pressure was increased. Residual stresses in the circumferential direction were significantly lower for subclavian arteries than for common iliac arteries at measurements after 30 min (p=0.05) and 16hrs (p=0.01). Investigations of the collagen microstructure revealed different collagen fiber orientations and dispersions in subclavian and iliac arteries. The difference in the collagen microstructure revealed further that the adventitia seems to contribute significantly to the passive mechanical response of the tested arteries at physiological loadings. Histological investigations indicated pronounced thickened intimal layers in subclavian and common iliac arteries, with a thickness comparable to the adventitial layer. In conclusion, we obtained biomechanical differences between subclavian and common iliac arteries, which possibly resulted from their different mechanical loadings/environments and respective in vivo movements caused by their anatomical locations. The biomechanical differences explored in this study are well reflected by the microstructure of the collagen and the histology of the investigated arteries, and the results can improve trauma patient care and endovascular implant design.

STATEMENT OF SIGNIFICANCE

During surgical interventions surgeons experienced that subclavian arteries (SAs) supplying the upper extremities, appear more fragile and prone to damage during surgical repair than common iliac arteries (CIAs), supplying the lower extremities. To investigate this difference in a systematic way the aim of this study was to compare the biomechanical properties of these two arteries from the same donors in terms of geometry, extension-inflation-torsion behavior, residual stresses, microstructure, and histology. In regard to cardiovascular medicine the material behavior of aged human arteries is of crucial interest. Moreover, the investigation of SA is important as it can help to improve surgical procedures at this challenging location. Over the long-term it might well be of value in the construction of artificial arteries for substituting native arteries. In addition, the analysis of mechanical stresses can improve design and material choice for endovascular implants to optimize long-term implant function.

Development of asymmetric stent for treatment of eccentric plaque

The selection of stent and balloon type is decisive in the stenting process. In the treatment of an eccentric plaque obstruction, a symmetric expansion from stent dilatation generates nonuniform stress distribution, which may aggravate fibrous cap prone to rupture. This paper developed a new stent design to treat eccentric plaque using structural transient dynamic analysis in ANSYS. A non-symmetric structural geometry of stent is generated to obtain reasonable stress distribution safe for the arterial layer surrounding the stent. To derive the novel structural geometry, a Sinusoidal stent type is modified by varying struts length and width, adding bridges, and varying curvature width of struts. An end ring of stent struts was also modified to eliminate dogboning phenomenon and to reduce the Ectropion angle. Two balloon types were used to deploy the stent, an ordinary cylindrical and offset balloon. Positive modification results were used to construct the final non-symmetric stent design, called an Asymmetric stent. Analyses of the deformation characteristics, changes in surface roughness and induced stresses within intact arterial layer were subsequently examined. Interaction between the stent and vessel wall was implemented by means of changes in surface roughness and stress distribution analyses. The Palmaz and the Sinusoidal stent were used for a comparative study. This study indicated that the Asymmetric stent types reduced the central radial recoiling and the dogboning phenomenon. In terms of changes in surface roughness and induced stresses, the Asymmetric stent has a comparable effect with that of the Sinusoidal stent. In addition, it could enhance the distribution of surface roughening as expanded by an offset balloon.

In vivo repeatability of the pulse wave inverse problem in human carotid arteries

Accurate arterial stiffness measurement would improve diagnosis and monitoring for many diseases. Atherosclerotic plaques and aneurysms are expected to involve focal changes in vessel wall properties; therefore, a method to image the stiffness variation would be a valuable clinical tool. The pulse wave inverse problem (PWIP) fits unknown parameters from a computational model of arterial pulse wave propagation to ultrasound-based measurements of vessel wall displacements by minimizing the difference between the model and measured displacements. The PWIP has been validated in phantoms, and this study presents the first in vivo demonstration. The common carotid arteries of five healthy volunteers were imaged five times in a single session with repositioning of the probe and subject between each scan. The 1D finite difference computational model used in the PWIP spanned from the start of the transducer to the carotid bifurcation, where a resistance outlet boundary condition was applied to approximately model the downstream reflection of the pulse wave. Unknown parameters that were estimated by the PWIP included a 10-segment linear piecewise compliance distribution and 16 discrete cosine transformation coefficients for each of the inlet boundary conditions. Input data was selected to include pulse waves resulting from the primary pulse and dicrotic notch. The recovered compliance maps indicate that the compliance increases close to the bifurcation, and the variability of the average pulse wave velocity estimated through the PWIP is on the order of 11%, which is similar to that of the conventional processing technique which tracks the wavefront arrival time (13%).

Uncertainty quantification and sensitivity analysis of an arterial wall mechanics model for evaluation of vascular drug therapies

Quantification of the uncertainty in constitutive model predictions describing arterial wall mechanics is vital towards non-invasive assessment of vascular drug therapies. Therefore, we perform uncertainty quantification to determine uncertainty in mechanical characteristics describing the vessel wall response upon loading. Furthermore, a global variance-based sensitivity analysis is performed to pinpoint measurements that are most rewarding to be measured more precisely. We used previously published carotid diameter–pressure and intima–media thickness (IMT) data (measured in triplicate), and Holzapfel–Gasser–Ogden models. A virtual data set containing 5000 diastolic and systolic diameter–pressure points, and IMT values was generated by adding measurement error to the average of the measured data. The model was fitted to single-exponential curves calculated from the data, obtaining distributions of constitutive parameters and constituent load bearing parameters. Additionally, we (1) simulated vascular drug treatment to assess the relevance of model uncertainty and (2) evaluated how increasing the number of measurement repetitions influences model uncertainty. We found substantial uncertainty in constitutive parameters. Simulating vascular drug treatment predicted a 6% point reduction in collagen load bearing (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_\mathrm {coll}$$\end{document}Lcoll), approximately 50% of its uncertainty. Sensitivity analysis indicated that the uncertainty in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_{\mathrm {coll}}$$\end{document}Lcoll was primarily caused by noise in distension and IMT measurements. Spread in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_{\mathrm {coll}}$$\end{document}Lcoll could be decreased by 50% when increasing the number of measurement repetitions from 3 to 10. Model uncertainty, notably that in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_{\mathrm {coll}}$$\end{document}Lcoll, could conceal effects of vascular drug therapy. However, this uncertainty could be reduced by increasing the number of measurement repetitions of distension and wall thickness measurements used for model parameterisation.

Comparative in vivo analysis of the role of the adventitia and the endothelium on arterial mechanical function: relevance for aortic counterpulsation

PURPOSE

The comparative effect of the intimal and adventitial layers on arterial biomechanics control, in basal and altered conditions, remains to be elucidated. This study aimed (1) to characterize the arterial conduit (CF) and buffering (distensibility) function of the iliac arteries in in vivo animals, in which the intimal and adventitial layers were removed; (2) to determine the effects of intra-aortic ballon pumping (IABP) on simultaneously de-adventitialized (DA) and de-endothelialized (DE) iliac arteries before and after induced heart failure.

METHODS

Pressure and diameter signals were measured in the iliac arteries of sheep (n = 7) in which the adventitial and intima layer were removed. Intra-aortic balloon pump (IABP) assistance was used in a control state and after heart failure induction.

RESULTS

Both DE and DA determined significant changes in arterial diameter, distensibility and CF. Changes were higher after DA than after DE in terms of distensibility and CF (p<0.05). DA followed by DE (DA + DE) showed significant increases in arterial diameter and CF, accompanied by a decrease in distensibility (p<0.05) with respect to intact arteries. Heart failure induction caused significant hemodynamic changes without modifying the already impaired local biomechanical parameters. Nonsignificant improvements in the biomechanical parameters of DA + DE iliac arteries were observed during IABP before and after heart failure induction.

CONCLUSIONS

Biomechanical changes caused by DA of iliac arteries were more important than those observed after DE. The DA + DE arteries showed significant differences with respect to intact arteries and with DA or DE arteries. IABP-related effects on arterial mechanics were absent in DA + DE arteries.

Mechanical behavior of porcine thoracic aorta in physiological and supra-physiological intraluminal pressures

Understanding the mechanical behavior of aorta under supra-physiological loadings is an important aspect of modeling tissue behavior in various applications that involve large deformations. Utilizing inflation–extension experiments, the mechanical behavior of porcine descending thoracic aortic segments under physiological and supra-physiological intraluminal pressures was investigated. The pressure was changed in the range of 0–70 kPa and the deformation of the segment was determined in three dimensions using a custom-made motion capture system. An orthotropic Fung-type constitutive model was characterized by implementing a novel computationally efficient framework that ensured material stability for numerical simulations. The nonlinear rising trend of circumferential stretch ratio [Formula: see text] from outer toward inner wall was significantly increased at higher pressures. The increase in [Formula: see text] from physiological pressure (13 kPa) to 70 kPa was 13% at the outer wall and 22% at the inner wall; in this pressure range, the longitudinal stretch ratio [Formula: see text] increased 20%. A significant nonlinearity in the material behavior was observed as in the same pressure range, and the circumferential and longitudinal Cauchy stresses at the inner wall were increased 16 and 18 times, respectively. The overall constitutive model was verified in several loading paths in the [Formula: see text] space to confirm its applicability in multi-axial loading conditions.