Hyperelastic modelling of arterial layers with distributed collagen fibre orientations

Effects of dispersion of fiber orientation on the mechanical property of the arterial wall

Effects of dispersion of the fiber orientation on the mechanical property of the arterial wall in health and disease subjected to the combined internal pressure and axial loading are examined within the framework of the finite deformation hyper-elasticity theory. Taking into account the residual stress, a two layer thick-walled circular cylindrical tube model with the fiber-reinforced incompressible composite hyper-elastic material is employed. Stress-radius curves and stress distributions of the arterial wall are given in health and disease considering dispersion of the fiber orientation. With instability/bifurcation analysis, it is found that not only the uniform inflation of the arterial wall in health, but also the instability or bifurcation problem for the arterial wall in disease may be described by this model.

Freeze-dried heart valve scaffolds

In this study, structure and biomechanical properties of freeze-dried decellularized porcine pulmonary heart valves were investigated. Heart valves were dissected from porcine hearts. The tissues were decellularized and separated in three groups: (1) without lyoprotectant, (2) with 5% sucrose, and (3) with a mixture of 2.5% sucrose and 2.5% hydroxyl ethylene starch (HES), and then underwent freeze-drying. Freeze-drying in the absence of lyoprotectants caused an overall more disintegrated appearance of the histological architecture of the porcine valves, especially between the fibrosa and the ventricularis layers. Freeze-dried tissues with lyoprotectants have a looser network of collagen and elastic fibers with bigger pore sizes. Tissue freeze-dried in the absence of lyoprotecants had the largest pore sizes, whereas the tissue freeze-dried in the presence of protectants showed pores of intermediate sizes between the decellularized tissue and the unprotected freeze-dried samples. Tissue freeze-dried with sucrose alone displayed less porosity than tissue freeze-dried with the sucrose/HES mixture, whereas no significant differences in biomechanical properties were observed. Decellularization decreased the elastic modulus of artery tissue. The elastic modulus of freeze-dried tissue without protectants resembled that of decellularized tissue. The elastic modulus values of freeze-dried tissue stabilized by lyoprotectants were greater than those of decellularized tissue, but similar to those of native tissue.

An essential requirement for β1 integrin in the assembly of extracellular matrix proteins within the vascular wall

β1 integrin has been shown to contribute to vascular smooth muscle cell differentiation, adhesion and mechanosensation in vitro. Here we showed that deletion of β1 integrin at the onset of smooth muscle differentiation resulted in interrupted aortic arch, aneurysms and failure to assemble extracellular matrix proteins. These defects result in lethality prior to birth. Our data indicates that β1 integrin is not required for the acquisition, but it is essential for the maintenance of the smooth muscle cell phenotype, as levels of critical smooth muscle proteins are gradually reduced in mutant mice. Furthermore, while deposition of extracellular matrix was not affected, its structure was disrupted. Interestingly, defects in extracellular matrix and vascular wall assembly, were restricted to the aortic arch and its branches, compromising the brachiocephalic and carotid arteries and to the exclusion of the descending aorta. Additional analysis of β1 integrin in the pharyngeal arch smooth muscle progenitors was performed using wnt1Cre. Neural crest cells deleted for β1 integrin were able to migrate to the pharyngeal arches and associate with endothelial lined arteries; but exhibited vascular remodeling defects and early lethality. This work demonstrates that β1 integrin is dispensable for migration and initiation of the smooth muscle differentiation program, however, it is essential for remodeling of the pharyngeal arch arteries and for the assembly of the vessel wall of their derivatives. It further establishes a critical role of β1 integrin in the protection against aneurysms that is particularly confined to the ascending aorta and its branches.

Impact of residual stretch and remodeling on collagen engagement in healthy and pulmonary hypertensive calf pulmonary arteries at physiological pressures

Understanding the mechanical behavior of proximal pulmonary arteries (PAs) is crucial to evaluating pulmonary vascular function and right ventricular afterload. Early and current efforts focus on these arteries' histological changes, in vivo pressure-diameter behavior and mechanical properties under in vitro mechanical testing. However, the in vivo stretch and stress states remain poorly characterized. To further understand the mechanical behavior of the proximal PAs under physiological conditions, this study computed the residual stretch and the in vivo circumferential stretch state in the main pulmonary arteries in both control and hypertensive calves by using in vitro and in vivo artery geometry data, and modeled the impact of residual stretch and arterial remodeling on the in vivo circumferential stretch distribution and collagen engagement in the main pulmonary artery. We found that the in vivo circumferential stretch distribution in both groups was nonuniform across the vessel wall with the largest stretch at the outer wall, suggesting that collagen at the outer wall would engage first. It was also found that the circumferential stretch was more uniform in the hypertensive group, partially due to arterial remodeling that occurred during their hypoxic treatment, and that their onset of collagen engagement occurred at a higher pressure. It is concluded that the residual stretch and arterial remodeling have strong impact on the in vivo stretch state and the collagen engagement and thus the mechanical behavior of the main pulmonary artery in calves.

A forward incremental prestressing method with application to inverse parameter estimations and eye-specific simulations of posterior scleral shells

Numerical simulations or inverse numerical analyses of individual eyes or eye segments are often based on an eye-specific geometry obtained from in vivo medical images such as computed tomography (CT) scans or from in vitro 3D digitiser scans. These eye-specific geometries are usually measured while the eye is subjected to internal pressure. Due to the nonlinear stiffening of the collagen fibril network in the eye, numerical incorporation of the pre-existing stress/strain state may be essential for realistic eye-specific computational simulations. Existing prestressing methods either compute accurate predictions of the prestressed state or guarantee a unique solution. In this contribution, a forward incremental prestressing method is presented which unifies the advantages of the existing approaches by providing accurate and unique predictions of the pre-existing stress/strain state at the true measured geometry. The impact of prestressing is investigated on (i) the inverse constitutive parameter identification of a synthetic sclera inflation test and (ii) an eye-specific simulation that estimates the realistic mechanical response of a pre-loaded posterior monkey scleral shell. Evaluation of the pre-existing stress/strain state in the inverse analysis had a significant impact on the reproducibility of the constitutive parameters but may be estimated based on an approximative approach. The eye-specific simulation of one monkey eye shows that prestressing is required for accurate displacement and stress/strain predictions. The numerical results revealed an increasing error in displacement, strain and stress predictions with increasing pre-existing pressure load when the pre-stress/strain state is disregarded. Disregarding the prestress may lead to a significant underestimation of the strain/stress environment in the sclera and overestimation in the lamina cribrosa.

