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

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

Review of cardiac-coronary interaction and insights from mathematical modeling

Cardiac-coronary interaction is fundamental to the function of the heart. As one of the highest metabolic organs in the body, the cardiac oxygen demand is met by blood perfusion through the coronary vasculature. The coronary vasculature is largely embedded within the myocardial tissue which is continually contracting and hence squeezing the blood vessels. The myocardium-coronary vessel interaction is two-ways and complex. Here, we review the different types of cardiac-coronary interactions with a focus on insights gained from mathematical models. Specifically, we will consider the following: (1) myocardial-vessel mechanical interaction; (2) metabolic-flow interaction and regulation; (3) perfusion-contraction matching, and (4) chronic interactions between the myocardium and coronary vasculature. We also provide a discussion of the relevant experimental and clinical studies of different types of cardiac-coronary interactions. Finally, we highlight knowledge gaps, key challenges, and limitations of existing mathematical models along with future research directions to understand the unique myocardium-coronary coupling in the heart. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

In Vivo Material Properties of Human Common Carotid Arteries: Trends and Sex Differences

INTRODUCTION

In vivo estimation of material properties of arterial tissue can provide essential insights into the development and progression of cardiovascular diseases. Furthermore, these properties can be used as an input to finite element simulations of potential medical treatments.

MATERIALS AND METHODS

This study uses non-invasively measured pressure, diameter and wall thickness of human common carotid arteries (CCAs) acquired in 103 healthy subjects. A non-linear optimization was performed to estimate material parameters of two different constitutive models: a phenomenological, isotropic model and a structural, anisotropic model. The effect of age, sex, body mass index and blood pressure on the parameters was investigated.

RESULTS AND CONCLUSION

Although both material models were able to model in vivo arterial behaviour, the structural model provided more realistic results in the supra-physiological domain. The phenomenological model predicted very high deformations for pressures above the systolic level. However, the phenomenological model has fewer parameters that were shown to be more robust. This is an advantage when only the physiological domain is of interest. The effect of stiffening with age, BMI and blood pressure was present for women, but not always for men. In general, sex had the biggest effect on the mechanical properties of CCAs. Stiffening trends with age, BMI and blood pressure were present but not very strong. The intersubject variability was high. Therefore, it can be concluded that finding a representative set of parameters for a certain age or BMI group would be very challenging. Instead, for purposes of patient-specific modelling of surgical procedures, we currently advise the use of patient-specific parameters.

Brain tissue mechanics is governed by microscale relations of the tissue constituents

Local mechanical tissue properties are a critical regulator of cell function in the central nervous system (CNS) during development and disorder. However, we still don't fully understand how the mechanical properties of individual tissue constituents, such as cell nuclei or myelin, determine tissue mechanics. Here we developed a model predicting local tissue mechanics, which induces non-affine deformations of the tissue components. Using the mouse hippocampus and cerebellum as model systems, we show that considering individual tissue components alone, as identified by immunohistochemistry, is not sufficient to reproduce the local mechanical properties of CNS tissue. Our results suggest that brain tissue shows a universal response to applied forces that depends not only on the amount and stiffness of the individual tissue constituents but also on the way how they assemble. Our model may unify current incongruences between the mechanics of soft biological tissues and the underlying constituents and facilitate the design of better biomedical materials and engineered tissues. To this end, we provide a freely-available platform to predict local tissue elasticity upon providing immunohistochemistry images and stiffness values for the constituents of the tissue.

A numerical investigation of the mechanics of intracranial aneurysms walls: Assessing the influence of tissue hyperelastic laws and heterogeneous properties on the stress and stretch fields

Numerical simulations have been extensively used in the past two decades for the study of intracranial aneurysms (IAs), a dangerous disease that occurs in the arteries that reach the brain and affect overall 3.2% of a population without comorbidity with up to 60% mortality rate, in case of rupture. The majority of those studies, though, assumed a rigid-wall model to simulate the blood flow. However, to also study the mechanics of IAs walls, it is important to assume a fluid-solid interaction (FSI) modeling. Progress towards more reliable FSI simulations is limited because FSI techniques pose severe numerical difficulties, but also due to scarce data on the mechanical behavior and material constants of IA tissue. Additionally, works that have investigated the impact of different wall modeling choices for patient-specific IAs geometries are a few and often with limited conclusions. Thus our present study investigated the effect of different modeling approaches to simulate the motion of an IA. We used three hyperelastic laws - the Yeoh law, the three-parameter Mooney-Rivlin law, and a Fung-like law with a single parameter - and two different ways of modeling the wall thickness and tissue mechanical properties - one assumed that both were uniform while the other accounted for the heterogeneity of the wall by using a "hemodynamics-driven" approach in which both thickness and material constants varied spatially with the cardiac-cycle-averaged hemodynamics. Pulsatile numerical simulations, with patient-specific vascular geometries harboring IAs, were carried out using the one-way fluid-solid interaction solution strategy implemented in solids4foam, an extension of OpenFOAM®, in which the blood flow is solved and applied as the driving force of the wall motion. We found that different wall morphology models yield smaller absolute differences in the mechanical response than different hyperelastic laws. Furthermore, the stretch levels of IAs walls were more sensitive to the hyperelastic and material constants than the stress. These findings could be used to guide modeling decisions on IA simulations, since the computational behavior of each law was different, for example, with the Yeoh law being the fastest to converge.

On the Single-Point Calculation of Stress-Strain Data under Large Deformations with Stress and Mixed Control

Stress-strain data with a given constitutive model of material can be calculated directly at a single material point. In this work, we propose a framework to perform single-point calculations under large deformations with stress and mixed control, to test and validate sophisticated constitutive models for materials. Inspired by Galerkin-FFT methods, a well-defined mask projector is used for stress and mixed control, and the derived nonlinear equations are solved in Newton iterations with Krylov solvers, simplifying implementation. One application example of the single-point calculator in developing sophisticated models for anisotropic single crystal rate-independent elastoplasticity is given, illustrating that the proposed algorithm can simulate asymmetrical deformation responses under uni-axial loading. Another example for artificial neural network models of the particle reinforced composite is also given, demonstrating that the commonly used machine learning or deep learning modeling frameworks can be directly incorporated into the proposed calculator. The central difference approximation of the tangent is validated so that derivative-free calculations for black-box constitutive models are possible. The proposed Python-coded single-point calculator is shown to be capable of quickly building, testing, and validating constitutive models with sophisticated or implicit structures, thus boosting the development of novel constitutive models for advanced solid materials.

Inverse problems in blood flow modeling: A review

Mathematical and computational modeling of the cardiovascular system is increasingly providing non-invasive alternatives to traditional invasive clinical procedures. Moreover, it has the potential for generating additional diagnostic markers. In blood flow computations, the personalization of spatially distributed (i.e., 3D) models is a key step which relies on the formulation and numerical solution of inverse problems using clinical data, typically medical images for measuring both anatomy and function of the vasculature. In the last years, the development and application of inverse methods has rapidly expanded most likely due to the increased availability of data in clinical centers and the growing interest of modelers and clinicians in collaborating. Therefore, this work aims to provide a wide and comparative overview of literature within the last decade. We review the current state of the art of inverse problems in blood flows, focusing on studies considering fully dimensional fluid and fluid-solid models. The relevant physical models and hemodynamic measurement techniques are introduced, followed by a survey of mathematical data assimilation approaches used to solve different kinds of inverse problems, namely state and parameter estimation. An exhaustive discussion of the literature of the last decade is presented, structured by types of problems, models and available data.

