A Comparative Analysis of Different Treatments for Distal Femur Fractures using the Finite Element Method

The main objective of this work is the evaluation, by means of the finite element method (FEM) of the mechanical stability and long-term microstructural modifications in bone induced to three different kinds of fractures of the distal femur by three types of implants: the Condyle Plate, the less invasive stabilization system plate (LISS) and the distal femur nail (DFN). The displacement and the stress distributions both in bone and implants and the internal bone remodelling process after fracture and fixation are obtained and analysed by computational simulation. The main conclusions of this work are that distal femoral fractures can be treated correctly with the Condyle Plate, the LISS plate and the DFN. The stresses both in LISS and DFN implant are high especially around the screws. When respect to remodelling, the LISS produces an important resorption in the fractured region, while the other two implants do not strongly modify bone tissue microstructure.     

Anisotropic adaptive finite element method for modelling blood flow

In this study, we present an adaptive anisotropic finite element method (FEM) and demonstrate how computational efficiency can be increased when applying the method to the simulation of blood flow in the cardiovascular system. We use the SUPG formulation for the transient 3D incompressible Navier–Stokes equations which are discretised by linear finite elements for both the pressure and the velocity field.

Given the pulsatile nature of the flow in blood vessels we have pursued adaptivity based on the average flow over a cardiac cycle. Error indicators are derived to define an anisotropic mesh metric field. Mesh modification algorithms are used to anisotropically adapt the mesh according to the desired size field. We demonstrate the efficiency of the method by first applying it to pulsatile flow in a straight cylindrical vessel and then to a porcine aorta with a stenosis bypassed by a graft. We demonstrate that the use of an anisotropic adaptive FEM can result in an order of magnitude reduction in computing time with no loss of accuracy compared to analyses obtained with uniform meshes.     

Search for critical loading condition of the spine–A meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2–L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

Tissue Aspiration and Inverse Finite Element Characterisation of Soft Tissues

In this work we examined the determination of soft tissue parameters via tissue aspiration experiments and inverse finite element characterisation. An aspiration tube was put against the target tissue. The deformation of the tissue inside the tube caused by weak suction was tracked with a video based system. A strain energy function was employed to model the elastic behaviour of soft tissue and viscoelasticity was accounted for by means of a quasi-linear viscoelastic formulation. Friction between the aspiration tube and the aspirated tissue was included in the model. Based on the assumed material model, an optimal set of material parameters was found, in order to best fit the experimental data obtained from ex-vivo experiments on pig kidney cortex. The inverse method resulted in robust determination of the unknown material parameters.     

Simulation and study of the geometric parameters in the inguinal area and the genesis of inguinal hernias

We are interested in studying the genesis of a very common pathology: the human inguinal hernia. How the human inguinal hernia appears is not definitively clear, but it is accepted that it is caused by a combination of mechanical and biochemical alterations, and that muscular simulation plays an important role in this. This study proposes a model to explain how some physical parameters affect the ability to simulate the region dynamically and how these parameters are involved in generating inguinal hernias. We are particularly interested in understanding the mechanical alterations in the inguinal region because little is known about them or how they behave dynamically. Our model corroborates the most important theories regarding the generation of inguinal hernias and is an initial approach to numerically evaluating this affection.     

Characterisation and simulation of an active microvalve for glaucoma

Glaucoma drainage device (GDD) has the potential to eliminate hypotony but still suffers from poor flow control and fibrosis. The ideal shunt should change its hydraulic resistance to achieve the desired intraocular pressure (IOP). In this study, the characterisation of a preliminary design of a new GDD is presented. This is activated by means of a diaphragm, which is actuated by conducting polymers. The valve can be manufactured employing microelectromechanical system technology by soft lithography. The characterisation process is performed by numerical simulation using the finite element method, considering the coupling between the fluid and the structure (diaphragm) obtaining the hydraulic resistance for several positions of the diaphragm. To analyse the hydraulic system of the microvalve implanted in a human eye, an equivalent circuit model was used. The parameters of the equivalent circuit model were obtained from numerical simulation. The hydraulic resistance of the designed GDD varies in the range of 13.08–0.36 mmHg min/μl compared with 3.38–0.43 mmHg min/μl for the Ahmed valve. The maximum displacement of the diaphragm in the vertical direction is 18.9 μm, and the strain in the plane is 2%. The proposed preliminary design allows to control the IOP by varying the hydraulic resistance in a greater range than the existing passive valves, and the numerical simulation facilitates the characterisation and the improvement of the design before its construction, reducing time and costs.     

Preliminary simulation model toward the study of the effects caused by different mandibular advancement devices in OSAS treatment

Abstract

The paper aims to evaluate the effects caused by a Mandibular Advancement Device (MAD) for Obstructive Sleep Apnoea Syndrome (OSAS) treatment. This study is based on Finite Element Method (FEM) for evaluating the load distribution on temporomandibular joint, especially on the mandibular condyle and disc, and on periodontal ligaments. The stress values on condyle and periodontal ligaments lead authors to consider MAD a safe procedure even for a long period. The obtained results also show the relationship between MAD material and load distribution at the periodontal ligaments. The paper is a step toward future analyses for studying and comparing the effects of MAD features, such as material, shape and dimensions, in order to allow the clinician prescribing the most fitting device.     

The use of preserved tissue in finite element modelling of fresh ovine vertebral behaviour

The aim of this study was to investigate whether the predicted finite element (FE) stiffness of vertebral bone is altered when using images of preserved rather than fresh tissue to generate specimen-specific FE models. Fresh ovine vertebrae were used to represent embalmed (n = 3) and macerated dry-bone (n = 3) specimens and treated accordingly. Specimens were scanned pre- and post-treatment using micro-computed tomography. A constant threshold level derived from these images was used to calculate the respective bone volume fraction (BV/TV) from which the conversion factor validated for fresh tissue was used to determine material properties that were assigned to corresponding FE models. Results showed a definite change in the BV/TV between the fresh and the preserved bone. However, the changes in the predicted FE stiffness were not generally greater than the variations expected from assignment of loading and boundary conditions. In conclusion, images of preserved tissue can be used to generate FE models that are representative of fresh tissue with a tolerable level of error.     

Implementation and validation of thoracic side impact injury prediction metrics in a human body model

This study's purpose was to implement injury metrics into the Total Human Model for Safety (THUMS) mirroring the spinal accelerometers, rib accelerometers and chest band instrumentation from two lateral post-mortem human subject sled test configurations. In both sled configurations, THUMS contacted a flat rigid surface (either a wall or beam) at 6.7 m/s. Sled A maximum simulated wall forces for the thorax, abdomen and pelvis were 7.1, 5.0 and 10.0 kN versus 5.7 ± 0.8, 3.4 ± 1.2 and 6.2 ± 2.7 kN experimentally. Sled B maximum simulated beam forces for the torso and pelvis were 8.0 and 7.6 kN versus 8.5 ± 0.2 and 7.9 ± 2.5 kN experimentally. Quantitatively, force magnitude contributed more to variation between simulated and experimental forces than phase shift. Acceleration-based injury metrics were within one standard deviation of experimental means except for the lower spine in the rigid wall sled test. These validated metrics will be useful for quantifying occupant loading conditions and calculating injury risks in various loading configurations.