Structure and torsion in the normal and situs inversus totalis cardiac left ventricle. II. Modeling cardiac adaptation to mechanical load

Mathematical models provide a suitable platform to test hypotheses on the relation between local mechanical stimuli and responses to cardiac structure and geometry. In the present model study, we tested hypothesized mechanical stimuli and responses in cardiac adaptation to mechanical load on their ability to estimate a realistic myocardial structure of the normal and situs inversus totalis (SIT) left ventricle (LV). In a cylindrical model of the LV, 1) mass was adapted in response to myofiber strain at the beginning of ejection and to global contractility (average systolic pressure), 2) cavity volume was adapted in response to fiber strain during ejection, and 3) myofiber orientations were adapted in response to myofiber strain during ejection and local misalignment between neighboring tissue parts. The model was able to generate a realistic normal LV geometry and structure. In addition, the model was also able to simulate the instigating situation in the rare SIT LV with opposite torsion and transmural courses in myofiber direction between the apex and base [Delhaas et al. (6)]. These results substantiate the importance of mechanical load in the formation and maintenance of cardiac structure and geometry. Furthermore, in the model, adapted myocardial architecture was found to be insensitive to fiber misalignment in the transmural direction, i.e., myofiber strain during ejection was sufficient to generate a realistic transmural variation in myofiber orientation. In addition, the model estimates that, despite differences in structure, global pump work and the mass of the normal and SIT LV are similar.     

Control of whole heart geometry by intramyocardial mechano-feedback: a model study

Geometry of the heart adapts to mechanical load, imposed by pressures and volumes of the cavities. We regarded preservation of cardiac geometry as a homeostatic control system. The control loop was simulated by a chain of models, starting with geometry of the cardiac walls, sequentially simulating circulation hemodynamics, myofiber stress and strain in the walls, transfer of mechano-sensed signals to structural changes of the myocardium, and finalized by calculation of resulting changes in cardiac wall geometry. Instead of modeling detailed mechano-transductive pathways and their interconnections, we used principles of control theory to find optimal transfer functions, representing the overall biological responses to mechanical signals. As biological responses we regarded tissue mass, extent of contractile myocyte structure and extent of the extra-cellular matrix. Mechano-structural stimulus-response characteristics were considered to be the same for atrial and ventricular tissue. Simulation of adaptation to self-generated hemodynamic load rendered physiologic geometry of all cardiac cavities automatically. Adaptation of geometry to chronic hypertension and volume load appeared also physiologic. Different combinations of mechano-sensors satisfied the condition that control of geometry is stable. Thus, we expect that for various species, evolution may have selected different solutions for mechano-adaptation.     

Relative crystallinity of plant biomass: studies on assembly, adaptation and acclimation

Plant biomechanical design is central to cell shape, morphogenesis, reproductive performance and protection against environmental and mechanical stress. The cell wall forms the central load bearing support structure for plant design, yet a mechanistic understanding of its synthesis is incomplete. A key tool for studying the structure of cellulose polymorphs has been x-ray diffraction and fourier transform infrared spectroscopy (FTIR). Relative crystallinity index (RCI) is based on the x-ray diffraction characteristics of two signature peaks and we used this technique to probe plant assembly, adaptation and acclimation. Confocal microscopy was used to visualize the dynamics of cellulose synthase in transgenic Arabidopsis plants expressing a homozygous YFP::CESA6. Assembly: RCI values for stems and roots were indistinguishable but leaves had 23.4 and 21.6% lower RCI than stems and roots respectively. Adaptation: over 3-fold variability in RCI was apparent in leaves from 35 plant species spanning Ordovician to Cretaceous periods. Within this study, RCI correlated positively with leaf geometric constraints and with mass per unit area, suggestive of allometry. Acclimation: biomass crystallinity was found to decrease under conditions of thigmomorphogenesis in Arabidopsis. Further, in etiolated pea hypocotyls, RCI values also decreased compared to plants that were grown in light, consistent with alterations in FTIR cellulose fingerprint peaks and live cell imaging experiments revealing rapid orientation of the YFP::cellulose synthase-6 array in response to light. Herein, results and technical challenges associated with the structure of the cell wall that gives rise to sample crystallinity are presented and examined with respect to adaptation, acclimation and assembly in ecosystem-level processes.     

Validation of a patient-specific one-dimensional model of the systemic arterial tree

