Computational Biomechanics: In-Silico Tools for the Investigation of Surgical Procedures and Devices

In silico assessment of the reliability and performance of artificial sphincter for urinary incontinence

BACKGROUND

The standard artificial urinary sphincter (AUS) is an implantable device for the treatment of urinary incontinence by applying a pressure loading around the urethra through an inflatable cuff, often inducing no-physiological stimulation up to tissue degenerative phenomena. A novel in silico approach is proposed to fill the gap of the traditional procedures by providing tools to quantitatively assess AUS reliability and performance based on AUS-urethra interaction.

METHODS

The approach requires the development of 3D numerical models of AUS and urethra, and experimental investigations to define their mechanical behaviors. Computational analyses are performed to simulate the urethral lumen occlusion by AUS inflation under different pressures, and the lumen opening by applying an intraluminal pressure progressively increased under the AUS action (Abaqus Explicit solver). The AUS reliability is evaluated in terms of tissue stimulation by the mechanical fields potentially responsible for vasoconstriction and tissue damage, while the performance by the intraluminal pressure that causes the lumen opening for a specific occlusive pressure, showing the maximum urethral pressure for which continence is guaranteed.

RESULTS

The present study implemented the procedure considering the gold standard AMS 800 and a novel patented AUS. Results provided the comparison between two sphincteric devices and the evaluation of the influence of different building materials and geometrical features on the AUS functionality.

CONCLUSIONS

The approach was developed for the AUS, but it could be adapted also to artificial sphincters for the treatment of other anatomical dysfunctions, widening the analyzable device configurations and reducing experimental and ethical efforts.

Ex vivo, in vivo and in silico studies of corneal biomechanics: a systematic review

Healthy cornea guarantees the refractive power of the eye and the protection of the inner components, but injury, trauma or pathology may impair the tissue shape and/or structural organization and therefore its material properties, compromising its functionality in the ocular visual process. It turns out that biomechanical research assumes an essential role in analysing the morphology and biomechanical response of the cornea, preventing pathology occurrence, and improving/optimising treatments. In this review, ex vivo, in vivo and in silico methods for the corneal mechanical characterization are reported. Experimental techniques are distinct in testing mode (e.g., tensile, inflation tests), samples' species (human or animal), shape and condition (e.g., healthy, treated), preservation methods, setup and test protocol (e.g., preconditioning, strain rate). The meaningful results reported in the pertinent literature are discussed, analysing differences, key features and weaknesses of the methodologies adopted. In addition, numerical techniques based on the finite element method are reported, incorporating the essential steps for the development of corneal models, such as geometry, material characterization and boundary conditions, and their application in the research field to extend the experimental results by including further relevant aspects and in the clinical field for diagnostic procedure, treatment and planning surgery. This review aims to analyse the state-of-art of the bioengineering techniques developed over the years to study the corneal biomechanics, highlighting their potentiality to improve diagnosis, treatment and healing process of the corneal tissue, and, at the same, pointing out the current limits in the experimental equipment and numerical tools that are not able to fully characterize in vivo corneal tissues non-invasively and discourage the use of finite element models in daily clinical practice for surgical planning.

Biomechanical Modelling of Porcine Kidney

In this study, the viscoelastic properties of porcine kidney in the upper, middle and lower poles were investigated using oscillatory shear tests. The viscoelastic properties were extracted in the form of the storage modulus and loss modulus in the frequency and time domain. Measurements were taken as a function of frequency from 0.1 Hz to 6.5 Hz at a shear strain amplitude of 0.01 and as function of strain amplitude from 0.001 to 0.1 at a frequency of 1 Hz. Measurements were also taken in the time domain in response to a step shear strain. Both the frequency and time domain data were fitted to a conventional Standard Linear Solid (SLS) model and a semi-fractional Kelvin-Voigt (SFKV) model with a comparable number of parameters. The SFKV model fitted the frequency and time domain data with a correlation coefficient of 0.99. Although the SLS model well fitted the time domain data and the storage modulus data in the frequency domain, it was not able to capture the variation in loss modulus with frequency with a correlation coefficient of 0.53. A five parameter Maxwell-Wiechert model was able to capture the frequency dependence in storage modulus and loss modulus better than the SLS model with a correlation of 0.85.

