The effects of impingement and dysplasia on stress distributions in the hip joint during sitting and walking: a finite element analysis

Soft tissue damage has been observed in hip joints with pathological geometries. Our primary goal was to study the relationship between morphological variations of the bony components of the hip and resultant stresses within the soft tissues of the joint during routine daily activities. The secondary goal was to find the range of morphological parameters in which stresses are minimized. Computational models of normal and pathological joints were developed based on variations of morphological parameters of the femoral head (Alpha angle) and acetabulum (CE angle). The Alpha angle was varied between 40 degrees (normal joint) and 80 degrees (cam joint). The CE angle was varied between 0 degrees (dysplastic joint) and 40 degrees (pincer joint). Dynamic loads and motions for walking and standing to sitting were applied to all joint configurations. Contact pressures and stresses were calculated and crosscompared to evaluate the influence of morphology. The stresses in the soft tissues depended strongly on the head and acetabular geometry. For the dysplastic joint, walking produced high acetabular rim stresses. Conversely, for impinging joints, standing-to-sitting activities that involved extensive motion were critical, inducing excessive distortion and shearing of the tissue-bone interface. Zones with high von Mises stresses corresponded with clinically observed damage zones in the acetabular cartilage and labrum. Hip joint morphological parameters that minimized were 20 degrees

Comparison of theoretical fixation stability of three devices employed in medial opening wedge high tibial osteotomy: a finite element analysis.

BACKGROUND Medial open wedge high tibial osteotomy is a well-established procedure for the treatment of unicompartmental osteoarthritis and symptomatic varus malalignment. We hypothesized that different fixation devices generate different fixation stability profiles for the various wedge sizes in a finite element (FE) analysis. METHODS Four types of fixation were compared: 1) first and 2) second generation Puddu plates, and 3) TomoFix plate with and 4) without bone graft. Cortical and cancellous bone was modelled and five different opening wedge sizes were studied for each model. Outcome measures included: 1) stresses in bone, 2) relative displacement of the proximal and distal tibial fragments, 3) stresses in the plates, 4) stresses on the upper and lower screw surfaces in the screw channels. RESULTS The highest load for all fixation types occurred in the plate axis. For the vast majority of the wedge sizes and fixation types the shear stress (von Mises stress) was dominating in the bone independent of fixation type. The relative displacements of the tibial fragments were low (in μm range). With an increasing wedge size this displacement tended to increase for both Puddu plates and the TomoFix plate with bone graft. For the TomoFix plate without bone graft a rather opposite trend was observed.For all fixation types the occurring stresses at the screw-bone contact areas pulled at the screws and exceeded the allowable threshold of 1.2 MPa for at least one screw surface. Of the six screw surfaces that were studied, the TomoFix plate with bone graft showed a stress excess of one out of twelve and without bone graft, five out of twelve. With the Puddu plates, an excess stress occurred in the majority of screw surfaces. CONCLUSIONS The different fixation devices generate different fixation stability profiles for different opening wedge sizes. Based on the computational simulations, none of the studied osteosynthesis fixation types warranted an intransigent full weight bearing per se. The highest fixation stability was observed for the TomoFix plates and the lowest for the first generation Puddu plate. These findings were revealed in theoretical models and need to be validated in controlled clinical settings.

Finite element analysis in asymptomatic, symptomatic, and ruptured abdominal aortic aneurysms: in search of new rupture risk predictors.

OBJECTIVES To compare biomechanical rupture risk parameters of asymptomatic, symptomatic and ruptured abdominal aortic aneurysms (AAA) using finite element analysis (FEA). STUDY DESIGN Retrospective biomechanical single center analysis of asymptomatic, symptomatic, and ruptured AAAs. Comparison of biomechanical parameters from FEA. MATERIALS AND METHODS From 2011 to 2013 computed tomography angiography (CTA) data from 30 asymptomatic, 15 symptomatic, and 15 ruptured AAAs were collected consecutively. FEA was performed according to the successive steps of AAA vessel reconstruction, segmentation and finite element computation. Biomechanical parameters Peak Wall Rupture Risk Index (PWRI), Peak Wall Stress (PWS), and Rupture Risk Equivalent Diameter (RRED) were compared among the three subgroups. RESULTS PWRI differentiated between asymptomatic and symptomatic AAAs (p < .0004) better than PWS (p < .1453). PWRI-dependent RRED was higher in the symptomatic subgroup compared with the asymptomatic subgroup (p < .0004). Maximum AAA external diameters were comparable between the two groups (p < .1355). Ruptured AAAs showed the highest values for external diameter, total intraluminal thrombus volume, PWS, RRED, and PWRI compared with asymptomatic and symptomatic AAAs. In contrast with symptomatic and ruptured AAAs, none of the asymptomatic patients had a PWRI value >1.0. This threshold value might identify patients at imminent risk of rupture. CONCLUSIONS From different FEA derived parameters, PWRI distinguishes most precisely between asymptomatic and symptomatic AAAs. If elevated, this value may represent a negative prognostic factor for asymptomatic AAAs.

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.

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.

Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup

Osteoporotic proximal femur fractures are caused by low energy trauma, typically when falling on the hip from standing height. Finite element simulations, widely used to predict the fracture load of femora in fall, usually include neither mass-related inertial effects, nor the viscous part of bone's material behavior. The aim of this study was to elucidate if quasi-static non-linear homogenized finite element analyses can predict in vitro mechanical properties of proximal femora assessed in dynamic drop tower experiments. The case-specific numerical models of thirteen femora predicted the strength (R2=0.84, SEE=540 N, 16.2%), stiffness (R2=0.82, SEE=233 N/mm, 18.0%) and fracture energy (R2=0.72, SEE=3.85 J, 39.6%); and provided fair qualitative matches with the fracture patterns. The influence of material anisotropy was negligible for all predictions. These results suggest that quasi-static homogenized finite element analysis may be used to predict mechanical properties of proximal femora in the dynamic sideways fall situation.

Mesh-based vs. Image-based Statistical Appearance Model of the Human Femur: a Preliminary Comparison Study for the Creation of Finite Element Meshes

Statistical models have been recently introduced in computational orthopaedics to investigate the bone mechanical properties across several populations. A fundamental aspect for the construction of statistical models concerns the establishment of accurate anatomical correspondences among the objects of the training dataset. Various methods have been proposed to solve this problem such as mesh morphing or image registration algorithms. The objective of this study is to compare a mesh-based and an image-based statistical appearance model approaches for the creation of nite element(FE) meshes. A computer tomography (CT) dataset of 157 human left femurs was used for the comparison. For each approach, 30 finite element meshes were generated with the models. The quality of the obtained FE meshes was evaluated in terms of volume, size and shape of the elements. Results showed that the quality of the meshes obtained with the image-based approach was higher than the quality of the mesh-based approach. Future studies are required to evaluate the impact of this finding on the final mechanical simulations.

Differences in endplate deformation of the adjacent and augmented vertebra following cement augmentation

Vertebral cement augmentation can restore the stiffness and strength of a fractured vertebra and relieve chronic pain. Previous finite element analysis, biomechanical tests and clinical studies have indirectly associated new adjacent vertebral fractures following augmentation to altered loading. The aim of this repeated measures in situ biomechanical study was to determine the changes in the adjacent and augmented endplate deformation following cement augmentation of human cadaveric functional spine units (FSU) using micro-computed tomography (micro-CT). The surrounding soft tissue and posterior elements of 22 cadaveric human FSU were removed. FSU were assigned to two groups, control (n = 8) (loaded on day 1 and day 2) and augmented (n = 14) (loaded on day 1, augmented 20% cement fill, and loaded on day 2). The augmented group was further subdivided into a prophylactic augmentation group (n = 9), and vertebrae which spontaneously fractured during loading on day 1 (n = 5). The FSU were axially loaded (200, 1,000, 1,500-2,000 N) within a custom made radiolucent, saline filled loading device. At each loading step, FSUs were scanned using the micro-CT. Endplate heights were determined using custom software. No significant increase in endplate deformation following cement augmentation was noted for the adjacent endplate (P > 0.05). The deformation of the augmented endplate was significantly reduced following cement augmentation for both the prophylactic and fracture group (P < 0.05, P < 0.01, respectively). Endplate deformation of the controls showed no statistically significant differences between loading on day 1 and day 2. A linear relationship was noted between the applied compressive load and endplate deflection (R (2) = 0.58). Evidence of significant endplate deformation differences between unaugmented and augmented FSU, while evident for the augmented endplate, was not present for the adjacent endplate. This non-invasive micro-CT method may also be useful to investigate endplate failure, and parameters that predict vertebral failure.

