Homogenized finite element simulations can predict the primary stability of dental implants in human jawbone.

Adequate primary stability is a pre-requisite for the osseointegration and long-term success of dental implants. Primary stability depends essentially on the bone mechanical integrity at the implantation site. Clinically, a qualitative evaluation can be made on medical images, but finite element (FE) simulations can assess the primary stability of a bone-implant construct quantitatively based on high-resolution CT images. However, FE models lack experimental validation on clinically relevant bone anatomy. The aim of this study is to validate such an FE model on human jawbones. Forty-seven bone biopsies were extracted from human cadaveric jawbones. Dental implants of two sizes (Ø3.5 mm and Ø4.0 mm) were inserted and the constructs were subjected to a quasi-static bending-compression loading protocol. Those mechanical tests were replicated with sample-specific non-linear homogenized FE models. Bone was modeled with an elastoplastic constitutive law that included damage. Density-based material properties were mapped based on µCT images of the bone samples. The experimental ultimate load was better predicted by FE (R2 = 0.83) than by peri-implant bone density (R2 = 0.54). Unlike bone density, the simulations were also able to capture the effect of implant diameter. The primary stability of a dental implant in human jawbones can be predicted quantitatively with FE simulations. This method may be used for improving the design and insertion protocols of dental implants.

Estimation of the effective yield properties of human trabecular bone using nonlinear micro-finite element analyses.

Micro-finite element ([Formula: see text]FE) analyses are often used to determine the apparent mechanical properties of trabecular bone volumes. Yet, these apparent properties depend strongly on the applied boundary conditions (BCs) for the limited size of volumes that can be obtained from human bones. To attenuate the influence of the BCs, we computed the yield properties of samples loaded via a surrounding layer of trabecular bone ("embedded configuration"). Thirteen cubic volumes (10.6 mm side length) were collected from [Formula: see text]CT reconstructions of human vertebrae and femora and converted into [Formula: see text]FE models. An isotropic elasto-plastic material model was chosen for bone tissue, and nonlinear [Formula: see text]FE analyses of six uniaxial, shear, and multi-axial load cases were simulated to determine the yield properties of a subregion (5.3 mm side length) of each volume. Three BCs were tested. Kinematic uniform BCs (KUBCs: each boundary node is constrained with uniform displacements) and periodicity-compatible mixed uniform BCs (PMUBCs: each boundary node is constrained with a uniform combination of displacements and tractions mimicking the periodic BCs for an orthotropic material) were directly applied to the subregions, while the embedded configuration was achieved by applying PMUBCs on the larger volumes instead. Yield stresses and strains, and element damage at yield were finally compared across BCs. Our findings indicate that yield strains do not depend on the BCs. However, KUBCs significantly overestimate yield stresses obtained in the embedded configuration (+43.1 ± 27.9%). PMUBCs underestimate (-10.0 ± 11.2%), but not significantly, yield stresses in the embedded situation. Similarly, KUBCs lead to higher damage levels than PMUBCs (+51.0 ± 16.9%) and embedded configurations (+48.4 ± 15.0%). PMUBCs are better suited for reproducing the loading conditions in subregions of the trabecular bone and deliver a fair estimation of their effective (asymptotic) yield properties.