Unveiling the effects of key factors in enhancing gastroesophageal reflux: A fluid-structure analysis before and after laparoscopic sleeve gastrectomy

BACKGROUND AND OBJECTIVES

Gastro-oesophageal reflux disease (GERD) consists in the passage of gastric acid content from the stomach to the oesophagus, causing burns and deteriorating the quality of life. Laparoscopic Sleeve Gastrectomy (LSG) could induce de novo GERD and worsen pre-existing GERD because of the higher gastric pressurisation, reduction of stomach volume and a wider His-angle. In the proposed work, various computational gastric 2D models were developed to understand the effects of variables such as the His-angle, the antral dimension, and the bolus viscosity on the reflux increase.

METHODS

Fluid-Structure Interaction (FSI) computational models which couple the solid mechanics of the gastric wall, and the fluid domain of the bolus, have been developed to shed light on biomechanical aspects of GERD after LSG. A closure was imposed to the lower oesophageal sphincter (LES) mimicking what happens physiologically after food intake.

RESULTS

Results showed that the configuration prone to higher reflux flow was the post-surgical 65° model with a staple line starting directly from the pylorus without antral preservation, for all considered viscosities. Increasing viscosity, reflux flow decreased. Post-surgical refluxes were higher than pre-ones and decreased with increasing antrum preservation.

CONCLUSIONS

These results could be a starting point for analysis of anatomical features, bariatric surgery and GERD occurrence. Further studies based on 3D geometries need to be performed.

Biomechanics of Hollow Organs: Experimental Testing and Computational Modeling

Hollow organs are visceral organs that are hollow tubes or pouches (such as the intestine or the stomach, respectively) or that include a cavity (such as the heart) and which subserve a vital function [...].

Open Knee(s): A Free and Open Source Library of Specimen-Specific Models and Related Digital Assets for Finite Element Analysis of the Knee Joint

There is a growing interest in the use of virtual representations of the knee for musculoskeletal research and clinical decision making, and to generate digital evidence for design and regulation of implants. Accessibility to previously developed models and related digital assets can dramatically reduce barriers to entry to conduct simulation-based studies of the knee joint and therefore help accelerate scientific discovery and clinical innovations. Development of models for finite element analysis is a demanding process that is both time consuming and resource intensive. It necessitates expertise to transform raw data to reliable virtual representations. Modeling and simulation workflow has many processes such as image segmentation, surface geometry generation, mesh generation and finally, creation of a finite element representation with relevant loading and boundary conditions. The outcome of the workflow is not only the end-point knee model but also many other digital by-products. When all of these data, derivate assets, and tools are freely and openly accessible, researchers can bypass some or all the steps required to build models and focus on using them to address their research goals. With provenance to specimen-specific anatomical and mechanical data and traceability of digital assets throughout the whole lifecycle of the model, reproducibility and credibility of the modeling practice can be established. The objective of this study is to disseminate Open Knee(s), a cohort of eight knee models (and relevant digital assets) for finite element analysis, that are based on comprehensive specimen-specific imaging data. In addition, the models and by-products of modeling workflows are described along with model development strategies and tools. Passive flexion served as a test simulation case, demonstrating an end-user application. Potential roadmaps for reuse of Open Knee(s) are also discussed.

Biomechanics Assist Measurement, Modeling, Engineering Applications, and Clinical Decision Making in Medicine

Biomechanical studies of surgeries and medical devices are usually performed with human or animal models [...].

Patient-specific stomach biomechanics before and after laparoscopic sleeve gastrectomy

BACKGROUND

Obesity has become a global epidemic. Bariatric surgery is considered the most effective therapeutic weapon in terms of weight loss and improvement of quality of life and comorbidities. Laparoscopic sleeve gastrectomy (LSG) is one of the most performed procedures worldwide, although patients carry a nonnegligible risk of developing post-operative GERD and BE.

OBJECTIVES

The aim of this work is the development of computational patient-specific models to analyze the changes induced by bariatric surgery, i.e., the volumetric gastric reduction, the mechanical response of the stomach during an inflation process, and the related elongation strain (ES) distribution at different intragastric pressures.

METHODS

