Convergence behavior of high resolution finite element models of trabecular bone

摘要

The convergence behavior of finite element models depends on the size of elements used, the element polynomial order, and on the complexity of the applied loads. For high-resolution models of trabecular bone, changes in architecture and densitymay also be important. The goal of this study was to investigate the influence of these factors on the convergence behavior of high-resolution models of trabecular bone. Two human vertebral and two bovine tibial trabecular bone specimens were modeled atfour resolutions ranging from 20 to 80μm and subjected to both compressive and shear loading. Results indicated that convergence behavior depended on both loading mode (axial versus shear) and volume fraction of the specimen. Compared to the 20μmresolution, the differences in apparent Young's modulus at 40μm resolution were less than 5 percent for all specimens, and for apparent shear modulus were less than 7 percent. By contrast, differences at 80μm resolution in apparent modulus were up to 41 percent, depending on the specimen tested and loading mode. Overall, differences in apparent properties were always less than 10 percent when the ratio of mean trabecular thickness to element size was greater than four. Use of higher order elements didnot improve the results. Tissue level parameters such as maximum principal strain did not converge. Tissue level strains converged when considered relative to a threshold value, but only if the strains were evaluated at Gauss points rather than elementcentroids. These findings indicate that good convergence can be obtained with this modeling technique, although element size should be chosen based on factors such as loading mode, mean trabecular thickness, and the particular output parameter of interest.