The mechanisms giving rise to vein tortuosity, which is often associated with varicosis, are poorly understood. Recent works suggest that significant biological changes in the wall of varicose veins ...may precede the mechanical aspects of the disease. To test the hypothesis of tortuosity being a consequence of these changes, a Finite Element model was developed based on previous experimental work on vein buckling. The model was then used to evaluate the effect of alterations of the mechanical behavior of the wall on tortuosity onset and severity. The results showed that increasing anisotropy toward the circumferential direction promotes tortuosity. An increase in wall stiffness tends to decrease the level of tortuosity but interestingly, if the vein segment is little or not pre-stretched such increase will not prevent, or it will even promote, the onset of tortuosity. These results provide additional arguments supporting the hypothesis of tortuosity being the consequence of biologically-induced changes in the varicose vein wall. Based on a 3D model of the leg and in vivo identification of the material properties of varicose veins, a clinical validation of these findings is being developed.
Display omitted
•We test the hypothesis that vein tortuosity results from wall mechanical degradation.•Finite Element model of venous tortuosity onset under internal pressure.•Identification of anisotropic mechanical behavior from previous inflation test data.•Alteration of the wall mechanical behavior has significant effects on tortuosity.
Several Finite Element (FE) models of the pelvis have been developed to comprehensively assess the onset of pathologies and for clinical and industrial applications. However, because of the ...difficulties associated with the creation of subject-specific FE mesh from CT scan and MR images, most of the existing models rely on the data of one given individual. Moreover, although several fast and robust methods have been developed for automatically generating tetrahedral meshes of arbitrary geometries, hexahedral meshes are still preferred today because of their distinct advantages but their generation remains an open challenge. Recently, approaches have been proposed for fast 3D reconstruction of bones based on X-ray imaging. In this study, we adapted such an approach for the fast and automatic generation of all-hexahedral subject-specific FE models of the pelvis based on the elastic registration of a generic mesh to the subject-specific target in conjunction with element regularity and quality correction. The technique was successfully tested on a database of 120 3D reconstructions of pelvises from biplanar X-ray images. For each patient, a full hexahedral subject-specific FE mesh was generated with an accurate surface representation.
BackgroundPressure-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.MethodsExperimental 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.FindingsThe 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).InterpretationThe 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.
Transmission of loads between the prosthetic socket and the residual limb is critical for the comfort and walking ability of people with transfemoral amputation. This transmission is mainly ...determined by the socket tightening, muscle forces, and socket ischial support. However, numerical investigations of the amputated gait, using modeling approaches such as MusculoSkeletal (MSK) modeling, ignore the weight-bearing role of the ischial support. This simplification may lead to errors in the muscle force estimation.
This study aims to propose a MSK model of the amputated gait that accounts for the interaction between the body and the ischial support for the estimation of the muscle forces of 13 subjects with unilateral transfemoral amputation.
Contrary to previous studies on the amputated gait which ignored the interaction with the ischial support, here, the contact on the ischial support was included in the external loads acting on the pelvis in a MSK model of the amputated gait.
Including the ischial support induced an increase in the activity of the main abductor muscles, while adductor muscles' activity was reduced. These results suggest that neglecting the interaction with the ischial support leads to erroneous muscle force distribution considering the gait of people with transfemoral amputation. Although subjects with various bone geometries, particularly femur lengths, were included in the study, similar results were obtained for all subjects.
Eventually, the estimation of muscle forces from MSK models could be used in combination with finite element models to provide quantitative data for the design of prosthetic sockets.
OBJECTIVE: 2D Ultrasound (US) imaging has been recently investigated as a more accessible alternative to 3D Magnetic Resonance Imaging (MRI) for the estimation of soft issue motion under external ...mechanical loading. In the context of pressure ulcer prevention, the aim of this pilot MRI study was to design an experiment to characterize the sacral soft tissue motion under a controlled mechanical loading. Such an experiment targeted the estimation of the discrepancy between tissue motion assessed using a 2D imaging modality (echography) versus tissue motion assessed using a (reference) 3D imaging modality (MRI).METHODS: One healthy male volunteer participated in the study. An MRI-compatible custom-made setup was designed and used to load the top region of the sacrum with a 3D-printed copy of the US transducer. Five MR images were collected, one in the unloaded and four in the different loaded configurations (400-1200 g). Then, a 3D displacement field for each loading configuration was extracted based on the results of digital volume correlation. Tissue motion was separated into the X, Y, Z directions of the MRI coordinate system and the ratios between the out-of-plane and in-plane components were assessed for each voxel of the selected region of interest.RESULTS: Ratios between the out-of-plane and in-plane displacement components were higher than 0.6 for more than half of the voxels in the region of interest for all load cases and higher than 1 for at least quarter of the voxels when loads of 400-800 g were used.CONCLUSION: The out-of-ultrasound-plane tissue displacement was not negligible, therefore 2D US imaging should be used with caution for the evaluation of the tissue motion in the sacrum region. The 3D US modality should be further investigated for this application.
The Finite Element Method (FEM) is a well-established procedure for computing approximate solutions to deterministic engineering problems described by partial differential equations. FEM produces ...discrete approximations of the solution with a discretisation error that can be an be quantified with \emph{a posteriori} error estimates. The practical relevance of error estimates for biomechanics problems, especially for soft tissue where the response is governed by large strains, is rarely addressed. In this contribution, we propose an implementation of \emph{a posteriori} error estimates targeting a user-defined quantity of interest, using the Dual Weighted Residual (DWR) technique tailored to biomechanics. The proposed method considers a general setting that encompasses three-dimensional geometries and model non-linearities, which appear in hyperelastic soft tissues. We take advantage of the automatic differentiation capabilities embedded in modern finite element software, which allows the error estimates to be computed generically for a large class of models and constitutive laws. First we validate our methodology using experimental measurements from silicone samples, and then illustrate its applicability for patient-specific computations of pressure ulcers on a human heel.
Background: Physiologic aging is associated with loss of mobility, sarcopenia, skin atrophy and loss of elasticity. These factors contribute, in the elderly, to the occurrence of a pressure ulcer ...(PU). Brightness mode ultrasound (US) and shear wave elastography (SWE) have been proposed as a patient-specific, bedside, and predictive tool for PU. However, reliability and clinical feasibility in application to the sacral region have not been clearly established.Method: The current study aimed to propose a simple bedside protocol combining US and SWE. The protocol was first tested on a group of 19 healthy young subjects by two operators. The measurements were repeated three times. Eight parameters were evaluated at the medial sacral crest. Intraclass Correlation Coefficient (ICC) was used for reliability assessment and the modified Bland Altman plot analysis for agreement assessment. The protocol was then evaluated for clinical feasibility on a healthy older group of 11 subjects with a mean age of 65 ± 2.4 yrs.Findings: ICC showed poor to good reliability except for skin SWE and hypodermis thickness with an ICC (reported as: mean (95%CI)) of 0.78 (0.50-0.91) and 0.98 (0.95-0.99) respectively. No significant differences were observed between the young and older group except for the muscle Shear Wave Speed (SWS) (respectively 2.11 ± 0.27 m/s vs 1.70 ± 0.17 m/s).Interpretation: This is the first protocol combining US and SWE that can be proposed on a large scale in nursing homes. Reliability, however, was unsatisfactory for most parameters despite efforts to standardize the protocol and measurement definitions. Further studies are needed to improve reliability.
Pressure ulcers are a severe disease affecting patients that are bedridden or in a wheelchair bound for long periods of time. These wounds can develop in the deep layers of the skin of specific parts ...of the body, mostly on heels or sacrum, making them hard to detect in their early stages. Strain levels have been identified as a direct danger indicator for triggering pressure ulcers. Prevention could be possible with the implementation of subject-specific Finite Element (FE) models. However, generation and validation of such FE models is a complex task, and the current implemented techniques offer only a partial solution of the entire problem considering only external displacements and pressures, or cadaveric samples. In this paper, we propose an in vivo solution based on the 3D non-rigid registration between two Magnetic Resonance (MR) images, one in an unloaded configuration and the other deformed by means of a plate or an indenter. From the results of the image registration, the displacement field and subsequent strain maps for the soft tissues were computed. An extensive study, considering different cases (on heel pad and sacrum regions) was performed to evaluate the reproducibility and accuracy of the results obtained with this methodology.The implemented technique can give insight for several applications. It adds a useful tool for better understanding the propagation of deformations in the heel soft tissues that could generate pressure ulcers. This methodology can be used to obtain data on the material properties of the soft tissues to define constitutive laws for FE simulations and finally it offers a promising technique for validating FE models.
Abstract A patient-specific finite-element (FE) model of the human leg is developed to model the stress distribution in and around a vein wall in order to determine the biomechanical response of ...varicose veins to compression treatment. The aim is to investigate the relationship between the local pressure on the soft tissues induced by wearing the compression garment and the development and evolution of varicose veins and various skin-related diseases such as varicose veins and ulcers. Because experimental data on the mechanical properties of healthy superficial veins and varicose veins are scarce in literature, ultrasound images of in vivo varicose veins are acquired and analysed to extract the material constants using Finite Element Model Updating. The decrease in trans-mural pressure, which conditions the effectiveness of compressive treatments, is computed from the simulation results. This constitutes the original added value of the developed model as decrease in trans-mural pressures cannot be assessed experimentally by any other means. Results show that external compression is effective in decreasing the trans-mural pressure, thereby having a positive effect in the control and treatment of vein-related diseases.
Soft biological tissues demonstrate strong time-dependent and strain-rate mechanical behavior, arising from their intrinsic visco-elasticity and fluid-solid interactions (especially at sufficiently ...large time scales). The time-dependent mechanical properties of soft tissues influence their physiological functions and are linked to several pathological processes. Poro-elastic modeling represents a promising approach because it allows the integration of multiscale/multiphysics data to probe biologically relevant phenomena at a smaller scale and embeds the relevant mechanisms at the larger scale. The implementation of multi-phasic flow poro-elastic models however is a complex undertaking, requiring extensive knowledge. The open-source software FEniCSx Project provides a novel tool for the automated solution of partial differential equations by the finite element method. This paper aims to provide the required tools to model the mixed formulation of poro-elasticity, from the theory to the implementation, within FEniCSx. Several benchmark cases are studied. A column under confined compression conditions is compared to the Terzaghi analytical solution, using the L2-norm. An implementation of poro-hyper-elasticity is proposed. A bi-compartment column is compared to previously published results (Cast3m implementation). For all cases, accurate results are obtained in terms of a normalized Root Mean Square Error (RMSE). Furthermore, the FEniCSx computation is found three times faster than the legacy FEniCS one. The benefits of parallel computation are also highlighted.