Development and evaluation of a new methodology for the fast generation of patient-specific Finite Element models of the buttock for sitting-acquired deep tissue injury prevention

Current poisson's ratio values of finite element models are too low to consider soft tissues nearly-incompressible: illustration on the human heel region

Tissues' nearly incompressibility was well reported in the literature but little effort has been made to compare volume variations computed by simulations with in vivo measurements. In this study, volume changes of the fat pad during controlled indentations of the human heel region were estimated from segmented medical images using digital volume correlation. The experiment was reproduced using finite element modelling with several values of Poisson's ratio for the fat pad, from 0.4500 to 0.4999. A single value of Poisson's ratio could not fit all the indentation cases. Estimated volume changes were between 0.9% - 11.7%.

Definition and evaluation of a finite element model of the human heel for diabetic foot ulcer prevention under shearing loads

Diabetic foot ulcers are triggered by mechanical loadings applied to the surface of the plantar skin. Strain is considered to play a crucial role in relation to ulcer etiology and can be assessed by Finite Element (FE) modeling. A difficulty in the generation of these models is the choice of the soft tissue material properties. In the literature, many studies attempt to model the behavior of the heel soft tissues by implementing constitutive laws that can differ significantly in terms of mechanical response. Moreover, current FE models lack of proper evaluation techniques that could estimate their ability to simulate realistic strains. In this article, we propose and evaluate a FE model of the human heel for diabetic foot ulcer prevention. Soft tissue constitutive laws are defined through the fitting of experimental stretch-stress curves published in the literature. The model is then evaluated through Digital Volume Correlation (DVC) based on non-rigid 3D Magnetic Resonance Image Registration. The results from FE analysis and DVC show similar strain locations in the fat pad and strain intensities according to the type of applied loads. For additional comparisons, different sets of constitutive models published in the literature are applied into the proposed FE mesh and simulated with the same boundary conditions. In this case, the results in terms of strains show great diversity in locations and intensities, suggesting that more research should be developed to gain insight into the mechanical properties of these tissues.

PIPER adult comfort: an open-source full body human body model for seating comfort assessment and its validation under static loading conditions

Introduction: In this paper we introduce an adult-sized FE full-body HBM for seating comfort assessments and present its validation in different static seating conditions in terms of pressure distribution and contact forces. Methods: We morphed the PIPER Child model into a male adult-sized model with the help of different target sources including his body surface scans, and spinal and pelvic bone surfaces and an open sourced full body skeleton. We also introduced soft tissue sliding under the ischial tuberosities (ITs). The initial model was adapted for seating applications with low modulus soft tissue material property and mesh refinements for buttock regions, etc. We compared the contact forces and pressure-related parameters simulated using the adult HBM with those obtained experimentally from the person whose data was used for the model development. Four seat configurations, with the seat pan angle varying from 0° to 15° and seat-to-back angle fixed at 100°, were tested. Results: The adult HBM could correctly simulate the contact forces on the backrest, seat pan, and foot support with an average error of less than 22.3 N and 15.5 N in the horizontal and vertical directions, which is small considering the body weight (785 N). In terms of contact area, peak, and mean pressure, the simulation matched well with the experiment for the seat pan. With soft tissue sliding, higher soft tissue compression was obtained in agreement with the observations from recent MRI studies. Discussion: The present adult model could be used as a reference using a morphing tool as proposed in PIPER. The model will be published openly online as part of the PIPER open-source project (www.PIPER-project.org) to facilitate its reuse and improvement as well as its specific adaptation for different applications.

Could Horse Gait and Induced Pelvic Dynamic Loads in Female Equestrians Be a Risk Factor in Pudendal Neuralgia?

Pudendal Neuralgia (PN) is a rare, debilitating disease caused by damage to the pudendal nerve, which innervates the anus, rectum, perineum, lower urinary tract, and genitalia. Although its etiology remains scientifically unknown, a number of sports practices, including horse-riding, are reported as triggering and/or aggravating factors. The present work summarizes the experimental measurements of the contact pressure at the interface between the rider and saddle, for a population of 12 experienced female riders. These tests reveal that dynamic horseback-riding leads to high levels of peak pressures in the perineal region, which confirms that the practice of equine sports may cause neuropathologies such as PN. All collected data will be used as boundary conditions in a future numerical 3D model aimed at locating the possible areas of pudendal nerve crushing.

