Computer simulating a trial of a load-bearing implant: example of an intramedullary prosthesis

Computational modelling is becoming ever more important for obtaining regulatory approval for new medical devices. An accepted approach is to infer performance in a population from an analysis conducted for an idealised or ‘average’ patient; we present here a method for predicting the performance of an orthopaedic implant when released into a population—effectively simulating a clinical trial. Specifically we hypothesise that an analysis based on a method for predicting the performance in a population will lead to different conclusions than an analysis based on an idealised or ‘average’ patient. To test this hypothesis we use a finite element model of an intramedullary implant in a bone whose size and remodelling activity is different for each individual in the population. We compare the performance of a low Young’s modulus implant (View the MathML source) to one with a higher Young’s modulus (200 GPa). Cyclic loading is applied and failure is assumed when the migration of the implant relative to the bone exceeds a threshold magnitude. The analysis for an idealised of ‘average’ patient predicts that the lower modulus device survives longer whereas the analysis simulating a clinical trial predicts no statistically-significant tendency (p=0.77) for the low modulus device to perform better. It is concluded that population-based simulations of implant performance–simulating a clinical trial–present a very valuable opportunity for more realistic computational pre-clinical testing of medical devices.

Numerical simulation of a patent technology for sealing of deep-sea oil wells using nonlinear finite element method

For over 50 years bridge plugs and cement have been used for well abandonment and work over and are still the material of choice. However the failures of cement abandonments using bridge plugs has been reported on many occasions, some of which have resulted in fatal consequences. A new patented product is designed to address the shortcomings associated with using bridge plugs and cement. The new developed tools use an alloy based on bismuth that is melted in situ using Thermite reaction. The tool uses the expansion properties of bismuth to seal the well. Testing the new technology in real field under more than 2 km deep sea water can be expensive. Virtual simulation of the new device under simulated thermal and mechanical environment can be achieved using nonlinear finite element method to validate the product and reduce cost. Experimental testing in the lab is performed to measure heat generated due to thermite reaction. Then, a sequential thermal mechanical explicit/implicit finite element solver is used to simulate the device under both testing lab and deep water conditions.

A comparison of stress intensity factors through crack closure integral and other approaches using eXtended element-free Galerkin method

In this paper, a new approach for extracting stress intensity factors (SIFs) by the extended element-free Galerkin method, through a crack closure integral (CCI) scheme, is proposed. The CCI calculation is used in conjunction with a local smoothing technique to improve the accuracy of the computed SIFs in a number of case studies of linear elastic fracture mechanics. The cases involve problems of mixed-mode, curved crack and thermo-mechanical loading. The SIFs by CCI, displacement and stress methods are compared with those based on the M-integral technique reported in the literature. The proposed CCI method involves very simple relations, and still gives good accuracy. The convergence of the results is also examined.

Predicting the ultimate load of a CFRP wingbox

The finite element method in conjunction with the Soutis-Fleck model is used to predict the residual strength after impact of a carbon-fibre reinforced plastic wingbox subjected to a cantilever type loading. The maximum stress failure criterion further validates the Soutis-Fleck model predictions. The Soutis-Fleck model predicts that the wingbox fails at a tip load of 99.2 kN, approximately 5.5% less than the experimental observation

Algorithmic method for scraper load-time optimization

Scrapers have established an important position in the earthmoving field as they are independently capable of accomplishing an earthmoving operation. Given that loading a scraper to its capacity does not entail its maximum production, optimizing the scraper’s loading time is an essential prerequisite for successful operations management. The relevant literature addresses the loading time optimization through a graphical method that is founded on the invalid assumption that the hauling time is independent of the load time. To correct this, a new algorithmic optimization method that incorporates the golden section search and the bisection algorithm is proposed. Comparison of the results derived from the proposed and the existing method demonstrates that the latter entails the systematic needless prolongation of the loading stage thus resulting in reduced hourly production and increased cost. Therefore, the proposed method achieves an improved modeling of scraper earthmoving operations and contributes toward a more efficient cost management.

Finite element modelling of the formation of pre-load damage in cement mantles of orthopaedic joint replacements

Finite element modeling of the formation of pre-loaded damage in cement mantles of orthopaedic joint replacements was presented. The existence of cracking suggested a high level of residual stress. The direction of maximum principal stress vectors corresponded well with the observed crack orientation. Results suggested that cracking depends upon a combination of residual stress, porosity and temperature rise during polymerization.

Simulation of three-dimension cutting force and tool deflection in the end milling operation based on finite element method

This paper presents a 3D simulation system which is employed in order to predict cutting forces and tool deflection during end-milling operation. In order to verify the accuracy of 3D simulation, results (cutting forces and tool deflection) were compared with those based on the theoretical relationships, in terms of agreement with experiments. The results obtained indicate that the simulation is capable of predicting the cutting forces and tool deflection.

Stiffness analysis of polymeric composites using the finite element method

The ability to predict the mechanical behavior of polymer composites is crucial for their design and manufacture. Extensive studies based on both macro- and micromechanical analyses are used to develop new insights into the behavior of composites. In this respect, finite element modeling has proved to be a particularly powerful tool. In this article, we present a Galerkin scheme in conjunction with the penalty method for elasticity analyses of different types of polymer composites. In this scheme, the application of Green's theorem to the model equation results in the appearance of interfacial flux terms along the boundary between the filler and polymer matrix. It is shown that for some types of composites these terms significantly affect the stress transfer between polymer and fillers. Thus, inclusion of these terms in the working equations of the scheme preserves the accuracy of the model predictions. The model is used to predict the most important bulk property of different types of composites. Composites filled with rigid or soft particles, and composites reinforced with short or continuous fibers are investigated. For each case, the results are compared with the available experimental results and data obtained from other models reported in the literature. Effects of assumptions made in the development of the model and the selection of the prescribed boundary conditions are discussed.