Fluid structure interaction versus rigid-wall approach in the study of the symptomatic stenosed carotid artery: Importance of wall compliance and resilience of loose connective tissue

The purpose of this paper is to demonstrate the importance of a compliant wall approach in modeling of non-Newtonian and non-physiological blood flows. A case study of a stenosed and symptomatic carotid bifurcation was considered to show the influence of the wall-resilience assumption on the flow parameters obtained with numerical simulations. Patient-specific data concerning the geometry and flow conditions were collected and used to carry out two-way coupled fluid structure interaction simulations of the pulsatile blood flow through carotid artery. The wall compliance was considered separately as related to the wall-elasticity and as associated with the reaction of the loose connective tissue surrounding the carotid bifurcation. The obtained hemodynamic parameters were compared to those which were found in rigid-wall simulations. The difference between the results obtained for rigid-wall and compliant-wall approaches for the peak-systolic area-averaged wall shear stress achieved 35%, whereas the difference between the time-averaged local vorticity and shear strain reached, respectively, 42% and 43%. The influence of the highly resilient wall on the monitored hemodynamic parameters was significant even if time-averaged values are compared, which suggests that these metrics are considerably overestimated if the wall compliance is not considered. Moreover, the findings show that the mechanical response of the loose connective tissue cannot be neglected in blood flow simulations. Additionally, this study indicates that stiffening of the arterial wall due to atherosclerosis significantly rises hemodynamic parameters. This explains the therapeutic benefits of surgical removal of plaque lesions formed in the carotid bifurcation (endarterectomy).

Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm

Abdominal aortic aneurysms are serious and difficult to detect, conditions can be deadly if they rupture. In this study, the heat transfer and flow physics of Abdominal Aortic Aneurysm (AAA) were discussed and associated with cardiac cycle to illustrate the cardiac thermal pulse (CTP) of AAA. A CTP and infrared thermography (IRT) evaluation-based on AAA and abdomen skin surface detection method was proposed, respectively. Infrared thermography (IRT) is a promising imaging technique that may detect AAA quicker and cheaper than other imaging techniques (as biomarker). From CFD rigid-wall and FSI Analysis, the transient bioheat transfer effect resulted in a distinct thermal signature (circular thermal elevation) on the temperature profile of midriff skin surface, at both regular body temperature and supine position, under normal clinical temperature. However, it is important to note that thermography is not a perfect technology, and it does have some limitations, such as lack of clinical trials. There is still work to be done to improve this imaging technique and make it a more viable and accurate method for detecting abdominal aortic aneurysms. However, thermography is currently one of the most convenient technologies in this field, and it has the potential to detect abdominal aortic aneurysms earlier than other techniques. CTP, on the other hand, was used to examine the thermal physics of AAA. In CFD rigid-wall Analysis, AAA had a CTP that only responded to systolic phase at regular body temperature. In contrast, a healthy abdominal aorta displayed a CTP that responded to the full cardiac cycle, including diastolic phase at all simulated cases. Besides, the findings from FSI Analysis suggest the influence of numerical simulation techniques on the prediction of thermal physics behaviours of AAA and abdominal skin surface. Lastly, this study correlated the relationship between natural convective heat transfer coefficient with AAA and provided reference for potential clinical diagnostic using IRT in clinical implications.

Probing the possibility of lesion formation/progression in vicinity of a primary atherosclerotic plaque: A fluid-solid interaction study and angiographic evidences

It is shown that certain locations in the arterial tree, such as coronary and cerebral arteries, are more prevalent to plaque formation. Endothelial activation and consequent plaque development are attributed to local hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT), and stress phase angle. After a certain level of plaque progression, these hemodynamic parameters are disturbed before and after the plaque. In the current study, it is hypothesized that the vicinity of a primary lesion is susceptible for further degeneration and second plaque formation. A fluid-solid interaction (FSI) model of the coronary artery with different levels of asymmetric constriction, is simulated and the trend of hemodynamic parameters were studied in both of the plaque side (PS) and the opposite wall (facing the plaque [PF]). Also, a novel factor is introduced that can identify the high-risk regions associated with WSS oscillations to negative values. Our results indicate that when more than half of the artery is constricted, the downstream of the plaque is highly exposed to endothelial pathogenesis the PS, such that negative WSS, and as well, critical values of OSI and RRT, that is, -1.2 Pa, 0.42 and 6.5 s, respectively arise in this region. PS endothelial cells in this region exposed to the highest risk of atherosclerosis based on the proposed index (3 out of 3). As well, three cases of angiographic images are provided that confirms existence of secondary lesion close to the primary one as predicted by our computational simulations.

Multiscale simulations suggest a protective role of neo-adventitia in abdominal aortic aneurysms

Abdominal aortic aneurysms (AAAs) are a dangerous cardiovascular disease, the pathogenesis of which is not yet fully understood. In the present work a recent mechanopathological theory, which correlates AAA progression with microstructural and mechanical alterations in the tissue, is investigated using multiscale models. The goal is to combine these changes, within the framework of mechanobiology, with possible mechanical cues that are sensed by vascular cells along the AAA pathogenesis. Particular attention is paid to the formation of a 'neo-adventitia' on the abluminal side of the aortic wall, which is characterized by a highly random (isotropic) distribution of collagen fibers. Macro- and micro-scale results suggest that the formation of an AAA, as expected, perturbs the micromechanical state of the aortic tissue and triggers a growth and remodeling (G&R) reaction by mechanosensing cells such as fibroblasts. This G&R then leads to the formation of a thick neo-adventitia that appears to bring the micromechanical state of the tissue closer to the original homeostatic level. In this context, this new layer could act like a protective sheath, similar to the tunica adventitia in healthy aortas. This potential 'attempt at healing' by vascular cells would have important implications on the stability of the AAA wall and thus on the risk of rupture. STATEMENT OF SIGNIFICANCE: Current clinical criteria for risk assessment in AAAs are still empirical, as the causes and mechanisms of the disease are not yet fully understood. The strength of the arterial tissue is closely related to its microstructure, which in turn is remodeled by mechanosensing cells in the course of the disease. In this study, multiscale simulations show a possible connection between mechanical cues at the microscopic level and collagen G&R in AAA tissue. It should be emphasized that these micromechanical cues cannot be visualized in vivo. Therefore, the results presented here will help to advance our current understanding of the disease and motivate future experimental studies, with important implications for AAA risk assessment.

Biomechanical constitutive modeling of the gastrointestinal tissues: a systematic review

The gastrointestinal (GI) tract is a continuous channel through the body that consists of the esophagus, the stomach, the small intestine, the large intestine, and the rectum. Its primary functions are to move the intake of food for digestion before storing and ultimately expulsion of feces. The mechanical behavior of GI tissues thus plays a crucial role for GI function in health and disease. The mechanical properties are characterized by a biomechanical constitutive model, which is a mathematical representation of the relation between load and deformation in a tissue. Hence, validated biomechanical constitutive models are essential to characterize and simulate the mechanical behavior of the GI tract. Here, a systematic review of these constitutive models is provided. This review is limited to studies where a model of the strain energy function is proposed to characterize the stress-strain relation of a GI tissue. Several needs are identified for more advanced modeling including: 1) Microstructural models that provide actual structure-function relations; 2) Validation of coupled electro-mechanical models accounting for active muscle contractions; 3) Human data to develop and validate models. The findings from this review provide guidelines for using existing constitutive models as well as perspective and directions for future studies.

Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases

Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.

