Cortical and trabecular mechanical properties in the femoral neck vary differently with changes in bone mineral density.

Graphical Abstract.

Critical loss of primary implant stability in osteosynthesis locking screws under cyclic overloading.

Primary implant stability, which refers to the stability of the implant during the initial healing period is a crucial factor in determining the long-term success of the implant and lays the foundation for secondary implant stability achieved through osseointegration. Factors affecting primary stability include implant design, surgical technique, and patient-specific factors like bone quality and morphology. In vivo, the cyclic nature of anatomical loading puts osteosynthesis locking screws under dynamic loads, which can lead to the formation of micro cracks and defects that slowly degrade the mechanical connection between the bone and screw, thus compromising the initial stability and secondary stability of the implant. Monotonic quasi-static loading used for testing the holding capacity of implanted screws is not well suited to capture this behavior since it cannot capture the progressive deterioration of peri­implant bone at small displacements. In order to address this issue, this study aims to determine a critical point of loss of primary implant stability in osteosynthesis locking screws under cyclic overloading by investigating the evolution of damage, dissipated energy, and permanent deformation. A custom-made test setup was used to test implanted 2.5 mm locking screws under cyclic overloading test. For each loading cycle, maximum forces and displacement were recorded as well as initial and final cycle displacements and used to calculate damage and energy dissipation evolution. The results of this study demonstrate that for axial, shear, and mixed loading significant damage and energy dissipation can be observed at approximately 20 % of the failure force. Additionally, at this load level, permanent deformations on the screw-bone interface were found to be in the range of 50 to 150 mm which promotes osseointegration and secondary implant stability. This research can assist surgeons in making informed preoperative decisions by providing a better understanding of the critical point of loss of primary implant stability, thus improving the long-term success of the implant and overall patient satisfaction.

Trabecular architecture of the distal femur in extant hominids.

Extant great apes are characterized by a wide range of locomotor, postural and manipulative behaviours that each require the limbs to be used in different ways. In addition to external bone morphology, comparative investigation of trabecular bone, which (re-)models to reflect loads incurred during life, can provide novel insights into bone functional adaptation. Here, we use canonical holistic morphometric analysis (cHMA) to analyse the trabecular morphology in the distal femoral epiphysis of Homo sapiens (n = 26), Gorilla gorilla (n = 14), Pan troglodytes (n = 15) and Pongo sp. (n = 9). We test two predictions: (1) that differing locomotor behaviours will be reflected in differing trabecular architecture of the distal femur across Homo, Pan, Gorilla and Pongo; (2) that trabecular architecture will significantly differ between male and female Gorilla due to their different levels of arboreality but not between male and female Pan or Homo based on previous studies of locomotor behaviours. Results indicate that trabecular architecture differs among extant great apes based on their locomotor repertoires. The relative bone volume and degree of anisotropy patterns found reflect habitual use of extended knee postures during bipedalism in Homo, and habitual use of flexed knee posture during terrestrial and arboreal locomotion in Pan and Gorilla. Trabecular architecture in Pongo is consistent with a highly mobile knee joint that may vary in posture from extension to full flexion. Within Gorilla, trabecular architecture suggests a different loading of knee in extension/flexion between females and males, but no sex differences were found in Pan or Homo, supporting our predictions. Inter- and intra-specific variation in trabecular architecture of distal femur provides a comparative context to interpret knee postures and, in turn, locomotor behaviours in fossil hominins.

High-resolution local trabecular strain within trabecular structure under cyclic loading.

