Influence of Smoothing on Voxel-Based Mesh Accuracy in Micro-Finite Element

The interest in automatic volume meshing for finite element analysis (FEA) has grown more since the appearance of microfocus CT (μCT), due to its high resolution, which allows for the assessment of mechanical behaviour at a high precision. Nevertheless, the basic meshing approach of generating one hexahedron per voxel produces jagged edges. To prevent this effect, smoothing algorithms have been introduced to enhance the topology of the mesh. However, whether smoothing also improves the accuracy of voxel-based meshes in clinical applications is still under question. There is a trade-off between smoothing and quality of elements in the mesh. Distorted elements may be produced by excessive smoothing and reduce accuracy of the mesh. In the present work, influence of smoothing on the accuracy of voxel-based meshes in micro-FE was assessed. An accurate 3D model of a trabecular structure with known apparent mechanical properties was used as a reference model. Virtual CT scans of this reference model (with resolutions of 16, 32 and 64 μm) were then created and used to build voxel-based meshes of the microarchitecture. Effects of smoothing on the apparent mechanical properties of the voxel-based meshes as compared to the reference model were evaluated. Apparent Young’s moduli of the smooth voxel-based mesh were significantly closer to those of the reference model for the 16 and 32 μm resolutions. Improvements were not significant for the 64 μm, due to loss of trabecular connectivity in the model. This study shows that smoothing offers a real benefit to voxel-based meshes used in micro-FE. It might also broaden voxel-based meshing to other biomechanical domains where it was not used previously due to lack of accuracy. As an example, this work will be used in the framework of the European project ContraCancrum, which aims at providing a patient-specific simulation of tumour development in brain and lungs for oncologists. For this type of clinical application, such a fast, automatic, and accurate generation of the mesh is of great benefit.

Hydrogels for nucleus replacement—Facing the biomechanical challenge

Hydrogels are considered promising for disc regeneration strategies. However, it is currently unknown whether the destruction of the natural interface between nucleus and surrounding structures caused by nucleotomy and an inadequate annulus closure diminishes the mechanical competence of the disc. This in vitro study aimed to clarify these mechanisms and to evaluate whether hydrogels are able to restore the biomechanical behaviour of the disc. Nucleus pressure in an ovine intervertebral disc was measured in vivo during day and night and adapted to an in vitro axial compressive diurnal (15min) and night (30min) load. Effects of different defects on disc height and nucleus pressure were subsequently measured in vitro using 30 ovine motion segments. Following cases were considered: intact; annulus incision repaired by suture and glue; annulus incision with removal and re-implantation of nucleus tissue; and two different hydrogels repaired by suture and glue. The intradiscal pressure in vivo was 0.75MPa during day and 0.5MPa during night corresponding to an in vitro axial compressive force of 130 and 58N, respectively. The compression test showed that neither the implantation of hydrogels nor the re-implantation of the natural nucleus, assumed as being the ideal implant, was able to restore the mechanical functionality of an intact disc. Results indicate the importance of the natural anchorage of the nucleus with its surrounding structures and the relevance of an appropriate annulus closure. Therefore, hydrogels that are able to mimic the mechanical behaviour of the native nucleus may fail in restoring the mechanical behaviour of the disc.

Statistical finite element model for bone shape and biomechanical properties

We present a framework for statistical finite element analysis combining shape and material properties, and allowing performing statistical statements of biomechanical performance across a given population. In this paper, we focus on the design of orthopaedic implants that fit a maximum percentage of the target population, both in terms of geometry and biomechanical stability. CT scans of the bone under consideration are registered non-rigidly to obtain correspondences in position and intensity between them. A statistical model of shape and intensity (bone density) is computed by means of principal component analysis. Afterwards, finite element analysis (FEA) is performed to analyse the biomechanical performance of the bones. Realistic forces are applied on the bones and the resulting displacement and bone stress distribution are calculated. The mechanical behaviour of different PCA bone instances is compared.

In situ micropillar compression reveals superior strength and ductility but an absence of damage in lamellar bone