Soft Tissue Hybrid Model for Real-Time Simulations

In this article, a recent formulation for real-time simulation is developed combining the strain energy density of the Spring Mass Model (SMM) with the equivalent representation of the Strain Energy Density Function (SEDF). The resulting Equivalent Energy Spring Model (EESM) is expected to provide information in real-time about the mechanical response of soft tissue when subjected to uniaxial deformations. The proposed model represents a variation of the SMM and can be used to predict the mechanical behavior of biological tissues not only during loading but also during unloading deformation states. To assess the accuracy achieved by the EESM, experimental data was collected from liver porcine samples via uniaxial loading and unloading tensile tests. Validation of the model through numerical predictions achieved a refresh rate of 31 fps (31.49 ms of computation time for each frame), achieving a coefficient of determination R² from 93.23% to 99.94% when compared to experimental data. The proposed hybrid formulation to characterize soft tissue mechanical behavior is fast enough for real-time simulation and captures the soft material nonlinear virgin and stress-softened effects with high accuracy.

Patient-specific computational modeling of endovascular aneurysm repair: State of the art and future directions

Endovascular aortic repair (EVAR) has become the preferred intervention option for aortic aneurysms and dissections. This is because EVAR is much less invasive than the alternative open surgery repair. While in-hospital mortality rates are smaller for EVAR than open repair (1%-2% vs. 3%-5%), the early benefits of EVAR are lost after 3 years due to larger rates of complications in the EVAR group. Clinicians follow instructions for use (IFU) when possible, but are left with personal experience on how to best proceed and what choices to make with respect to stent-graft (SG) model choice, sizing, procedural options, and their implications on long-term outcomes. Computational modeling of SG deployment in EVAR and tissue remodeling after intervention offers an alternative way of testing SG designs in silico, in a personalized way before intervention, to ultimately select the strategies leading to better outcomes. Further, computational modeling can be used in the optimal design of SGs in cases of complex geometries. In this review, we address some of the difficulties and successes associated with computational modeling of EVAR procedures. There is still work to be done in all areas of EVAR in silico modeling, including model validation, before models can be applied in the clinic, but much progress has already been made. Critical to clinical implementation are current efforts focusing on developing fast algorithms that can achieve (near) real-time solutions, as well as ways of dealing with inherent uncertainties related to patient aortic wall degradation on an individualized basis. We are optimistic that EVAR modeling in the clinic will soon become a reality to help clinicians optimize EVAR interventions and ultimately reduce EVAR-associated complications.

Membrane curvature and connective fiber alignment in guinea pig round window membrane

The round window membrane (RWM) covers an opening between the perilymph fluid-filled inner ear space and the air-filled middle ear space. As the only non-osseous barrier between these two spaces, the RWM is an ideal candidate for aspiration of perilymph for diagnostics purposes and delivery of medication for treatment of inner ear disorders. Routine access across the RWM requires the development of new surgical tools whose design can only be optimized with a thorough understanding of the RWM's structure and properties. The RWM possesses a layer of collagen and elastic fibers so characterization of the distribution and orientation of these fibers is essential. Confocal and two-photon microscopy were conducted on intact RWMs in a guinea pig model to characterize the distribution of collagen and elastic fibers. The fibers were imaged via second-harmonic-generation, autofluorescence, and Rhodamine B staining. Quantitative analyses of both fiber orientation and geometrical properties of the RWM uncovered a significant correlation between mean fiber orientations and directions of zero curvature in some portions of the RWM, with an even more significant correlation between the mean fiber orientations and linear distance along the RWM in a direction approximately parallel to the cochlear axis. The measured mean fiber directions and dispersions can be incorporated into a generalized structure tensor for use in the development of continuum anisotropic mechanical constitutive models that in turn will enable optimization of surgical tools to access the cochlea. STATEMENT OF SIGNIFICANCE: The Round Window Membrane (RWM) is the only non-osseous barrier separating the middle and inner ear spaces, and thus is an ideal portal for medical access to the cochlea. An understanding of RWM structure and mechanical response is necessary to optimize the design of surgical tools for this purpose. The RWM geometry and the connective fiber orientation and dispersion are measured via confocal and 2-photon microscopy. A region of the RWM geometry is characterized as a hyperbolic paraboloid and another region as a tapered parabolic cylinder. Predominant fiber directions correlate well with directions of zero curvature in the hyperbolic paraboloid region. Overall fiber directions correlate well with position along a line approximately parallel to the central axis of the cochlea's spiral.

Uncertainty quantification in subject-specific estimation of local vessel mechanical properties

Quantitative estimation of local mechanical properties remains critically important in the ongoing effort to elucidate how blood vessels establish, maintain, or lose mechanical homeostasis. Recent advances based on panoramic digital image correlation (pDIC) have made high-fidelity 3D reconstructions of small-animal (e.g., murine) vessels possible when imaged in a variety of quasi-statically loaded configurations. While we have previously developed and validated inverse modeling approaches to translate pDIC-measured surface deformations into biomechanical metrics of interest, our workflow did not heretofore include a methodology to quantify uncertainties associated with local point estimates of mechanical properties. This limitation has compromised our ability to infer biomechanical properties on a subject-specific basis, such as whether stiffness differs significantly between multiple material locations on the same vessel or whether stiffness differs significantly between multiple vessels at a corresponding material location. In the present study, we have integrated a novel uncertainty quantification and propagation pipeline within our inverse modeling approach, relying on empirical and analytic Bayesian techniques. To demonstrate the approach, we present illustrative results for the ascending thoracic aorta from three mouse models, quantifying uncertainties in constitutive model parameters as well as circumferential and axial tangent stiffness. Our extended workflow not only allows parameter uncertainties to be systematically reported, but also facilitates both subject-specific and group-level statistical analyses of the mechanics of the vessel wall.

Nonlinear Free and Forced Vibrations of a Hyperelastic Micro/Nanobeam Considering Strain Stiffening Effect

In recent years, the static and dynamic response of micro/nanobeams made of hyperelasticity materials received great attention. In the majority of studies in this area, the strain-stiffing effect that plays a major role in many hyperelastic materials has not been investigated deeply. Moreover, the influence of the size effect and large rotation for such a beam that is important for the large deformation was not addressed. This paper attempts to explore the free and forced vibrations of a micro/nanobeam made of a hyperelastic material incorporating strain-stiffening, size effect, and moderate rotation. The beam is modelled based on the Euler–Bernoulli beam theory, and strains are obtained via an extended von Kármán theory. Boundary conditions and governing equations are derived by way of Hamilton’s principle. The multiple scales method is applied to obtain the frequency response equation, and Hamilton’s technique is utilized to obtain the free undamped nonlinear frequency. The influence of important system parameters such as the stiffening parameter, damping coefficient, length of the beam, length-scale parameter, and forcing amplitude on the frequency response, force response, and nonlinear frequency is analyzed. Results show that the hyperelastic microbeam shows a nonlinear hardening behavior, which this type of nonlinearity gets stronger by increasing the strain-stiffening effect. Conversely, as the strain-stiffening effect is decreased, the nonlinear frequency is decreased accordingly. The evidence from this study suggests that incorporating strain-stiffening in hyperelastic beams could improve their vibrational performance. The model proposed in this paper is mathematically simple and can be utilized for other kinds of micro/nanobeams with different boundary conditions.

