Numerical simulation of the failure of ventricular tissue due to deep penetration: The impact of constitutive properties

Specimen width affects vascular tissue integrity for in-vitro characterisation

The preparation of slender specimens for in-vitro tissue characterisation could potentially alter mechanical tissue properties. To investigate this factor, rectangular specimens were prepared from the wall of the porcine aorta for uniaxial tensile loading. Varying strip widths of 16 mm, 8 mm, and 4 mm were achieved by excising zero, one, and three cuts within the specimen along the loading direction, respectively. While specimens loaded along the vessel's circumferential direction acquired consistent tissue properties, the width of test specimens influenced the results of axially loaded tissue; vascular wall stiffness was reduced by approximately 40% in specimens with strips 4 mm wide. In addition, the cross-loading stretch was strongly influenced by specimen strip width, and fiber sliding contributed to the softening of slender tensile specimens, an outcome from finite element analysis of test specimens. We may, therefore, conclude that cutting orthogonal to the main direction of collagen fibers introduces mechanical trauma that weakens slender tensile specimens, compromising the determination of representative mechanical vessel wall properties.

Fracture of the porcine aorta. Part 1: symconCT fracture testing and DIC

Tissue failure and damage are inherent parts of vascular diseases and tightly linked to clinical events. Additionally, experimental set-ups designed to study classical engineering materials are suboptimal in the exploration of vessel wall fracture properties. The classical Compact Tension (CT) test was augmented to enable stable fracture propagation, resulting in the symmetry-constraint Compact Tension (symconCT) test, a suitable set-up for fracture testing of vascular tissue. The test was combined with Digital Image Correlation (DIC) to study tissue fracture in 45 porcine aorta specimens. Test specimens were loaded in axial and circumferential directions in a physiological solution at 37∘ C. Loading the aortic vessel wall in the axial direction resulted in mode I tissue failure and a fracture path aligned with the circumferential vessel direction. Circumferential loading resulted in mode I-dominated failure with multiple deflections of the fracture path. The aorta ruptured at a principal Green-Lagrange strain of approximately 0.7, and strain rate peaks that develop ahead of the crack tip reached nearly 400 times the strain rate on average over the test specimen. It required approximately 70% more external work to fracture the aorta by circumferential than axial load; normalised with the fracture surface, similar energy levels are, however, observed. The symconCT test resulted in a stable fracture propagation, which, combined with DIC, provided a set-up for the in-depth analysis of vascular tissue failure. The high strain rates ahead of the crack tip indicate the significance of rate effects in the constitutive description of vascular tissue fracture. STATEMENT OF SIGNIFICANCE: This paper represents a significant step forward in understanding the fracture properties of porcine aorta. Inspired by the Compact Tension test, we developed an ad hoc experimental protocol to investigate stable crack propagation in soft materials, providing new insights into the mechanical processes that lead to the rupture of vascular tissue. The set-up enables the assessment of strains and strain rates ahead of the crack tip, and our findings could improve the clinical risk assessment of vascular pathologies as well as optimize medical device design.

Microstructural deformation observed by Mueller polarimetry during traction assay on myocardium samples

Despite recent advances, the myocardial microstructure remains imperfectly understood. In particular, bundles of cardiomyocytes have been observed but their three-dimensional organisation remains debated and the associated mechanical consequences unknown. One of the major challenges remains to perform multiscale observations of the mechanical response of the heart wall. For this purpose, in this study, a full-field Mueller polarimetric imager (MPI) was combined, for the first time, with an in-situ traction device. The full-field MPI enables to obtain a macroscopic image of the explored tissue, while providing detailed information about its structure on a microscopic scale. Specifically it exploits the polarization of the light to determine various biophysical quantities related to the tissue scattering or anisotropy properties. Combined with a mechanical traction device, the full-field MPI allows to measure the evolution of such biophysical quantities during tissue stretch. We observe separation lines on the tissue, which are associated with a fast variation of the fiber orientation, and have the size of cardiomyocyte bundles. Thus, we hypothesize that these lines are the perimysium, the collagen layer surrounding these bundles. During the mechanical traction, we observe two mechanisms simultaneously. On one hand, the azimuth shows an affine behavior, meaning the orientation changes according to the tissue deformation, and showing coherence in the tissue. On the other hand, the separation lines appear to be resistant in shear and compression but weak against traction, with a forming of gaps in the tissue.

Feasibility and Safety of Percutaneous Lead Revision for Subacute and Delayed Cardiac Device Lead Perforation

OBJECTIVES

This study assessed the management approach and outcome of subacute (1 to 30 days post-implantation) and delayed (>30 days) cardiac perforation by pacemaker or implantable cardioverter-defibrillator (ICD) leads.