The aim of this study is to develop and validate a patient-specific distributed model of the systemic arterial tree. This model is built using geometric and hemodynamic data measured on a specific person and validated with noninvasive measurements of flow and pressure on the same person, providing thus a patient-specific model and validation. The systemic arterial tree geometry was obtained from MR angiographic measurements. A nonlinear viscoelastic constitutive law for the arterial wall is considered. Arterial wall distensibility is based on literature data and adapted to match the wave propagation velocity of the main arteries of the specific subject, which were estimated by pressure waves traveling time. The intimal shear stress is modeled using the Witzig-Womersley theory. Blood pressure is measured using applanation tonometry and flow rate using transcranial ultrasound and phase-contrast-MRI. The model predicts pressure and flow waveforms in good qualitative and quantitative agreement with the in vivo measurements, in terms of wave shape and specific wave features. Comparison with a generic one-dimensional model shows that the patient-specific model better predicts pressure and flow at specific arterial sites. These results obtained let us conclude that a patient-specific one-dimensional model of the arterial tree is able to predict well pressure and flow waveforms in the main systemic circulation, whereas this is not always the case for a generic one-dimensional model.     

A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse

We propose a multiscale mechanobiological model of bone remodelling to investigate the site-specific evolution of bone volume fraction across the midshaft of a femur. The model includes hormonal regulation and biochemical coupling of bone cell populations, the influence of the microstructure on bone turnover rate, and mechanical adaptation of the tissue. Both microscopic and tissue-scale stress/strain states of the tissue are calculated from macroscopic loads by a combination of beam theory and micromechanical homogenisation. This model is applied to simulate the spatio-temporal evolution of a human midshaft femur scan subjected to two deregulating circumstances: (i) osteoporosis and (ii) mechanical disuse. Both simulated deregulations led to endocortical bone loss, cortical wall thinning and expansion of the medullary cavity, in accordance with experimental findings. Our model suggests that these observations are attributable to a large extent to the influence of the microstructure on bone turnover rate. Mechanical adaptation is found to help preserve intracortical bone matrix near the periosteum. Moreover, it leads to non-uniform cortical wall thickness due to the asymmetry of macroscopic loads introduced by the bending moment. The effect of mechanical adaptation near the endosteum can be greatly affected by whether the mechanical stimulus includes stress concentration effects or not.     

Myocardial contraction in the reverse development of adaptation to short-term exposure to stress

Contractile function of an isolated right atrium was studied in short-term stressor effects-adapted male Wistar rats at different times after adaptation was completed. Adaptation to short-term stressor effects was shown to produce a restricted decrease of myocardial contractility shortly after adaptation was completed. At the 3d day another decrease of contractile function was noted. However, contractile function returned to the control level by the 5th day. At the same time adaptation completely prevented the impairment of myocardial contractility, induced by prolonged stress. The protective effect was seen immediately after adaptation, by days 3 and 5 after it, being reduced by day 10. It is assumed that at the 5th day after adaptation, the animals experience the post-adaptation state marked by disappearance of the negative adaptation effect and by remarkable protective effect of adaptation. As a result, all the characteristics of myocardial contractility evaluated after prolonged stress experienced by the animals at the 5th day following short-term stressor effects do not differ from control parameters.     

Effect of preliminary adaptation to short-term stress exposure on the resistance of the spontaneously contracting myocardium to a lipid peroxidation inducer

Contractile function of the isolated right atrium was studied in male Wistar rats adapted to short-term stressor exposures at varying times after adaptation was completed. Adaptation to short-term stressor exposures induced a limited decrease in myocardial contractility immediately after adaptation was over. On the 3d day an additional reduction in the characteristics of contractile function was still observed. However, by the 5th day the characteristics recovered to the control level. At the same time adaptation completely prevented the derangement of myocardial contractility, induced by exposure to a prolonged stress. That protective effect was observed as early as adaptation was completed, on days 3 and 5 after adaptation, and became lessened on the 10th day. It is assumed that on the 5th day after adaptation the animals are in a postadaptation state where the untoward effect of adaptation disappears whereas the protective effect is demonstrable to a full extent. As a result all the characteristics of myocardial contractility following a prolonged stress on the 5th day after completion of short-term stressor exposures differed in no way from the control parameters.     

Adaptation of cardiac structure by mechanical feedback in the environment of the cell: a model study.