Trabecular bone structure is a complex microstructure consisting of rods and plates, which poses challenges for its mechanical characterization. Digital image correlation (DIC) offers the possibility to characterize the strain response on the surface of trabecular bone. This study employed DIC equipped with a telecentric lens to investigate the strain state of individual trabeculae within their trabecular structure by assessing the longitudinal strain of the trabeculae at both the middle and near the edges of the trabeculae. Due to the high-resolution of the used DIC system, local surface strain of trabeculae was analyzed too. Lastly, the correlation between longitudinal trabecular strain and the orientation and slenderness of the trabeculae was investigated. The results showed that the strain magnification close to the edge of the trabeculae was higher and reached up to 8-folds the strain along the middle of the trabeculae. On the contrary, no strain magnification was found for most of the trabeculae between the longitudinal trabecular strain along the middle of the trabeculae and the globally applied strain. High-resolution full-field strain maps were obtained on the surface of trabeculae showing heterogeneous strain distribution with increasing load. No significant correlation was found between longitudinal trabecular strain and its orientation or slenderness. These findings and the applied methodology can be used to broaden our understanding of the deformation mechanisms of trabeculae within the trabecular network.

The deep trabecular structure of first metacarpals in extant hominids.

OBJECTIVES: Recent studies have associated subarticular trabecular bone distribution in the extant hominid first metacarpal (Mc1) with observed thumb use, to infer fossil hominin thumb use. Here, we analyze the entire Mc1 to test for interspecific differences in: (1) the absolute volume of trabecular volume fraction, (2) the distribution of the deeper trabecular network, and (3) the distribution of trabeculae in the medullary cavity, especially beneath the Mc1 disto-radial flange. MATERIALS AND METHODS: Trabecular bone was imaged using micro-computed tomography in a sample of Homo sapiens (n = 11), Pan paniscus (n = 10), Pan troglodytes (n = 11), Gorilla gorilla (n = 10) and Pongo sp., (n = 7). Using Canonical Holistic Morphometric Analysis (cHMA), we tested for interspecific differences in the trabecular bone volume fraction (BV/TV) and its relative distribution (rBV/TV) throughout the Mc1, including within the head, medullary cavity, and base. RESULTS: P. paniscus had the highest, and H. sapiens the lowest, BV/TV relative to other species. rBV/TV distribution statistically distinguished the radial concentrations and lack of medullary trabecular bone in the H. sapiens Mc1 from all other hominids. H. sapiens and, to a lesser extent, G. gorilla also had a significantly higher trabecular volume beneath the disto-radial flange relative to other hominids. DISCUSSION: These results are consistent with differences in observed thumb use in these species and may also reflect systemic differences in bone volume fraction. The trabecular bone extension into the medullary cavity and concentrations beneath the disto-radial flange may represent crucial biomechanical signals that will aid in the inference of fossil hominin thumb use.

The 2-layer elasto-visco-plastic rheological model for the material parameter identification of bone tissue extended by a damage law.

The response of bone tissue to mechanical load is complex and includes plastic hardening, viscosity and damage. The quantification of these effects plays a mayor role in bone research and in biomechanical clinical trials as to better understand related diseases. In this study, the damage growth in individual wet human trabeculae subjected to cyclic overloading is quantified by inverse rheological modeling. Therefore, an already published rheological material model, that includes linear elasticity, plasticity and viscosity is extended by a damage law. The model is utilized in an optimization process to identify the corresponding material parameters and damage growth in single human trabeculae under tensile load. Results show that the damage model is leading to a better fit of the test data with an average root-mean-square-error (RMSE) of 2.52 MPa compared to the non-damage model with a RMSE of 3.03 MPa. Although this improvement is not significant, the damage model qualitatively better represents the data as it accounts for the visible stiffness reduction along the load history. It returns realistic stiffness values of 11.92 GPa for the instantaneous modulus and 5.73 GPa for the long term modulus of wet trabecular human bone. Further, the growth of damage in the tissue along the load history is substantial, with values above 0.8 close to failure. The relative loss of stiffness per cycle is in good agreement with comparable literature. Inverse rheological modeling proves to be a valuable tool for quantifying complex constitutive behavior from a single mechanical measurement. The evolution of damage in the tissue can be identified continuously over the load history and separated from other effects.

Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy.

BACKGROUND AND OBJECTIVE: Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure. METHODS: hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases. RESULTS: The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role. CONCLUSIONS: hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.

Accuracy of osseointegrated screw-bone construct stiffness and peri-implant loading predicted by homogenized FE models relative to micro-FE models.