BACKGROUND

Implantation of pacemaker and ICD leads is associated with a small but serious risk of cardiac perforation. Appropriate management remains uncertain.

METHODS

The study population included all patients referred to a single institution for subacute or delayed lead perforation after pacemaker or ICD implantation (identified after hospital discharge) during the period from 2007 to 2020. The approach and outcome of lead management were retrospectively assessed.

RESULTS

Fifty-four cases of cardiac perforation were identified (35 females; mean age: 75.5 ± 9.7 years). Cardiac perforation was related to a pacemaker lead in 36 patients, and the perforating leads were originally placed in the right ventricular apex in 41 patients. The average time from lead implantation to first presentation of symptoms of perforation was 60.8 ± 89.1 days (range 2 to 412 days). Symptoms suggestive of cardiac perforation were reported by 31 patients (57.4%). Twenty three patients were asymptomatic, in whom lead perforation was discovered incidentally on radiographic imaging, suggesting lead migration or anomalous electrical data on device interrogation. In all patients, the leads were removed or repositioned by the percutaneous approach, with no major periprocedural complications and without surgical intervention.

CONCLUSIONS

In this largest series to date of subacute or delayed cardiac device lead perforation, percutaneous repositioning or replacement of the perforating lead was found to be a safe and effective management approach.

Structural modelling of the cardiovascular system

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

Calibration of interface properties and application to a finite element model for predicting vena cava filter-induced vein wall failure

We present a computational framework that integrates experimental techniques and finite element modeling to calibrate material fracture parameters of the vena cava and the interaction properties between a retrievable filter (Günther Tulip) and the vena cava wall. The fitted parameters were then used to analyze the interaction of the inferior vena cava filter with the vena cava during the deployment process. An idealized cava finite element model was then developed including residual stresses and physiological pressure conditions. Filter deployment was simulated, and a comprehensive study of tissue-filter interaction was performed by cohesive surface modeling. Simulations predict that there are no fracture areas for either model, so we can conclude that there is no penetration of the anchor into the vena cava. This suggests there are other physiological situations, such as the Valsalva maneuver, which could produce this penetration observed on some patients.

Contribution of computational model for assessment of heart tissue local stress caused by suture in LVAD implantation

STUDY

Implantation of a Left Ventricular Assist Device (LVAD) may produce both excessive local tissue stress and resulting strain-induced tissue rupture that are potential iatrogenic factors influencing the success of the surgical attachment of the LVAD into the myocardium. By using a computational simulation compared to mechanical tests, we sought to investigate the characteristics of stress-induced suture material on porcine myocardium.

METHODS

Tensile strength experiments (n = 8) were performed on bulk left myocardium to establish a hyperelastic reduced polynomial constitutive law. Simultaneously, suture strength tests on left myocardium (n = 6) were performed with a standard tensile test setup. Experiments were made on bulk ventricular wall with a single U-suture (polypropylene 3-0) and a PTFE pledget. Then, a Finite Element simulation of a LVAD suture case was performed. Strength versus displacement behavior was compared between mechanical and numerical experiments. Local stress fields in the model were thus analyzed.

RESULTS

A strong correlation between the experimental and the numerical responses was observed, validating the relevance of the numerical model. A secure damage limit of 100 kPa on heart tissue was defined from mechanical suture testing and used to describe numerical results. The impact of suture on heart tissue could be accurately determined through new parameters of numerical data (stress diffusion, triaxiality stress). Finally, an ideal spacing between sutures of 2 mm was proposed.

CONCLUSION

Our computational model showed a reliable ability to provide and predict various local tissue stresses created by suture penetration into the myocardium. In addition, this model contributed to providing valuable information useful to design less traumatic sutures for LVAD implantation. Therefore, our computational model is a promising tool to predict and optimize LVAD myocardial suture.

An Invariant-Based Damage Model for Human and Animal Skins

