Finite element modelling approaches for well-ordered porous metallic materials for orthopaedic applications: cost effectiveness and geometrical considerations

The mechanical properties of well-ordered porous materials are related to their geometrical parameters at the mesoscale. Finite element (FE) analysis is a powerful tool to design well-ordered porous materials by analysing the mechanical behaviour. However, FE models are often computationally expensive. This article aims to develop a cost-effective FE model to simulate well-ordered porous metallic materials for orthopaedic applications. Solid and beam FE modelling approaches are compared, using finite size and infinite media models considering cubic unit cell geometry. The model is then applied to compare two unit cell geometries: cubic and diamond. Models having finite size provide similar results than the infinite media model approach for large sample sizes. In addition, these finite size models also capture the influence of the boundary conditions on the mechanical response for small sample sizes. The beam FE modelling approach showed little computational cost and similar results to the solid FE modelling approach. Diamond unit cell geometry appeared to be more suitable for orthopaedic applications than the cubic unit cell geometry.     

Nanoconfined polymers: modelling and simulation approaches

In this article, application of molecular simulation methods for studying molecular pictures of nanoconfined polymers is reviewed and discussed. The simulation methods, covering a range from atomistic to systematically parameterised coarse-grained models, employed in the literature to study nanoconfined polymers are reviewed and their results are compared together. The effect of polymer–surface interactions, surface curvature and surface area on the alteration of polymer structure and dynamics from the unperturbed (bulk) polymer properties are discussed. The length scales over which the surface influences the polymer structure and dynamics and the magnitude of surface effect on dynamics deceleration in the interphase are addressed in terms of different local and global chain properties.     

Comparing parametric solid modelling/reconfiguration,global shape modelling and free-form deformation for the generation of 3D digital models of femurs from X-ray images

At present, computer assisted surgery systems help orthopaedic surgeons both plan and perform surgical procedures. To enable these systems to function, it is crucial to have at one's disposal 3D models of anatomical structures, surgical tools and prostheses (if required). This paper analyses and compares three methods for generating 3D digital models of anatomical structures starting from X-ray images: parametric solid modelling/reconfiguration, global shape modelling and free-form deformation. Seven experiences involving the generation of a femur model were conducted by software developers and different skilled users. These experiences are described in detail and compared at different stages and from different points of view.     

Planning and simulation of minimally-invasive robotic heart surgery

Due to their numerous advantages, mainly in terms of patient benefit, mini-invasive robotically assisted interventions are gaining in importance in various surgical fields. However, this conversion has its own challenges that stem from both its novelty and complexity. In this paper we propose to accompany the surgeons in their transition, by offering an integrated environment that enables them to make better use of this new technology. The proposed system is patient-dependent, and enables the planning, validation, simulation, teaching and archiving of robotically assisted interventions. The approach is illustrated for a coronary bypass graft using the daVinci tele-operated robot.     

Advanced material modelling in numerical simulation of primary acetabular press-fit cup stability

Primary stability of artificial acetabular cups, used for total hip arthroplasty, is required for the subsequent osteointegration and good long-term clinical results of the implant. Although closed-cell polymer foams represent an adequate bone substitute in experimental studies investigating primary stability, correct numerical modelling of this material depends on the parameter selection.

Material parameters necessary for crushable foam plasticity behaviour were originated from numerical simulations matched with experimental tests of the polymethacrylimide raw material. Experimental primary stability tests of acetabular press-fit cups consisting of static shell assembly with consecutively pull-out and lever-out testing were subsequently simulated using finite element analysis.

Identified and optimised parameters allowed the accurate numerical reproduction of the raw material tests. Correlation between experimental tests and the numerical simulation of primary implant stability depended on the value of interference fit. However, the validated material model provides the opportunity for subsequent parametric numerical studies.     

Standards for molecular dynamics modelling and simulation of relaxation

An attempt is made to formulate a set of requirements for simulation and modelling of relaxation in dense media. Each requirement is illustrated by examples of numerical simulation of particles with different types of interaction given by soft-sphere, Lennard–Jones, embedded atom method or Coulomb potential. The approaches developed are expected to be universal for some classes of relaxation processes in liquids, fluids, crystals and plasmas.     

Soft tissue modelling for applications in virtual surgery and surgical robotics

Soft tissue modelling has gained a great deal of importance, for a large part due to its application in surgical training simulators for minimally invasive surgery (MIS). This article provides a structured overview of different continuum-mechanical models that have been developed over the years. It aims at facilitating model choice for specific soft tissue modelling applications. According to the complexity of the model, different features of soft biological tissue will be incorporated, i.e. nonlinearity, viscoelasticity, anisotropy, heterogeneity and finally, tissue damage during deformation. A brief summary of experimental methods for material characterisation and an introduction to methods for geometric modelling are also provided.

