Application of design space optimization to bone remodeling simulation of trabecular architecture in human proximal femur for higher computational efficiency

Since the 1990s, topology optimization has been used to computationally investigate the bone remodeling under the assumption that bone remodeling progresses such that the bone material is used in a structurally optimal way. The foremost concerns in applying topology optimization to bone remodeling include the difficulty of handling large-scale problems and associated huge computational cost. In this paper, we applied a recently developed topology optimization algorithm, design space optimization (DSO), to bone remodeling simulation in order to determine trabecular architecture in human proximal femur with higher computational efficiency. We represented the full trabecular architecture in human proximal femur using a two dimensional micro-FE model with 50@mm pixel resolution and performed simulation under three load cases in daily activities. From the quantitative comparison with conventional topology optimization results as well as the actual trabecular architecture, it was shown that DSO produced structurally equivalent trabecular architecture with shorter computing time and smaller memory requirement. As future research, a three-dimensional bone remodeling simulation will require the preparation and management of tens of millions of FE elements, and therefore DSO would be essential to handle such a massive problem as a ''more efficient'' algorithm.

[1]  Hwj Rik Huiskes,et al.  Trabecular Bone Tissue Strains in the Healthy and Osteoporotic Human Femur , 2003, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[2]  In Gwun Jang,et al.  Computational study of Wolff's law with trabecular architecture in the human proximal femur using topology optimization. , 2008, Journal of biomechanics.

[3]  M. Bendsøe Optimal shape design as a material distribution problem , 1989 .

[4]  L. S. Matthews,et al.  Trabecular bone remodeling: an experimental model. , 1991, Journal of biomechanics.

[5]  J. Petersson,et al.  Large-scale topology optimization in 3D using parallel computing , 2001 .

[6]  J. C. Simo,et al.  Numerical instabilities in bone remodeling simulations: the advantages of a node-based finite element approach. , 1995, Journal of biomechanics.

[7]  M. Srámek,et al.  Magnetic resonance imaging of trabecular and cortical bone in mice: comparison of high resolution in vivo and ex vivo MR images with corresponding histology. , 2004, European journal of radiology.

[8]  P Rüegsegger,et al.  Tissue stresses and strain in trabeculae of a canine proximal femur can be quantified from computer reconstructions. , 1999, Journal of biomechanics.

[9]  S. Goldstein,et al.  Variations in Three‐Dimensional Cancellous Bone Architecture of the Proximal Femur in Female Hip Fractures and in Controls , 2000, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[10]  M. Bendsøe,et al.  Generating optimal topologies in structural design using a homogenization method , 1988 .

[11]  K. Svanberg The method of moving asymptotes—a new method for structural optimization , 1987 .

[12]  H Weinans,et al.  A physiological approach to the simulation of bone remodeling as a self-organizational control process. , 1994, Journal of biomechanics.

[13]  N. Kikuchi,et al.  A homogenization sampling procedure for calculating trabecular bone effective stiffness and tissue level stress. , 1994, Journal of biomechanics.

[14]  D M Nunamaker,et al.  Trabecular bone remodeling around smooth and porous implants in an equine patellar model. , 1987, Journal of biomechanics.

[15]  G. Niebur,et al.  High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone. , 2000, Journal of biomechanics.

[16]  Byung Man Kwak,et al.  Design space optimization using a numerical design continuation method , 2002 .

[17]  B. Kwak,et al.  Evolutionary topology optimization using design space adjustment based on fixed grid , 2006 .

[18]  G. Beaupré,et al.  An approach for time‐dependent bone modeling and remodeling—application: A preliminary remodeling simulation , 1990, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[19]  Taiji Adachi,et al.  Functional adaptation of cancellous bone in human proximal femur predicted by trabecular surface remodeling simulation toward uniform stress state. , 2002, Journal of biomechanics.

[20]  Sharmila Majumdar,et al.  Micro-Computed Tomography Evaluation of Trabecular Bone Structure on Loaded Mice Tail Vertebrae , 2003, Spine.

[21]  J. Petersson Some convergence results in perimeter-controlled topology optimization , 1999 .

[22]  Edward J. Haug,et al.  Design Sensitivity Analysis of Structural Systems , 1986 .

[23]  B. Kwak,et al.  Design space optimization using design space adjustment and refinement , 2007 .

[24]  W. J. Whitehouse,et al.  Scanning electron microscope studies of trabecular bone in the proximal end of the human femur. , 1974, Journal of anatomy.

[25]  M. Ashby,et al.  Cellular solids: Structure & properties , 1988 .

[26]  Maarten Moesen,et al.  A symmetry invariant formulation of the relationship between the elasticity tensor and the fabric tensor. , 1985, Mechanics of materials : an international journal.

[27]  Y. Kim,et al.  Parallelized structural topology optimization for eigenvalue problems , 2004 .

[28]  R Huiskes,et al.  Comparison of micro-level and continuum-level voxel models of the proximal femur. , 2006, Journal of biomechanics.

[29]  H. Grootenboer,et al.  The behavior of adaptive bone-remodeling simulation models. , 1992, Journal of biomechanics.

[30]  K. Bathe Finite Element Procedures , 1995 .

[31]  R E Guldberg,et al.  Trabecular bone adaptation to variations in porous-coated implant topology. , 1997, Journal of biomechanics.

[32]  M Evans,et al.  Induction of bone formation in rat tail vertebrae by mechanical loading. , 1993, Bone and mineral.

[33]  S. Cowin,et al.  Wolff's law of trabecular architecture at remodeling equilibrium. , 1986, Journal of biomechanical engineering.

[34]  Glen L Niebur,et al.  Biomechanical effects of intraspecimen variations in tissue modulus for trabecular bone. , 2002, Journal of biomechanics.

[35]  S J Hollister,et al.  Trabecular surface remodeling simulation for cancellous bone using microstructural voxel finite element models. , 2001, Journal of biomechanical engineering.

[36]  M Bagge,et al.  A model of bone adaptation as an optimization process. , 2000, Journal of biomechanics.

[37]  Ole Sigmund,et al.  On the Design of Compliant Mechanisms Using Topology Optimization , 1997 .

[38]  Il Yong Kim,et al.  Analogy of strain energy density based bone-remodeling algorithm and structural topology optimization. , 2009, Journal of biomechanical engineering.

[39]  R. Huiskes,et al.  A new method to determine trabecular bone elastic properties and loading using micromechanical finite-element models. , 1995, Journal of biomechanics.