Daniel Simões Lopes

Finite element simulations of a hip joint with femoroacetabular impingement.

In this study, a three-dimensional finite element (FE) model based on the specific anatomy of a patient presenting a femoroacetabular impingement of the ‘cam’-type is developed. The FE meshes of the structures of interest are obtained from arthrographic magnetic resonance images. All soft tissues are considered linear elastic and isotropic, and the bones were assumed rigid. A compression of the femur on the acetabular cavity as well as flexural movements and internal rotations are applied. Stresses and contact pressures are evaluated in this patient-specific model in order to better interpret the mechanism of aggression of the femoral and acetabular cartilages. The corresponding results are presented and discussed. The values obtained for the contact pressures are similar to those reported by other models based on idealised geometries. An FE analysis of a non-cam hip is also performed for comparison with the pathological case.

Finite element studies of the mechanical behaviour of the diaphragm in normal and pathological cases

The diaphragm is a muscular membrane separating the abdominal and thoracic cavities, and its motion is directly linked to respiration. In this study, using data from a 59-year-old female cadaver obtained from the Visible Human Project, the diaphragm is reconstructed and, from the corresponding solid object, a shell finite element mesh is generated and used in several analyses performed with the ABAQUS 6.7 software. These analyses consider the direction of the muscle fibres and the incompressibility of the tissue. The constitutive model for the isotropic strain energy as well as the passive and active strain energy stored in the fibres is adapted from Humphreys model for cardiac muscles. Furthermore, numerical results for the diaphragmatic floor under pressure and active contraction in normal and pathological cases are presented.

A superellipsoid-plane model for simulating foot-ground contact during human gait.

Abstract Musculoskeletal models and forward dynamics simulations of human movement often include foot–ground interactions, with the foot–ground contact forces often determined using a constitutive model that depends on material properties and contact kinematics. When using soft constraints to model the foot–ground interactions, the kinematics of the minimum distance between the foot and planar ground needs to be computed. Due to their geometric simplicity, a considerable number of studies have used point–plane elements to represent these interacting bodies, but few studies have provided comparisons between point contact elements and other geometrically based analytical solutions. The objective of this work was to develop a more general-purpose superellipsoid–plane contact model that can be used to determine the three-dimensional foot–ground contact forces. As an example application, the model was used in a forward dynamics simulation of human walking. Simulation results and execution times were compared with a point-like viscoelastic contact model. Both models produced realistic ground reaction forces and kinematics with similar computational efficiency. However, solving the equations of motion with the surface contact model was found to be more efficient (~18% faster), and on average numerically ~37% less stiff. The superellipsoid–plane elements are also more versatile than point-like elements in that they allow for volumetric contact during three-dimensional motions (e.g. rotating, rolling, and sliding). In addition, the superellipsoid–plane element is geometrically accurate and easily integrated within multibody simulation code. These advantages make the use of superellipsoid–plane contact models in musculoskeletal simulations an appealing alternative to point-like elements.

Geometric Detection Algorithms for Cavities on Protein Surfaces in Molecular Graphics: A Survey

Detecting and analysing protein cavities provides significant information about active sites for biological processes (e.g. protein–protein or protein–ligand binding) in molecular graphics and modelling. Using the three‐dimensional (3D) structure of a given protein (i.e. atom types and their locations in 3D) as retrieved from a PDB (Protein Data Bank) file, it is now computationally viable to determine a description of these cavities. Such cavities correspond to pockets, clefts, invaginations, voids, tunnels, channels and grooves on the surface of a given protein. In this work, we survey the literature on protein cavity computation and classify algorithmic approaches into three categories: evolution‐based, energy‐based and geometry‐based. Our survey focuses on geometric algorithms, whose taxonomy is extended to include not only sphere‐, grid‐ and tessellation‐based methods, but also surface‐based, hybrid geometric, consensus and time‐varying methods. Finally, we detail those techniques that have been customized for GPU (graphics processing unit) computing.

A Survey on Implicit Surface Polygonization

Implicit surfaces (IS) are commonly used in image creation, modeling environments, modeling objects, and scientific data visualization. In this article, we present a survey of different techniques for fast visualization of IS. The main classes of visualization algorithms are identified along with the advantages of each in the context of the different types of IS commonly used in computer graphics. We focus closely on polygonization methods, as they are the most suited to fast visualization. Classification and comparison of existing approaches are presented using criteria extracted from current research. This enables the identification of the best strategies according to the number of specific requirements, such as speed, accuracy, quality, or stylization.

A finite element model to simulate femoral fractures in calves: testing different polymers for intramedullary interlocking nails.