Computational predictions of stiffness and peri-implant loading of screw-bone constructs are highly relevant to investigate and improve bone fracture fixations. Homogenized finite element (hFE) models have been used for this purpose in the past, but their accuracy has been questioned given the numerous simplifications, such as neglecting screw threads and modelling the trabecular bone structure as a continuum. This study aimed to investigate the accuracy of hFE models of an osseointegrated screw-bone construct when compared to micro-FE models considering the simplified screw geometry and different trabecular bone material models. Micro-FE and hFE models were created from 15 cylindrical bone samples with a virtually inserted, osseointegrated screw (fully bonded interface). Micro-FE models were created including the screw with threads (=reference models) and without threads to quantify the error due to screw geometry simplification. In the hFE models, the screws were modelled without threads and four different trabecular bone material models were used, including orthotropic and isotropic material derived from homogenization with kinematic uniform boundary conditions (KUBC), as well as from periodicity-compatible mixed uniform boundary conditions (PMUBC). Three load cases were simulated (pullout, shear in two directions) and errors in the construct stiffness and the volume average strain energy density (SED) in the peri-implant region were evaluated relative to the micro-FE model with a threaded screw. The pooled error caused by only omitting screw threads was low (max: 8.0%) compared to the pooled error additionally including homogenized trabecular bone material (max: 92.2%). Stiffness was predicted most accurately using PMUBC-derived orthotropic material (error: -0.7 ± 8.0%) and least accurately using KUBC-derived isotropic material (error: +23.1 ± 24.4%). Peri-implant SED averages were generally well correlated (R2 ≥ 0.76), but slightly over- or underestimated by the hFE models and SED distributions were qualitatively different between hFE and micro-FE models. This study suggests that osseointegrated screw-bone construct stiffness can be predicted accurately using hFE models when compared to micro-FE models and that volume average peri-implant SEDs are well correlated. However, the hFE models are highly sensitive to the choice of trabecular bone material properties. PMUBC-derived isotropic material properties represented the best trade-off between model accuracy and complexity in this study.

Comparison of linear and nonlinear stepwise µFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies.

Digital volume correlation (DVC) enables to evaluate the ability of µFE models in predicting experimental results on the mesoscale. In this study predicted displacement fields of three different linear and materially nonlinear µFE simulation methods were compared to DVC measured displacement fields at specific load steps in the elastic regime (StepEl) and after yield (StepUlt). Five human trabecular bone biopsies from a previous study were compressed in several displacement steps until failure. At every compression step, µCT images (resolution: 36 µm) were recorded. A global DVC algorithm was applied to compute the displacement fields at all loading steps. The unloaded 3D images were then used to generate homogeneous, isotropic, linear and materially nonlinear µFE models. Three different µFE simulation methods were used: linear (L), nonlinear (NL), and nonlinear stepwise (NLS). Regarding L and NL, the boundary conditions were derived from the interpolated displacement fields at StepEl and StepUlt, while for the NLS method nonlinear changes of the boundary conditions of the experiments were captured using the DVC displacement field of every available load step until StepEl and StepUlt. The predicted displacement fields of all µFE simulation methods were in good agreement with the DVC measured displacement fields (individual specimens: R2>0.83 at StepEl and R2>0.59 at StepUlt; pooled data: R2>0.97 at StepEl and R2>0.92 at StepUlt). At StepEl, all three simulation methods showed similar intercepts, slopes, and coefficients of determination while the nonlinear µFE models improved the prediction of the displacement fields slightly in all Cartesian directions at StepUlt (individual specimens: L: R2>0.59 and NL, NLS: R2>0.68; pooled data: L: R2>0.92 and NL, NLS: R2>0.94). Damaged/overstrained elements in L, NL, and NLS occurred at similar locations but the number of overstrained elements was overestimated when using the L simulation method. Considering the increased solving time of the nonlinear µFE models as well as the acceptable performance in displacement prediction of the linear µFE models, one can conclude that for similar use cases linear µFE models represent the best compromise between computational effort and accuracy of the displacement field predictions.

Feasibility of aluminum phantom radiography for osteoporosis detection in postmenopausal women with a fragility fracture of the distal radius compared to DXA and HR-pQCT.

