Enhanced Algorithm for Desktop PC-Based Micro-Finite Element Modeling of Whole-Section Stiffness from In Vivo MR Images

Ning Zhang¹, Jeremy F. Magland¹, Chamith S. Rajapakse¹, Yusuf A. Bhagat¹, Felix W. Wehrli¹

High resolution MR image based FE simulations have been shown to provide valuable information on changes in bone mechanical properties during aging and in response to intervention. However, the scales of the FE models are currently restricted to computing resources. In this work we describe an optimized algorithm for MR image-based finite-element simulation of trabecular bone, with the objective of running these simulations within the limitations of standard workstations. Simulations based on in-vivo MR images and ex-vivo micro-CT images were performed. Memory usage and CPU times were compared. The optimized algorithm enables simulations of systems with as many as 30 million elements on desktop computers.