The overview is non-exhaustive, focusing on the most important general models and models with specific biological applications. A trade-off in complexity must be made for enabling real-time simulation, but still maintaining realistic representation of the organ deformation. Depending on the organ and tissue types, different models with emphasis on certain features will prove to be more appropriate, meaning the optimal model choice is organ and tissue-dependent.     

Tribology of polymer brush: microscale modelling and simulation

In order to examine the tribological behaviour of the polymer brush, we constructed a mesoscale particle model of a polymer brush system, based on a Brownian dynamics scheme. The polymer model consists of coarse-grained beads connected with harmonic springs. The Lennard-Jones type interaction is assumed between beads. The flow velocity is affected by the local packing fraction of beads. With this model, we executed a series of molecular dynamics simulation to investigate the mechanical responses of brush to shear flow.     

Non-linear elastic deformable models for real-time surgery simulation

In this article, we describe the latest developments of the minimally invasive hepatic surgery simulator prototype developed at INRIA. A key problem with such a simulator is the physical modelling of soft tissues. We propose a new deformable model based on non-linear elasticity and the finite element method. This model is valid for large displacements, which means in particular that it is invariant with respect to rotations. This property improves the realism of the deformations and solves the problems related to the shortcomings of linear elasticity, which is only valid for small displacements. We also address the problem of volume variations by adding to our model incompressibility constraints. Finally, we demonstrate the relevance of this approach for the real-time simulation of laparoscopic surgical gestures on the liver.     

An Efficient,Box Shape Independent Non-Bonded Force and Virial Algorithm for Molecular Dynamics

Abstract

A notation is introduced and used to transform a conventional specification of the non-bonded force and virial algorithm in the case of periodic boundary conditions into an alternative specification. The implementation of the transformed specification is simpler and typically a factor of 1.5 faster than a conventional implementation. Moreover, it is generic with respect to the shape of the simulated system, i.e. the same routines can be used to handle triclinic boxes, truncated octahedron boxes etc. An implementation of this method is presented, and the speed achieved on various machines is given. Essence of the new method is that the number of calculations of image particle positions is strongly reduced during non-bonded force calculations.     

Patient-specific modelling of pulmonary airflow using GPU cluster for the application in medical practice

In this paper, we propose a novel patient-specific method of modelling pulmonary airflow using graphics processing unit (GPU) computation that can be applied in medical practice. To overcome the barriers imposed by computation speed, installation price and footprint to the application of computational fluid dynamics, we focused on GPU computation and the lattice Boltzmann method (LBM). The GPU computation and LBM are compatible due to the characteristics of the GPU. As the optimisation of data access is essential for the performance of the GPU computation, we developed an adaptive meshing method, in which an airway model is covered by isotropic subdomains consisting of a uniform Cartesian mesh. We found that 4³ size subdomains gave the best performance. The code was also tested on a small GPU cluster to confirm its performance and applicability, as the price and footprint are reasonable for medical applications.     

Computational modelling and optimisation of soft tissue outcome in cranio-maxillofacial surgery planning

Cranio-maxillofacial (CMF) surgery operations are associated with rearrangement of facial hard and soft tissues, leading to dramatic changes in facial geometry. Often, correction of the aesthetical patient's appearance is the primary objective of the surgical intervention. Due to the complexity of the facial anatomy and the biomechanical behaviour of soft tissues, the result of the surgical impact cannot always be predicted on the basis of surgeon's intuition and experience alone. Computational modelling of soft tissue outcome using individual tomographic data and consistent numerical simulation of soft tissue mechanics can provide valuable information for surgeons during the planning stage. In this article, we present a general framework for computer-assisted planning of CMF surgery interventions that is based on the reconstruction of patient's anatomy from 3D computer tomography images and finite element analysis of soft tissue deformations. Examples from our clinical case studies that deal with the solution of direct and inverse surgical problems (i.e. soft tissue prediction, inverse implant shape design) demonstrate that the developed approach provides a useful tool for accurate prediction and optimisation of aesthetic surgery outcome.     

石油钻井噪声与振动对鲤鱼(Cyprinus carpio)能量收支的影响