Recently, promising results have been reported for detection of osteoporosis with use of an aluminum phantom. Therefore, the aim of this study was to evaluate the feasibility of radiography-based bone mineral density (BMD) measurement using a graded aluminum phantom. This study included 27 postmenopausal women with a distal radius fracture. Aluminum phantom radiography of the healthy radius was conducted as well as high-resolution peripheral quantitative computed tomography (HR-pQCT) measurement of the ultradistal radius and dual energy X-ray absorptiometry (DXA) of the radius, spine, and hip. A strong correlation was observed between aluminum phantom radiography-based mean gray value (mGV) and DXA-derived BMD, especially for the ultradistal radius (ρ = 0.75; p < 0.001). A moderate correlation for the femoral neck (ρ = 0.61 and p < 0.001) between modalities was found. Radius mGV and HR-pQCT-derived BMD only showed a moderate correlation (ρ = 0.48; p < 0.09). Aluminum phantom radiography might serve as a cost efficient, highly available, low-radiation dose screening, and diagnostic method for osteoporosis additively to DXA measurements. Especially, an application in areas with constrained DXA availability and such as preoperative trauma settings would be beneficial. However, further investigation and assessment of specificity and sensitivity is needed.

A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models.

Inverse bone (re)modeling (IBR) can infer physiological loading conditions from the bone microstructure. IBR scales unit loads, imposed on finite element (FE) models of a bone, such that the trabecular microstructure is homogeneously loaded and the difference to a target stimulus is minimized. Micro-FE (µFE) analyses are typically used to model the microstructure, but computationally more efficient, homogenized FE (hFE) models, where the microstructure is replaced by an equivalent continuum, could be used instead. However, also the target stimulus has to be translated from the tissue to the continuum level. In this study, a new continuum-level target stimulus relating relative bone density and strain energy density is proposed. It was applied using different types of hFE models to predict the physiological loading of 21 distal radii sections, which was subsequently compared to µFE-based IBR. The hFE models were able to correctly identify the dominant load direction and showed a high correlation of the predicted forces, but mean magnitude errors ranged from - 14.7 to 26.6% even for the best models. While µFE-based IBR can still be regarded as a gold standard, hFE-based IBR enables faster predictions, the usage of more sophisticated boundary conditions, and the usage of clinical images.

The biomechanical behavior of 3D printed human femoral bones based on generic and patient-specific geometries.

BACKGROUND: Bone is a highly complex composite material which makes it hard to find appropriate artificial surrogates for patient-specific biomechanical testing. Despite various options of commercially available bones with generic geometries, these are either biomechanically not very realistic or rather expensive. METHODS: In this work, additive manufacturing was used for the fabrication of artificial femoral bones. These were based on CT images of four different commercially available femoral bone surrogates and three human bones with varying bone density. The models were 3D printed using a low-budget fused deposition modeling (FDM) 3D printer and PLA filament. The infill density was mechanically calibrated and varying cortical thickness was used. Compression tests of proximal femora simulating stance were performed and the biomechanical behavior concerning ultimate force, spring stiffness, and fracture pattern were evaluated as well as compared to the results of commercial and cadaveric bones. RESULTS: Regarding the ultimate forces and spring stiffness, the 3D printed analogs showed mechanical behavior closer to their real counterparts than the commercially available polyurethan-based surrogates. Furthermore, the increase in ultimate force with increasing bone density observed in human femoral bones could be reproduced well. Also, the fracture patterns observed match well with fracture patterns observed in human hip injuries. CONCLUSION: Consequently, the methods presented here show to be a promising alternative for artificial generic surrogates concerning femoral strength testing. The manufacturing is straightforward, cheap, and patient-specific geometries are possible.

Fracture toughness determination of porcine muscle tissue based on AQLV model derived viscous dissipated energy.