Objective To verify if the finite element method can correctly estimate the performance of polyacetal and polyamide 6 intramedullary nails in stabilizing a femoral fracture in calves and to estimate the performance of a polypropylene nail in same conditions. Study Design Computational and experimental study. Sample Population Finite element models (FEMs). Methods Based on a 3-dimensional finite element method (FEM) of the femoral diaphysis, 3 models were constructed to simulate an oblique simple fracture stabilized by an intramedullary nail composed of 1 of 3 distinct polymers. Models were tested under 6 loading conditions that simulated a static calf or a calf in different walking phases. Maximum bone and implant stresses were compared to yield and rupture stresses of specific materials. Results Under static conditions, all polymers were resistant to critical deformation and rupture because maximum von Mises stresses were lower than the respective yield and rupture stresses. However, during walking, maximum stresses exceeded the yield and rupture limits of the polymers, in agreement with a previous in vivo study, which used polyacetal and polyamide nails. Conclusions FEM correctly estimated that polyacetal and polyamide 6 nails would fail to immobilize an oblique femoral diaphyseal fracture in calves that were allowed to walk freely during the early postoperative period. FEM can be useful in the development of new bovine orthopedic devices.OBJECTIVE To verify if the finite element method can correctly estimate the performance of polyacetal and polyamide 6 intramedullary nails in stabilizing a femoral fracture in calves and to estimate the performance of a polypropylene nail in same conditions. STUDY DESIGN Computational and experimental study. SAMPLE POPULATION Finite element models (FEMs). METHODS Based on a 3-dimensional finite element method (FEM) of the femoral diaphysis, 3 models were constructed to simulate an oblique simple fracture stabilized by an intramedullary nail composed of 1 of 3 distinct polymers. Models were tested under 6 loading conditions that simulated a static calf or a calf in different walking phases. Maximum bone and implant stresses were compared to yield and rupture stresses of specific materials. RESULTS Under static conditions, all polymers were resistant to critical deformation and rupture because maximum von Mises stresses were lower than the respective yield and rupture stresses. However, during walking, maximum stresses exceeded the yield and rupture limits of the polymers, in agreement with a previous in vivo study, which used polyacetal and polyamide nails. CONCLUSIONS FEM correctly estimated that polyacetal and polyamide 6 nails would fail to immobilize an oblique femoral diaphyseal fracture in calves that were allowed to walk freely during the early postoperative period. FEM can be useful in the development of new bovine orthopedic devices.

Bone remodelling analysis of a bovine femur for a veterinary implant design

The response of bovine bone to the presence of an implant is analysed with the aim of simulating bone remodelling in a developing model of a polymeric intramedullary interlocking nail for veterinary use. A 3-D finite element model of the femur diaphysis is built based on computed tomography images and using a CAD-based modelling pipeline. The bone remodelling process after the surgery is analysed and compared with the healthy bone. The remodelling law assumes that bone adapts to the mechanical environment. For the analyses a consistent set of loads is determined for the bovine walk cycle. The remodelling results reproduce the morphologic features of bone and provide evidence of the difference on the bone behaviour when comparing metallic and polymeric nails. Our findings indicate that an intramedullary polymeric nail has the advantage over the metallic one of improving long-term bone healing and possibly avoiding the need of the implant removal.

Geometric Detection Algorithms for Cavities on Protein Surfaces in Molecular Graphics: A Survey: Detection Algorithms for Cavities

Detecting and analysing protein cavities provides significant information about active sites for biological processes (e.g. protein–protein or protein–ligand binding) in molecular graphics and modelling. Using the three‐dimensional (3D) structure of a given protein (i.e. atom types and their locations in 3D) as retrieved from a PDB (Protein Data Bank) file, it is now computationally viable to determine a description of these cavities. Such cavities correspond to pockets, clefts, invaginations, voids, tunnels, channels and grooves on the surface of a given protein. In this work, we survey the literature on protein cavity computation and classify algorithmic approaches into three categories: evolution‐based, energy‐based and geometry‐based. Our survey focuses on geometric algorithms, whose taxonomy is extended to include not only sphere‐, grid‐ and tessellation‐based methods, but also surface‐based, hybrid geometric, consensus and time‐varying methods. Finally, we detail those techniques that have been customized for GPU (graphics processing unit) computing.

VRRRRoom: Virtual Reality for Radiologists in the Reading Room

Reading room conditions such as illumination, ambient light, human factors and display luminance, play an important role on how radiologists analyze and interpret images. Indeed, serious diagnostic errors can appear when observing images through everyday monitors. Typically, these occur whenever professionals are ill-positioned with respect to the display or visualize images under improper light and luminance conditions. In this work, we show that virtual reality can assist radiodiagnostics by considerably diminishing or cancel out the effects of unsuitable ambient conditions. Our approach combines immersive head-mounted displays with interactive surfaces to support professional radiologists in analyzing medical images and formulating diagnostics. We evaluated our prototype with two senior medical doctors and four seasoned radiology fellows. Results indicate that our approach constitutes a viable, flexible, portable and cost-efficient option to traditional radiology reading rooms.

Shape Analysis of the Femoral Head: A Comparative Study Between Spherical, (Super)Ellipsoidal, and (Super)Ovoidal Shapes

In this work, MacConaills classification that the articular surface of the femoral head is better represented by ovoidal shapes rather than purely spherical shapes is computationally tested. To test MacConaills classification, a surface fitting framework was developed to fit spheres, ellipsoids, superellipsoids, ovoids, and superovoids to computed tomography (CT) data of the femoral proximal epiphysis. The framework includes several image processing and computational geometry techniques, such as active contour segmentation and mesh smoothing, where implicit surface fitting is performed with genetic algorithms. By comparing the surface fitting error statistics, the results indicate that (super)ovoids fit femoral articular surfaces better than spherical or (super)ellipsoidal shapes.