Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

Protocol for constructing subject-specific biomechanical models of knee joint

A robust protocol for building subject-specific biomechanical models of the human knee joint is proposed which uses magnetic resonance imaging, motion analysis and force platform data in conjunction with detailed 3D finite element models. The proposed protocol can be used for determining stress and strain distributions and contact kinetics in different knee elements at different body postures during various physical activities. Several examples are provided to highlight the capabilities and potential applications of the proposed protocol. This includes preliminary results on the role of body weight on the stresses and strains induced in the knee articular cartilages and meniscus during single-leg stance and calculations of the induced stresses and ligament forces during the gait cycle.     

Comparison between mechanical stress and bone mineral density in the femur after total hip arthroplasty by using subject-specific finite element analyses

The mechanism underling bone mineral density (BMD) loss that occurs in the femur after total hip arthroplasty (THA) remains unknown. We compared the equivalent stress and strain energy density (SED) to BMD in the femur after THA using subject-specific finite element analyses. Twenty-four patients who had undergone primary cementless THA were analysed. BMD was measured using dual-energy X-ray absorptiometry (DEXA) at 1 week and 3, 6 and 12 months after THA. Seven regions of interest (ROIs) were defined in accordance with Gruen's system (ROIs 1–7). Computed tomography images of the femurs were acquired pre- and postoperatively, and the images were converted into three-dimensional finite element (FE) models. Equivalent stress and SED were analysed and compared with DEXA data. BMD was maintained 1 year after THA in ROIs 3, 4, 5 and 6, whereas BMD decreased in ROIs 1, 2 and 7. FE analysis revealed that equivalent stress in ROIs 3, 4, 5 and 6 was much higher than that in ROIs 1, 2 and 7. A significant correlation was observed between the rate of changes in BMD and equivalent stress. Reduction of equivalent stress may contribute to decrease in BMD in the femur after THA.     

Prognosis of implant longevity in terms of annual bone loss: a methodological finite element study

Dental implant failure is mainly the consequence of bone loss at peri-implant area. It usually begins in crestal bone. Due to this gradual loss, implants cannot withstand functional force without bone overload, which promotes complementary loss. As a result, implant lifetime is significantly decreased. To estimate implant success prognosis, taking into account 0.2 mm annual bone loss for successful implantation, ultimate occlusal forces for the range of commercial cylindrical implants were determined and changes of the force value for each implant due to gradual bone loss were studied. For this purpose, finite element method was applied and von Mises stresses in implant–bone interface under 118.2 N functional occlusal load were calculated. Geometrical models of mandible segment, which corresponded to Type II bone (Lekholm & Zarb classification), were generated from computed tomography images. The models were analyzed both for completely and partially osseointegrated implants (bone loss simulation). The ultimate value of occlusal load, which generated 100 MPa von Mises stresses in the critical point of adjacent bone, was calculated for each implant. To estimate longevity of implants, ultimate occlusal loads were correlated with an experimentally measured 275 N occlusal load (Mericske-Stern & Zarb). These findings generally provide prediction of dental implants success.     

Mechanical significance of femoral head trabecular bone structure in Loris and Galago evaluated using micromechanical finite element models

Work on the interspecific and intraspecific variation of trabecular bone in the proximal femur of primates demonstrates important architectural variation between animals with different locomotor behaviors. This variation is thought to be related to the processes of bone adaptation whereby bone structure is optimized to the mechanical environment. Micromechanical finite element models were created for the proximal femur of the leaping Galago senegalensis and the climbing and quadrupedal Loris tardigradus by converting bone voxels from high-resolution X-ray computed tomography scans of the femoral head to eight-noded brick elements. The resulting models had approximately 1.8 million elements each. Loading conditions representing takeoff phase of a leap and more generalized load orientations were applied to the models, and the models were solved using the iterative "row-by-row" matrix-vector multiplication algorithm. The principal strain and Von Mises stress results for the leaping model were similar for both species at each load orientation. Similar hip joint reaction forces in the range of 4.9 x to 12 x body weight were calculated for both species under each loading condition, but the hip reaction values estimated for Loris were higher than predicted based on locomotor behavior. These results suggest that functional adaptation to hip joint loading may not fully explain the differences in femoral head trabecular bone structure in Galago and Loris. The finite element method represents a unique and useful tool for analyzing the functional adaptation of trabecular bone in a diversity of animals and for reconstructing locomotor behavior in extinct taxa.     

Prediction of the structural response of the femoral shaft under dynamic loading using subject-specific finite element models