Regional annulus fibre orientations used as a tool for the calibration of lumbar intervertebral disc finite element models

The collagen network of the annulus fibrosus largely controls the functional biomechanics of the lumbar intervertebral discs (IVDs). Quantitative anatomical examinations have shown bundle orientation patterns, possibly coming from regional adaptations of the annulus mechanics. This study aimed to show that the regional differences in annulus mechanical behaviour could be reproduced by considering only fibre orientation changes. Using the finite element method, a lumbar annulus was modelled as a poro-hyperelastic material in which fibres were represented by a direction-dependent strain energy density term. Fibre orientations were calibrated to reproduce the annulus tensile behaviours measured for four different regions: posterior outer, anterior outer, posterior inner and anterior inner. The back-calculated fibre angles and regional patterns as well as the global disc behaviour were comparable with anatomical descriptions reported in the literature. It was concluded that annulus fibre variations might be an effective tool to calibrate lumbar spine IVD and segment models.

The effects of angiotensin II on the coupled microstructural and biomechanical response of C57BL/6 mouse aorta

RATIONALE

Abdominal aortic aneurysm (AAA) is a complex disease that leads to a localized dilation of the infrarenal aorta, the rupture of which is associated with significant morbidity and mortality. Animal models of AAA can be used to study how changes in the microstructural and biomechanical behavior of aortic tissues develop as disease progresses in these animals. We chose here to investigate the effect of angiotensin II (AngII) in C57BL/6 mice as a first step towards understanding how such changes occur in the established ApoE(-/-) AngII infused mouse model of AAA.

OBJECTIVE

The objective of this study was to utilize a recently developed device in our laboratory to determine how the microstructural and biomechanical properties of AngII-infused C57BL/6 wildtype mouse aorta change following 14 days of AngII infusion.

METHODS

C57BL/6 wildtype mice were infused with either saline or AngII for 14 day. Aortas were excised and tested using a device capable of simultaneously characterizing the biaxial mechanical response and load-dependent (unfixed, unfrozen) extracellular matrix organization of mouse aorta (using multiphoton microscopy). Peak strains and stiffness values were compared across experimental groups, and both datasets were fit to a Fung-type constitutive model. The mean mode and full width at half maximum (FWHM) of fiber histograms from two photon microscopy were quantified in order to assess the preferred fiber distribution and degree of fiber splay, respectively.

RESULTS

The axial stiffness of all mouse aorta was found to be an order of magnitude larger than the circumferential stiffness. The aortic diameter was found to be significantly increased for the AngII infused mice as compared to saline infused control (p=0.026). Aneurysm, defined as a percent increase in maximum diameter of 30% (defined with respect to saline control), was found in 3 of the 6 AngII infused mice. These three mice displayed adventitial collagen that lacked characteristic fiber crimp. The biomechanical response in the AngII infused mice showed significantly reduced circumferential compliance. We also noticed that the ability of the adventitial collagen fibers in AngII infused mice to disperse in reaction to circumferential loading was suppressed.

CONCLUSIONS

Collagen remodeling is present following 14 days of AngII infusion in C57BL/6 mice. Aneurysmal development occurred in 50% of our AngII infused mice, and these dilatations were accompanied with adventitial collagen remodeling and decreased circumferential compliance.

Determination of the layer-specific distributed collagen fibre orientations in human thoracic and abdominal aortas and common iliac arteries

The established method of polarized microscopy in combination with a universal stage is used to determine the layer-specific distributed collagen fibre orientations in 11 human non-atherosclerotic thoracic and abdominal aortas and common iliac arteries (63 ± 15.3 years, mean ± s.d.). A dispersion model is used to quantify over 37 000 recorded fibre angles from tissue samples. The study resulted in distinct fibre families, fibre directions, dispersion and thickness data for each layer and all vessels investigated. Two fibre families were present for the intima, media and adventitia in the aortas, with often a third and sometimes a fourth family in the intima in the respective axial and circumferential directions. In all aortas, the two families were almost symmetrically arranged with respect to the cylinder axis, closer to the axial direction in the adventitia, closer to the circumferential direction in the media and in between in the intima. The same trend was found for the intima and adventitia of the common iliac arteries; however, there was only one preferred fibre alignment present in the media. In all locations and layers, the observed fibre orientations were always in the tangential plane of the walls, with no radial components and very small dispersion through the wall thickness. A wider range of in-plane fibre orientations was present in the intima than in the media and adventitia. The mean total wall thickness for the aortas and the common iliac artery was 1.39 and 1.05 mm, respectively. For the aortas, a slight thickening of the intima and a thinning of the media in increasingly distal regions were observed. A clear intimal thickening was present distal to the branching of the celiac arteries. All data, except for the media of the common iliac arteries, showed two prominent collagen fibre families for all layers so that two-fibre family models seem most appropriate.

A novel approach to modelling and simulating the contact behaviour between a human hand model and a deformable object

A deeper understanding of biomechanical behaviour of human hands becomes fundamental for any human hand-operated activities. The integration of biomechanical knowledge of human hands into product design process starts to play an increasingly important role in developing an ergonomic product-to-user interface for products and systems requiring high level of comfortable and responsive interactions. Generation of such precise and dynamic models can provide scientific evaluation tools to support product and system development through simulation. This type of support is urgently required in many applications such as hand skill training for surgical operations, ergonomic study of a product or system developed and so forth. The aim of this work is to study the contact behaviour between the operators' hand and a hand-held tool or other similar contacts, by developing a novel and precise nonlinear 3D finite element model of the hand and by investigating the contact behaviour through simulation. The contact behaviour is externalised by solving the problem using the bi-potential method. The human body's biomechanical characteristics, such as hand deformity and structural behaviour, have been fully modelled by implementing anisotropic hyperelastic laws. A case study is given to illustrate the effectiveness of the approach.