Simulating the time evolving geometry, mechanical properties, and fibrous structure of bioprosthetic heart valve leaflets under cyclic loading

Currently, the most common replacement heart valve design is the 'bioprosthetic' heart valve (BHV), which has important advantages in that it does not require permanent anti-coagulation therapy, operates noiselessly, and has blood flow characteristics similar to the native valve. BHVs are typically fabricated from glutaraldehyde-crosslinked pericardial xenograft tissue biomaterials (XTBs) attached to a rigid, semi-flexible, or fully collapsible stent in the case of the increasingly popular transcutaneous aortic valve replacement (TAVR). While current TAVR assessments are positive, clinical results to date are generally limited to <2 years. Since TAVR leaflets are constructed using thinner XTBs, their mechanical demands are substantially greater than surgical BHV due to the increased stresses during in vivo operation, potentially resulting in decreased durability. Given the functional complexity of heart valve operation, in-silico predictive simulations clearly have potential to greatly improve the TAVR development process. As such simulations must start with accurate material models, we have developed a novel time-evolving constitutive model for pericardial xenograft tissue biomaterials (XTB) utilized in BHV (doi: 10.1016/j.jmbbm.2017.07.013). This model was able to simulate the observed tissue plasticity effects that occur in approximately in the first two years of in vivo function (50 million cycles). In the present work, we implemented this model into a complete simulation pipeline to predict the BHV time evolving geometry to 50 million cycles. The pipeline was implemented within an isogeometric finite element formulation that directly integrated our established BHV NURBS-based geometry (doi: 10.1007/s00466-015-1166-x). Simulations of successive loading cycles indicated continual changes in leaflet shape, as indicated by spatially varying increases in leaflet curvature. While the simulation model assumed an initial uniform fiber orientation distribution, anisotropic regional changes in leaflet tissue plastic strain induced a complex changes in regional fiber orientation. We have previously noted in our time-evolving constitutive model that the increases in collagen fiber recruitment with cyclic loading placed an upper bound on plastic strain levels. This effect was manifested by restricting further changes in leaflet geometry past 50 million cycles. Such phenomena was accurately captured in the valve-level simulations due to the use of a tissue-level structural-based modeling approach. Changes in basic leaflet dimensions agreed well with extant experimental studies. As a whole, the results of the present study indicate the complexity of BHV responses to cyclic loading, including changes in leaflet shape and internal fibrous structure. It should be noted that the later effect also influences changes in local mechanical behavior (i.e. changes in leaflet anisotropic tissue stress-strain relationship) due to internal fibrous structure resulting from plastic strains. Such mechanism-based simulations can help pave the way towards the application of sophisticated simulation technologies in the development of replacement heart valve technology.

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.

Microstructural and mechanical characterization of the layers of human descending thoracic aortas

The mechanical properties of human aortas are linked to the layered tissue and its microstructure at different length scales. Each layer has specific mechanical and structural properties. While the ground substance and the elastin play an important role in tissue stiffness at small strain, collagen fibers carry most of the load at larger strains, which corresponds to the physiological conditions of the aorta at maximum pulsatile blood pressure. In fact, collagen fibers are crimped in the unloaded state. Collagen fibers show different orientation distributions when they are observed in a plane that is tangent to the aortic wall (in-plane section) or along a direction orthogonal to it (out-of-plane section). This was systematically investigated using large images (2500 × 2500 µm) with high resolution obtained by second harmonic generation (SHG) in order to homogenize tissue heterogeneity after a convergence analysis, which is a main goal of the study. In addition, collagen fibers show lateral interactions due to entanglements and the presence of transverse elastin fibers, observed on varying length scales using atomic force microscopy and a three-dimensional rendering obtained by stacking a sequence of SHG and two-photon fluorescence images; this is another important contribution. Human descending thoracic aortas from 13 heartbeat donors aged 28 to 66 years were examined. Uniaxial tensile tests were carried out on the longitudinal and circumferential strips of the aortic wall and the three separated layers (intima, media and adventitia). A structurally-motivated material model with (i) a term to describe the combined response of ground substance and elastin and (ii) terms to consider four families of collagen fibers with different directions was applied. The exclusion of compressed fibers was implemented in the fitting process of the experimental data, which was optimized by a genetic algorithm. The results show that a single fiber family with directional and dispersion parameters measured from SHG images can describe the mechanical response of all 39 layers (3 layers for each of the 13 aortas) with very good accuracy when a second (auxiliary) family of aligned fibers is introduced in the orthogonal direction to account for lateral fiber interaction. Indeed, all observed distributions of collagen directions can be accurately fitted by a single bivariate von Mises distribution. Statistical analysis of in-plane and out-of-plane dispersion of fiber orientations reveals structural differences between the three layers and a change of collagen dispersion parameters with age. STATEMENT OF SIGNIFICANCE: The stiffness of healthy young aortas is adjusted so that a diameter expansion of about 10 % is possible during the heartbeat. This creates the Windkessel effect, which smooths out the pulsating nature of blood flow and benefits organ perfusion. The specific elastic properties of the aorta that are required to achieve this effect are related to the microstructure of the aortic tissue at different length scales. An increase in the aortic stiffness, in addition to reducing cyclic expansion and worsening perfusion, is a risk factor for clinical hypertension. The present study relates the microstructure of healthy human aortas to the mechanical response and examines the changes in microstructural parameters with age, which is a key factor in increasing stiffness.

Numerical Models Can Assist Choice of an Aortic Phantom for In Vitro Testing

(1) Background: The realization of appropriate aortic replicas for in vitro experiments requires a suitable choice of both the material and geometry. The matching between the grade of details of the geometry and the mechanical response of the materials is an open issue that deserves attention. (2) Methods: To explore this issue, we performed a series of Fluid–Structure Interaction simulations, which compared the dynamics of three aortic models. Specifically, we reproduced a patient-specific geometry with a wall of biological tissue or silicone, and a parametric geometry based on in vivo data made in silicone. The biological tissue and the silicone were modeled with a fiber-oriented anisotropic and isotropic hyperelastic model, respectively. (3) Results: Clearly, both the aorta’s geometry and its constitutive material contribute to the determination of the aortic arch deformation; specifically, the parametric aorta exhibits a strain field similar to the patient-specific model with biological tissue. On the contrary, the local geometry affects the flow velocity distribution quite a lot, although it plays a minor role in the helicity along the arch. (4) Conclusions: The use of a patient-specific prototype in silicone does not a priori ensure a satisfactory reproducibility of the real aorta dynamics. Furthermore, the present simulations suggest that the realization of a simplified replica with the same compliance of the real aorta is able to mimic the overall behavior of the vessel.

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.

Biomechanical characterization of the passive response of the thoracic aorta in chronic hypoxic newborn lambs using an evolutionary strategy

