Simulation of arterial dissection by a penetrating external body using cohesive zone modelling

A Review on Damage and Rupture Modelling for Soft Tissues

Computational modelling of damage and rupture of non-connective and connective soft tissues due to pathological and supra-physiological mechanisms is vital in the fundamental understanding of failures. Recent advancements in soft tissue damage models play an essential role in developing artificial tissues, medical devices/implants, and surgical intervention practices. The current article reviews the recently developed damage models and rupture models that considered the microstructure of the tissues. Earlier review works presented damage and rupture separately, wherein this work reviews both damage and rupture in soft tissues. Wherein the present article provides a detailed review of various models on the damage evolution and tear in soft tissues focusing on key conceptual ideas, advantages, limitations, and challenges. Some key challenges of damage and rupture models are outlined in the article, which helps extend the present damage and rupture models to various soft tissues.

On a phase-field approach to model fracture of small intestine walls

We address anisotropic elasticity and fracture in small intestine walls (SIWs) with both experimental and computational methods. Uniaxial tension experiments are performed on porcine SIW samples with varying alignments and quantify their nonlinear elastic anisotropic behavior. Fracture experiments on notched SIW strips reveal a high sensitivity of the crack propagation direction and the failure stress on the tissue orientation. From a modeling point of view, the observed anisotropic elastic response is studied with a continuum mechanical model stemming from a strain energy density with a neo-Hookean component and an anisotropic component with four families of fibers. Fracture is addressed with the phase-field approach, featuring two-fold anisotropy in the fracture toughness. Elastic and fracture model parameters are calibrated based on the experimental data, using the maximum and minimum limits of the experimental stress-stretch data set. A very good agreement between experimental data and computational results is obtained, the role of anisotropy being effectively captured by the proposed model in both the elastic and the fracture behavior. STATEMENT OF SIGNIFICANCE: This article reports a comprehensive experimental data set on the mechanical failure behavior of small intestinal tissue, and presents the corresponding protocols for preparing and testing the samples. On the one hand, the results of this study contribute to the understanding of small intestine mechanics and thus to understanding of load transfer mechanisms inside the tissue. On the other hand, these results are used as input for a phase-field modelling approach, presented in this article. The presented model can reproduce the mechanical failure behavior of the small intestine in an excellent way and is thus a promising tool for the future mechanical description of diseased small intestinal tissue.

Development of a test method to investigate mode II fracture and dissection of arteries

The current study presents the development and implementation of a bespoke experimental technique to generate and characterise mode II crack initiation and propagation in arterial tissue. The current study begins with a demonstration that lap-shear testing of arterial tissue results in mixed mode fracture, rather than mode II. We perform a detailed computational design of a bespoke experimental method (which we refer to as a shear fracture ring test (SFRT)) to robustly and repeatably generate mode II crack initiation and propagation in arteries. This method is based on generating a localised region of high shear adjacent to a cylindrical loading bar. Placement of a radial notch in this region of high shear stress is predicted to result in a kinking of the crack during a mode II initiation and propagation of the crack over a long distance in the circumferential (c)-direction along the circumferential-axial (c-a) plane. Fabrication and experimental implementation of the SFRT on excised ovine aorta specimens confirms that the bespoke test method results in pure mode II initiation and propagation. We demonstrate that the mode II fracture strength along the c-a plane is eight times higher than the corresponding mode I strength determined from a standard peel test. We also calibrate the mode II fracture energy based on our measurement of crack propagation rates. The mechanisms of fracture uncovered in the current study, along with our quantification of mode II fracture properties have significant implications for current understanding of the biomechanical conditions underlying aortic dissection.

Computational modeling reveals the relationship between intrinsic failure properties and uniaxial biomechanical behavior of arterial tissue

Biomechanical failure of the artery wall can lead to rupture, a catastrophic event with a high rate of mortality. Thus, there is a pressing need to understand failure behavior of the arterial wall. Uniaxial testing remains the most common experimental technique to assess tissue failure properties. However, the relationship between intrinsic failure parameters of the tissue and measured uniaxial failure properties is not fully established. Furthermore, the effect of the experimental variables, such as specimen shape and boundary conditions, on the measured failure properties is not well understood. We developed a finite element model capable of recapitulating pre-failure and post-failure uniaxial biomechanical response of the arterial tissue specimen. Intrinsic stiffness, strength and fracture toughness of the vessel wall tissue were used as the input material parameters to the model. Two uniaxial testing protocols were considered: a conventional setup with a rectangular specimen held at the grips by cardboard inserts, and the other used a dogbone specimen with soft foam inserts at the grips. Our computational study indicated negligible differences in the peak stress and post-peak mechanical behavior between these two testing protocols. It was also found that the tissue experienced only modest localized failure until higher levels of applied stretch beyond the peak stress. A robust cohesive model was capable of modeling the post-peak biomechanical response, which was primarily governed by tissue fracture toughness. Our results suggest that the post-peak region, in conjunction with the peak stress, must be considered to evaluate the complete biomechanical failure behavior of the soft tissue.

