Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7-72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient [Formula: see text]. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.
Mechanical properties of human trabecular bone play an important role in age-related bone fragility and implant stability. Microfinite element (lFE) analysis allows computing the apparent elastic properties of trabecular bone for use in homogenized FE (hFE) analysis,but the results depend unfortunately on the type of applied boundary conditions(BCs). In this study, 167 human femoral trabecular cubic regions with a side length of 5.3mm were extracted from three proximal femora and analyzed using lFE analysis to compare systematically their stiffness with kinematic uniform BCs (KUBCs) and periodicity-compatible mixed uniform BCs (PMUBCs). The obtained elastic constants were then used in the volume fraction and fabric-based orthotropic Zysset–Curnier model to identify their respective model parameters. As expected, PMUBCs lead to more compliant apparent elastic properties than KUBCs, especially in shear. The differences in stiffness decreased with bone volume fraction and mean intercept length (MIL). Unlike KUBCs, PMUBCs were sensitive to heterogeneity of the biopsies. The Zysset–Curnier model fitted the apparent elastic constants successfully in both cases with adjusted coefficients of determination (r2adj) of 0.986 for KUBCs and 0.975 for PMUBCs. The proper use of these BCs for hFE analysis of whole bones will need to be investigated in future work.
The micro-architecture of cancellous bone is considered a major determinant of the fracture risk. Yet, if morphometry tells about alterations of the trabecular network, its elastic behaviour is best described by bone volume fraction (BV/TV) and the fabric tensor, which gives the anisotropy of the trabecular structure. This remains to be proven for yield strength, the onset of bone failure. The microstructure of 126 samples extracted from femoral heads of two female subjects was evaluated on micro-computed tomography scans via 25 structural indices. Parameters such as plate and rod decomposition via ITS and textural analyses by ISV, similar to the trabecular bone score, were also examined. The degree of collinearity between indices was assessed. The indices considered sufficiently independent were included in multi-linear regression models predicting stiffness or yield strength measured via nonlinear micro finite element analyses. The models' accuracy was checked and the contributions of all explanatory variables to the prediction were compared. Our results show that BV/TV alone explained most of the predicted yield strength (76%) and stiffness (89%). BV/TV together with the fabric tensor explained more than 98% of both measures! The fabric tensor also had a larger impact on yield strength (23%) than on the stiffness predictions (9%). On the other hand, the predictive value of the other independent factors (Tb.Th.SD, Tb.Sp.SD, rTb.Th, RR.Junc.D, ISV) was negligible (<1%). In conclusion, just as stiffness, yield strength of femoral trabecular bone is also best explained by BV/TV and trabecular anisotropy, the latter being even more relevant in its post-elastic behaviour.
Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r=0.81, femur-specific fabric: r=0.88, mean fabric: r=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.
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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.