Abstract
While physiological loading on lower long bones changes during bone development, the bone cross section either remains circular or slowly changes from nearly circular to other shapes such as oval and roughly triangular. Bone is said to be an optimal structure, where strength is maximized using the optimal distribution of bone mass (also called Wolff’s law). One of the most appropriate mathematical validations of this law would be a structural optimization-based formulation where total strain energy is minimized against a mass and a space constraint. Assuming that the change in cross section during bone development and homeostasis after adulthood is direct result of the change in physiological loading, this work investigates what optimization problem formulation (collectively, design variables, objective function, constraints, loading conditions, etc.) results in mathematically optimal solutions that resemble bones under actual physiological loading. For this purpose, an advanced structural optimization-based computational model for cortical bone development and defect repair is presented. In the optimization problem, overall bone stiffness is maximized first against a mass constraint, and then also against a polar first moment of area constraint that simultaneously constrains both mass and space. The investigation is completed in two stages. The first stage is developmental stage when physiological loading on lower long bones (tibia) is a random combination of axial, bending and torsion. The topology optimization applied to this case with the area moment constraint results into circular and elliptical cross sections similar to that found in growing mouse or human. The second investigation stage is bone homeostasis reached in adulthood when the physiological loading has a fixed pattern. A drill hole defect is applied to the adult mouse bone, which would disrupt the homeostasis. The optimization applied after the defect interestingly brings the damaged section back to the original intact geometry. The results, however, show that cortical bone geometry is optimal for the physiological loading only when there is also a constraint on polar moment of area. Further numerical experiments show that application of torsion along with the gait-analysis-based physiological loading improves the results, which seems to indicate that the cortical bone geometry is optimal for some amount of torsion in addition to the gait-based physiological loading. This work has a potential to be extended to bone growth/development models and fracture healing models, where topology optimization and polar moment of area constraint have not been introduced earlier.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Andrade-Campos A, Ramos A, Simões JA (2012) A model of bone adaptation as a topology optimization process with contact. J Biomed Sci Eng 05:229–244. https://doi.org/10.4236/jbise.2012.55030
Andreassen E, Clausen A, Schevenels M, Lazarov BS, Sigmund O (2011) Efficient topology optimization in MATLAB using 88 lines of code. Struct Multidiscip Optim 43:1–16. https://doi.org/10.1007/s00158-010-0594-7
Bagge M (2000) A model of bone adaptation as an optimization process. J Biomech 33:1349–1357. https://doi.org/10.1016/S0021-9290(00)00124-X
Bendsoe MP, Sigmund O (2004) Topology optimization: theory, methods, and applications, 2nd edn. Springer, Berlin. https://doi.org/10.1007/978-3-662-05086-6
Bendsøe MP, Sigmund O (2003) Topology optimization: theory, methods and applications. Springer, Berlin
Bertram JEA, Swartz SM (1991) The ‘law of bone transformation’: a case of crying Wolff? Biol Rev 66:245–273. https://doi.org/10.1111/j.1469-185X.1991.tb01142.x
Campbell TM, Wong WT, Mackie EJ (2003) Establishment of a model of cortical bone repair in mice. Calcif Tissue Int 73:49–55. https://doi.org/10.1007/s00223-002-2120-4
Capozza RF et al (2010) Structural analysis of the human tibia by tomographic (pQCT) serial scans. J Anat 216:470–481. https://doi.org/10.1111/j.1469-7580.2009.01201.x
Carter DR, Van der Meulen MCH, Beaupré GS (1996) Mechanical factors in bone growth and development. Bone 18:S5–S10. https://doi.org/10.1016/8756-3282(95)00373-8
Chao EY, Inoue N (2003) Biophysical stimulation of bone fracture repair, regeneration and remodelling. Eur Cells Mater 6:72–84 (discussion 84-75)
Clarke KA, Still J (2001) Development and consistency of gait in the mouse. Physiol Behav 73:159–164
Colloca M (2010) Bone remodeling process as an optimal structural design. Dessertation. Sapienza University of Rome
Cowgill LW, Warrener A, Pontzer H, Ocobock C (2010) Waddling and toddling: the biomechanical effects of an immature gait. Am J Phys Anthropol 143:52–61. https://doi.org/10.1002/ajpa.21289
Diaz A, Sigmund O (1995) Checkerboard patterns in layout optimization. Struct Optim 10:40–45. https://doi.org/10.1007/BF01743693
Dini P, David A (2009) Repeatability of spatiotemporal gait parameters: comparison between normal children and children with hemiplegic spastic cerebral palsy. Braz J Phys Ther 13:215–222
Diop M, Rahmani A, Belli A, Gautheron V, Geyssant A, Cottalorda J (2005) Influence of speed variation and age on ground reaction forces and stride parameters of children’s normal gait. Int J Sports Med 26:682–687. https://doi.org/10.1055/s-2004-830382
Ethier CR, Simmons CA (2007) Introductory biomechanics: from cells to organisms. Cambridge Texts in Biomedical Engineering. Cambridge University Press, Cambridge. https://doi.org/10.1017/cbo9780511809217
Feather-Schussler DN, Ferguson TS (2016) A battery of motor tests in a neonatal mouse model of cerebral palsy. J Vis Exp. https://doi.org/10.3791/53569
Fernandes P, Rodrigues H, Jacobs C (1999) A model of bone adaptation using a global optimisation criterion based on the trajectorial theory of Wolff. Comput Methods Biomech Biomed Eng 2:125–138. https://doi.org/10.1080/10255849908907982
Fowler SC, Birkestrand BR, Chen R, Moss SJ, Vorontsova E, Wang G, Zarcone TJ (2001) A force-plate actometer for quantitating rodent behaviors: illustrative data on locomotion, rotation, spatial patterning, stereotypies, and tremor. J Neurosci Methods 107:107–124. https://doi.org/10.1016/S0165-0270(01)00359-4
Francu J, Novácková P (2012) Torsion of a non-circular bar. Eng Mech 19:45–60
Franklyn M, Oakes B, Field B, Wells P, Morgan D (2008) Section modulus is the optimum geometric predictor for stress fractures and medial tibial stress syndrome in both male and female athletes. Am J Sports Med 36:1179–1189. https://doi.org/10.1177/0363546508314408
Frost HM (1994) Wolff’s Law and bone’s structural adaptations to mechanical usage: an overview for clinicians. Angle Orthod 64:175–188. https://doi.org/10.1043/0003-3219(1994)064%3c0175:WLABSA%3e2.0.CO;2
Frost HM (2001) From Wolff’s law to the Utah paradigm: insights about bone physiology and its clinical applications. Anat Rec 262:398–419. https://doi.org/10.1002/ar.1049
Frost HM (2004) A 2003 update of bone physiology and Wolff’s law for clinicians. Angle Orthod 74:3–15. https://doi.org/10.1043/0003-3219(2004)074%3c0003:auobpa%3e2.0.co;2
Galea GL, Hannuna S, Meakin LB, Delisser PJ, Lanyon LE, Price JS (2015) Quantification of alterations in cortical bone geometry using site specificity software in mouse models of aging and the responses to ovariectomy and altered loading. Front Endocrinol 6:52. https://doi.org/10.3389/fendo.2015.00052
Goldman HM, McFarlin SC, Cooper DM, Thomas CD, Clement JG (2009) Ontogenetic patterning of cortical bone microstructure and geometry at the human mid-shaft femur. Anat Rec 292:48–64. https://doi.org/10.1002/ar.20778
Gosman JH, Hubbell ZR, Shaw CN, Ryan TM (2013) Development of cortical bone geometry in the human femoral and tibial diaphysis. Anat Rec 296:774–787. https://doi.org/10.1002/ar.22688
Harrigan TP, Hamilton JJ (1994) Bone remodeling and structural optimization. J Biomech 27:323–328. https://doi.org/10.1016/0021-9290(94)90008-6
Hart NH, Nimphius S, Rantalainen T, Ireland A, Siafarikas A, Newton R (2017) Mechanical basis of bone strength: influence of bone material, bone structure and muscle action. J Musculoskelet Neuronal Interact 17:114
Hibbeler RC (1994) Mechanics of materials. Macmillan College Pub., Kerala
Hollister SJ, Kikuchi N, Goldstein SA (1993) Do bone ingrowth processes produce a globally optimized structure? J Biomech 26:391–407. https://doi.org/10.1016/0021-9290(93)90003-W
Huang X, Xie YM (2011) Evolutionary topology optimization of continuum structures including design-dependent self-weight loads. Finite Elem Anal Des 47:942–948. https://doi.org/10.1016/j.finel.2011.03.008
Huiskes R (1997) Simulation of self-organization and functional adaptation in bone Hefte zu der. Unfallchirurg 261:299–320
Huiskes R (2000) If bone is the answer, then what is the question? J Anat 197:145–156. https://doi.org/10.1046/j.1469-7580.2000.19720145.x
Jacobs BY, Kloefkorn HE, Allen KD (2014) Gait analysis methods for rodent models of osteoarthritis. Curr Pain Headache Rep 18:456. https://doi.org/10.1007/s11916-014-0456-x
Jang IG, Kim IY (2008) Computational study of Wolff’s law with trabecular architecture in the human proximal femur using topology optimization. J Biomech 41:2353–2361. https://doi.org/10.1016/j.jbiomech.2008.05.037
Jang IG, Kim IY (2009) Computational simulation of trabecular adaptation progress in human proximal femur during growth. J Biomech 42:573–580. https://doi.org/10.1016/j.jbiomech.2008.12.009
Jang IG, Kim IY, Kwak BM (2008) Analogy of strain energy density based bone-remodeling algorithm and structural topology optimization. J Biomech Eng 131:011012–011001–011012–011017. https://doi.org/10.1115/1.3005202
Jang IG, Kim IY, Kwak BM (2009) Analogy of strain energy density based bone-remodeling algorithm and structural topology optimization. J Biomech Eng 131:7. https://doi.org/10.1115/1.3005202
Karihaloo BL, Hemp WS (1987) Optimum sections for given torsional and flexural rigidity. Proc R Soc Lond A Math Phys Sci 409:67
Land C, Schoenau E (2008) Fetal and postnatal bone development: reviewing the role of mechanical stimuli and nutrition. Best Pract Res Clin Endocrinol Metab 22:107–118. https://doi.org/10.1016/j.beem.2007.09.005
Lee JW, Kawahara K, Nakano T (2010) Bone regeneration based on orientation of biological apatite (bap) c-axis in osteopetrotic (op/op) mice. Mater Sci Forum 638–642:588–593. https://doi.org/10.4028/www.scientific.net/MSF.638-642.588
Lieberman DE, Polk JD, Demes B (2004) Predicting long bone loading from cross-sectional geometry. Am J Phys Anthropol 123:156–171. https://doi.org/10.1002/ajpa.10316
Lim SS, Lee B, Kim IS, Hwang SJ (2017) Differential modulation of zoledronate and etidronate in osseous healing of an extracted socket and tibia defect. Oral Surg Oral Med Oral Pathol Oral Radiol 123:8–19. https://doi.org/10.1016/j.oooo.2016.08.009
Liu L et al (2012) Rosiglitazone inhibits bone regeneration and causes significant accumulation of fat at sites of new bone formation. Calcif Tissue Int 91:139–148. https://doi.org/10.1007/s00223-012-9623-4
Liu C et al (2018) Effects of mechanical loading on cortical defect repair using a novel mechanobiological model of bone healing. Bone 108:145–155. https://doi.org/10.1016/j.bone.2017.12.027
Matsumoto T, Goto D, Sato S (2013) Subtraction micro-computed tomography of angiogenesis and osteogenesis during bone repair using synchrotron radiation with a novel contrast agent. Lab Investig J Tech Methods Pathol 93:1054–1063. https://doi.org/10.1038/labinvest.2013.87
Matsumoto T, Sato D, Hashimoto Y (2016) Individual and combined effects of noise-like whole-body vibration and parathyroid hormone treatment on bone defect repair in ovariectomized mice. Proc Inst Mech Eng [H] 230:30–38. https://doi.org/10.1177/0954411915616987
Meadows TH, Bronk JT, Chao YS, Kelly PJ (1990) Effect of weight-bearing on healing of cortical defects in the canine tibia. J Bone Joint Surg Am 72:1074–1080
Milgrom C et al (1989) The area moment of inertia of the tibia: a risk factor for stress fractures. J Biomech 22:1243–1248. https://doi.org/10.1016/0021-9290(89)90226-1
Monfoulet L, Rabier B, Chassande O, Fricain J-C (2010) Drilled hole defects in mouse femur as models of intramembranous cortical and cancellous bone regeneration. Calcif Tissue Int 86:72–81. https://doi.org/10.1007/s00223-009-9314-y
Piccinini M, Cugnoni J, Botsis J, Ammann P, Wiskott A (2014) Influence of gait loads on implant integration in rat tibiae: experimental and numerical analysis. J Biomech 47:3255–3263. https://doi.org/10.1016/j.jbiomech.2014.08.023
Prasad J, Wiater BP, Nork SE, Bain SD, Gross TS (2010) Characterizing gait induced normal strains in a murine tibia cortical bone defect model. J Biomech 43:2765–2770. https://doi.org/10.1016/j.jbiomech.2010.06.030
Rietz A (2001) Sufficiency of a finite exponent in SIMP (power law) methods. Struct Multidiscip Optim 21:159–163. https://doi.org/10.1007/s001580050180
Rossi J-M, Wendling-Mansuy S (2007) A topology optimization based model of bone adaptation. Comput Methods Biomech Biomed Eng 10:419–427. https://doi.org/10.1080/10255840701550303
Rubin CT, Lanyon LE (1985) Regulation of bone mass by mechanical strain magnitude. Calcif Tissue Int 37:411–417
Serradj N, Jamon M (2007) Age-related changes in the motricity of the inbred mice strains 129/sv and C57BL/6j. Behav Brain Res 177:80–89. https://doi.org/10.1016/j.bbr.2006.11.001
Shaw CN, Stock JT (2009) Intensity, repetitiveness, and directionality of habitual adolescent mobility patterns influence the tibial diaphysis morphology of athletes. Am J Phys Anthropol 140:149–159. https://doi.org/10.1002/ajpa.21064
Sigmund O, Petersson J (1998) Numerical instabilities in topology optimization: a survey on procedures dealing with checkerboards, mesh-dependencies and local minima. Struct Optim 16:68–75. https://doi.org/10.1007/bf01214002
Tiwari AK, Prasad J (2016) Computer modelling of bone’s adaptation: the role of normal strain, shear strain and fluid flow. Biomech Model Mechanobiol. https://doi.org/10.1007/s10237-016-0824-z
Tsubota K, Suzuki Y, Yamada T, Hojo M, Makinouchi A, Adachi T (2009) Computer simulation of trabecular remodeling in human proximal femur using large-scale voxel FE models: approach to understanding Wolff’s law. J Biomech 42:1088–1094. https://doi.org/10.1016/j.jbiomech.2009.02.030
Uusitalo H et al (2001) A metaphyseal defect model of the femur for studies of murine bone healing. Bone 28:423–429. https://doi.org/10.1016/S8756-3282(01)00406-9
Weidauer LA, Binkley T, Vukovich M, Specker B (2014) Greater polar moment of inertia at the tibia in athletes who develop stress fractures. Orthop J Sports Med 2:2325967114541411. https://doi.org/10.1177/2325967114541411
Wolff J (1893) Das gesetz der transformation der knochen. DMW-Deutsche Medizinische Wochenschrift 19:1222–1224
Xinghua Z, He G, Bingzhao G (2005) The application of topology optimization on the quantitative description of the external shape of bone structure. J Biomech 38:1612–1620. https://doi.org/10.1016/j.jbiomech.2004.06.029
Yang P-F, Sanno M, Ganse B, Koy T, Brüggemann G-P, Müller LP, Rittweger J (2014) Torsion and antero-posterior bending in the in vivo human tibia loading regimes during walking and running. PLoS ONE 9:e94525. https://doi.org/10.1371/journal.pone.0094525
Zumwalt AC, Hamrick M, Schmitt D (2006) Force plate for measuring the ground reaction forces in small animal locomotion. J Biomech 39:2877–2881. https://doi.org/10.1016/j.jbiomech.2005.10.006
Acknowledgements
Authors acknowledge IIT Ropar for the Institute Fellowship. The authors also acknowledge the computational facilities provided by Science and Engineering Board (SERB), Department of Science and Technology (DST), Government of India via Grant No. SB/S3/MMER/046/2013 and also by Defense Research and Development Organization (DRDO) via Grant No. TBRL/CARS/75/2016.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Sen, C., Prasad, J. Exploring conditions that make cortical bone geometry optimal for physiological loading. Biomech Model Mechanobiol 18, 1335–1349 (2019). https://doi.org/10.1007/s10237-019-01147-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10237-019-01147-z