Finite element based estimation of contact areas and pressures of the human scaphoid in various functional positions of the hand

The scaphoid is the most frequently fractured carpal bone. When investigating fixation stability, which may influence healing, knowledge of forces and moments acting on the scaphoid is essential. The aim of this study was to evaluate cartilage contact forces acting on the intact scaphoid in various functional wrist positions using finite element modeling. A novel methodology was utilized as an attempt to overcome some limitations of earlier studies, namely, relatively coarse imaging resolution to assess geometry, assumption of idealized cartilage thicknesses and neglected cartilage pre-stresses in the unloaded joint. Carpal bone positions and articular cartilage geometry were obtained independently by means of high resolution CT imaging and incorporated into finite element (FE) models of the human wrist in eight functional positions. Displacement driven FE analyses were used to resolve inter-penetration of cartilage layers, and provided contact areas, forces and pressure distribution for the scaphoid bone. The results were in the range reported by previous studies. Novel findings of this study were: (i) cartilage thickness was found to be heterogeneous for each bone and vary considerably between carpal bones; (ii) this heterogeneity largely influenced the FE results and (iii) the forces acting on the scaphoid in the unloaded wrist were found to be significant. As major limitations, accuracy of the method was found to be relatively low, and the results could not be compared to independent experiments. The obtained results will be used in a following study to evaluate existing and recently developed screws used to fix scaphoid fractures.

A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads

Due to the inherent limitations of DXA, assessment of the biomechanical properties of vertebral bodies relies increasingly on CT-based finite element (FE) models, but these often use simplistic material behaviour and/or single loading cases. In this study, we applied a novel constitutive law for bone elasticity, plasticity and damage to FE models created from coarsened pQCT images of human vertebrae, and compared vertebral stiffness, strength and damage accumulation for axial compression, anterior flexion and a combination of these two cases. FE axial stiffness and strength correlated with experiments and were linearly related to flexion properties. In all loading modes, damage localised preferentially in the trabecular compartment. Damage for the combined loading was higher than cumulated damage produced by individual compression and flexion. In conclusion, this FE method predicts stiffness and strength of vertebral bodies from CT images with clinical resolution and provides insight into damage accumulation in various loading modes.

Cement distribution, volume, and compliance in vertebroplasty: some answers from an anatomy-based nonlinear finite element study

STUDY DESIGN: The biomechanics of vertebral bodies augmented with real distributions of cement were investigated using nonlinear finite element (FE) analysis. OBJECTIVES: To compare stiffness, strength, and stress transfer of augmented versus nonaugmented osteoporotic vertebral bodies under compressive loading. Specifically, to examine how cement distribution, volume, and compliance affect these biomechanical variables. SUMMARY OF BACKGROUND DATA: Previous FE studies suggested that vertebroplasty might alter vertebral stress transfer, leading to adjacent vertebral failure. However, no FE study so far accounted for real cement distributions and bone damage accumulation. METHODS: Twelve vertebral bodies scanned with high-resolution pQCT and tested in compression were augmented with various volumes of cements and scanned again. Nonaugmented and augmented pQCT datasets were converted to FE models, with bone properties modeled with an elastic, plastic and damage constitutive law that was previously calibrated for the nonaugmented models. The cement-bone composite was modeled with a rule of mixture. The nonaugmented and augmented FE models were subjected to compression and their stiffness, strength, and stress map calculated for different cement compliances. RESULTS: Cement distribution dominated the stiffening and strengthening effects of augmentation. Models with cement connecting either the superior or inferior endplate (S/I fillings) were only up to 2 times stiffer than the nonaugmented models with minimal strengthening, whereas those with cement connecting both endplates (S + I fillings) were 1 to 8 times stiffer and 1 to 12 times stronger. Stress increases above and below the cement, which was higher for the S + I cases and was significantly reduced by increasing cement compliance. CONCLUSION: The developed FE approach, which accounts for real cement distributions and bone damage accumulation, provides a refined insight into the mechanics of augmented vertebral bodies. In particular, augmentation with compliant cement bridging both endplates would reduce stress transfer while providing sufficient strengthening.