In the cardiac left ventricle during systole mechanical load of the myocardial fibers is distributed uniformly. A mechanism is proposed by which control of mechanical load is distributed over many individual control units acting in the environment of the cell. The mechanics of the equatorial region of the left ventricle was modeled by a thick-walled cylinder composed of 6-1500 shells of myocardial fiber material. In each shell a separate control unit was simulated. The direction of the cells was varied so that systolic fiber shortening approached a given optimum of 15%. End-diastolic sarcomere length was maintained at 2.1 microns. Regional early-systolic stretch and global contractility stimulated growth of cellular mass. If systolic shortening was more than normal the passive extracellular matrix stretched. The design of the load-controlling mechanism was derived from biological experiments showing that cellular processes are sensitive to mechanical deformation. After simulating a few hundred adaptation cycles, the macroscopic anatomical arrangement of helical pathways of the myocardial fibers formed automatically. If pump load of the ventricle was changed, wall thickness and cavity volume adapted physiologically. We propose that the cardiac anatomy may be defined and maintained by a multitude of control units for mechanical load, each acting in the cellular environment. Interestingly, feedback through fiber stress is not a compelling condition for such control.     

A volumetric model for growth of arterial walls with arbitrary geometry and loads

Stress and deformation in arterial wall tissue are factors which may influence significantly its response and evolution. In this work we develop models based on nonlinear elasticity and finite element numerical solutions for the mechanical behaviour and the remodelling of the soft tissue of arteries, including anisotropy induced by the presence of collagen fibres. Remodelling and growth in particular constitute important features in order to interpret stenosis and atherosclerosis. The main object of this work is to model accurately volumetric growth, induced by fluid shear stress in the intima and local wall stress in arteries with patient-specific geometry and loads. The model is implemented in a nonlinear finite element setting which may be applied to realistic 3D geometries obtained from in vivo measurements. The capabilities of this method are demonstrated in several examples. Firstly a stenotic process on an idealised geometry induced by a non-uniform shear stress distribution is considered. Following the growth of a right coronary artery from an in vivo reconstructed geometry is presented. Finally, experimental measurements for growth under hypertension for rat carotid arteries are modelled.     

Fast Simulation of Mechanical Heterogeneity in the Electrically Asynchronous Heart Using the MultiPatch Module

Cardiac electrical asynchrony occurs as a result of cardiac pacing or conduction disorders such as left bundle-branch block (LBBB). Electrically asynchronous activation causes myocardial contraction heterogeneity that can be detrimental for cardiac function. Computational models provide a tool for understanding pathological consequences of dyssynchronous contraction. Simulations of mechanical dyssynchrony within the heart are typically performed using the finite element method, whose computational intensity may present an obstacle to clinical deployment of patient-specific models. We present an alternative based on the CircAdapt lumped-parameter model of the heart and circulatory system, called the MultiPatch module. Cardiac walls are subdivided into an arbitrary number of patches of homogeneous tissue. Tissue properties and activation time can differ between patches. All patches within a wall share a common wall tension and curvature. Consequently, spatial location within the wall is not required to calculate deformation in a patch. We test the hypothesis that activation time is more important than tissue location for determining mechanical deformation in asynchronous hearts. We perform simulations representing an experimental study of myocardial deformation induced by ventricular pacing, and a patient with LBBB and heart failure using endocardial recordings of electrical activation, wall volumes, and end-diastolic volumes. Direct comparison between simulated and experimental strain patterns shows both qualitative and quantitative agreement between model fibre strain and experimental circumferential strain in terms of shortening and rebound stretch during ejection. Local myofibre strain in the patient simulation shows qualitative agreement with circumferential strain patterns observed in the patient using tagged MRI. We conclude that the MultiPatch module produces realistic regional deformation patterns in the asynchronous heart and that activation time is more important than tissue location within a wall for determining myocardial deformation. The CircAdapt model is therefore capable of fast and realistic simulations of dyssynchronous myocardial deformation embedded within the closed-loop cardiovascular system.     

A 3-D model of tumor progression based on complex automata driven by particle dynamics