The present study involves experiments and modelling aimed at characterizing the passive structural mechanical behavior of the chronic hypoxic lamb thoracic aorta, whose gestation, birth and postnatal period were carried at high altitude (3600 masl). To this end, the mechanical response was studied via tensile and pressurization tests. The tensile and pressurization tests measurements were used simultaneously to calibrate the material parameters of the Gasser-Holzapfel-Ogden (GHO) hyperelasctic anisotropic constitutive model through an analytical-numerical optimization procedure solved with an evolutionary strategy that guarantees a stable response of the model. The model and procedure of calibration adequately adjust to the material behavior in a wide deformation range with an appropriate physical description. The results of this study predict the mechanical response of the lamb thoracic aorta under generalized loading states like those that can occur in physiological conditions and/or in systemic arterial hypertension. Finally, the novel use of the evolutionary strategy, together with the set of experiments and tools used in this study, provide a robust alternative to validate biomechanical characterizations.

Incorporating fiber recruitment in hyperelastic modeling of vascular tissues by means of kinematic average

We present a framework for considering the gradual recruitment of collagen fibers in hyperelastic constitutive modeling. An effective stretch, which is a response variable representing the true stretch at the tissue-scale, is introduced. Properties of the effective stretch are discussed in detail. The effective stretch and strain invariants derived from it are used in selected hyperelastic constitutive models to describe the tissue response. This construction is investigated in conjunction with Holzapfel-Gasser-Ogden family strain energy functions. The ensuing models are validated against a large body of uniaxial and bi-axial stress-strain response data from human aortic aneurysm tissues. Both the descriptive and the predictive capabilities are examined. The former is evaluated by the quality of constitutive fitting, and the latter is assessed using finite element simulation. The models significantly improve the quality of fitting, and reproduce the experiment displacement, stress, and strain distributions with high fidelity in the finite element simulation.

Reduced biomechanical models for precision-cut lung-slice stretching experiments

Precision-cut lung-slices (PCLS), in which viable airways embedded within lung parenchyma are stretched or induced to contract, are a widely used ex vivo assay to investigate bronchoconstriction and, more recently, mechanical activation of pro-remodelling cytokines in asthmatic airways. We develop a nonlinear fibre-reinforced biomechanical model accounting for smooth muscle contraction and extracellular matrix strain-stiffening. Through numerical simulation, we describe the stresses and contractile responses of an airway within a PCLS of finite thickness, exposing the importance of smooth muscle contraction on the local stress state within the airway. We then consider two simplifying limits of the model (a membrane representation and an asymptotic reduction in the thin-PCLS-limit), that permit analytical progress. Comparison against numerical solution of the full problem shows that the asymptotic reduction successfully captures the key elements of the full model behaviour. The more tractable reduced model that we develop is suitable to be employed in investigations to elucidate the time-dependent feedback mechanisms linking airway mechanics and cytokine activation in asthma.

Ex Vivo Models to Decipher the Molecular Mechanisms of Genetic Notch Cardiovascular Disorders

Notch is an evolutionary, conserved, cell-cell signaling pathway that is central to several biological processes, from tissue morphogenesis to homeostasis. It is therefore not surprising that several genetic mutations of Notch components cause inherited human diseases, especially cardiovascular disorders. Despite numerous efforts, current in vivo models are still insufficient to unravel the underlying mechanisms of these pathologies, hindering the development of utmost needed medical therapies. In this perspective review, we discuss the limitations of current murine models and outline how the combination of microphysiological systems (MPSs) and targeted computational models can lead to breakthroughs in this field. In particular, while MPSs enable the experimentation on human cells in controlled and physiological environments, in silico models can provide a versatile tool to translate the in vitro findings to the more complex in vivo setting. As a showcase example, we focus on Notch-related cardiovascular diseases, such as Alagille syndrome, Adams-Oliver syndrome, and cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL). Impact statement In this review, a comprehensive overview of the limitations of current in vivo models of genetic Notch cardiovascular diseases is provided, followed by a discussion over the potential of microphysiological systems and computational models in overcoming these limitations and in potentiating drug testing and modeling of these pathologies.

An integrated set-up for ex vivo characterisation of biaxial murine artery biomechanics under pulsatile conditions

Ex vivo characterisation of arterial biomechanics enables detailed discrimination of the various cellular and extracellular contributions to arterial stiffness. However, ex vivo biomechanical studies are commonly performed under quasi-static conditions, whereas dynamic biomechanical behaviour (as relevant in vivo) may differ substantially. Hence, we aim to (1) develop an integrated set-up for quasi-static and dynamic biaxial biomechanical testing, (2) quantify set-up reproducibility, and (3) illustrate the differences in measured arterial stiffness between quasi-static and dynamic conditions. Twenty-two mouse carotid arteries were mounted between glass micropipettes and kept fully vasodilated. While recording pressure, axial force (F), and inner diameter, arteries were exposed to (1) quasi-static pressure inflation from 0 to 200 mmHg; (2) 300 bpm dynamic pressure inflation (peaking at 80/120/160 mmHg); and (3) axial stretch (λ_z) variation at constant pressures of 10/60/100/140/200 mmHg. Measurements were performed in duplicate. Single-point pulse wave velocities (PWV; Bramwell-Hill) and axial stiffness coefficients (c_ax = dF/dλ_z) were calculated at the in vivo value of λ_z. Within-subject coefficients of variation were ~ 20%. Dynamic PWVs were consistently higher than quasi-static PWVs (p < 0.001); c_ax increased with increasing pressure. We demonstrated the feasibility of ex vivo biomechanical characterisation of biaxially-loaded murine carotid arteries under pulsatile conditions, and quantified reproducibility allowing for well-powered future study design.

Synthetic photoplethysmography (PPG) of the radial artery through parallelized Monte Carlo and its correlation to body mass index (BMI)

Cardiovascular disease is one of the leading causes of death in the United States and obesity significantly increases the risk of cardiovascular disease. The measurement of blood pressure (BP) is critical in monitoring and managing cardiovascular disease hence new wearable devices are being developed to make BP more accessible to physicians and patients. Several wearables utilize photoplethysmography from the wrist vasculature to derive BP assessment although many of these devices are still at the experimental stage. With the ultimate goal of supporting instrument development, we have developed a model of the photoplethysmographic waveform derived from the radial artery at the volar surface of the wrist. To do so we have utilized the relation between vessel biomechanics through Finite Element Method and Monte Carlo light transport model. The model shows similar features to that seen in PPG waveform captured using an off the shelf device. We observe the influence of body mass index on the PPG signal. A degradation the PPG signal of up to 40% in AC to DC signal ratio was thus observed.

Numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation

Background

The effects of arterial wall compliance on blood flow have been revealed using fluid-structure interaction in last decades. However, microcirculation is not considered in previous researches. In fact, microcirculation plays a key role in regulating blood flow. Therefore, it is very necessary to involve microcirculation in arterial hemodynamics.

Objective

The main purpose of the present study is to investigate how wall compliance affects the flow characteristics and to establish the comparisons of these flow variables with rigid wall when microcirculation is considered.

Methods

We present numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation. A novel outlet boundary condition is employed to prescribe microcirculation in an idealised model.

Results

The novel finding in this work is that wall compliance under the consideration of microcirculation leads to the increase of wall shear stress in contrast to rigid wall, contrary to the traditional result that wall compliance makes wall shear stress decrease when a constant or time dependent pressure is specified at an outlet.