The ability of soft collagenous tissue (SCT) to withstand propagation of a defect in the presence of a macroscopic crack is termed the 'fracture toughness parameter'. In soft tissues not undergoing significant plastic deformation, it is purported that a considerable amount of additional energy is dissipated during failure processes, due to viscoelasticity. Hence the total work, measured experimentally during failure, is the sum of fracture and viscoelastic energies. Previous authors have aimed to apply constitutive modeling to describe viscoelastic hysteresis for fracture toughness determination with a tendency of models to either over or underestimate the viscous energy. In this study, the fracture toughness of porcine muscle tissue is determined using two strategies. Firstly, it was determined experimentally by calculation of the difference in dissipated energy of notched and unnotched tissue specimens undergoing cyclic 'triangular wave' excitation with increasing strain levels in uniaxial tension. The second strategy involved the extension and use of the adaptive quasi-linear viscoelastic model (AQLV) to model cyclic loading (model parameters were obtained from a previous study) and sequentially the dissipated energy was calculated. The mean value of the dissipated energy based on the AQLV approach was then subtracted from the total dissipated energy of notched porcine muscle tissue samples to determine the fracture toughness. The mean experimental viscous dissipated energy ratio was 0.24 ± 0.04 in the experimental approach, compared to 0.28 ± 0.03 for the AQLV model. Fracture toughness determined experimentally yielded 0.84 ± 0.80 kJ/m2, and 0.71 ± 0.76 kJ/m2 for the AQLV model, without a significant difference (p = 0.87). Hence, the AQLV model enables a reasonable estimation of viscous dissipated energy in porcine muscle tissue with the advantage to perform tests only on notched specimens, instead of testing additional unnotched samples. Moreover, the AQLV model will help to better understand the constitutive viscoelastic behaviour of SCTs and might also serve as a basis for future fracture toughness determination with constitutive model simulations.

Predicting the trabecular bone apparent stiffness tensor with spherical convolutional neural networks.

The apparent stiffness tensor is relevant for characterizing trabecular bone quality. Previous studies have used morphology-stiffness relationships for estimating the apparent stiffness tensor. In this paper, we propose to train spherical convolutional neural networks (SphCNNs) to estimate this tensor. Information of the edges, trabecular thickness, and spacing are summarized in functions on the unitary sphere used as inputs for the SphCNNs. The concomitant dimensionality reduction makes it possible to train neural networks on relatively small datasets. The predicted tensors were compared to the stiffness tensors computed by using the micro-finite element method (µFE), which was considered as the gold standard, and models based on fourth-order fabric tensors. Combining edges and trabecular thickness yields significant improvements in the accuracy compared to the methods based on fourth-order fabric tensors. From the results, SphCNNs are promising for replacing the more expensive µFE stiffness estimations.

A computational framework for canonical holistic morphometric analysis of trabecular bone.

Bone is a remarkable, living tissue that functionally adapts to external loading. Therefore, bone shape and internal structure carry information relevant to many disciplines, including medicine, forensic science, and anthropology. However, morphometric comparisons of homologous regions across different individuals or groups are still challenging. In this study, two methods were combined to quantify such differences: (1) Holistic morphometric analysis (HMA) was used to quantify morphometric values in each bone, (2) which could then be mapped to a volumetric mesh of a canonical bone created by a statistical free-form deformation model (SDM). Required parameters for this canonical holistic morphometric analysis (cHMA) method were identified and the robustness of the method was evaluated. The robustness studies showed that the SDM converged after one to two iterations, had only a marginal bias towards the chosen starting image, and could handle large shape differences seen in bones of different species. Case studies were performed on metacarpal bones and proximal femora of different primate species to confirm prior study results. The differences between species could be visualised and statistically analysed in both case studies. cHMA provides a framework for performing quantitative comparisons of different morphometric quantities across individuals or groups. These comparisons facilitate investigation of the relationship between spatial morphometric variations and function or pathology, or both.

A parameter reduced adaptive quasi-linear viscoelastic model for soft biological tissue in uniaxial tension.

