Finite element analysis modelling of proximal femoral fractures, including post-fixation periprosthetic fractures.

Biomechanical testing has been a cornerstone for the development of surgical implants used in fracture stabilisation. In a multi-disciplinary collaboration complex at the University of Wales, Swansea, novel computerised clinically relevant models were developed using advanced computational engineering. In-house software (developed initially for commercial aerospace engineering), allowed accurate finite element analysis (FEA) models of the whole femur to be created, including the internal architecture of the bone, by means of linear interpolation of greyscale images from multiaxial CT scans. This allowed for modelling the changing trabecular structure and bone mineral density as seen in progressive osteoporosis. Falls from standing were modelled in a variety of directions (with and without muscle action) using analysis programmes which resulted in fractures consistent with those seen in clinical practice. By meshing implants into these models and repeating the mechanism of injury in simulation, periprosthetic fractures were also recreated. Further development with simulated physiological activities (e.g. walking and rising from sitting) along with attrition in the bone (in the boundary zones where stress concentration occurs) will allow further known modes of failure in implants to be reproduced. Robust simulation of macro and micro-scale events will allow the testing of novel new designs in simulations far more complex than conventional biomechanical testing will allow.

Unstructured mesh methods for the solution of the unsteady compressible flow equations with moving boundary components.

A review of a procedure for the simulation of time-dependent, inviscid and turbulent viscous, compressible flows involving geometries that change in time is presented. The adopted discretization technique employs unstructured meshes and both explicit and implicit time-stepping schemes. A dual time-stepping procedure and an ALE formulation enable flows involving moving boundary components to be included. Techniques that have been developed to maintain the validity of the unstructured mesh and to allow for the capture of moving flow features are also reviewed. Using the in-house developed techniques, some examples are included to demonstrate the use of the approach for the simulation of a number of flows of practical industrial interest.

A low-order unstructured-mesh approach for computational electromagnetics in the time domain.

Maxwell's curl equations in the time domain are solved using an explicit linear finite-element approach implemented on unstructured tetrahedral meshes. For the simulation of scattering problems, a perfectly matched layer is added at the artificial far-field boundary, created by the truncation of the physical domain prior to the numerical solution. The complete solution procedure is parallelized. The computational challenges that are encountered when attempting simulations at higher frequencies suggest that the implementation of a hybrid algorithm could have certain advantages. The hybrid approach adopted uses a combination of the finite-element procedure and the well-known low operation count/low storage finite-difference time-domain method. Examples are included to demonstrate the numerical performance of the techniques that are described.