The dynamics of a growing tumor involving mechanical remodeling of healthy tissue and vasculature is neglected in most of the existing tumor models. This is due to the lack of efficient computational framework allowing for simulation of mechanical interactions. Meanwhile, just these interactions trigger critical changes in tumor growth dynamics and are responsible for its volumetric and directional progression. We describe here a novel 3-D model of tumor growth, which combines particle dynamics with cellular automata concept. The particles represent both tissue cells and fragments of the vascular network. They interact with their closest neighbors via semi-harmonic central forces simulating mechanical resistance of the cell walls. The particle dynamics is governed by both the Newtonian laws of motion and the cellular automata rules. These rules can represent cell life-cycle and other biological interactions involving smaller spatio-temporal scales. We show that our complex automata, particle based model can reproduce realistic 3-D dynamics of the entire system consisting of the tumor, normal tissue cells, blood vessels and blood flow. It can explain phenomena such as the inward cell motion in avascular tumor, stabilization of tumor growth by the external pressure, tumor vascularization due to the process of angiogenesis, trapping of healthy cells by invading tumor, and influence of external (boundary) conditions on the direction of tumor progression. We conclude that the particle model can serve as a general framework for designing advanced multiscale models of tumor dynamics and it is very competitive to the modeling approaches presented before.     

A lumped parameter model of mechanically mediated acute and long-term adaptations of contractility and geometry in lymphatics for characterization of lymphedema

A primary purpose of the lymphatic system is to transport fluid from peripheral tissues to the central venous system in order to maintain tissue–fluid balance. Failure to perform this task results in lymphedema marked by swelling of the affected limb as well as geometric remodeling and reduced contractility of the affected lymphatic vessels. The mechanical environment has been implicated in the regulation of lymphatic contractility, but it is unknown how changes in the mechanical environment are related to loss of contractile function and remodeling of the tissue. The purpose of this paper was to introduce a new theoretical framework for acute and long-term adaptations of lymphatic vessels to changes in mechanical loading. This theoretical framework combines a simplified version of a published lumped parameter model for lymphangion function and lymph transport, a published microstructurally motivated constitutive model for the active and passive mechanical behavior of isolated rat thoracic ducts, and novel models for acute mechanically mediated vasoreactive adaptations and long-term volumetric growth to simulate changes in muscle contractility and geometry of a single isolated rat thoracic duct in response to a sustained elevation in afterload. The illustrative examples highlight the potential role of the mechanical environment in the acute maintenance of contractility and long-term geometric remodeling, presumably aimed at meeting fluid flow demands while also maintaining mechanical homeostasis. Results demonstrate that contractility may adapt in response to shear stress to meet fluid flow demands and show that pressure-induced long-term geometric remodeling may attenuate these adaptations and reduce fluid flow. The modeling framework and illustrative simulations help suggest relevant experiments that are necessary to accurately quantify and predict the acute and long-term adaptations of lymphangions to altered mechanical loading.     

A versatile hybrid agent-based,particle and partial differential equations method to analyze vascular adaptation

Peripheral arterial occlusive disease is a chronic pathology affecting at least 8–12 million people in the USA, typically treated with a vein graft bypass or through the deployment of a stent in order to restore the physiological circulation. Failure of peripheral endovascular interventions occurs at the intersection of vascular biology, biomechanics, and clinical decision making. It is our hypothesis that the majority of endovascular treatment approaches share the same driving mechanisms and that a deep understanding of the adaptation process is pivotal in order to improve the current outcome of the procedure. The postsurgical adaptation of vein graft bypasses offers the perfect example of how the balance between intimal hyperplasia and wall remodeling determines the failure or the success of the intervention. Accordingly, this work presents a versatile computational model able to capture the feedback loop that describes the interaction between events at cellular/tissue level and mechano-environmental conditions. The work here presented is a generalization and an improvement of a previous work by our group of investigators, where an agent-based model uses a cellular automata principle on a fixed hexagonal grid to reproduce the leading events of the graft’s restenosis. The new hybrid model here presented allows a more realistic simulation both of the biological laws that drive the cellular behavior and of the active role of the membranes that separate the various layers of the vein. The novel feature is to use an immersed boundary implementation of a highly viscous flow to represent SMC motility and matrix reorganization in response to graft adaptation. Our implementation is modular, and this makes us able to choose the right compromise between closeness to the physiological reality and complexity of the model. The focus of this paper is to offer a new modular implementation that combines the best features of an agent-based model, continuum mechanics, and particle-tracking methods to cope with the multiscale nature of the adaptation phenomena. This hybrid method allows us to quickly test various hypotheses with a particular attention to cellular motility, a process that we demonstrated should be driven by mechanical homeostasis in order to maintain the right balance between cells and extracellular matrix in order to reproduce a distribution similar to histological experimental data from vein grafts.

     

Nonlinear anisotropic stress analysis of anatomically realistic cerebral aneurysms