The goal of this study was to predict the structural response of the femoral shaft under dynamic loading conditions using subject-specific finite element (SS-FE) models and to evaluate the prediction accuracy of the models in relation to the model complexity. In total, SS-FE models of 31 femur specimens were developed. Using those models, dynamic three-point bending and combined loading tests (bending with four different levels of axial compression) of bare femurs were simulated, and the prediction capabilities of five different levels of model complexity were evaluated based on the impact force time histories: baseline, mass-based scaled, structure-based scaled, geometric SS-FE, and heterogenized SS-FE models. Among the five levels of model complexity, the geometric SS-FE and the heterogenized SS-FE models showed statistically significant improvement on response prediction capability compared to the other model formulations whereas the difference between two SS-FE models was negligible. This result indicated the geometric SS-FE models, containing detailed geometric information from CT images with homogeneous linear isotropic elastic material properties, would be an optimal model complexity for prediction of structural response of the femoral shafts under the dynamic loading conditions. The average and the standard deviation of the RMS errors of the geometric SS-FE models for all the 31 cases was 0.46 kN and 0.66 kN, respectively. This study highlights the contribution of geometric variability on the structural response variation of the femoral shafts subjected to dynamic loading condition and the potential of geometric SS-FE models to capture the structural response variation of the femoral shafts.     

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction

In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2–5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R² = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R² = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem.     

Biomechanical behaviour of cancellous bone on patellofemoral arthroplasty with Journey prosthesis: a finite element study

Isolated patellofemoral (PF) arthritis of the knee is a common cause of anterior knee pain and disability. Patellofemoral arthroplasty (PFA) is a bone conserving solution for patients with PF degeneration. Failure mechanisms of PFA include growing tibiofemoral arthritis and loosening of components. The implant loosening can be associated with bone resorption or fatigue-failure of bone by overload. This research work aims at determining the structural effects of the implantation of PF prosthesis Journey PFJ (Smith & Nephew, Inc., Memphis, TN, USA) on femoral cancellous bone. For this purpose, the finite element method is considered to perform computational simulations for different conditions, such as well-fixed and loosening scenarios. From the global results obtained, in the well-fixed scenario, a decrease in strain on cancellous bone was noticed, which can be related to bone resorption. In the loosening scenario, when the cement layer becomes inefficient, a significant increase in cancellous bone strain was observed, which can be associated with bone fatigue-failure.These strain changes suggest a weakness of the femur after PFA.     

Predicting the permeability of trabecular bone by micro-computed tomography and finite element modeling

The intrinsic permeability of bone plays an important role in the transport of nutrients and minerals within the tissue, and affects the mechanical stimuli that are related to the fate of the stem cells. The objective of this study was to establish a method to assess trabecular bone permeability using experimental and finite element (FE) modeling approaches based on micro computed tomography (µCT) images. Human cadaveric tibia cube specimens (N=23) were scanned with µCT. The permeability was measured experimentally using a custom-developed constant-head permeameter, and computationally by a poroelastic formulation to simulate the fluid flow within the discretized bone matrix and pore phase. The average of the experimentally measured permeability was 4.84×10⁻¹⁰ m² with a standard deviation of 3.70×10⁻¹⁰ m². A regression model of the µCT determined that the maximum bone area to total area ratio (maxBA/TA) for all slices that are perpendicular to the direction of fluid flow explained 84% of the variability of the natural logarithm of the experimentally measured permeability. The 2D measure of maxBA/TA performed better than 3D measures in general, although some parameters were reasonably well associated with permeability such as bone volume ratio (BV/TV, r=−0.71), the bone surface/bone volume (BS/BV, r=0.73), and the trabecular thickness (TbTh, r=−0.71). The correlation between the permeability predicted with FE models and experimentally measured permeability was reasonable (r=0.69), but the FE approach did not accurately represent the wide variability of permeability measured experimentally. The results of this study suggest that the changes in the trabecular bone microarchitecture have an exponential relationship with permeability, and the use of µCT-based 2D measurement of maxBA/TA performs well at predicting permeability, thus providing a convenient approach to measure this important aspect affecting biomechanical functions in the tissue.     

Maxillary expansion treatment using bone anchors: development and validation of a 3D finite element model

Objective: Develop a finite element (FE) model of a skull to perform biomechanical studies of maxillary expansion using bone anchors (BA).

Materials and methods: A skull model was developed and assigned material properties based on Hounsfield unit (HU) values of cone-beam computerized tomography (CBCT) images. A 3 mm diameter cylindrical BA was modelled and inserted in the palatal bone. A 4 mm transverse displacement was applied on the anchor. An evaluation on the effect on local stresses of BA implantation inclination angle was performed.

Results: Proper displacement results and strain–stress trends for the expansion process were present. Stress distribution patterns were similar as reported in the literature. No significant difference between BA inclination angles was found.

Conclusion: This work leads to a better understanding and prediction of craniofacial and maxillary bone remodelling during ME with BA treatments and is a first step towards the development of patient specific treatments.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

Hexahedral meshing of subject-specific anatomic structures using mapped building blocks

