Abdominal aortic aneurysms (AAAs) are frequently characterized by the development of an intra-luminal thrombus (ILT), which is known to have multiple biochemical and biomechanical implications. Development of the ILT is not well understood, and shear–stress-triggered activation of platelets could be the first step in its evolution. Vortical structures (VSs) in the flow affect platelet dynamics, which motivated the present study of a possible correlation between VS and ILT formation in AAAs. VSs educed by the λ2-method using computational fluid dynamics simulations of the backward-facing step problem, normal aorta, fusiform AAA and saccular AAA were investigated. Patient-specific luminal geometries were reconstructed from computed tomography scans, and Newtonian and Carreau–Yasuda models were used to capture salient rheological features of blood flow. Particularly in complex flow domains, results depended on the constitutive model. VSs developed all along the normal aorta, showing that a clear correlation between VSs and high wall shear stress (WSS) existed, and that VSs started to break up during late systole. In contrast, in the fusiform AAA, large VSs developed at sites of tortuous geometry and high WSS, occupying the entire lumen, and lasting over the entire cardiac cycle. Downward motion of VSs in the AAA was in the range of a few centimetres per cardiac cycle, and with a VS burst at that location, the release (from VSs) of shear-stress-activated platelets and their deposition to the wall was within the lower part of the diseased artery, i.e. where the thickest ILT layer is typically observed. In the saccular AAA, only one VS was found near the healthy portion of the aorta, while in the aneurysmatic bulge, no VSs occurred. We present a fluid-dynamics-motivated mechanism for platelet activation, convection and deposition in AAAs that has the potential of improving our current understanding of the pathophysiology of fluid-driven ILT growth.

No abstract

The elastic strain energy potential for nonlinear fibre-reinforced materials is customarily obtained by superposition of the potentials of the matrix and of each family of fibres. Composites with statistically oriented fibres, such as biological tissues, can be seen as being reinforced by a continuous infinity of fibre families, the orientation of which can be represented by means of a probability density function defined on the unit sphere (i.e. the solid angle). In this case, the superposition procedure gives rise to an integral form of the elastic potential such that the deformation features in the integral, which therefore cannot be calculated a priori. As a consequence, an analytical use of this potential is impossible. In this paper, we implemented this integral form of the elastic potential into a numerical procedure that evaluates the potential, the stress and the elasticity tensor at each deformation step. The numerical integration over the unit sphere is performed by means of the method of spherical designs, in which the result of the integral is approximated by a suitable sum over a discrete subset of the unit sphere. As an example of application, we modelled the collagen fibre distribution in articular cartilage, and used it in simulating displacementcontrolled tests: the unconfined compression of a cylindrical sample and the contact problem in the hip joint.

No abstract

Biomechanical studies on abdominal aortic aneurysms (AAA) seek to provide for better decision criteria to undergo surgical intervention for AAA repair. More accurate results can be obtained by using appropriate material models for the tissues along with accurate geometric models and more realistic boundary conditions for the lesion. However, patient-specific AAA models are generated from gated medical images in which the artery is under pressure. Therefore, identification of the AAA zero pressure geometry would allow for a more realistic estimate of the aneurysmal wall mechanics. This study proposes a novel iterative algorithm to find the zero pressure geometry of patient-specific AAA models. The methodology allows considering the anisotropic hyperelastic behavior of the aortic wall, its thickness and accounts for the presence of the intraluminal thrombus. Results on 12 patient-specific AAA geometric models indicate that the procedure is computational tractable and efficient, and preserves the global volume of the model. In addition, a comparison of the peak wall stress computed with the zero pressure and CT-based geometries during systole indicates that computations using CT-based geometric models underestimate the peak wall stress by 59 ± 64 and 47 ± 64 kPa for the isotropic and anisotropic material models of the arterial wall, respectively.

No abstract

The paper deals with the modeling of balloon angioplasty by considering the balloon-induced overstretch of remnant non-diseased tissues in atherosclerotic arteries. A stenotic artery is modeled as a heterogenous structure composed of adventitia, media and a model plaque, and residual stresses are considered. The constitutive models are able to capture the anisotropic elastic tissue response in addition to the inelastic phenomena associated with tissue stretches beyond the physiological domain. The inelastic model describes the experimentally-observed changes of the wall during balloon inflation, i.e. non-recoverable deformation, and tissue weakening. The contact of the artery with a balloon catheter is simulated by a point-to-surface strategy. The states of deformations and stresses within the artery before, during and after balloon inflation are computed, compared and discussed. The 3D stress states at physiological loading conditions before and after balloon inflation differ significantly, and even compressive normal stresses may occur in the media after dilation.

Remodelling of soft biological tissue is characterized by interacting biochemical and biomechanical events, which change the tissue's microstructure, and, consequently, its macroscopic mechanical properties. Remodelling is a well-defined stage of the healing process, and aims at recovering or repairing the injured extracellular matrix. Like other physiological processes, remodelling is thought to be driven by homeostasis, i.e. it tends to re-establish the properties of the uninjured tissue. However, homeostasis may never be reached, such that remodelling may also appear as a continuous pathological transformation of diseased tissues during aneurysm expansion, for example. A simple constitutive model for soft biological tissues that regards remodelling as homeostatic-driven turnover is developed. Specifically, the recoverable effective tissue damage, whose rate is the sum of a mechanical damage rate and a healing rate, serves as a scalar internal thermodynamic variable. In order to integrate the biochemical and biomechanical aspects of remodelling, the healing rate is, on the one hand, driven by mechanical stimuli, but, on the other hand, subjected to simple metabolic constraints. The proposed model is formulated in accordance with continuum damage mechanics within an open-system thermodynamics framework. The numerical implementation in an in-house finite-element code is described, particularized for Ogden hyperelasticity. Numerical examples illustrate the basic constitutive characteristics of the model and demonstrate its potential in representing aspects of remodelling of soft tissues. Simulation results are verified for their plausibility, but also validated against reported experimental data.