Ageing societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the hierarchical organization of bone. Here, we investigate the mechanical response under monotonic and cyclic compression of both single osteonal lamellae and macroscopic samples containing numerous osteons. Micropillar compression tests in a scanning electron microscope, microindentation and macroscopic compression tests were performed on dry ovine bone to identify the elastic modulus, yield stress, plastic deformation, damage accumulation and failure mechanisms. We found that isolated lamellae exhibit a plastic behaviour, with higher yield stress and ductility but no damage. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behaviour of bone at the microscale to a quasi-brittle response driven by the growth of cracks along interfaces or in the vicinity of pores at the macroscale.

Continuum damage interactions between tension and compression in osteonal bone

Skeletal diseases such as osteoporosis impose a severe socio-economic burden to ageing societies. Decreasing mechanical competence causes a rise in bone fracture incidence and mortality especially after the age of 65 y. The mechanisms of how bone damage is accumulated under different loading modes and its impact on bone strength are unclear. We hypothesise that damage accumulated in one loading mode increases the fracture risk in another. This study aimed at identifying continuum damage interactions between tensile and compressive loading modes. We propose and identify the material constants of a novel piecewise 1D constitutive model capable of describing the mechanical response of bone in combined tensile and compressive loading histories. We performed several sets of loading–reloading experiments to compute stiffness, plastic strains, and stress-strain curves. For tensile overloading, a stiffness reduction (damage) of 60% at 0.65% accumulated plastic strain was detectable as stiffness reduction of 20% under compression. For compressive overloading, 60% damage at 0.75% plastic strain was detectable as a stiffness reduction of 50% in tension. Plastic strain at ultimate stress was the same in tension and compression. Compression showed softening and tension exponential hardening in the post-yield regime. The hardening behaviour in compression is unaffected by a previous overload in tension but the hardening behaviour in tension is affected by a previous overload in compression as tensile reloading strength is significantly reduced. This paper demonstrates how damage accumulated under one loading mode affects the mechanical behaviour in another loading mode. To explain this and to illustrate a possible implementation we proposed a theoretical model. Including such loading mode dependent damage and plasticity behaviour in finite element models will help to improve fracture risk analysis of whole bones and bone implant structures.

The long-term mechanical integrity of non-reinforced PEEK-OPTIMA polymer for demanding spinal applications: experimental and finite-element analysis

Polyetheretherketone (PEEK) is a novel polymer with potential advantages for its use in demanding orthopaedic applications (e.g. intervertebral cages). However, the influence of a physiological environment on the mechanical stability of PEEK has not been reported. Furthermore, the suitability of the polymer for use in highly stressed spinal implants such as intervertebral cages has not been investigated. Therefore, a combined experimental and analytical study was performed to address these open questions. A quasi-static mechanical compression test was performed to compare the initial mechanical properties of PEEK-OPTIMA polymer in a dry, room-temperature and in an aqueous, 37 degrees C environment (n=10 per group). The creep behaviour of cylindrical PEEK polymer specimens (n=6) was measured in a simulated physiological environment at an applied stress level of 10 MPa for a loading duration of 2000 hours (12 weeks). To compare the biomechanical performance of different intervertebral cage types made from PEEK and titanium under complex loading conditions, a three-dimensional finite element model of a functional spinal unit was created. The elastic modulus of PEEK polymer specimens in a physiological environment was 1.8% lower than that of specimens tested at dry, room temperature conditions (P<0.001). The results from the creep test showed an average creep strain of less than 0.1% after 2000 hours of loading. The finite element analysis demonstrated high strain and stress concentrations at the bone/implant interface, emphasizing the importance of cage geometry for load distribution. The stress and strain maxima in the implants were well below the material strength limits of PEEK. In summary, the experimental results verified the mechanical stability of the PEEK-OPTIMA polymer in a simulated physiological environment, and over extended loading periods. Finite element analysis supported the use of PEEK-OPTIMA for load-bearing intervertebral implants.