To extend the use of computational techniques like finite element analysis to clinical settings, it would be beneficial to have the ability to generate a unique model for every subject quickly and efficiently. This work is an extension of two previously developed mapped meshing tools that utilised force and displacement control to map a template mesh to a subject-specific surface. The objective of this study was to map a template block structure, common to multiblock meshing techniques, to a subject-specific surface. The rationale is that the blocks are considerably less refined and may be readily edited after mapping, thereby yielding a mesh of high quality in less time than mapping the mesh itself. In this paper, the versatility and robustness of the method was verified by processing four data-sets. The method was found to be robust enough to cope with the variability of bony surface size, spatial position and geometry, producing building block structures (BBSs) that generated meshes comparable to those produced using BBSs that were created manually.     

An accurate estimation of bone density improves the accuracy of subject-specific finite element models

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.     

Calibration of crushable foam plasticity models for synthetic bone material for use in finite element analysis of acetabular cup deformation and primary stability

Polyurethane (PU) foam is a material often used in biomechanical experiments and demands for the definition of crushable foam plasticity (CFP) in numerical simulations of the primary stability and deformation of implants, to describe the crushing behaviour appropriately. Material data of PU foams with five different densities (10–40 pounds per cubic foot were ascertained experimentally in uniaxial compression test and used to calibrate CFP models for finite element modelling. Additionally, experimental and numerical deformation, push-out and lever-out tests of press-fit acetabular cups were carried out to assess the influence of the chosen material definition (linear elastic and CFP) on the numerical results. Comparison of the experimentally and numerically determined force–displacement curves of the uniaxial compression test showed a mean deviation of less than 3%. In primary stability testing, the deviation between the experimental and numerical results was in a range of 0%–27% for CFP modelling and 64%–341% for the linear elastic model. The material definition selected, highly influenced the numerical results in the current study. The use of a CFP model is recommended for further numerical simulations, when a deformation of the foam beyond the yield strength is likely to occur.     

Achilles tendon stress is more sensitive to subject-specific geometry than subject-specific material properties: A finite element analysis

This study used subject-specific measures of three-dimensional (3D) free Achilles tendon geometry in conjunction with a finite element method to investigate the effect of variation in subject-specific geometry and subject-specific material properties on tendon stress during submaximal isometric loading. Achilles tendons of eight participants (Aged 25–35 years) were scanned with freehand 3D ultrasound at rest and during a 70% maximum voluntary isometric contraction. Ultrasound images were segmented, volume rendered and transformed into subject-specific 3D finite element meshes. The mean (±SD) lengths, volumes and cross-sectional areas of the tendons at rest were 62 ± 13 mm, 3617 ± 984 mm³ and 58 ± 11 mm² respectively. The measured tendon strain at 70% MVIC was 5.9 ± 1.3%. Subject-specific material properties were obtained using an optimisation approach that minimised the difference between measured and modelled longitudinal free tendon strain. Generic geometry was represented by the average mesh and generic material properties were taken from the literature. Local stresses were subsequently computed for combinations of subject-specific and generic geometry and material properties. For a given geometry, changing from generic to subject-specific material properties had little effect on the stress distribution in the tendon. In contrast, changing from generic to subject-specific geometry had a 26-fold greater effect on tendon stress distribution. Overall, these findings indicate that the stress distribution experienced by the living free Achilles tendon of a young and healthy population during voluntary loading are more sensitive to variation in tendon geometry than variation in tendon material properties.     

Influence of anisotropic bone properties on the biomechanical behavior of the acetabular cup implant: a multiscale finite element study

Although the biomechanical behavior of the acetabular cup (AC) implant is determinant for the surgical success, it remains difficult to be assessed due to the multiscale and anisotropic nature of bone tissue. The aim of the present study was to investigate the influence of the anisotropic properties of peri-implant trabecular bone tissue on the biomechanical behavior of the AC implant at the macroscopic scale. Thirteen bovine trabecular bone samples were imaged using micro-computed tomography (μCT) with a resolution of 18 μm. The anisotropic biomechanical properties of each sample were determined at the scale of the centimeter based on a dedicated method using asymptotic homogenization. The material properties obtained with this multiscale approach were used as input data in a 3D finite element model to simulate the macroscopic mechanical behavior of the AC implant under different loading conditions. The largest stress and strain magnitudes were found around the equatorial rim and in the polar area of the AC implant. All macroscopic stiffness quantities were significantly correlated (R² > 0.85, p < 6.5 e-6) with BV/TV (bone volume/total volume). Moreover, the maximum value of the von Mises stress field was significantly correlated with BV/TV (R² > 0.61, p < 1.6 e-3) and was always found at the bone-implant interface. However, the mean value of the microscopic stress (at the scale of the trabeculae) decrease as a function of BV/TV for vertical and torsional loading and do not depend on BV/TV for horizontal loading. These results highlight the importance of the anisotropic properties of bone tissue.