Effect of cyclic deformation on xenogeneic heart valve biomaterials

Glutaraldehyde-fixed bovine pericardium is currently the most popular biomaterial utilized in the creation of bioprosthetic heart valves. However, recent studies indicate that glutaraldehyde fixation results in calcification and structural valve deterioration, limiting the longevity of bioprosthetic heart valves. Additionally, glutaraldehyde fixation renders the tissue incompatible with constructive recipient cellular repopulation, remodeling and growth. Use of unfixed xenogeneic biomaterials devoid of antigenic burden has potential to overcome the limitations of current glutaraldehyde-fixed biomaterials. Heart valves undergo billion cycles of opening and closing throughout the patient’s lifetime. Therefore, understanding the response of unfixed tissues to cyclic loading is crucial to these in a heart valve leaflet configuration. In this manuscript we quantify the effect of cyclic deformation on cycle dependent strain, structural, compositional and mechanical properties of fixed and unfixed tissues. Glutaraldehyde-fixed bovine pericardium underwent marked cyclic dependent strain, resulting from significant changes in structure, composition and mechanical function of the material. Conversely, unfixed bovine pericardium underwent minimal strain and maintained its structure, composition and mechanical integrity. This manuscript demonstrates that unfixed bovine pericardium can withstand cyclic deformations equivalent to 6 months of in vivo heart valve leaflet performance.

Computational modeling of progressive damage and rupture in fibrous biological tissues: application to aortic dissection

This study analyzes the lethal clinical condition of aortic dissections from a numerical point of view. On the basis of previous contributions by Gültekin et al. (Comput Methods Appl Mech Eng 312:542-566, 2016 and 331:23-52, 2018), we apply a holistic geometrical approach to fracture, namely the crack phase-field, which inherits the intrinsic features of gradient damage and variational fracture mechanics. The continuum framework captures anisotropy, is thermodynamically consistent and is based on finite strains. The balance of linear momentum and the crack evolution equation govern the coupled mechanical and phase-field problem. The solution scheme features the robust one-pass operator-splitting algorithm upon temporal and spatial discretizations. Based on experimental data of diseased human thoracic aortic samples, the elastic material parameters are identified followed by a sensitivity analysis of the anisotropic phase-field model. Finally, we simulate an incipient propagation of an aortic dissection within a multi-layered segment of a thoracic aorta that involves a prescribed initial tear. The finite element results demonstrate a severe damage zone around the initial tear and exhibit a rather helical crack pattern, which aligns with the fiber orientation. It is hoped that the current contribution can provide some directions for further investigations of this disease.

3D patient-specific numerical modeling of the soft palate considering adhesion from the tongue

Collapse of the soft palate in the upper airway contributes to obstructive sleeping apnea (OSA). In this study, we investigate the influence of the adhesion from the tongue on the soft palate global response. This is achieved using a cohesive zone finite element approach. A traction-separation law is determined to describe the adhesion effect from the surface tension of the lining liquid between the soft palate and the tongue. According to pull-off experimental tests of human lining liquid from the oral surface of the soft palate, the corresponding cohesive properties, including the critical normal traction stress and the failure separation displacement, are obtained. The 3D patient-specific soft palate geometry is accounted for, based on one specific patient's computed tomography (CT) images. The calculation results show that influence of the adhesion from the tongue surface on the global response of the soft palate depends on the length ratio between the cohesive length and the soft palate length. When the length of the cohesive zone is smaller than half of the soft palate length, the adhesion's influence is negligible. When the adhesion length is larger than 70 percent of soft palate length, the adhesion force contributes to preventing the soft palate from collapsing towards to the pharynx wall, i.e. the closing pressure is more negative than in the no adhesion case. These results may provide useful information to the clinical treatment of OSA patients.