Determining frictional properties of pants and cushion cover materials using human soft tissue and a rigid sled and how they affect seated shear forces

Shear forces on the buttocks while seated are directly linked to friction, yet the frictional properties at the seat interface are unknown. Shear forces are one of the factors related to increase risk of pressure injury formation. The goals of this study included determining coefficients of friction between three cushion covers and two clothing fabrics using a mechanical system as well as human participants and to evaluate the impact of the cushion covers on shear loading on the buttocks while seated. A chair with separate seat pan tilt and back recline movements was built and instrumented with reflective markers and a load cell. A motion capture system and load cell were used to determine the angles of seat pan tilt at which the sled and participants started sliding, as well as shear forces at three recline angles for three cushion covers (vinyl, one-layer nylon, and two-layer nylon). Results showed the vinyl and two-layer nylon cushion covers respectively had the largest and smallest coefficients of friction for both pants materials. The coefficients of friction calculated with the human participants and rigid sled were within 10% of each other, demonstrating similar results. Further, increasing back recline increased shear load on the buttocks, while the two-layer nylon cover reduced shear forces seen on the buttocks. This work furthers the understanding of shear loading on the buttocks, will aid in the protocols for reducing pressure injuries, and suggests that coefficients of friction found using rigid bodies may be applied to deformable bodies.

Mechanical analysis of deep tissue injury during sitting in patients with spinal cord injury via parametric finite element model

Spinal cord injury patients are prone to develop deep tissue injury because of long-term mechanical load. However, there is a lack of statistical research on the influence of tissue characteristics on the internal mechanical state of soft tissue. This study aimed to investigate the influence of tissue characteristics on the internal mechanical state of buttock in spinal cord injury patients. A three-dimensional reference buttock model was established and a visualization program was generated to modify the parameter values. Through changing the muscle atrophy, body mass index and the radius of curvature of the ischial tuberosity, 96 different model variants were simulated and validated in this study. With body mass index increasing from 16 to 40, the principal shear stress was 10.4 times principal compressive stress, the maximum shear strain and the max cluster volume increased by 1.2 (P < 0.001) and 8.8 times (P < 0.001), respectively. The interaction between BMI and muscle atrophy was significant when BMI was greater than or equal to 22.5 kg/m². In all BMI stages, when the radius of curvature of the ischial tuberosity was 19 mm, the internal stress of the tissue always occupies the highest value. The results demonstrate that body mass index is the most important factor affecting the risk of buttock deep tissue injury. This study provides insights into investigation of inter-individual factors influencing the soft tissue response and assessment of deep tissue injury risk during sitting.

Société de Biomécanique Young Investigator Award 2021: Numerical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue in the context of pressure ulcer prevention

BACKGROUND

Pressure-induced tissue strain is one major pathway for Pressure Ulcer development and, especially, Deep Tissue Injury. Biomechanical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue is therefore essential. In the literature, a viscoelastic formulation is generally assumed for the experimental characterization of skeletal muscles, with the limitation that the underlying physical mechanisms that give rise to the time dependent stress-strain behaviour are not known. The objective of this study is to explore the capability of poroelasticity to reproduce the apparent viscoelastic behaviour of passive muscle tissue under confined compression.

METHODS

Experimental stress-relaxation response of 31 cylindrical porcine samples tested under fast and slow confined compression by Vaidya and collaborators were used. An axisymmetric Finite Element model was developed in ABAQUS and, for each sample a one-to-one inverse analysis was performed to calibrate the specimen-specific constitutive parameters, namely, the drained Young's modulus, the void ratio, hydraulic permeability, the Poisson's ratio, the solid grain's and fluid's bulk moduli.

FINDINGS

The peak stress and consolidation were recovered for most of the samples (N=25) by the poroelastic model (normalised root-mean-square error ≤0.03 for fast and slow confined compression conditions).

INTERPRETATION