Mechanical characterisation of soft viscous materials is essential for many applications including aerospace industries, material models for surgical simulation, and tissue mimicking materials for anatomical models. Constitutive material models are, therefore, necessary to describe soft biological tissues in physiologically relevant strain ranges. Hereby, the adaptive quasi-linear viscoelastic (AQLV) model enables accurate modelling of the strain-dependent non-linear viscoelastic behaviour of soft tissues with a high flexibility. However, the higher flexibility produces a large number of model parameters. In this study, porcine muscle and liver tissue samples were modelled in the framework of the originally published AQLV (3-layers of Maxwell elements) model using four incremental ramp-hold experiments in uniaxial tension. AQLV model parameters were reduced by decreasing model layers (M) as well as the number of experimental ramp-hold steps (N). Leave One out cross validation tests show that the original AQLV model (3M4N) with 19 parameters, accurately describes porcine muscle tissue with an average R2 of 0.90 and porcine liver tissue, R2 of 0.86. Reducing the number of layers (N) in the model produced acceptable model fits for 1-layer (R2 of 0.83) and 2-layer models (R2 of 0.89) for porcine muscle tissue and 1-layer (R2 of 0.84) and 2-layer model (R2 of 0.85) for porcine liver tissue. Additionally, a 2 step (2N) ramp-hold experiment was performed on additional samples of porcine muscle tissue only to further reduce model parameters. Calibrated spring constant values for 2N ramp-hold tests parameters k1 and k2 had a 16.8% and 38.0% deviation from those calibrated for a 4 step (4N) ramp hold experiment. This enables further reduction of material parameters by means of step reduction, effectively reducing the number of parameters required to calibrate the AQLV model from 19 for a 3M4N model to 8 for a 2M2N model, with the added advantage of reducing the time per experiment by 50%. This study proposes a 'reduced-parameter' AQLV model (2M2N) for the modelling of soft biological tissues at finite strain ranges. Sequentially, the comparison of model parameters of soft tissues is easier and the experimental burden is reduced.

Material design of soft biological tissue replicas using viscoelastic micromechanical modelling.

Anatomical models for research and education are often made of artificial materials that attempt to mimic biological tissues in terms of their mechanical properties. Recent developments in additive manufacturing allow tuning mechanical properties with microstructural designs. We propose a strategy for designing material microstructures to mimic soft tissue viscoelastic behaviour, based on a micromechanical Mori-Tanaka model. The model was applied to predict homogenised viscoelastic properties of materials, exhibiting a matrix-inclusion microstructure with varying inclusion volume fractions. The input properties were thereby obtained from compression relaxation tests on silicone elastomers. Validation of the model was done with experimental results for composite samples. Finally, different combinations of silicones were compared to mechanical properties of soft tissues (hepatic, myocardial, adipose, cervical, and prostate tissue), found in literature, in order to design microstructures for replicating these tissues in terms of viscoelasticity. The viscoelastic Mori-Tanaka model showed good agreement with the corresponding experimental results for low inclusion volume fractions, while high fractions lead to underestimation of the complex modulus by the model. Predictions for the loss tangent were reasonably accurate, even for higher inclusion volume fractions. Based on the model, designs for 3D printed microstructures can be extracted in order to replicate the viscoelastic properties of soft tissues.

Experimental validation of a micro-CT finite element model of a human cadaveric mandible rehabilitated with short-implant-supported partial dentures.

PURPOSE: This study aimed to address the predictive value of a micro-computed tomography (µCT)-based finite element (µFE) model of a human cadaveric edentulous posterior mandible, rehabilitated by short dental implants. Hereby, three different prosthetic/implant configurations of fixed partial dentures ("Sp"-3 splinted crowns on 3 implants, "Br" - Bridge: 3 splinted crowns on 2 implants, and "Si"- 3 single crowns) were analysed by comparing the computational predictions of the global stiffness with experimental data. METHODS: Experimental displacement of the bone/implant/prosthesis system was measured under axial and oblique loads of 100 N using an optical deformation system (GOM Aramis) and the overall movement of the testing machine (Zwick Z030). Together with the measured machine force, an "Aramis" (optical markers) and "Zwick" (test machine) stiffness were calculated. FE models were created based on µCT-scans of the cadaveric mandible sample (n = 1) before and after implantation and using stl-files of the crowns. The same load tests and boundary conditions were simulated on the models and the µFE-results were compared to experimental data using linear regression analysis. RESULTS: The regression line through a plot of pooled stiffness values (N/mm) for the optical displacement recording (true local displacement) and the test machine (machine compliance included) had a slope of 0.57 and a correlation coefficient R2 of 0.82. The average pooled correlation of global stiffness between the experiment and FE-analysis (FEA) showed a R2 of 0.80, but the FEA-stiffness was 7.2 times higher. The factor was highly dependent on the test configuration. Sp-configuration showed the largest stiffness followed by Br-configuration (17% difference in experiment and 21% in FEA). CONCLUSIONS: The current study showed good qualitative agreement between the experimental and predicted global stiffness of different short implant configurations. It could be deduced that 1:1 splinting of the short implants by the crowns is most favorable for the stiffness of the implant/prosthesis system. However, in the clinical context, the absolute in silico readings must be interpreted cautiously, as the FEA showed a considerable overestimation of the values.