Remodeling of the collagen fiber architecture due to compaction in small vessels under tissue engineered conditions

Mechanical loading protocols in tissue engineering (TE) aim to improve the deposition of a properly organized collagen fiber network. In addition to collagen remodeling, these conditioning protocols can result in tissue compaction. Tissue compaction is beneficial to tissue collagen alignment, yet it may lead to a loss of functionality of the TE construct due to changes in geometry after culture. Here, a mathematical model is presented to relate the changes in collagen architecture to the local compaction within a TE small blood vessel, assuming that under static conditions, compaction is the main factor responsible for collagen fiber organization. An existing structurally based model is extended to incorporate volumetric tissue compaction. Subsequently, the model is applied to describe the collagen architecture of TE constructs under either strain based or stress based stimulus functions. Our computations indicate that stress based simulations result in a helical collagen fiber distribution along the vessel wall. The helix pitch angle increases from a circumferential direction in the inner wall, over about 45 deg in the middle vessel layer, to a longitudinal direction in the outer wall. These results are consistent with experimental data from TE small diameter blood vessels. In addition, our results suggest a stress dependent remodeling of the collagen, suggesting that cell traction is responsible for collagen orientation. These findings may be of value to design improved mechanical conditioning protocols to optimize the collagen architecture in engineered tissues.

A passive strain-energy function for elastic and muscular arteries: correlation of material parameters with histological data

A plethora of phenomenological and structure-motivated constitutive models have thus far been used as pseudoelastic descriptors in arterial biomechanics, but their parameters have not been explicitly correlated with histology. This study associated biaxial histological data with strain-energy function (SEF) parameters derived from uniaxial tension data of arteries from different topographical sites (carotid artery vs. thoracic aorta vs. femoral artery). A two-term SEF fitted the passive stress-strain data of healthy porcine tissue, justified by the biphasic response characterizing elastin-rich tissues. Selection of a quadratic (orthotropic) over the neo-Hookean (isotropic) term was dictated by the directional dissimilarities in low-stress mechanical response, consistent with our histological data indicating orthotropic symmetry for unstressed elastin. Use of the exponential term was dictated by mechanical dissimilarities at high stresses and variations in unstressed collagen composition and orientation. Accurate fits were attained; topographical variations and anisotropy in material parameters were accounted by respective variations in histomorphometrical data.

Comparative analysis of the biaxial mechanical behavior of carotid wall tissue and biological and synthetic materials used for carotid patch angioplasty

Patch angioplasty is the most common technique used for the performance of carotid endarterectomy. A large number of patching materials are available for use while new materials are being continuously developed. Surprisingly little is known about the mechanical properties of these materials and how these properties compare with those of the carotid artery wall. Mismatch of the mechanical properties can produce mechanical and hemodynamic effects that may compromise the long-term patency of the endarterectomized arterial segment. The aim of this paper was to systematically evaluate and compare the biaxial mechanical behavior of the most commonly used patching materials. We compared PTFE (n  =  1), Dacron (n  =  2), bovine pericardium (n  =  10), autogenous greater saphenous vein (n  =  10), and autogenous external jugular vein (n  =  9) with the wall of the common carotid artery (n  =  18). All patching materials were found to be significantly stiffer than the carotid wall in both the longitudinal and circumferential directions. Synthetic patches demonstrated the most mismatch in stiffness values and vein patches the least mismatch in stiffness values compared to those of the native carotid artery. All biological materials, including the carotid artery, demonstrated substantial nonlinearity, anisotropy, and variability; however, the behavior of biological and biologically-derived patches was both qualitatively and quantitatively different from the behavior of the carotid wall. The majority of carotid arteries tested were stiffer in the circumferential direction, while the opposite anisotropy was observed for all types of vein patches and bovine pericardium. The rates of increase in the nonlinear stiffness over the physiological stress range were also different for the carotid and patching materials. Several carotid wall samples exhibited reverse anisotropy compared to the average behavior of the carotid tissue. A similar characteristic was observed for two of 19 vein patches. The obtained results quantify, for the first time, significant mechanical dissimilarity of the currently available patching materials and the carotid artery. The results can be used as guidance for designing more efficient patches with mechanical properties resembling those of the carotid wall. The presented systematic comparative mechanical analysis of the existing patching materials provides valuable information for patch selection in the daily practice of carotid surgery and can be used in future clinical studies comparing the efficacy of different patches in the performance of carotid endarterectomy.

The capsule's contribution to total hip construct stability--a finite element analysis

Instability is a significant concern in total hip arthroplasty (THA), particularly when there is structural compromise of the capsule due to pre-existing pathology or due to necessities of surgical approach. An experimentally grounded fiber-direction-based finite element model of the hip capsule was developed, and was integrated with an established three-dimensional model of impingement/dislocation. Model validity was established by close similarity to results from a cadaveric experiment in a servohydraulic hip simulator. Parametric computational runs explored effects of graded levels of capsule thickness, of regional detachment from the capsule's femoral or acetabular insertions, of surgical incisions of capsule substance, and of capsule defect repairs. Depending strongly upon the specific site, localized capsule defects caused varying degrees of construct stability compromise, with several specific situations involving over 60% decrement in dislocation resistance. Construct stability was returned substantially toward intact-capsule levels following well-conceived repairs, although the suture sites involved were often at substantial risk of failure. These parametric model results underscore the importance of retaining or robustly repairing capsular structures in THA, in order to maximize overall construct stability. 

Evolving biaxial mechanical properties of mouse carotid arteries in hypertension