Patient-specific pre- and post-surgical models were extracted from Magnetic Resonance Imaging (MRI) scans of patients with morbid obesity submitted to LSG. Twenty-three patients were analyzed, resulting in forty-six 3D-geometries and related computational analyses.

RESULTS

A significant difference between the mechanical behavior of pre- and post-surgical stomach subjected to the same internal gastric pressure was observed, that can be correlated to a change in the global stomach stiffness and a minor gastric wall tension, resulting in unusual activations of mechanoreceptors following food intake and satiety variation after LSG.

CONCLUSIONS

Computational patient-specific models may contribute to improve the current knowledge about anatomical and physiological changes induced by LSG, aiming at reducing post-operative complications and improving quality of life in the long run.

Mechanical Behavior of Subcutaneous and Visceral Abdominal Adipose Tissue in Patients with Obesity

The mechanical characterization of adipose tissues is important for various medical purposes, including plastic surgery and biomechanical applications, such as computational human body models for the simulation of surgical procedures or injury prediction, for example, in the evaluation of vehicle crashworthiness. In this context, the measurement of human subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) mechanical properties in relation to subject characteristics may be really relevant. The aim of this work was to properly characterize the mechanical response of adipose tissues in patients with obesity. Then, the data were exploited to develop a reliable finite element model of the adipose tissues characterized by a constitutive material model that accounted for nonlinear elasticity and time dependence. Mechanical tests have been performed on both SAT and VAT specimens, which have been harvested from patients with severe obesity during standard laparoscopic sleeve gastrectomy intervention. The experimental campaign included indentation tests, which permitted us to obtain the initial/final indentation stiffnesses for each specimen. Statistical results revealed a higher statistical stiffness in SAT than in VAT, with an initial/final indentation stiffness of 1.65 (SD ± 0.29) N/30.30 (SD ± 20) N compared to 1.29 (SD ± 0.30) N/21.00 (SD ± 16) N. Moreover, the results showed that gender, BMI, and age did not significantly affect the stiffness. The experimental results were used in the identification of the constitutive parameters to be inserted in the constitutive material model. Such constitutive characterization of VAT and SAT mechanics can be the starting point for the future development of more accurate computational models of the human adipose tissue and, in general, of the human body for the optimization of numerous medical and biomechanical procedures and applications.

Mechanical behavior of infrapatellar fat pad of patients affected by osteoarthritis

The infrapatellar fat pad (IFP) is an adipose tissue present in the knee that lies between the patella, femur, meniscus and tibia, filling the space between these structures. IFP facilitates the distribution of the synovial fluid and may act to absorb impulsive actions generated through the joint. IFP in osteoarthritis (OA) pathology undergoes structural changes characterized by inflammation, hypertrophy and fibrosis. The aim of the present study is to analyze the mechanical behavior of the IFP in patients affected by end-stage OA. A specific test fixture was designed and indentation tests were performed on IFP specimens harvested from OA patients who underwent total knee arthroplasty. Experiments allowed to assess the typical features of mechanical response, such as non-linear stress-strain behavior and time-dependent effects. Results from mechanical experimentations were implemented within the framework of a visco-hyperelastic constitutive theory, with the aim to provide data for computational modelling of OA IFP role in knee mechanics. Initial and final indentation stiffness were calculated for all subjects and statistical results reveled that OA IFP mechanics was not significantly influenced by gender, BMI and sample preparation. OA IFP mechanical behavior was also compared to that of other adipose tissues. OA IFP appeared to be a stiffer adipose tissue compared to subcutaneous, visceral adipose tissues and heel fat pads. It is reasonable that fibrosis induces a modification of the tissue destabilizing the normal distribution of forces in the joint during movement, causing a worsening of the disease.

Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle

Interventional neuroradiology is characterized by engineering- and experience-driven device development with design improvements every few months. However, clinical validation of these new devices requires lengthy and expensive randomized controlled trials. This contribution proposes a machine learning-based in silico study design to evaluate new devices more quickly with a small sample size. Acute diffusion- and perfusion-weighted MRI, segmented one-week follow-up imaging, and clinical variables were available for 90 acute ischemic stroke patients. Three treatment option-specific random forest models were trained to predict the one-week follow-up lesion segmentation for (1) patients successfully recanalized using intra-arterial mechanical thrombectomy, (2) patients successfully recanalized using intravenous thrombolysis, and (3) non-recanalizing patients as an analogue for conservative treatment for each patient in the sample, independent of the true group membership. A repeated-measures analysis of the three predicted follow-up lesions for each patient revealed significantly larger lesions for the non-recanalizing group compared to the successful intravenous thrombolysis treatment group, which in turn showed significantly larger lesions compared to the successful mechanical thrombectomy treatment group (p < 0.001). A groupwise comparison of the true follow-up lesions for the three treatment options showed the same trend but did not reach statistical significance (p = 0.19). We conclude that the proposed machine learning-based in silico trial design leads to clinically feasible results and can support new efficacy studies by providing additional power and potential early intermediate results.

Computational methods for the investigation of ski boots ergonomics

Ski boots are known to cause vasoconstriction in the wearer’s lower limbs and, thus, cause a “cold leg” phenomenon. To address this problem, this work provides a computational framework for analysing interactions between the ski boot and the lower limb. The geometry of the lower limb was derived from magnetic resonance imaging and computed tomography techniques and anthropometric data. The geometry of the ski boot shell was obtained by means of three-dimensional computer aided design models from a manufacturer. Concerning the ski boot liner, laser scanning techniques were implemented to capture the geometry of each layer. The mechanical models of the ski boot and the lower limb were identified and validated by means of coupled experimental investigations and computational analyses. The computational models were exploited to simulate the buckling process and to investigate interaction phenomena between the boot and the lower limb. Similarly, experimental activities were performed to further analyse the buckling phenomena. The obtained computational and experimental results were compared regarding both interaction pressure and displacements between the buckle and the corresponding buckle hooks. These comparisons provided reasonable agreement (mean value of discrepancy between the model and mean experimental results in the tibial region: 20%), underlining the model’s capability to correctly interpret results from experimental measurements. Results identified the critical areas of the leg, such as the tibial region, the calcaneal region of the foot and the anterior sole, which may suffer the most due to the hydrostatic pressure and compressive strain exerted on them. The results highlight that computational methods allow investigation of the interaction phenomena between the lower leg and ski boot, potentially providing an effective framework for a more comfortable and ergonomic design of ski boots.

Design of a force-measuring setup for colorectal compression anastomosis and first ex-vivo results

Purpose

The introduction of novel endoscopic instruments is essential to reduce trauma in visceral surgery. However, endoscopic device development is hampered by challenges in respecting the dimensional restrictions, due to the narrow access route, and by achieving adequate force transmission. As the overall goal of our research is the development of a patient adaptable, endoscopic anastomosis manipulator, biomechanical and size-related characterization of gastrointestinal organs are needed to determine technical requirements and thresholds to define functional design and load-compatible dimensioning of devices.

Methods

We built an experimental setup to measure colon tissue compression piercing forces. We tested 54 parameter sets, including variations of three tissue fixation configurations, three piercing body configurations (four, eight, twelve spikes) and insertion trajectories of constant velocities (5 mms⁻¹, 10 mms⁻¹,15 mms⁻¹) and constant accelerations (5 mms⁻², 10 mms⁻², 15 mms⁻²) each in 5 samples. Furthermore, anatomical parameters (lumen diameter, tissue thickness) were recorded.

Results

There was no statistically significant difference in insertion forces neither between the trajectory groups, nor for variation of tissue fixation configurations. However, we observed a statistically significant increase in insertion forces for increasing number of spikes. The maximum mean peak forces for four, eight and twelve spikes were 6.4 ± 1.5 N, 13.6 ± 1.4 N and 21.7 ± 5.8 N, respectively. The 5th percentile of specimen lumen diameters and pierced tissue thickness were 24.1 mm and 2.8 mm, and the 95th percentiles 40.1 mm and 4.8 mm, respectively.

Conclusion

The setup enabled reliable biomechanical characterization of colon material, on the base of which design specifications for an endoscopic anastomosis device were derived. The axial implant closure unit must enable axial force transmission of at least 28 N (22 ± 6 N). Implant and applicator diameters must cover a range between 24 and 40 mm, and the implant gap, compressing anastomosed tissue, between 2 and 5 mm.