Engineering, 2014, 6, 113-126 Published Online March 2014 in SciRes. http://www.scirp.org/journal/eng http://dx.doi.org/10.4236/eng.2014.63015 How to cite this paper: Latifi, H., et al. (2014) Computational Simulations of Bone Remodeling under Natural Mechanical Loading or Muscle Malfunction Using Evolutionary Structural Optimization Method. Engineering, 6, 113-126. http://dx.doi.org/10.4236/eng.2014.63015 Computational Simulations of Bone Remodeling under Natural Mechanical Loading or Muscle Malfunction Using Evolutionary Structural Optimization Method Hadi Latifi 1 , Yi Min Xie 1 , Xiaodong Huang 1 , Mehmet Bilgen 2 1 Centre for Innovative Structures and Materials, School of Civil, Environmental and Chemical Engineering, RMIT University, Melbourne, Australia 2 Biophysics Department, Faculty of Medicine, Erciyes University, Kayseri, Turkey Email: Mehmet.Bilgen@yahoo.com Received 25 December 2013; revised 24 January 2014; accepted 3 February 2014 Copyright © 2014 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/ Abstract Live bone inherently responds to applied mechanical stimulus by altering its internal tissue com- position and ultimately biomechanical properties, structure and function. The final formation may structurally appear inferior by design but complete by function. To understand the loading re- sponse, this paper numerically investigated structural remodeling of mature sheep femur using evolutionary structural optimization method (ESO). Femur images from Computed Tomography scanner were used to determine the elastic modulus variation and subsequently construct finite element model of the femur with stiffest elasticity measured. Major muscle forces on dominant phases of healthy sheep gait were imposed on the femur under static mode. ESO was applied to progressively alter the remodeling of numerically simulated femur from its initial to final design by iteratively removing elements with low strain energy density (SED). The computations were repeated with two different mesh sizes to test the convergence. The elements within the medul- lary canal had low SEDs and therefore were removed during the optimization. The SEDs in the re- maining elements varied with angle around the circumference of the shaft. Those elements with low SED were inefficient in supporting the load and thus fundamentally explained how bone re- models itself with less stiff inferior tissue to meet load demand. This was in line with the Wolff’s law of transformation of bone. Tissue growth and remodeling process was found to shape the sheep femur to a mechanically optimized structure and this was initiated by SED in macro-scale according to traditional principle of Wolff’s law.