Conclusions

This work provides the valuable study of hemodynamics under physiological and realistic boundary conditions and proves that wall compliance may have a positive impact on wall shear stress based on this model. This methodology in this paper could be used in real model simulations.

Functional Grading of a Transversely Isotropic Hyperelastic Model with Applications in Modeling Tricuspid and Mitral Valve Transition Regions

Surgical simulators and injury-prediction human models require a combination of representative tissue geometry and accurate tissue material properties to predict realistic tool-tissue interaction forces and injury mechanisms, respectively. While biological tissues have been individually characterized, the transition regions between tissues have received limited research attention, potentially resulting in inaccuracies within simulations. In this work, an approach to characterize the transition regions in transversely isotropic (TI) soft tissues using functionally graded material (FGM) modeling is presented. The effect of nonlinearities and multi-regime nature of the TI model on the functional grading process is discussed. The proposed approach has been implemented to characterize the transition regions in the leaflet (LL), chordae tendinae (CT) and the papillary muscle (PM) of porcine tricuspid valve (TV) and mitral valve (MV). The FGM model is informed using high resolution morphological measurements of the collagen fiber orientation and tissue composition in the transition regions, and deformation characteristics predicted by the FGM model are numerically validated to experimental data using X-ray diffraction imaging. The results indicate feasibility of using the FGM approach in modeling soft-tissue transitions and has implications in improving physical representation of tissue deformation throughout the body using a scalable version of the proposed approach.

A full-field 3D digital image correlation and modelling technique to characterise anterior cruciate ligament mechanics ex vivo

It is limiting to use conventional methods when characterising material properties of complex biological tissues with inhomogeneous and anisotropic structure, such as the anterior cruciate ligament (ACL) in the knee joint. This study aims to develop and utilise a three-dimensional digital image correlation method (3D DIC) for the purpose of determining material properties of femur-ACL-tibia complex across the surface without any contact between the tissue and the loading equipment. A full-field (360° view) 3D DIC test setup consisting of six digital single-lens reflex cameras was developed and ACL specimens from skeletally mature dog knee joints were tested. The six cameras were arranged into three pairs and the cameras within each pair were positioned with 25° in between to obtain the desired stereovision output. The test setup was calibrated twice: first to obtain the intrinsic and extrinsic parameters within camera pairs, and second to align the 3D surfaces from each camera pair in order to generate the full view of the ACLs. Using the undeformed 3D surfaces of the ligaments, ACL-specific finite element models were generated. Longitudinal deformation of ligaments under tensile loads obtained from the 3D DIC, and this was analysed to serve as input for the inverse finite element analysis. As a result, hyperelastic coefficients from the first-order Ogden model that characterise ACL behaviour were determined with a marginal error of <1.5%. This test setup and methodology provides a means to accurately determine inhomogeneous and anisotropic material properties of ACL. The methodology described in this study could be adopted to investigate other biological and cultured tissues with complex structure. STATEMENT OF SIGNIFICANCE: Determining the material properties of soft tissues with complex anatomical structure, such as the anterior cruciate ligament (ACL), is important to better understand their contribution to musculoskeletal biomechanics. Current conventional methods for characterising material properties of the ACL are often limited to a contact measurement approach, however an improved understanding of the mechanics of this complex tissue is vital in terms of preventing injury and developing novel therapies. This article reports the development and utilisation of non-contact optical methodology involving full-field three-dimensional digital image correlation and finite element analysis to accurately investigate material properties of the ACL, in a controlled environment. This technique reduces inaccuracies due to specimen clamping and more importantly considers the inhomogeneous nature of the examined tissue.

Tissue Ingrowth Markedly Reduces Mechanical Anisotropy and Stiffness in Fibre Direction of Highly Aligned Electrospun Polyurethane Scaffolds

PURPOSE

The lack of long-term patency of synthetic vascular grafts currently available on the market has directed research towards improving the performance of small diameter grafts. Improved radial compliance matching and tissue ingrowth into the graft scaffold are amongst the main goals for an ideal vascular graft.

METHODS

Biostable polyurethane scaffolds were manufactured by electrospinning and implanted in subcutaneous and circulatory positions in the rat for 7, 14 and 28 days. Scaffold morphology, tissue ingrowth, and mechanical properties of the scaffolds were assessed before implantation and after retrieval.

RESULTS

Tissue ingrowth after 24 days was 96.5 ± 2.3% in the subcutaneous implants and 77.8 ± 5.4% in the circulatory implants. Over the 24 days implantation, the elastic modulus at 12% strain decreased by 59% in direction of the fibre alignment whereas it increased by 1379% transverse to the fibre alignment of the highly aligned scaffold of the subcutaneous implants. The lesser aligned scaffold of the circulatory graft implants exhibited an increase of the elastic modulus at 12% strain by 77% in circumferential direction.

CONCLUSION

Based on the observations, it is proposed that the mechanism underlying the softening of the highly aligned scaffold in the predominant fibre direction is associated with scaffold compaction and local displacement of fibres by the newly formed tissue. The stiffening of the scaffold, observed transverse to highly aligned fibres and for more a random fibre distribution, represents the actual mechanical contribution of the tissue that developed in the scaffold.

A refined dynamic finite-strain shell theory for incompressible hyperelastic materials: equations and two-dimensional shell virtual work principle

Based on previous work for the static problem, in this paper, we first derive one form of dynamic finite-strain shell equations for incompressible hyperelastic materials that involve three shell constitutive relations. In order to single out the bending effect as well as to reduce the number of shell constitutive relations, a further refinement is performed, which leads to a refined dynamic finite-strain shell theory with only two shell constitutive relations (deducible from the given three-dimensional (3D) strain energy function) and some new insights are also deduced. By using the weak formulation of the shell equations and the variation of the 3D Lagrange functional, boundary conditions and the two-dimensional shell virtual work principle are derived. As a benchmark problem, we consider the extension and inflation of an arterial segment. The good agreement between the asymptotic solution based on the shell equations and that from the 3D exact one gives verification of the former. The refined shell theory is also applied to study the plane-strain vibrations of a pressurized artery, and the effects of the axial pre-stretch, pressure and fibre angle on the vibration frequencies are investigated in detail.

A Data-Driven Learning Method for Constitutive Modeling: Application to Vascular Hyperelastic Soft Tissues

We address the problem of machine learning of constitutive laws when large experimental deviations are present. This is particularly important in soft living tissue modeling, for instance, where large patient-dependent data is found. We focus on two aspects that complicate the problem, namely, the presence of an important dispersion in the experimental results and the need for a rigorous compliance to thermodynamic settings. To address these difficulties, we propose to use, respectively, Topological Data Analysis techniques and a regression over the so-called General Equation for the Nonequilibrium Reversible-Irreversible Coupling (GENERIC) formalism (M. Grmela and H. Ch. Oettinger, Dynamics and thermodynamics of complex fluids. I. Development of a general formalism. Phys. Rev. E 56, 6620, 1997). This allows us, on one hand, to unveil the true “shape” of the data and, on the other, to guarantee the fulfillment of basic principles such as the conservation of energy and the production of entropy as a consequence of viscous dissipation. Examples are provided over pseudo-experimental and experimental data that demonstrate the feasibility of the proposed approach.