BACKGROUND: Static deformation analysis and estimation of wall stress distribution of patient-specific cerebral aneurysms can provide useful insights into the disease process and rupture. METHOD OF APPROACH: The three-dimensional geometry of saccular cerebral aneurysms from 27 patients (18 unruptured and nine ruptured) was reconstructed based on computer tomography angiography images. The aneurysm wall tissue was modeled using a nonlinear, anisotropic, hyperelastic material model (Fung-type) which was incorporated in a user subroutine in ABAQUS. Effective material fiber orientations were assumed to align with principal surface curvatures. Static deformation of the aneurysm models were simulated assuming uniform wall thickness and internal pressure load of 100 mm Hg. RESULTS: The numerical analysis technique was validated by quantitative comparisons to results in the literature. For the patient-specific models, in-plane stresses in the aneurysm wall along both the stiff and weak fiber directions showed significant regional variations with the former being higher. The spatial maximum of stress ranged from as low as 0.30 MPa in a small aneurysm to as high as 1.06 MPa in a giant aneurysm. The patterns of distribution of stress, strain, and surface curvature were found to be similar. Sensitivity analyses showed that the computed stress is mesh independent and not very sensitive to reasonable perturbations in model parameters, and the curvature-based criteria for fiber orientations tend to minimize the total elastic strain energy in the aneurysms wall. Within this small study population, there were no statistically significant differences in the spatial means and maximums of stress and strain values between the ruptured and unruptured groups. However, the ratios between the stress components in the stiff and weak fiber directions were significantly higher in the ruptured group than those in the unruptured group. CONCLUSIONS: A methodology for nonlinear, anisotropic static deformation analysis of geometrically realistic aneurysms was developed, which can be used for a more accurate estimation of the stresses and strains than previous methods and to facilitate prospective studies on the role of stress in aneurysm rupture.     

The effect of mating system on invasiveness: some genetic load may be advantageous when invading new environments

The role of adaptation in determining invasion success has been acknowledged recently, notably through the accumulation of case studies of rapid evolution during bioinvasions. Despite this growing body of empirical evidence, there is still a need to develop the theoretical background of invasions with adaptation. Specifically, the impact of mating system on the dynamics of adaptation during invasion of a new environment remains only partially understood. Here, we analyze a simulation demo-genetic model of bioinvasion accounting for partial asexuality rates. We simulate two levels of recurrent immigration from a source population at mutation–drift–selection equilibrium to a new empty environment with a different adaptive landscape (black-hole sink). Adaptation relies on a quantitative trait coded explicitly by 10 loci under mutation, selection and genetic drift. Using this model, we confirm previous results on the positive effects on invasiveness of migration, mutation and similarity of local phenotypic optima. We further show how the invasion dynamics of the introduced population is affected by the rate of asexuality. Purely asexual species have lower invasion success in terms of probability and time to invasion than species with other mating systems. Among species with mixed mating systems, the greatest invasiveness is observed for species with high asexual rates. We suggest that this pattern is due to inflated genetic variance in the source population through the Hill-Robertson effect (i.e., clonal interference). An interesting consequence is that species with the highest genetic load in their source environment have greatest invasiveness in the new environment.     

Computer modeling of deployment and mechanical expansion of neurovascular flow diverter in patient-specific intracranial aneurysms

Flow diverter (FD) is an emerging neurovascular device based on self-expandable braided stent for treating intracranial aneurysms. Variability in FD outcome has underscored a need for investigating the hemodynamic effect of fully deployed FD in patient-specific aneurysms. Image-based computational fluid dynamics, which can provide important hemodynamic insight, requires accurate representation of FD in deployed states. We developed a finite element analysis (FEA) based workflow for simulating mechanical deployment of FD in patient-specific aneurysms. We constructed FD models of interlaced wires emulating the Pipeline Embolization Device, using 3D finite beam elements to account for interactions between stent strands, and between the stent and other components. The FEA analysis encompasses all steps that affect the final deployed configuration including stent crimping, delivery and expansion. Besides the stent, modeling also includes key components of the FD delivery system such as microcatheter, pusher, and distal coil. Coordinated maneuver of these components allowed the workflow to mimic clinical operation of FD deployment and to explore clinical strategies. The workflow was applied to two patient-specific aneurysms. Parametric study indicated consistency of the deployment result against different friction conditions, but excessive intra-stent friction should be avoided. This study demonstrates for the first time mechanical modeling of braided FD stent deployment in cerebral vasculature to produce realistic deployed configuration, thus paving the way for accurate CFD analysis of flow diverters for reliable prediction and optimization of treatment outcome.     