Quantifying the time course of load-induced changes in arterial wall geometry, microstructure, and properties is fundamental to developing mathematical models of growth and remodeling. Arteries adapt to altered pressure and flow by modifying wall thickness, inner diameter, and axial length via marked cell and matrix turnover. To estimate particular biomaterial implications of such adaptations, we used a 4-fiber family constitutive relation to quantify passive biaxial mechanical behaviors of mouse carotid arteries 0 (control), 7-10, 10-14, or 35-56 days after an aortic arch banding surgery that increased pulse pressure and pulsatile flow in the right carotid artery. In vivo circumferential and axial stretches at mean arterial pressure were, for example, 11% and 26% lower, respectively, in hypertensive carotids 35-56 days after banding than in normotensive controls; this finding is consistent with observations that hypertension decreases distensibility. Interestingly, the strain energy W stored in the carotids at individual in vivo conditions was also less in hypertensive compared with normotensive carotids. For example, at 35-56 days after banding, W was 24%, 39%, and 47% of normal values at diastolic, mean, and systolic pressures, respectively. The energy stored during the cardiac cycle, W(sys)-W(dias), also tended to be less, but this reduction did not reach significance. When computed at normal in vivo values of biaxial stretch, however, W was well above normal for the hypertensive carotids. This net increase resulted from an overall increase in the collagen-related anisotropic contribution to W despite a decrease in the elastin-related isotropic contribution. The latter was consistent with observed decreases in the mass fraction of elastin.

Modeling the collagen fibril network of biological tissues as a nonlinearly elastic material using a continuous volume fraction distribution function

Despite distinct mechanical functions, biological soft tissues have a common microstructure in which a ground matrix is reinforced by a collagen fibril network. The microstructural properties of the collagen network contribute to continuum mechanical tissue properties that are strongly anisotropic with tensile-compressive asymmetry. In this study, a novel approach based on a continuous distribution of collagen fibril volume fractions is developed to model fibril reinforced soft tissues as a nonlinearly elastic and anisotropic material. Compared with other approaches that use a normalized number of fibrils for the definition of the distribution function, this representation is based on a distribution parameter (i.e. volume fraction) that is commonly measured experimentally while also incorporating pre-stress of the collagen fibril network in a tissue natural configuration. After motivating the form of the collagen strain energy function, examples are provided for two volume fraction distribution functions. Consequently, collagen second-Piola Kirchhoff stress and elasticity tensors are derived, first in general form and then specifically for a model that may be used for immature bovine articular cartilage. It is shown that the proposed strain energy is a convex function of the deformation gradient tensor and, thus, is suitable for the formation of a polyconvex tissue strain energy function.

Structural analysis of articular cartilage using multiphoton microscopy: input for biomechanical modeling

The 3-D morphology of chicken articular cartilage was quantified using multiphoton microscopy (MPM) for use in continuum-mechanical modeling. To motivate this morphological study we propose aspects of a new, 3-D finite strain constitutive model for articular cartilage focusing on the essential load-bearing morphology: an inhomogeneous, poro-(visco)elastic solid matrix reinforced by an anisotropic, (visco)elastic dispersed fiber fabric which is saturated by an incompressible fluid residing in strain-dependent pores. Samples of fresh chicken cartilage were sectioned in three orthogonal planes and imaged using MPM, specifically imaging the collagen fibers using second harmonic generation. Employing image analysis techniques based on Fourier analysis, we derived the principal directionality and dispersion of the collagen fiber fabric in the superficial layer. In the middle layer, objective thresholding techniques were used to extract the volume fraction occupied by extracellular collagen matrix. In conjunction with information available in the literature, or additional experimental testing, we show how this data can be used to derive a 3-D map of the initial solid volume fraction and Darcy permeability.

3D constitutive modeling of the biaxial mechanical response of intact and layer-dissected human carotid arteries

Human arteries with non-atherosclerotic intimal thickening consist of three distinct layers: adventitia, media and intima. From a series of axial extension and inflation tests on intact and layer-dissected human carotid arteries (adventitia and media-intima composite), a 3D structurally-based strain-energy function (SEF) is calibrated, and a set of five material parameters is identified which is not yet available in the literature. The zero-stress states of the artery tubes investigated are considered in the calibration process, and the related kinematics for the finite deformation of the individual continuum are described in detail. The SEF employed is capable of describing the different mechanical properties of the intact and layer-dissected tissue tubes (arterial segments) investigated at different pressure domains and axial stretches. The correlation coefficients and error measures determined indicate good correlation between the model and the experimental data for all tested tubes. Mean values of each individual material parameter provide a kind of 'master model' that characterizes the mean response of all mechanical data obtained from the human carotid arteries. The material parameters and the 3D constitutive model serve as a foundation for finite element simulations, and hence the analysis of more complex patient-specific boundary-value problems in the human carotid physiology and pathology.

Characterization of the anisotropic mechanical properties of excised human skin

The mechanical properties of skin are important for a number of applications including surgery, dermatology, impact biomechanics and forensic science. In this study, we have investigated the influence of location and orientation on the deformation characteristics of 56 samples of excised human skin. Uniaxial tensile tests were carried out at a strain rate of 0.012 s(-1) on skin from the back. Digital Image Correlation was used for 2D strain measurement and a histological examination of the dermis was also performed. The mean ultimate tensile strength (UTS) was 21.6±8.4 MPa, the mean failure strain 54%±17%, the mean initial slope 1.18±0.88 MPa, the mean elastic modulus 83.3±34.9 MPa and the mean strain energy was 3.6±1.6 MJ/m(3). A multivariate analysis of variance has shown that these mechanical properties of skin are dependent upon the orientation of the Langer lines (P<0.0001-P=0.046). The location of specimens on the back was also found to have a significant effect on the UTS (P=0.0002), the elastic modulus (P=0.001) and the strain energy (P=0.0052). The histological investigation concluded that there is a definite correlation between the orientation of the Langer lines and the preferred orientation of collagen fibres in the dermis (P<0.001). The data obtained in this study will provide essential information for those wishing to model the skin using a structural constitutive model.

Experimental and modeling study of collagen scaffolds with the effects of crosslinking and fiber alignment