A damage model for collagen fibres with an application to collagenous soft tissues

We propose a mechanical model to account for progressive damage in collagen fibres within fibrous soft tissues. The model has a similar basis to the pseudoelastic model that describes the Mullins effect in rubber but it also accounts for the effect of cross-links between collagen fibres. We show that the model is able to capture experimental data obtained from rat tail tendon fibres, and the combined effect of damage and collagen cross-links is illustrated for a simple shear test. The proposed three-dimensional framework allows a straightforward implementation in finite-element codes, which are needed to analyse more complex boundary-value problems for soft tissues under supra-physiological loading or tissues weakened by disease.

Towards microstructure-informed material models for human brain tissue

Emerging evidence suggests that the mechanical behavior of the brain plays a critical role in development, disease, and aging. Recent studies have begun to characterize the mechanical behavior of gray and white matter tissue and to identify sets of material models that best reproduce the stress-strain behavior of different brain regions. Yet, these models are mainly phenomenological in nature, their parameters often lack clear physical interpretation, and they fail to correlate the mechanical behavior to the underlying microstructural composition. Here we make a first attempt towards identifying general relations between microstructure and mechanics with the ultimate goal to develop microstructurally motivated constitutive equations for human brain tissue. Using histological staining, we analyze the microstructure of brain specimens from different anatomical regions, the cortex, basal ganglia, corona radiata, and corpus callosum, and identify the regional stiffness and viscosity under multiple loading conditions, simple shear, compression, and tension. Strikingly, our study reveals a negative correlation between cell count and stiffness, a positive correlation between myelin content and stiffness, and a negative correlation between proteoglycan content and stiffness. Additionally, our analysis shows a positive correlation between lipid and proteoglycan content and viscosity. We demonstrate how understanding the microstructural origin of the macroscopic behavior of the brain can help us design microstructure-informed material models for human brain tissue that inherently capture regional heterogeneities. This study represents an important step towards using brain tissue stiffness and viscosity as early diagnostic markers for clinical conditions including chronic traumatic encephalopathy, Alzheimer's and Parkinson's disease, or multiple sclerosis. STATEMENT OF SIGNIFICANCE: The complex and heterogeneous mechanical properties of brain tissue play a critical role for brain function. To understand and predict how brain tissue properties vary in space and time, it will be key to link the mechanical behavior to the underlying microstructural composition. Here we use histological staining to quantify area fractions of microstructural components of mechanically tested specimens and evaluate their individual contributions to the nonlinear macroscopic mechanical response. We further propose a microstructure-informed material model for human brain tissue that inherently captures regional heterogeneities. The current work provides unprecedented insights into the biomechanics of human brain tissue, which are highly relevant to develop refined computational models for brain tissue behavior or to advance neural tissue engineering.

Mechanical and structural changes in human thoracic aortas with age

Aortic mechanical and structural characteristics have profound effects on pathophysiology, but many aspects of physiologic stress-stretch state and intramural changes due to aging remain poorly understood in human tissues. While difficult to assess in vivo due to residual stresses and pre-stretch, physiologic stress-stretch characteristics can be calculated using experimentally-measured mechanical properties and constitutive modeling. Mechanical properties of 76 human descending thoracic aortas (TA) from 13 to 78-year-old donors (mean age 51±18 years) were measured using multi-ratio planar biaxial extension. Constitutive parameters were derived for aortas in 7 age groups, and the physiologic stress-stretch state was calculated. Intramural characteristics were quantified from histological images and related to aortic morphometry and mechanics. TA stiffness increased with age, and aortas became more nonlinear and anisotropic. Systolic and diastolic elastic energy available for pulsation decreased with age from 30 to 8 kPa and from 18 to 5 kPa, respectively. Cardiac cycle circumferential stretch dropped from 1.14 to 1.04, and circumferential and longitudinal physiologic stresses decreased with age from 90 to 72 kPa and from 90 to 17 kPa, respectively. Aortic wall thickness and radii increased with age, while the density of elastin in the tunica media decreased. The number of elastic lamellae and circumferential physiologic stress per lamellae unit remained constant with age at 102±10 and 0.85±0.04 kPa, respectively. Characterization of mechanical, physiological, and structural features in human aortas of different ages can help understand aortic pathology, inform the development of animal models that simulate human aging, and assist with designing devices for open and endovascular aortic repairs. STATEMENT OF SIGNIFICANCE: This manuscript describes mechanical and structural changes occurring in human thoracic aortas with age, and presents material parameters for 4 commonly used constitutive models. Presented data can help better understand aortic pathology, inform the development of animal models that simulate human aging, and assist with designing devices for open and endovascular aortic repairs.

Flow-Structure Interaction Simulations of the Aortic Heart Valve at Physiologic Conditions: The Role of Tissue Constitutive Model

The blood flow patterns in the region around the aortic valve depend on the geometry of the aorta and on the complex flow-structure interaction between the pulsatile flow and the valve leaflets. Consequently, the flow depends strongly on the constitutive properties of the tissue, which can be expected to vary between healthy and diseased heart valves or native and prosthetic valves. The main goal of this work is to qualitatively demonstrate that the choice of the constitutive model of the aortic valve is critical in analysis of heart hemodynamics. To accomplish that two different constitutive models were used in curvilinear immersed boundary-finite element-fluid-structure interaction (CURVIB-FE-FSI) method developed by Gilmanov et al. (2015, "A Numerical Approach for Simulating Fluid Structure Interaction of Flexible Thin Shells Undergoing Arbitrarily Large Deformations in Complex Domains," J. Comput. Phys., 300, pp. 814-843.) to simulate an aortic valve in an anatomic aorta at physiologic conditions. The two constitutive models are: (1) the Saint-Venant (StV) model and (2) the modified May-Newman&Yin (MNY) model. The MNY model is more general and includes nonlinear, anisotropic effects. It is appropriate to model the behavior of both prosthetic and biological tissue including native valves. Both models are employed to carry out FSI simulations of the same valve in the same aorta anatomy. The computed results reveal dramatic differences in both the vorticity dynamics in the aortic sinus and the wall shear-stress patterns on the aortic valve leaflets and underscore the importance of tissue constitutive models for clinically relevant simulations of aortic valves.

Arterial Stiffness: Different Metrics, Different Meanings

Findings from basic science and clinical studies agree that arterial stiffness is fundamental to both the mechanobiology and the biomechanics that dictate vascular health and disease. There is, therefore, an appropriately growing literature on arterial stiffness. Perusal of the literature reveals, however, that many different methods and metrics are used to quantify arterial stiffness, and reported values often differ by orders of magnitude and have different meanings. Without clear definitions and an understanding of possible inter-relations therein, it is increasingly difficult to integrate results from the literature to glean true understanding. In this paper, we briefly review methods that are used to infer values of arterial stiffness that span studies on isolated cells, excised intact vessels, and clinical assessments. We highlight similarities and differences and identify a single theoretical approach that can be used across scales and applications and thus could help to unify future results. We conclude by emphasizing the need to move toward a synthesis of many disparate reports, for only in this way will we be able to move from our current fragmented understanding to a true appreciation of how vascular cells maintain, remodel, or repair the arteries that are fundamental to cardiovascular properties and function.

An investigation of the glycosaminoglycan contribution to biaxial mechanical behaviours of porcine atrioventricular heart valve leaflets