Constitutive modelling of skins that account for damage effects is important to provide insight for various clinical applications, such as skin trauma and injury, artificial skin design, skin aging, disease diagnosis, surgery, as well as comparative studies of skin biomechanics between species. In this study, a new damage model for human and animal skins is proposed for the first time. The model is nonlinear, anisotropic, invariant-based, and is based on the Gasser-Ogden-Holzapfel constitutive law initially developed for arteries. Taking account of the mean collagen fibre orientation and its dispersion, the new model can describe a wide range of skins with damage. The model is first tested on the uniaxial test data of human skin and then applied to nine groups of uniaxial test data for the human, swine, rabbit, bovine and rhino skins. The material parameters can be inversely estimated based on uniaxial tests using the optimization method in MATLAB with a root mean square error ranged between 2.15% and 12.18%. A sensitivity study confirms that the fibre orientation dispersion and the mean fibre angle are among the most important factors that influence the behaviour of the damage model. In addition, these two parameters can only be reliably estimated if some histological information is provided. We also found that depending on the location of skins, the tissue damage may be brittle controlled by the fibre breaking limit (i.e., when the fibre stretch is greater than 1.13-1.32, depending on the species), or ductile (due to both the fibre and the matrix damages). The brittle damages seem to occur mostly in the back, and the ductile damages are seen from samples taken from the belly. The proposed constitutive model may be applied to various clinical applications that require knowledge of the mechanical response of human and animal skins.

Failure properties of vena cava tissue due to deep penetration during filter insertion

In this work, we use an in-vitro mechanical test to explore the resistance of biaxially stretched vena cava tissue against deep perforation and a methodology which integrates experimental and numerical modeling to identify constitutive fracture properties of the vena cava. Six sheep vena cava were harvested just after killing, and cyclic uniaxial tension tests in longitudinal and circumferential directions and biaxial deep penetration tests were performed. After that, we use a nonlinear finite element model to simulate in vitro penetration of the cava tissue in order to fit the fracture properties under penetration of the vena cava by defining a cohesive fracture zone. An iterative process was developed in order to fit the fracture properties of the vena cava using the previously obtained experimental results. The proposed solutions were obtained with fracture energy of 0.22 or 0.33 N/mm. In comparison with the experimental data, the simulation using [Formula: see text], [Formula: see text], and [Formula: see text] parameters ([Formula: see text]) is in good agreement with results from penetration experiments of cava tissue. It is noticeable that the parameter estimation process of the fracture behavior is more accurate than the estimation process of the elastic behavior for the toe region of the curve.

Computational Investigation of Fibrin Mechanical and Damage Properties at the Interface Between Native Cartilage and Implant

Scaffold-based tissue-engineered constructs as well as cell-free implants offer promising solutions to focal cartilage lesions. However, adequate mechanical stability of these implants in the lesion is required for successful repair. Fibrin is the most common clinically available adhesive for cartilage implant fixation, but fixation quality using fibrin is not well understood. The objectives of this study were to investigate the conditions leading to damage in the fibrin adhesive and to determine which adhesive properties are important in preventing delamination at the interface. An idealized finite element model of the medial compartment of the knee was created, including a circular defect and an osteochondral implant. Damage and failure of fibrin at the interface was represented by a cohesive zone model with coefficients determined from an inverse finite element method and previously published experimental data. Our results demonstrated that fibrin glue alone may not be strong enough to withstand physiologic loads in vivo while fibrin glue combined with chondrocytes more effectively prevents damage at the interface. The results of this study suggest that fibrin fails mainly in shear during off-axis loading and that adhesive materials that are stronger or more compliant than fibrin may be good alternatives due to decreased failure at the interface. The present model may be used to improve design and testing protocols of bioadhesives and give insight into the failure mechanisms of cartilage implant fixation in the knee joint.

Development of a quantitative mechanical test of atherosclerotic plaque stability

Atherosclerotic plaque rupture is the main cause of myocardial infarction and stroke. Both clinical and computational studies indicate that the shoulder region, where a plaque joins the vessel wall, is rupture-prone. Previous mechanistic studies focused on mechanical properties of the fibrous cap and tensile stresses, which could lead to tearing of the cap. Based on clinical observations of "mobile floating plaques," we postulate that de-adhesion between the fibrous cap and the underlying vessel wall may also play a role in plaque failure. Thus, measuring adhesive strength of the bond between plaque and vascular wall may provide useful new insights into plaque stability. Delamination experiments, widely used in examining inter-laminar adhesive strength of biological materials, were used to measure adhesive strength of advanced plaques in apolipoprotein E-knockout (apoE-KO) mice after 8 months on Western diet. We measured adhesive strength in terms of local energy release rate, G, during controlled plaque delamination. As a measure of the fracture energy required to delaminate a unit area of plaque from the underlying internal elastic lamina (IEL), G provides a quantitative measure of local adhesive strength of the plaque-IEL interface. The values for G acquired from 16 plaques from nine apoE-KO mouse aortas formed a positively skewed distribution with a mean of 24.5 J/m(2), median of 19.3 J/m(2), first quartile of 10.8 J/m(2), and third quartile of 34.1 J/m(2). These measurements are in the lower range of values reported for soft tissues. Histological studies confirmed delamination occurred at the interface between plaque and IEL.

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.