Effects of Osteoporosis on Bone Morphometry and Material Properties of Individual Human Trabeculae in the Femoral Head.

Osteoporosis is the most common bone disease and is conventionally classified as a decrease of total bone mass. Current diagnosis of osteoporosis is based on clinical risk factors and dual energy X-ray absorptiometry (DEXA) scans, but changes in bone quantity (bone mass) and quality (trabecular structure, material properties, and tissue composition) are not distinguished. Yet, osteoporosis is known to cause a deterioration of the trabecular network, which might be related to changes at the tissue scale-the material properties. The goal of the current study was to use a previously established test method to perform a thorough characterization of the material properties of individual human trabeculae from femoral heads in cyclic tensile tests in a close to physiologic, wet environment. A previously developed rheological model was used to extract elastic, viscous, and plastic aspects of material behavior. Bone morphometry and tissue mineralization were determined with a density calibrated micro-computed tomography (µCT) set-up. Osteoporotic trabeculae neither showed a significantly changed material or mechanical behavior nor changes in tissue mineralization, compared with age-matched healthy controls. However, donors with osteopenia indicated significantly reduced apparent yield strain and elastic work with respect to osteoporosis, suggesting possible initial differences at disease onset. Bone morphometry indicated a lower bone volume to total volume for osteoporotic donors, caused by a smaller trabecular number and a larger trabecular separation. A correlation of age with tissue properties and bone morphometry revealed a similar behavior as in osteoporotic bone. In the range studied, age does affect morphometry but not material properties, except for moderately increased tissue strength in healthy donors and moderately increased hardening exponent in osteoporotic donors. Taken together, the distinct changes of trabecular bone quality in the femoral head caused by osteoporosis and aging could not be linked to suspected relevant changes in material properties or tissue mineralization. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Femoral strength can be predicted from 2D projections using a 3D statistical deformation and texture model with finite element analysis.

Ultimate force of the proximal human femur can be predicted using Finite Element Analysis (FEA), but the models rely on 3D computed tomography images. Landmark-based statistical appearance models (SAM) and B-Spline transformation-based statistical deformation models (SDM) have been used to estimate 3D images from 2D projections, which facilitates model generation and reduces the radiation dose. However, there is no literature on the accuracy of SDM-based FEA models of bones with respect to experimental results. In this study, a methodology for an enhanced SDM with textural information is presented. The statistical deformation and texture models (SDTMs) are based on a set of 37 quantitative CT (QCT) images. They were used to estimate 3D images from two or one projections of the set in a leave-one-out setup. These estimations where then used to create FEA models. The ultimate force predicted by FEA models estimated from two or one projection using the SDTMs were compared to the experimental ultimate force from a previous study on the same femora and to the results of standard QCT-based FEA models. High correlations between predictions and experimental measurements were found for FEA models reconstructed from 2D projections with R2=0.835 when based on two projections and R2=0.724 when using one projection. The correlations were comparable to those reached with standard QCT-based FE-models with the experimental results (R2=0.795). This study shows the high potential of SDTM-based 3D image reconstruction and FEA modelling from 2D projections to predict femoral ultimate force.