The atrioventricular heart valve (AHV) leaflets have a complex microstructure composed of four distinct layers: atrialis, ventricularis, fibrosa and spongiosa. Specifically, the spongiosa layer is primarily proteoglycans and glycosaminoglycans (GAGs). Quantification of the GAGs' mechanical contribution to the overall leaflet function has been of recent focus for aortic valve leaflets, but this characterization has not been reported for the AHV leaflets. This study seeks to expand current GAG literature through novel mechanical characterizations of GAGs in AHV leaflets. For this characterization, mitral and tricuspid valve anterior leaflets (MVAL and TVAL, respectively) were: (i) tested by biaxial mechanical loading at varying loading ratios and by stress-relaxation procedures, (ii) enzymatically treated for removal of the GAGs and (iii) biaxially mechanically tested again under the same protocols as in step (i). Removal of the GAG contents from the leaflet was conducted using a 100 min enzyme treatment to achieve approximate 74.87% and 61.24% reductions of all GAGs from the MVAL and TVAL, respectively. Our main findings demonstrated that biaxial mechanical testing yielded a statistically significant difference in tissue extensibility after GAG removal and that stress-relaxation testing revealed a statistically significant smaller stress decay of the enzyme-treated tissue than untreated tissues. These novel findings illustrate the importance of GAGs in AHV leaflet behaviour, which can be employed to better inform heart valve therapeutics and computational models.

Determination of the Material Parameters in the Holzapfel-Gasser-Ogden Constitutive Model for Simulation of Age-Dependent Material Nonlinear Behavior for Aortic Wall Tissue under Uniaxial Tension

In this study, computational simulations and experiments were performed to investigate the mechanical behavior of the aorta wall because of the increasing occurrences of aorta-related diseases. The study focused on the deformation and strength of porcine and healthy human abdominal aortic tissues under uniaxial tensile loading. The experiments for the mechanical behavior of the arterial tissue were conducted using a uniaxial tensile test apparatus to validate the simulation results. In addition, the strength and stretching of the tissues in the abdominal aorta of a healthy human as a function of age were investigated based on the uniaxial tensile tests. Moreover, computational simulations using the ABAQUS finite element analysis program were conducted on the experimental scenarios based on age, and the Holzapfel–Gasser–Ogden (HGO) model was applied during the simulation. The material parameters and formulae to be used in the HGO model were proposed to identify the failure stress and stretch correlation with age.

Image-Guided Fluid-Structure Interaction Simulation of Transvalvular Hemodynamics: Quantifying the Effects of Varying Aortic Valve Leaflet Thickness

When flow-induced forces are altered at the blood vessel, maladaptive remodeling can occur. One reason such remodeling may occur has to do with the abnormal functioning of the aortic heart valve due to disease, calcification, injury, or an improperly-designed prosthetic valve, which restricts the opening of the valve leaflets and drastically alters the hemodynamics in the ascending aorta. While the specifics underlying the fundamental mechanisms leading to changes in heart valve function may differ from one cause to another, one common and important change is in leaflet stiffness and/or mass. Here, we examine the link between valve stiffness and mass and the hemodynamic environment in aorta by coupling magnetic resonance imaging (MRI) with high-resolution fluid–structure interaction (FSI) computational fluid dynamics to simulate blood flow in a patient-specific model. The thoracic aorta and a native aortic valve were re-constructed in the FSI model from the MRI data and used for the simulations. The effect of valve stiffness and mass is parametrically investigated by varying the thickness (h) of the leaflets (h = 0.6, 2, 4 mm). The FSI simulations were designed to investigate systematically progressively higher levels of valve stiffness by increasing valve thickness and quantifying hemodynamic parameters known to be linked to aortopathy and valve disease. The computed results reveal dramatic differences in all hemodynamic parameters: (1) the geometric orifice area (GOA), (2) the maximum velocity V max of the jet passing through the aortic orifice area, (3) the rate of energy dissipation E ˙ diss ( t ) , (4) the total loss of energy E diss , (5) the kinetic energy of the blood flow E kin ( t ) , and (6) the average magnitude of vorticity Ω a ( t ) , illustrating the change in hemodynamics that occur due to the presence of aortic valve stenosis.

A multiscale active structural model of the arterial wall accounting for smooth muscle dynamics

Arterial wall dynamics arise from the synergy of passive mechano-elastic properties of the vascular tissue and the active contractile behaviour of smooth muscle cells (SMCs) that form the media layer of vessels. We have developed a computational framework that incorporates both these components to account for vascular responses to mechanical and pharmacological stimuli. To validate the proposed framework and demonstrate its potential for testing hypotheses on the pathogenesis of vascular disease, we have employed a number of pharmacological probes that modulate the arterial wall contractile machinery by selectively inhibiting a range of intracellular signalling pathways. Experimental probes used on ring segments from the rabbit central ear artery are: phenylephrine, a selective α1-adrenergic receptor agonist that induces vasoconstriction; cyclopiazonic acid (CPA), a specific inhibitor of sarcoplasmic/endoplasmic reticulum Ca²⁺-ATPase; and ryanodine, a diterpenoid that modulates Ca²⁺ release from the sarcoplasmic reticulum. These interventions were able to delineate the role of membrane versus intracellular signalling, previously identified as main factors in smooth muscle contraction and the generation of vessel tone. Each SMC was modelled by a system of nonlinear differential equations that account for intracellular ionic signalling, and in particular Ca²⁺ dynamics. Cytosolic Ca²⁺ concentrations formed the catalytic input to a cross-bridge kinetics model. Contractile output from these cellular components forms the input to the finite-element model of the arterial rings under isometric conditions that reproduces the experimental conditions. The model does not account for the role of the endothelium, as the nitric oxide production was suppressed by the action of L-NAME, and also due to the absence of shear stress on the arterial ring, as the experimental set-up did not involve flow. Simulations generated by the integrated model closely matched experimental observations qualitatively, as well as quantitatively within a range of physiological parametric values. The model also illustrated how increased intercellular coupling led to smooth muscle coordination and the genesis of vascular tone.

A discrete fibre dispersion method for excluding fibres under compression in the modelling of fibrous tissues

Recently, micro-sphere-based methods derived from the angular integration approach have been used for excluding fibres under compression in the modelling of soft biological tissues. However, recent studies have revealed that many of the widely used numerical integration schemes over the unit sphere are inaccurate for large deformation problems even without excluding fibres under compression. Thus, in this study, we propose a discrete fibre dispersion model based on a systematic method for discretizing a unit hemisphere into a finite number of elementary areas, such as spherical triangles. Over each elementary area, we define a representative fibre direction and a discrete fibre density. Then, the strain energy of all the fibres distributed over each elementary area is approximated based on the deformation of the representative fibre direction weighted by the corresponding discrete fibre density. A summation of fibre contributions over all elementary areas then yields the resultant fibre strain energy. This treatment allows us to exclude fibres under compression in a discrete manner by evaluating the tension-compression status of the representative fibre directions only. We have implemented this model in a finite-element programme and illustrate it with three representative examples, including simple tension and simple shear of a unit cube, and non-homogeneous uniaxial extension of a rectangular strip. The results of all three examples are consistent and accurate compared with the previously developed continuous fibre dispersion model, and that is achieved with a substantial reduction of computational cost.