Collagen type I scaffolds are commonly used due to its abundance, biocompatibility, and ubiquity. Most applications require the scaffolds to operate under mechanical stresses. Therefore understanding and being able to control the structural-functional integrity of collagen scaffolds becomes crucial. Using a combined experimental and modeling approach, we studied the structure and function of Type I collagen gel with the effects of spatial fiber alignment and crosslinking. Aligned collagen scaffolds were created through the flow of magnetic particles enmeshed in collagen fibrils to mimic the anisotropy seen in native tissue. Inter- and intra- molecular crosslinking was modified chemically with Genipin to further improve the stiffness of collagen scaffolds. The anisotropic mechanical properties of collagen scaffolds were characterized using a planar biaxial tensile tester and parallel plate rheometer. The tangent stiffness from biaxial tensile test is two to three orders of magnitude higher than the storage moduli from rheological measurements. The biphasic nature of collagen gel was discussed and used to explain the mechanical behavior of collagen scaffolds under different types of mechanical tests. An anisotropic hyperelastic constitutive model was used to capture the characteristics of the stress-strain behavior exhibited by collagen scaffolds.

The effect of sustained compression on oxygen metabolic transport in the intervertebral disc decreases with degenerative changes

Intervertebral disc metabolic transport is essential to the functional spine and provides the cells with the nutrients necessary to tissue maintenance. Disc degenerative changes alter the tissue mechanics, but interactions between mechanical loading and disc transport are still an open issue. A poromechanical finite element model of the human disc was coupled with oxygen and lactate transport models. Deformations and fluid flow were linked to transport predictions by including strain-dependent diffusion and advection. The two solute transport models were also coupled to account for cell metabolism. With this approach, the relevance of metabolic and mechano-transport couplings were assessed in the healthy disc under loading-recovery daily compression. Disc height, cell density and material degenerative changes were parametrically simulated to study their influence on the calculated solute concentrations. The effects of load frequency and amplitude were also studied in the healthy disc by considering short periods of cyclic compression. Results indicate that external loads influence the oxygen and lactate regional distributions within the disc when large volume changes modify diffusion distances and diffusivities, especially when healthy disc properties are simulated. Advection was negligible under both sustained and cyclic compression. Simulating degeneration, mechanical changes inhibited the mechanical effect on transport while disc height, fluid content, nucleus pressure and overall cell density reductions affected significantly transport predictions. For the healthy disc, nutrient concentration patterns depended mostly on the time of sustained compression and recovery. The relevant effect of cell density on the metabolic transport indicates the disturbance of cell number as a possible onset for disc degeneration via alteration of the metabolic balance. Results also suggest that healthy disc properties have a positive effect of loading on metabolic transport. Such relation, relevant to the maintenance of the tissue functional composition, would therefore link disc function with disc nutrition.

Identifying a minimal rheological configuration: a tool for effective and efficient constitutive modeling of soft tissues

We describe a modeling methodology intended as a preliminary step in the identification of appropriate constitutive frameworks for the time-dependent response of biological tissues. The modeling approach comprises a customizable rheological network of viscous and elastic elements governed by user-defined 1D constitutive relationships. The model parameters are identified by iterative nonlinear optimization, minimizing the error between experimental and model-predicted structural (load-displacement) tissue response under a specific mode of deformation. We demonstrate the use of this methodology by determining the minimal rheological arrangement, constitutive relationships, and model parameters for the structural response of various soft tissues, including ex vivo perfused porcine liver in indentation, ex vivo porcine brain cortical tissue in indentation, and ex vivo human cervical tissue in unconfined compression. Our results indicate that the identified rheological configurations provide good agreement with experimental data, including multiple constant strain rate load/unload tests and stress relaxation tests. Our experience suggests that the described modeling framework is an efficient tool for exploring a wide array of constitutive relationships and rheological arrangements, which can subsequently serve as a basis for 3D constitutive model development and finite-element implementations. The proposed approach can also be employed as a self-contained tool to obtain simplified 1D phenomenological models of the structural response of biological tissue to single-axis manipulations for applications in haptic technologies.

A mechano-acoustic indentor system for in vivo measurement of nonlinear elastic properties of soft tissue

OBJECTIVES

Soft tissue exhibits nonlinear stress-strain behavior under compression. Characterizing its nonlinear elasticity may aid detection, diagnosis, and treatment of soft tissue abnormality. The purposes of this study were to develop a rate-controlled Mechano-Acoustic Indentor System and a corresponding finite element optimization method to extract nonlinear elastic parameters of soft tissue and evaluate its test-retest reliability.

METHODS

An indentor system using a linear actuator to drive a force-sensitive probe with a tip-mounted ultrasound transducer was developed. Twenty independent sites at the upper lateral quadrant of the buttock from 11 asymptomatic subjects (7 men and 4 women from a chiropractic college) were indented at 6% per second for 3 sessions, each consisting of 5 trials. Tissue thickness, force at 25% deformation, and area under the load-deformation curve from 0% to 25% deformation were calculated. Optimized hyperelastic parameters of the soft tissue were calculated with a finite element model using a first-order Ogden material model. Load-deformation response on a standardized block was then simulated, and the corresponding area and force parameters were calculated. Between-trials repeatability and test-retest reliability of each parameter were evaluated using coefficients of variation and intraclass correlation coefficients, respectively.

RESULTS

Load-deformation responses were highly reproducible under repeated measurements. Coefficients of variation of tissue thickness, area under the load-deformation curve from 0% to 25% deformation, and force at 25% deformation averaged 0.51%, 2.31%, and 2.23%, respectively. Intraclass correlation coefficients ranged between 0.959 and 0.999, indicating excellent test-retest reliability.

CONCLUSIONS

The automated Mechano-Acoustic Indentor System and its corresponding optimization technique offers a viable technology to make in vivo measurement of the nonlinear elastic properties of soft tissue. This technology showed excellent between-trials repeatability and test-retest reliability with potential to quantify the effects of a wide variety of manual therapy techniques on the soft tissue elastic properties.

The fiber orientation in the coronary arterial wall at physiological loading evaluated with a two-fiber constitutive model