The impact of glenoid labrum thickness and modulus on labrum and glenohumeral capsule function

The glenoid labrum is an integral component of the glenohumeral capsule's insertion into the glenoid, and changes in labrum geometry and mechanical properties may lead to the development of glenohumeral joint pathology. The objective of this research was to determine the effect that changes in labrum thickness and modulus have on strains in the labrum and glenohumeral capsule during a simulated physical examination for anterior instability. A labrum was incorporated into a validated, subject-specific finite element model of the glenohumeral joint, and experimental kinematics were applied simulating application of an anterior load at 0 deg, 30 deg, and 60 deg of external rotation and 60 deg of glenohumeral abduction. The radial thickness of the labrum was varied to simulate thinning tissue, and the tensile modulus of the labrum was varied to simulate degenerating tissue. At 60 deg of external rotation, a thinning labrum increased the average and peak strains in the labrum, particularly in the labrum regions of the axillary pouch (increased 10.5% average strain) and anterior band (increased 7.5% average strain). These results suggest a cause-and-effect relationship between age-related decreases in labrum thickness and increases in labrum pathology. A degenerating labrum also increased the average and peak strains in the labrum, particularly in the labrum regions of the axillary pouch (increased 15.5% strain) and anterior band (increased 10.4% strain). This supports the concept that age-related labrum pathology may result from tissue degeneration. This work suggests that a shift in capsule reparative techniques may be needed in order to include the labrum, especially as activity levels in the aging population continue to increase. In the future validated, finite element models of the glenohumeral joint can be used to explore the efficacy of new repair techniques for glenoid labrum pathology.     

Vinculin Force-Sensitive Dynamics at Focal Adhesions Enable Effective Directed Cell Migration

Cell migration is a complex process, requiring coordination of many subcellular processes including membrane protrusion, adhesion, and contractility. For efficient cell migration, cells must concurrently control both transmission of large forces through adhesion structures and translocation of the cell body via adhesion turnover. Although mechanical regulation of protein dynamics has been proposed to play a major role in force transmission during cell migration, the key proteins and their exact roles are not completely understood. Vinculin is an adhesion protein that mediates force-sensitive processes, such as adhesion assembly under cytoskeletal load. Here, we elucidate the mechanical regulation of vinculin dynamics. Specifically, we paired measurements of vinculin loads using a Förster resonance energy transfer-based tension sensor and vinculin dynamics using fluorescence recovery after photobleaching to measure force-sensitive protein dynamics in living cells. We find that vinculin adopts a variety of mechanical states at adhesions, and the relationship between vinculin load and vinculin dynamics can be altered by the inhibition of vinculin binding to talin or actin or reduction of cytoskeletal contractility. Furthermore, the force-stabilized state of vinculin required for the stabilization of membrane protrusions is unnecessary for random migration, but is required for directional migration along a substrate-bound cue. These data show that the force-sensitive dynamics of vinculin impact force transmission and enable the mechanical integration of subcellular processes. These results suggest that the regulation of force-sensitive protein dynamics may have an underappreciated role in many cellular processes.     

3D reconstruction and manufacture of real abdominal aortic aneurysms: from CT scan to silicone model

Abdominal aortic aneurysm (AAA) can be defined as a permanent and irreversible dilation of the infrarenal aorta. AAAs are often considered to be an aorta with a diameter 1.5 times the normal infrarenal aorta diameter. This paper describes a technique to manufacture realistic silicone AAA models for use with experimental studies. This paper is concerned with the reconstruction and manufacturing process of patient-specific AAAs. 3D reconstruction from computed tomography scan data allows the AAA to be created. Mould sets are then designed for these AAA models utilizing computer aided designcomputer aided manufacture techniques and combined with the injection-moulding method. Silicone rubber forms the basis of the resulting AAA model. Assessment of wall thickness and overall percentage difference from the final silicone model to that of the computer-generated model was performed. In these realistic AAA models, wall thickness was found to vary by an average of 9.21%. The percentage difference in wall thickness recorded can be attributed to the contraction of the casting wax and the expansion of the silicone during model manufacture. This method may be used in conjunction with wall stress studies using the photoelastic method or in fluid dynamic studies using a laser-Doppler anemometry. In conclusion, these patient-specific rubber AAA models can be used in experimental investigations, but should be assessed for wall thickness variability once manufactured.