Computationally-Optimized Bone Mechanical Modeling from High-Resolution Structural Images

Image-based mechanical modeling of the complex micro-structure of human bone has shown promise as a non-invasive method for characterizing bone strength and fracture risk in vivo. In particular, elastic moduli obtained from image-derived micro-finite element (μFE) simulations have been shown to correlate well with results obtained by mechanical testing of cadaveric bone. However, most existing large-scale finite-element simulation programs require significant computing resources, which hamper their use in common laboratory and clinical environments. In this work, we theoretically derive and computationally evaluate the resources needed to perform such simulations (in terms of computer memory and computation time), which are dependent on the number of finite elements in the image-derived bone model. A detailed description of our approach is provided, which is specifically optimized for μFE modeling of the complex three-dimensional architecture of trabecular bone. Our implementation includes domain decomposition for parallel computing, a novel stopping criterion, and a system for speeding up convergence by pre-iterating on coarser grids. The performance of the system is demonstrated on a dual quad-core Xeon 3.16 GHz CPUs equipped with 40 GB of RAM. Models of distal tibia derived from 3D in-vivo MR images in a patient comprising 200,000 elements required less than 30 seconds to converge (and 40 MB RAM). To illustrate the system's potential for large-scale μFE simulations, axial stiffness was estimated from high-resolution micro-CT images of a voxel array of 90 million elements comprising the human proximal femur in seven hours CPU time. In conclusion, the system described should enable image-based finite-element bone simulations in practical computation times on high-end desktop computers with applications to laboratory studies and clinical imaging

New Suggestions for the Mechanical Control of Bone Remodeling

Bone is constantly renewed over our lifetime through the process of bone (re)modeling. This process is important for bone to allow it to adapt to its mechanical environment and to repair damage from everyday life. Adaptation is thought to occur through the mechanosensitive response controlling the bone-forming and -resorbing cells. This report shows a way to extract quantitative information about the way remodeling is controlled using computer simulations. Bone resorption and deposition are described as two separate stochastic processes, during which a discrete bone packet is removed or deposited from the bone surface. The responses of the bone-forming and -resorbing cells to local mechanical stimuli are described by phenomenological remodeling rules. Our strategy was to test different remodeling rules and to evaluate the time evolution of the trabecular architecture in comparison to what is known from μ-CT measurements of real bone. In particular, we tested the reaction of virtual bone to standard therapeutic strategies for the prevention of bone deterioration, i.e., physical activity and medications to reduce bone resorption. Insensitivity of the bone volume fraction to reductions in bone resorption was observed in the simulations only for a remodeling rule including an activation barrier for the mechanical stimulus above which bone deposition is switched on. This is in disagreement with the commonly used rules having a so-called lazy zone