A patient-specific mechanical description of the coronary arterial wall is indispensable for individualized diagnosis and treatment of coronary artery disease. A way to determine the artery's mechanical properties is to fit the parameters of a constitutive model to patient-specific experimental data. Clinical data, however, essentially lack information about the stress-free geometry of an artery, which is necessary for constitutive modeling. In previous research, it has been shown that a way to circumvent this problem is to impose extra modeling constraints on the parameter estimation procedure. In this study, we propose a new modeling constraint concerning the in-situ fiber orientation (β (phys)). β (phys), which is a major contributor to the arterial stress-strain behavior, was determined for porcine and human coronary arteries using a mixed numerical-experimental method. The in-situ situation was mimicked using in-vitro experiments at a physiological axial pre-stretch, in which pressure-radius and pressure-axial force were measured. A single-layered, hyperelastic, thick-walled, two-fiber material model was accurately fitted to the experimental data, enabling the computation of stress, strain, and fiber orientation. β (phys) was found to be almost equal for all vessels measured (36.4 ± 0.3)°, which theoretically can be explained using netting analysis. In further research, this finding can be used as an extra modeling constraint in parameter estimation from clinical data.

A constrained von Mises distribution to describe fiber organization in thin soft tissues

The semi-circular von Mises distribution is widely used to describe the unimodal planar organization of fibers in thin soft tissues. However, it cannot accurately describe the isotropic subpopulation of fibers present in such tissues, and therefore an improved mathematical description is needed. We present a modified distribution, formed as a weighted mixture of the semi-circular uniform distribution and the semi-circular von Mises distribution. It is described by three parameters: β, which weights the contribution from each mixture component; k, the fiber concentration factor; and θ ( p ), the preferred fiber orientation. This distribution was used to fit data obtained by small-angle light scattering experiments from various thin soft tissues. Initial use showed that satisfactory fits of fiber distributions could be obtained (error generally < 1%), but at the cost of non-physically meaningful values of k and β. To address this issue, an empirical constraint between the parameters k and β was introduced, resulting in a constrained 2-parameter fiber distribution. Compared to the 3-parameter distribution, the constrained 2-parameter distribution fits experimental data well (error generally < 2%) and had the advantage of producing physically meaningful parameter values. In addition, the constrained 2-parameter approach was more robust to experimental noise. The constrained 2-parameter fiber distribution can replace the semi-circular von Mises distribution to describe unimodal planar organization of fibers in thin soft tissues. Inclusion of such a function in constitutive models for finite element simulations should provide better quantitative estimates of soft tissue biomechanics under normal and pathological conditions.

Pulmonary vascular stiffness: measurement, modeling, and implications in normal and hypertensive pulmonary circulations

This article introduces the concept of pulmonary vascular stiffness, discusses its increasingly recognized importance as a diagnostic marker in the evaluation of pulmonary vascular disease, and describes methods to measure and model it clinically, experimentally, and computationally. It begins with a description of systems-level methods to evaluate pulmonary vascular compliance and recent clinical efforts in applying such techniques to better predict patient outcomes in pulmonary arterial hypertension. It then progresses from the systems-level to the local level, discusses proposed methods by which upstream pulmonary vessels increase in stiffness, introduces concepts around vascular mechanics, and concludes by describing recent work incorporating advanced numerical methods to more thoroughly evaluate changes in local mechanical properties of pulmonary arteries.

A constituent-based model of age-related changes in conduit arteries

In the present report, a constituent-based theoretical model of age-related changes in geometry and mechanical properties of conduit arteries is proposed. The model was based on the premise that given the time course of the load on an artery and the accumulation of advanced glycation end-products in the arterial tissue, the initial geometric dimensions and properties of the arterial tissue can be predicted by a solution of a boundary value problem for the governing equations that follow from finite elasticity, structure-based constitutive modeling within the constrained mixture theory, continuum damage theory, and global growth approach for stress-induced structure-based remodeling. An illustrative example of the age-related changes in geometry, structure, composition, and mechanical properties of a human thoracic aorta is considered. Model predictions were in good qualitative agreement with available experimental data in the literature. Limitations and perspectives for refining the model are discussed.

Constitutive modeling of coronary arterial media--comparison of three model classes

Accurate modeling of arterial elasticity is imperative for predicting pulsatile blood flow and transport to the periphery, and for evaluating the mechanical microenvironment of the vessel wall. The goal of the present study is to compare a recently developed structural model of porcine left anterior descending artery media to two commonly used typical representatives of phenomenological and structure-motivated invariant-based models, in terms of the number of model parameters, model descriptive and predictive powers, and requisite different test protocols for reliable parameter estimation. The three models were compared against 3D data of radial inflation, axial extension, and twist tests. Also checked are the models predictive capabilities to response data not used for estimation, including both tests outside the range of estimation database, as well as protocols of a different nature. The results show that the descriptive estimation error (model fit to estimation database), measured by the sum of squared residuals (SSE) between full 3D data and model predictions, was about twice as low for the structural (4.58%) model compared to the other two (9.71 and 8.99% for the phenomenological and structure-motivated models, respectively). Similar SSE ratios were obtained for the predictive capabilities. Prediction SSE at high stretch based on estimation of two low stretches yielded an SSE value of 2.81% for the structural model, and 10.54% and 7.87% for the phenomenological and structure-motivated models, respectively. For the prediction of twist from inflation-extension data, SSE values for the torsional stiffness was 1.76% for the structural model and 39.62 and 2.77% for the phenomenological and structure-motivated models. The required number of model parameters for the structural model is four, whereas the phenomenological model requires six to nine and the structure-motivated has four parameters. These results suggest that modeling based on the tissue structural features improves model reliability in describing given data and in predicting the tissue general response.

An irreversible constitutive model for fibrous soft biological tissue: a 3-D microfiber approach with demonstrative application to abdominal aortic aneurysms

Understanding the failure and damage mechanisms of soft biological tissue is critical to a sensitive and specific characterization of tissue injury tolerance and its relation to biological responses. Despite increasing experimental and analytical efforts, failure-related irreversible effects of soft biological tissue are still poorly understood. There is still no clear definition of what "damage" of a soft biological material is, and conventional macroscopic indicators, as known from damage of engineering materials for example, may not identify the tissue's tolerance to injury appropriately. To account for the complex three-dimensional arrangement of collagen, a microfiber model approach is applied, where constitutive relations for collagen fibers are integrated over the unit sphere, which in turn defines the tissue's macroscopic properties. A collagen fiber is represented by a bundle of proteoglycan cross-linked collagen fibrils that undergoes irreversible deformations when exceeding its elastic tensile limit. The proposed constitutive model is able to predict strain stiffening at physiological strain levels and does not exhibit a clear macroscopic elastic limit, two typical features known from soft biological tissue testing. An elastic-predictor/plastic-corrector implementation of the model is followed and constitutive parameters are estimated from in vitro test data from a particular abdominal aortic aneurysm (AAA). Damage-based structural instabilities of the AAA under different inflation conditions are investigated, where the collagen orientation density has been estimated from its in vivo stress state.

