A model of bone adaptation as a topology optimization process with contact — Oak Academic Publishing
Research ArticleOpen AccessGoogle Scholar indexed
A model of bone adaptation as a topology optimization process with contact
Department of Mechanical Engineering, Centre for Mechanical Technology & Automation, University of Aveiro, Campus Universitário de Santiago, Aveiro, Portugal
,
Department of Mechanical Engineering, Centre for Mechanical Technology & Automation, University of Aveiro, Campus Universitário de Santiago, Aveiro, Portugal
,
Department of Mechanical Engineering, Centre for Mechanical Technology & Automation, University of Aveiro, Campus Universitário de Santiago, Aveiro, Portugal
1 Department of Mechanical Engineering, Centre for Mechanical Technology & Automation, University of Aveiro, Campus Universitário de Santiago, Aveiro, Portugal
2 Department of Mechanical Engineering, Centre for Mechanical Technology & Automation, University of Aveiro, Campus Universitário de Santiago, Aveiro, Portugal
3 Department of Mechanical Engineering, Centre for Mechanical Technology & Automation, University of Aveiro, Campus Universitário de Santiago, Aveiro, Portugal
Topology optimization is presently used in most diverse scientific, technologic and industrial areas, including biomechanics. Bone remodelling models and structural optimization has mutually provided inspiration for new developments in biomechanics and biomedicine. Considering that bone has the ability to adapt its internal structure to mechanical loading (Wolff’s law and Roux’s paradigm), it is possible to model the behaviour of the bone structure by the use of a topology optimization methodology whose optimization variables can be the relative densities and the orthotropic directions. In this work, the internal bone adaptation of a proximal femur is considered. The bone-remodelling scheme is numerically described by a time-dependent evolutionary procedure with anisotropic material parameters. The remodelling rate equation is obtained from the structural optimization task of maximizing the stiffness subject to a biological cost associated with metabolic maintenance of bone tissue in time. The situation of multiple load conditions is considered for a three-dimensional finite element model of the proximal femur. The bone density distribution of a real femur is used as the initial design for the onset of the remodelling mechanism. Examples of bone adaptation resulting from load changes are presented. The three-dimensional finite element model of the proximal femur with initial bone density distribution was adapted to implant a cementless stem. A remeshing technique is used to assign the bone relative density distribution to the new geometry and mesh. The time adaptation of the bone is assessed considering contact with friction at the bone-stem interface. Results of bone density evolution and osteointegration distribution are obtained.
Bens?e, M.P. and Sigmund, O. (2003) Topology optimization—Theory, methods and applications. Springer-Verlag, BH.
Huiskes, R. (2000) If bone is the answer, then what is the question? Journal of Anatomy, 197, 145-156. doi:10.1046/j.1469-7580.2000.19720145.x
Bagge, M. (2000) A model of bone adaptation as an optimization process. Journal of Biomechanics, 33, 1349-1357. doi:10.1016/S0021-9290(00)00124-X
Fernandes, P., Rodrigues, H. and Jacobs, C. (1999) A model of bone adaptation using a global optimisation criterion based on the trajectorial theory of wolff. Computer Methods in Biomechanics and Biomedical Engineering, 2, 125-138. doi:10.1080/10255849908907982
Fernandes, P., Guedes, J.M. and Rodrigues, H. (1999) Topology optimization of 3D linear elastic structures with a constraint on perimeter. Computers & Structures, 73, 583-594. doi:10.1016/S0045-7949(98)00312-5
Fernandes, P.R., Folgado, J., Jacobs, C. and Pellegrini, V. (2002) A contact model with ingrowth control for bone remodelling around cementless stems. Journal of Biomechanics, 35, 167-176. doi:10.1016/S0021-9290(01)00204-4
Folgado, J., Fernandes, P.R., Guedes, J.M. and Rodrigues, H. (2004) Evaluation of osteoporotic bone quality by a computational model for bone remodelling. Computers & Structures, 82, 1381-1388. doi:10.1016/j.compstruc.2004.03.033
Stolk, J., Verdonschot, N. and Huiskes, R. (2001) Hip joint and abductor-muscle forces adequately represent in vivo loading of a cemented total hip reconstruction. Journal of Biomechanics, 34, 917-926. doi:10.1016/S0021-9290(00)00225-6
Waide, V., Cristofolini, L., Stolk, J., Verdonschot, N. and Toni, A. (2003) Experimental investigation of bone remodeling using composite femurs. Clinical Biomechanics, 18, 523-536. doi:10.1016/S0268-0033(03)00072-X
Viceconti, M., Cavalloti, G., Andrisano, A.O. and Toni, A. (1995) Discussion on the design of a hip joint simulator. Medical Engineering & Physics, 18, 234-240. doi:10.1016/1350-4533(95)00026-7
Viceconti, M., Muccini, R., Bernakiewicz, M., Baleani, M. and Cristofolini, L. (2000) Large-sliding contact elements accurately predict levels of bone-implant micromotion relevant to osseointegration. Journal of Biomechanics, 33, 1611-1618. doi:10.1016/S0021-9290(00)00140-8
Viceconti, M., Pancanti, A., Dotti, M., Traina, F., Cristofolini, L. (2004) Effect of the initial implant fitting on the predicted secondary stability of a cementless stem. Medical and Biological Engineering and Computing, 42, 222-229. doi:10.1007/BF02344635
Viceconti, M., Ricci, S., Pancanti, A. and Cappello, A. (2004b) Numerical model to predict the long-term mechanical stability of cementless orthopaedic implants. Medical and Biological Engineering and Computing, 42, 747-753. doi:10.1007/BF02345207
Moreo, P., Péreza, M.A., García-Aznar, J.M. and Doblaré, M. (2007) Modelling the mechanical behaviour of living bony interfaces. Computer Methods in Applied Mechanics and Engineering, 196, 3300-3314. doi:10.1016/j.cma.2007.03.020
Reina, E., Gómez-Benito, M.J., Garcia-Aznar, J.M., Dominguez, J. and Doblaré, M. (2007) Modelo computacional de distracción ósea. Proceeding of CEMNI/ CILAMCE Conference, Sociedad Espa?ola de Métodos Numéricos en Ingeniería (SEMNI), CD-ROM paper 132, 1-10.
Reina, E., Gómez-Benito, M.J., Garcia-Aznar, J.M., Dominguez, J. and Doblaré, M. (2008) Modeling distraction osteogenesis: Analysis of the distraction rate. Biomechanics and Modeling in Mechanobiology, 8, 323-335.
Coelho, P., Fernandes, P.R., Guedes, J.M. and Rodrigues, H.C. (2008) A hierarchical model for concurrent material and topology optimization of three-dimensional structures. Structural and Multidisciplinary Optimization, 35, 107-115. doi:10.1007/s00158-007-0141-3
Coelho, P.G., Fernandes, P.R., Rodrigues, H.C., Cardoso, J.B. and Guedes, J.M. (2009) Numerical modeling of bone tissue adaptation—A hierarchical approach for bone apparent density and trabecular structure. Journal of Biomechanics, 42, 830-837. doi:10.1016/j.jbiomech.2009.01.020
Hardy, D.C., Frayssinet, P., Guilhem, A., Lafontaine, M.A. and Delince, P.E. (1991) Bonding of hydroxyapatite-coated femoral prostheses. Histopathology of specimens from four cases. Journal of Bone & Joint Surgery, 73, 732-740.
Tonino, A.J., Therin, M. and Doyle, C. (1999) Hydroxyapatitecoated femoral stems. Histology and histomorphometry around five components retrieved at post mortem. Journal of Bone & Joint Surgery (British Volume), 81, 148-154. doi:10.1302/0301-620X.81B1.8948
Fernandes, P. (1998) Optimiza??o da Topologia de Estruturas Tridimensionais. PhD Thesis, Instituto Superior Técnico, UTL.
Oliveira, J.A., Pinho-da-Cruz, J.A., Andrade-Campos, A. and Teixeira-Dias, F. (2006) Modelling the elastic behaviour of composites materials with asymptotic expansion homogenisation. Proceeding of COMPTEST, 151-152.
Pinho-da-Cruz, J., Oliveira, J.A., Teixeira-Dias, F. (2009) Asymptotic homogenisation in linear elasticity. Part I: Mathematical formulation and finite element modelling. Computational Materials Science, 45, 1073-1080. doi:10.1016/j.commatsci.2009.02.025
Oliveira, J.A., Pinho-da-Cruz, J. and Teixeira-Dias, F. (2009) Asymptotic homogenisation in linear elasticity. Part II: Finite element procedures and multiscale applications. Computational Materials Science, 45, 1081-1096. doi:10.1016/j.commatsci.2009.01.027
ABAQUS version 6.7 user manual (2006) Hibbitt. Karlsson & Sorenson. Inc.
MSC.MARC (2005) Theory and User Information and User subroutines and special Routines. Msc. Software Corporation.
Svanberg, K. (1987) The method of moving asymptotes – A new method for structural optimization. International Journal for Numerical Methods in Engineering, 24, 358-373. doi:10.1002/nme.1620240207
Medical Multimedia group (1999) A Patient’s Guide to Artificial Hip Replacement. http://www.tri-countyortho.com/mmg/pated/joints/hip/hip_replacement.html. Accessed 15 Jan 2008
Desmorat, B. (2007) Structural rigidity optimization with frictionless unilateral contact. International Journal of Solids and Structures, 44, 1132-1144. doi:10.1016/j.ijsolstr.2006.06.010
Fancello, E.A. (2006) Topology optimization for minimum mass design considering local failure constraints and contact boundary conditions. Structural and Multidisciplinary Optimization, 32, 229-240. doi:10.1007/s00158-006-0019-9
Folgado, J. (2004) Modelos computacionais para análise e projecto e próteses ortopédicas. PhD Thesis, IST-UTL.
Ramos, A. (2006) Estudo numérico e experimental de uma nova componente Femural da prótese de Anca Cimentada. PhD Thesis, Universidade de Aveiro.
Ramos, A., Fonseca, F. and Sim?es, J.A. (2006) Simulation of physiological loading in total hip replacements. Journal of Biomechanical Engineering, 128, 579-588. doi:10.1115/1.2205864
Ramos, A. and Sim?es, J.A. (2006) Tetrahedral versus hexahedral finite elements in numerical modelling of the proximal femur. Medical Engineering & Physics, 28, 916-924. doi:10.1016/j.medengphy.2005.12.006
Egan, J.M. and Marsden, D.C. (2001) A spring network model for the analysis of load transfer and tissue reactions in intramedullary fixation. Clinical Biomechanics, 16, 71-79. doi:10.1016/S0268-0033(00)00070-X
Sadegh, A.M., Luo, G.M. and Cowin, S.C. (1993) Bone ingrowth: An application of the boundary element method to bone remodeling at the implant interface. Journal of Biomechanics, 26, 167-182. doi:10.1016/0021-9290(93)90047-I
Spears, I.R., Pfleiderer, M., Schneider, E., Hille, E., Bergmannn, G. and Morlock, M.M. (2000) Interfacial conditions between a press-fit acetabulax cup and bone during daily activities: Implications for achieving bone in-growth. Journal of Biomechanics, 33, 1471-1477. doi:10.1016/S0021-9290(00)00096-8
Keaveny, T. and Bartel, D. (1993) Mechanical consequences of bone ingrowth in a hip prosthesis inserted without cement. Journal of Bone & Joint Surgery, 77, 911-923.