Inverse material parameter estimation of patient-specific finite element models at the carotid bifurcation: The impact of excluding the zero-pressure configuration and residual stress

Comparative analysis of Zero Pressure Geometry and prestress methods in cardiovascular Fluid-Structure Interaction

BACKGROUND AND OBJECTIVE

Modelling patient-specific aortic biomechanics with advanced computational techniques, such as Fluid-Structure Interaction (FSI), can be crucial to provide effective decision-making indices to enhance current clinical practices. To effectively simulate Ascending Thoracic Aortic Aneurysms (ATAA), the stress-free configuration must be defined. The Zero Pressure Geometry (ZPG) and the Prestress Tensor (PT) are two of the main approaches to tackle this issue. However, their impact on the numerical results is yet to be analysed. Computed Tomography Angiography (CTA) and Magnetic Resonance Imaging (MRI) data were used to develop patient-specific 2-way FSI frameworks.

METHODS

Three models were developed considering different tissue prestressing approaches to account for the reference configuration and their numerical results were compared. The selected approaches were: (i) ZPG, (ii) PT and (iii) a combination of the PT approach with a regional mapping of material properties (PTCAL).

RESULTS

The pressure fields estimated by all models were equivalent. The estimation of Wall Shear Stress (WSS) based metrics revealed good correspondence between all models except the Relative Residence Time (RRT). Regarding ATAA wall mechanics, the proposed extension to the PT approach presented a closer agreement with the ZPG model than its counterpart. Additionally, the PT and PTCAL approaches required around 60% fewer iterations to achieve cycle-to-cycle convergence than the ZPG algorithm.

CONCLUSION

Using a regional mapping of material properties in combination with the PT method presented a better correspondence with the ZPG approach. The outcomes of this study can pave the way for advancing the accuracy and convergence of ATAA numerical models using the PT methodology.

Characterisation of human penile tissue properties using experimental testing combined with multi-target inverse finite element modelling

This paper presents an inverse finite element (FE) approach aimed at estimating multi-layered human penile tissues. The inverse FE approach integrates experimental force-displacement and boundary deformation data of penile tissues with a developed FE model and uses new experimental data on human penile tissue. The experimental study encompasses whole organ plate-compression tests and individual layer tensile and compression tests, providing comprehensive insights into the tissue's mechanical behaviour. The biomechanical characterisation of penile tissue is of crucial significance for understanding its mechanical behaviour under various physiological and pathological conditions. The FE model is constructed using the realistic geometry of the penile segment and appropriate constitutive models for each tissue layer to leverage the accuracy and consistency of the model. Through systematic variation of tissue parameters in the inverse FE algorithm, simulations achieve the best match with both force-displacement and deformed boundary results obtained from the whole organ plate-compression tests. Test results from individual tissue layers are also utilised to assess the estimated parameters. The proposed inverse FE approach allows for the estimation of penile tissue parameters with high precision and reliability, shedding light on the mechanical properties of this complex biological organ. This work has applications not only in urology but also for researchers in various disciplines of biomechanics. As a result, our study contributes to advancing the understanding of human penile tissue mechanics whilst the methodology could also be applied to a range of other soft biological tissues. STATEMENT OF SIGNIFICANCE: This research uses a multi-target inverse finite element (FE) approach for estimating the material parameters of human penile tissues. By integrating experimental data and a realistic FE model, this study achieves high-precision constitutive model parameter estimation, offering key insights into penile tissue mechanics under various loading conditions. The significance of this work lies in the use of this inverse FE approach for fresh-frozen human penile tissues, to identify the mechanical properties and constitutive models for both segregated tunica albuginea and corpus cavernosum as well as intact penile tissue segments. The study's scientific impact lies in its advancement of the understanding of human urological tissue mechanics, impacting researchers and clinicians alike.

Experimental testing combined with inverse-FE for mechanical characterisation of penile tissues

Erectile dysfunction (ED) predominantly affects men in their 40-70s and can lead to poor quality of life. One option for ED treatment is surgical implantation of an inflatable penile prosthesis (IPP). However, they can be associated with negative outcomes including infection, migration or fibrosis. To improve outcomes, the interaction between the IPP device and surrounding tissues needs further investigation and this could be achieved using pre-clinical testbeds, but they need to be informed by extensive tissue testing. In this study, an experimental approach is adopted to characterise the mechanics of horse penile tissue and establish a testing protocol for penile tissue. The whole penis segments were tested in plate compression tests to obtain whole penis behaviour which is necessary for validation of a pre-clinical testbed, whilst tensile and compression tests were performed on individual penile tissues, namely corpus cavernosa and tunica albuginea. The second part of the paper deals with the development of a computational model employing an inverse finite element approach to estimate the material parameters of each tissue layer. These material parameters are in good agreement with the experimental results obtained from the individual tissue layers and whole organ tissue tests. This paper presents the first study proposing realistic nonlinear elastic material parameters for penile tissues and offers a validated testbed for IPPs. STATEMENT OF SIGNIFICANCE: Erectile Dysfunction (ED) affects over half the male population aged 40-70 potentially leading to poor quality of life. Patients not responding to conventional treatments of ED, are advised to use penile prostheses which can create an erection using implanted inflatable cylinders. A significant drawback of such prostheses, however, is the substantial tissue damage they can induce during their usage. Preclinical testbeds, including computational and bench-top models, could offer an efficient means of improving device designs to mitigate this damage but such testbeds require extensive knowledge of penile tissue properties. In this study, the authors determine penile tissue mechanics and apply an inverse FE approach to characterise the penile material properties required to validate preclinical models of the penis.

Computational Study of Abdominal Aortic Aneurysm Walls Accounting for Patient-Specific Non-Uniform Intraluminal Thrombus Thickness and Distinct Material Models: A Pre- and Post-Rupture Case

An intraluminal thrombus (ILT) is present in the majority of abdominal aortic aneurysms, playing a crucial role in their growth and rupture. Although most computational studies do not include the ILT, in the present study, this is taken into account, laying out the whole simulation procedure, namely, from computed tomography scans to medical image segmentation, geometry reconstruction, mesh generation, biomaterial modeling, finite element analysis, and post-processing, all carried out in open software. By processing the tomography scans of a patient's aneurysm before and after rupture, digital twins are reconstructed assuming a uniform aortic wall thickness. The ILT and the aortic wall are assigned different biomaterial models; namely, the first is modeled as an isotropic linear elastic material, and the second is modeled as the Mooney-Rivlin hyperelastic material as well as the transversely isotropic hyperelastic Holzapfel-Gasser-Ogden nonlinear material. The implementation of the latter requires the designation of local Cartesian coordinate systems in the aortic wall, suitably oriented in space, for the proper orientation of the collagen fibers. The composite aneurysm geometries (ILT and aortic wall structures) are loaded with normal and hypertensive static intraluminal pressure. Based on the calculated stress and strain distributions, ILT seems to be protecting the aneurysm from a structural point of view, as the highest stresses appear in the thrombus-free areas of the aneurysmal wall.