Ning Zhang1, Jeremy F. Magland1, Chamith S. Rajapakse1, Yusuf A. Bhagat1, Felix W. Wehrli1
1Department of Radiology,
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.