On a phenomenological model for active smooth muscle contraction

This paper presents a three-dimensional phenomenological model for the description of smooth muscle activation. A strain energy function is proposed as sum of the strain energy stored in the passive tissue, consisting of elastin and collagen, and an active calcium-driven energy related to the chemical contraction of the smooth muscle cells. Further, the proposed model includes the dispersions of the orientations of smooth muscle cells and collagen. These dispersions, measured in experiments, can be directly inserted into the model. The approach is implemented into the framework of the finite element method. Consequently, beside a validation with experiments the modelling concept is used for a three-dimensional numerical study.

Mechanical characterisation of the human thoracic descending aorta: experiments and modelling

This work presents experiments and modelling aimed at characterising the passive mechanical behaviour of the human thoracic descending aorta. To this end, uniaxial tension and pressurisation tests on healthy samples corresponding to newborn, young and adult arteries are performed. Then, the tensile measurements are used to calibrate the material parameters of the Holzapfel constitutive model. This model is found to adequately adjust the material behaviour in a wide deformation range; in particular, it captures the progressive stiffness increase and the anisotropy due to the stretching of the collagen fibres. Finally, the assessment of these material parameters in the modelling of the pressurisation test is addressed. The implication of this study is the possibility to predict the mechanical response of the human thoracic descending aorta under generalised loading states like those that can occur in physiological conditions and/or in medical device applications.

A constrained mixture approach to mechano-sensing and force generation in contractile cells

Biological tissues are very particular types of materials that have the ability to change their structure, properties and chemistry in response to external cues. Contractile cells, i.e. fibroblasts, are key players of tissue adaptivity as they are capable of reorganizing their surrounding extra-cellular matrix (ECM) by contracting and generating mechanical forces. This contractile behavior is attributed to the development of a stress-fiber (SF) network within the cell's cytoskeleton, a process that is known to be highly dependent of the nature of the mechanical environment (such as ECM stiffness or the presence of stress and strain). To describe these processes in a consistent manner, the present paper introduces a mutiphasic formulation (fluid/solid/solute mixture) that accounts for four major elements of cell contraction: cytoskeleton, cytosol, SF and actin monomers, as well as their interactions. The model represents the cross-talks between mechanics and chemistry through various means: (a) a mechano-sensitive formation and dissociation of an anisotropic SF network described by mass exchange between actin monomer and polymers, (b) a bio-mechanical model for SF contraction that captures the well-known length-tension and velocity-tension relation for muscles cells and (c) a convection/diffusion description for the transport of fluid and monomers within the cell. Numerical investigations show that the multiphasic model is able to capture the dependency of cell contraction on the stiffness of the mechanical environment and accurately describes the development of an oriented SF network observed in contracting fibroblasts.

A microstructurally driven model for pulmonary artery tissue

A new constitutive model for elastic, proximal pulmonary artery tissue is presented here, called the total crimped fiber model. This model is based on the material and microstructural properties of the two main, passive, load-bearing components of the artery wall, elastin, and collagen. Elastin matrix proteins are modeled with an orthotropic neo-Hookean material. High stretch behavior is governed by an orthotropic crimped fiber material modeled as a planar sinusoidal linear elastic beam, which represents collagen fiber deformations. Collagen-dependent artery orthotropy is defined by a structure tensor representing the effective orientation distribution of collagen fiber bundles. Therefore, every parameter of the total crimped fiber model is correlated with either a physiologic structure or geometry or is a mechanically measured material property of the composite tissue. Further, by incorporating elastin orthotropy, this model better represents the mechanics of arterial tissue deformation. These advancements result in a microstructural total crimped fiber model of pulmonary artery tissue mechanics, which demonstrates good quality of fit and flexibility for modeling varied mechanical behaviors encountered in disease states.

Effects of intima stiffness and plaque morphology on peak cap stress

BACKGROUND

Rupture of the cap of a vulnerable plaque present in a coronary vessel may cause myocardial infarction and death. Cap rupture occurs when the peak cap stress exceeds the cap strength. The mechanical stress within a cap depends on the plaque morphology and the material characteristics of the plaque components. A parametric study was conducted to assess the effect of intima stiffness and plaque morphology on peak cap stress.

METHODS

Models with idealized geometries based on histology images of human coronary arteries were generated by varying geometric plaque features. The constructed multi-layer models contained adventitia, media, intima, and necrotic core sections. For adventitia and media layers, anisotropic hyperelastic material models were used. For necrotic core and intima sections, isotropic hyperelastic material models were employed. Three different intima stiffness values were used to cover the wide range reported in literature. According to the intima stiffness, the models were classified as stiff, intermediate and soft intima models. Finite element method was used to compute peak cap stress.

RESULTS

The intima stiffness was an essential determinant of cap stresses. The computed peak cap stresses for the soft intima models were much lower than for stiff and intermediate intima models. Intima stiffness also affected the influence of morphological parameters on cap stresses. For the stiff and intermediate intima models, the cap thickness and necrotic core thickness were the most important determinants of cap stresses. The peak cap stress increased three-fold when the cap thickness was reduced from 0.25 mm to 0.05 mm for both stiff and intermediate intima models. Doubling the thickness of the necrotic core elevated the peak cap stress by 60% for the stiff intima models and by 90% for the intermediate intima models. Two-fold increase in the intima thickness behind the necrotic core reduced the peak cap stress by approximately 25% for both intima models. For the soft intima models, cap thickness was less critical and changed the peak cap stress by 55%. However, the necrotic core thickness was more influential and changed the peak cap stress by 100%. The necrotic core angle emerged as a critical determinant of cap stresses where a larger angle lowered the cap stresses. Contrary to the stiff and intermediate intima models, a thicker intima behind the necrotic core increased the peak cap stress by approximately 25% for the soft intima models. Adventitia thickness and local media regression had limited effects for all three intima models.