The strength of the proposed model is its fewer number of variables (N=6 for the proposed poroelastic model versus N=18 for the viscohyperelastic model proposed by Vaidya and collaborators). The incorporation of poroelasticity to clinical models of Pessure Ulcer formation could lead to more precise and mechanistic explorations of soft tissue injury risk factors.

A review of foot finite element modelling for pressure ulcer prevention in bedrest: Current perspectives and future recommendations

Pressure ulcers (PUs) are a major public health challenge, having a significant impact on healthcare service and patient quality of life. Computational biomechanical modelling has enhanced PU research by facilitating the investigation of pressure responses in subcutaneous tissue and skeletal muscle. Extensive work has been undertaken on PUs on patients in the seated posture, but research into heel ulcers has been relatively neglected. The aim of this review was to address the key challenges that exist in developing an effective FE foot model for PU prevention and the confusion surrounding the wide range of outputs reported. Nine FE foot studies investigating heel ulcers in bedrest were identified and reviewed. Six studies modelled the posterior part of the heel, two included the calf and foot, and one modelled the whole body. Due to the complexity of the foot anatomy, all studies involved simplification or assumptions regarding parts of the foot structure, boundary conditions and material parameters. Simulations aimed to understand better the stresses and strains exhibited in the heel soft tissues of the healthy foot. The biomechanical properties of soft tissue derived from experimental measurements are critical for developing a realistic model and consequently guiding clinical decisions. Yet, little to no validation was reported in each of the studies. If FE models are to address future research questions and clinical applications, then sound verification and validation of these models is required to ensure accurate conclusions and prediction of patient outcomes. Recommendations and considerations for future FE studies are therefore proposed.

Combining Freehand Ultrasound-Based Indentation and Inverse Finite Element Modeling for the Identification of Hyperelastic Material Properties of Thigh Soft Tissues

Finite element analysis (FEA) is a numerical modeling tool vastly employed in research facilities to analyze and predict load transmission between the human body and a medical device, such as a prosthesis or an exoskeleton. Yet, the use of finite element modeling (FEM) in a framework compatible with clinical constraints is hindered by, among others, heavy and time-consuming assessments of material properties. Ultrasound (U.S.) imaging opens new and unique opportunities for the assessment of in vivo material properties of soft tissues. Confident of these advances, a method combining a freehand U.S. probe and a force sensor was developed in order to compute the hyperelastic constitutive parameters of the soft tissues of the thigh in both relaxed (R) and contracted (C) muscles' configurations. Seven asymptomatic subjects were included for the experiment. Two operators in each configuration performed the acquisitions. Inverse FEM allowed for the optimization of an Ogden's hyperelastic constitutive model of soft tissues of the thigh in large displacement. The mean shear modulus identified for configurations R and C was, respectively, 3.2 ± 1.3 kPa and 13.7 ± 6.5 kPa. The mean alpha parameter identified for configurations R and C was, respectively, 10 ± 1 and 9 ± 4. An analysis of variance showed that the configuration had an effect on constitutive parameters but not on the operator.

Is a simplified Finite Element model of the gluteus region able to capture the mechanical response of the internal soft tissues under compression?

BACKGROUND

Internal soft tissue strains have been shown to be one of the main factors responsible for the onset of Pressure Ulcers and to be representative of its risk of development. However, the estimation of this parameter using Finite Element (FE) analysis in clinical setups is currently hindered by costly acquisition, reconstruction and computation times. Ultrasound (US) imaging is a promising candidate for the clinical assessment of both morphological and material parameters.

METHOD

The aim of this study was to investigate the ability of a local FE model of the region beneath the ischium with a limited number of parameters to capture the internal response of the gluteus region predicted by a complete 3D FE model. 26 local FE models were developed, and their predictions were compared to those of the patient-specific reference FE models in sitting position.

FINDINGS

A high correlation was observed (R = 0.90, p-value < 0.01). A sensitivity analysis showed that the most influent parameters were the mechanical behaviour of the muscle tissues, the ischium morphology and the external mechanical loading.

INTERPRETATION

Given the progress of US for capturing both morphological and material parameters, these results are promising because they open up the possibility to use personalised simplified FE models for risk estimation in daily clinical routine.