应用现场模拟实验方法测定了石油钻井噪声与振动对鲤鱼能量收支的影响.结果表明,该污染对鲤鱼摄食、生长、代谢和排泄等各能量收支组分均有显著影响,但其影响程度和趋势却不尽相同.鲤鱼摄食、生长和生长转换效率均随噪声与振动强度增大呈减小趋势;其中摄食较为敏感,但生长的受影响程度却显著大于摄食.与摄食和生长的变化趋势不同,在受到显著影响范围内,排泄能和代谢能随噪声与振动强度增大呈U型变化趋势;其敏感程度与摄食相同,受影响程度也显著低于生长.钻井噪声与振动对鲤鱼影响可分两个阶段;第一阶段为,鲤鱼能量摄入量降低导致各能量支出水平不同程度的下降;第二阶段为,噪声与振动进一步增强,除使摄食量继续减小外,同时造成体内代谢量和排泄量的增大,这两种影响的协同作用使鲤鱼生长大幅度降低.通过所建立的不同钻井噪声与振动强度下鲤鱼的能量收支模型,可以看出,鲤鱼生长大幅度降低总是与代谢量和排泄量增大相对应.     

GINsim: a software suite for the qualitative modelling, simulation and analysis of regulatory networks

This paper presents GINsim, a Java software suite devoted to the qualitative modelling, analysis and simulation of genetic regulatory networks. Formally, our approach leans on discrete mathematical and graph-theoretical concepts. GINsim encompasses an intuitive graph editor, enabling the definition and the parameterisation of a regulatory graph, as well as a simulation engine to compute the corresponding qualitative dynamical behaviour. Our computational approach is illustrated by a preliminary model analysis of the inter-cellular regulatory network activating Notch at the dorsal-ventral boundary in the wing imaginal disc of Drosophila. We focus on the cross-regulations between five genes (within and between two cells), which implements the dorsal-ventral border in the developing imaginal disc. Our simulations qualitatively reproduce the wild-type developmental pathway, as well as the outcome of various types of experimental perturbations, such as loss-of-function mutations or ectopically induced gene expression.     

Identification of Plasmodium falciparum apicoplast-targeted tRNA-guanine transglycosylase and its potential inhibitors using comparative genomics,molecular modelling,docking and simulation studies

tRNA modifications play an important role in the proper folding of tRNA and thereby determine its functionality as an adaptor molecule. Notwithstanding the centrality of this basic process in translation, a major gap in the genomics of Plasmodium falciparum is unambiguous identification of enzymes catalysing the various tRNA modifications. In this study, tRNA-modifying enzymes of P. falciparum were annotated using homology-based approach. Based on the presence of these identified enzymes, the modifications were compared with those of prokaryotic and eukaryotic organisms. Through sequence comparison and phylogenetic analysis, we have identified P. falciparum apicoplast tRNA-guanine 34 transglycosylase (TGT, EC: 2.4.2.29), which shows evidence of its prokaryotic origin. The docking analysis of the modelled TGT structures revealed that binding of quinazolinone derivatives is more favourable with P. falciparum apicoplast TGT as compared to human TGT. Molecular dynamic simulation and molecular mechanics/generalized Born surface area analysis of the complex confirmed the greater binding affinity of the ligand in the binding pocket of P. falciparum TGT protein. Further, evolutionary patterning analysis identified the amino acids of P. falciparum apicoplast TGT that are under purifying selection pressure and hence can be good inhibitor-targeting sites. Based on these computational studies, we suggest that P. falciparum apicoplast tRNA-guanine 34 transglycosylase can be a promising drug target.     

A novel approach for extracting viscoelastic parameters of living cells through combination of inverse finite element simulation and Atomic Force Microscopy

Dynamic mechanical behaviour of living cells has been described by viscoelasticity. However, quantitation of the viscoelastic parameters for living cells is far from sophisticated. In this paper, combining inverse finite element (FE) simulation with Atomic Force Microscope characterization, we attempt to develop a new method to evaluate and acquire trustworthy viscoelastic index of living cells. First, influence of the experiment parameters on stress relaxation process is assessed using FE simulation. As suggested by the simulations, cell height has negligible impact on shape of the force–time curve, i.e. the characteristic relaxation time; and the effect originates from substrate can be totally eliminated when stiff substrate (Young’s modulus larger than 3 GPa) is used. Then, so as to develop an effective optimization strategy for the inverse FE simulation, the parameters sensitivity evaluation is performed for Young’s modulus, Poisson’s ratio, and characteristic relaxation time. With the experiment data obtained through typical stress relaxation measurement, viscoelastic parameters are extracted through the inverse FE simulation by comparing the simulation results and experimental measurements. Finally, reliability of the acquired mechanical parameters is verified with different load experiments performed on the same cell.     

Modelling and simulation techniques for membrane biology

One of the most important aspects of Computational Cell Biology is the understanding of the complicated dynamical processes that take place on plasma membranes. These processes are often so complicated that purely temporal models cannot always adequately capture the dynamics. On the other hand, spatial models can have large computational overheads. In this article, we review some of these issues with respect to chemistry, membrane microdomains and anomalous diffusion and discuss how to select appropriate modelling and simulation paradigms based on some or all the following aspects: discrete, continuous, stochastic, delayed and complex spatial processes.