CONCLUSIONS

For the stiff and intermediate intima models, the cap thickness was the most important morphological risk factor. However for soft intima models, the necrotic core thickness and necrotic core angle had a bigger impact on the peak cap stress. We therefore need to enhance our knowledge of intima material properties if we want to derive critical morphological plaque features for risk evaluation.

Wall stress of the cervical carotid artery in patients with carotid dissection: a case-control study

Spontaneous internal carotid artery (ICA) dissection (sICAD) results from an intimal tear located around the distal carotid sinus. The mechanisms causing the tear are unknown. This case-control study tested the hypotheses that head movements increase the wall stress in the cervical ICA and that the stress increase is greater in patients with sICAD than in controls. Five patients with unilateral, recanalized, left sICAD and five matched controls were investigated before and after maximal head rotation to the left and neck hyperextension after 45° head rotation to the left. The anatomy of the extracranial carotid arteries was assessed by magnetic resonance imaging and used to create finite element models of the right ICA. Wall stress increased after head movements. Increases above the 80th and 90th percentile were located at the intimal side of the artery wall from 7.4 mm below to 10 mm above the cranial edge of the carotid sinus, i.e., at the same location as histologically confirmed tears in patients with sICAD. Wall stress increase did not differ between patients and controls. The present findings suggest that wall stress increases at the intimal side of the artery wall surrounding the distal edge of the carotid bulb after head movements may be important for the development of carotid dissection. The lack of wall stress difference between the two groups indicates that the carotid arteries of patients with carotid dissection have either distinct functional or anatomical properties or endured unusually heavy wall stresses to initiate dissection.

A structural constitutive model considering angular dispersion and waviness of collagen fibres of rabbit facial veins

Background

Structural constitutive models of vascular wall integrate information on composition and structural arrangements of tissue. In blood vessels, collagen fibres are arranged in coiled and wavy bundles and the individual collagen fibres have a deviation from their mean orientation. A complete structural constitutive model for vascular wall should incorporate both waviness and orientational distribution of fibres. We have previously developed a model, for passive properties of vascular wall, which considers the waviness of collagen fibres. However, to our knowledge there is no structural model of vascular wall which integrates both these features.

Methods

In this study, we have suggested a structural strain energy function that incorporates not only the waviness but also the angular dispersion of fibres. We studied the effect of parameters related to the orientational distribution on macro-mechanical behaviour of tissue during inflation-extension tests. The model was further applied on experimental data from rabbit facial veins.

Results

Our parametric study showed that the model is less sensitive to the orientational dispersion when fibres are mainly oriented circumferentially. The macro-mechanical response is less sensitive to changes in the mean orientation when fibres are more dispersed. The model accurately fitted the experimental data of veins, while not improving the quality of the fit compared to the model without dispersion. Our results showed that the orientational dispersion of collagen fibres could be compensated by a less abrupt and shifted to higher strain collagen engagement pattern. This should be considered when the model is fitted to experimental data and model parameters are used to study structural modifications of collagen fibre network in physiology and disease.

Conclusions

The presented model incorporates structural features related to waviness and orientational distribution of collagen fibres and thus offers possibilities to better understand the relation between structure and function in the vascular wall. Also, the model can be used to further study mechanically induced collagen remodelling in vascular tissue in health and disease.

Experimental study and constitutive modelling of the passive mechanical properties of the porcine carotid artery and its relation to histological analysis: Implications in animal cardiovascular device trials

The present study focusses on the determination, comparison and constitutive modelling of the passive mechanical properties of the swine carotid artery over very long stretches in both proximal and distal regions. Special attention is paid to the histological and mechanical variations of these properties depending on the proximity to the heart. The results can have clinical relevance, especially in the research field of intravascular device design. Before the final clinical trials on humans, research in the vascular area is conducted on animal models, swine being the most common due to the similarities between the human and swine cardiovascular systems as well as the fact that the swine size is suitable for testing devices, in this case endovascular carotid systems. The design of devices usually involves numerical techniques, and an important feature is the appropriate modelling of the mechanical properties of the vessel. Fourteen carotid swine arteries were harvested just after sacrifice and cyclic uniaxial tension tests in longitudinal and circumferential directions were performed for distal and proximal samples. The stress-stretch curves obtained were fitted with a hyperelastic anisotropic model. Stress-free configuration states were also analyzed. Finally, human and swine samples were processed in a histological laboratory and images were used to quantify their microconstituents. The statistical analysis revealed significant differences between the mechanical behavior of proximal and distal locations in the circumferential but not in the longitudinal direction. Circumferential direction samples show clear differences both in residual stretches and tensile curves between the two locations, while the features of longitudinal specimens are independent of the axial position. The statistical analysis provides significant evidence of changes depending on the position of the sample, mainly in elastin and SMC quantification.

Experimentally validated microstructural 3D constitutive model of coronary arterial media

Accurate modeling of arterial response to physiological or pathological loads may shed light on the processes leading to initiation and progression of a number of vascular diseases and may serve as a tool for prediction and diagnosis. In this study, a microstructure based hyperelastic constitutive model is developed for passive media of porcine coronary arteries. The most general model contains 12 independent parameters representing the three-dimensional inner fibrous structure of the media and includes the effects of residual stresses and osmotic swelling. Parameter estimation and model validation were based on mechanical data of porcine left anterior descending (LAD) media under radial inflation, axial extension, and twist tests. The results show that a reduced four parameter model is sufficient to reliably predict the passive mechanical properties. These parameters represent the stiffness and the helical orientation of each lamellae fiber and the stiffness of the interlamellar struts interconnecting these lamellae. Other structural features, such as orientational distribution of helical fibers and anisotropy of the interlamellar network, as well as possible transmural distribution of structural features, were found to have little effect on the global media mechanical response. It is shown that the model provides good predictions of the LAD media twist response based on parameters estimated from only biaxial tests of inflation and extension. In addition, good predictive capabilities are demonstrated for the model behavior at high axial stretch ratio based on data of law stretches.