Research Papers

Physically Realizable Three-Dimensional Bone Prosthesis Design With Interpolated Microstructures

[+] Author and Article Information
Andrew D. Cramer

School of Mathematics and Physics,
The University of Queensland,
Brisbane QLD 4072, Australia
e-mail: a.cramer@uq.edu.au

Vivien J. Challis

School of Mathematics and Physics,
The University of Queensland,
Brisbane QLD 4072, Australia
e-mail: vchallis@maths.uq.edu.au

Anthony P. Roberts

School of Mathematics and Physics,
The University of Queensland,
Brisbane QLD 4072, Australia
e-mail: apr@maths.uq.edu.au

Manuscript received June 13, 2016; final manuscript received November 29, 2016; published online January 24, 2017. Assoc. Editor: Kristen Billiar.

J Biomech Eng 139(3), 031013 (Jan 24, 2017) (8 pages) Paper No: BIO-16-1254; doi: 10.1115/1.4035481 History: Received June 13, 2016; Revised November 29, 2016

We present a new approach to designing three-dimensional, physically realizable porous femoral implants with spatially varying microstructures and effective material properties. We optimize over a simplified design domain to reduce shear stress at the bone-prosthetic interface with a constraint on the bone resorption measured using strain energy. This combination of objective and constraint aims to reduce implant failure and allows a detailed study of the implant designs obtained with a range of microstructure sets and parameters. The microstructure sets are either specified directly or constructed using shape interpolation between a finite number of microstructures optimized for multifunctional characteristics. We demonstrate that designs using varying microstructures outperform designs with a homogeneous microstructure for this femoral implant problem. Further, the choice of microstructure set has an impact on the objective values achieved and on the optimized implant designs. A proof-of-concept metal prototype fabricated via selective laser melting (SLM) demonstrates the manufacturability of designs obtained with our approach.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Zheng, X. , Lee, H. , Weisgraber, T. H. , Shusteff, M. , Deotte, J. , Duoss, E. B. , Kuntz, J. D. , Biener, M. M. , Ge, Q. , Jackson, J. A. , Kucheyev, S. O. , Fang, N. X. , and Spadaccini, C. M. , 2014, “ Ultralight, Ultrastiff Mechanical Metamaterials,” Science, 344(6190), pp. 1373–1377. [CrossRef] [PubMed]
Bendsøe, M. P. , and Sigmund, O. , 2003, Topology Optimization: Theory, Methods, and Applications, Springer Verlag, Berlin.
Rodrigues, H. , Guedes, J. M. , and Bendsoe, M. , 2002, “ Hierarchical Optimization of Material and Structure,” Struct. Multidiscip. Optim., 24(1), pp. 1–10. [CrossRef]
Coelho, P. , Fernandes, P. , Guedes, J. , and Rodrigues, H. , 2008, “ A Hierarchical Model for Concurrent Material and Topology Optimisation of Three-Dimensional Structures,” Struct. Multidiscip. Optim., 35(2), pp. 107–115. [CrossRef]
Cramer, A. D. , Challis, V. J. , and Roberts, A. P. , 2016, “ Microstructure Interpolation for Macroscopic Design,” Struct. Multidiscip. Optim., 53(3), pp. 489–500. [CrossRef]
Khanoki, S. A. , and Pasini, D. , 2012, “ Multiscale Design and Multiobjective Optimization of Orthopedic Hip Implants With Functionally Graded Cellular Material,” ASME J. Biomech. Eng., 134(3), p. 031004. [CrossRef]
Malchau, H. , Herberts, P. , Eisler, T. , Garellick, G. , and Söderman, P. , 2002, “ The Swedish Total Hip Replacement Register,” J. Bone Joint Surg., 84(Suppl 2), pp. S2–S20. [CrossRef]
Kuiper, J. , and Huiskes, R. , 1997, “ Mathematical Optimization of Elastic Properties: Application to Cementless Hip Stem Design,” ASME J. Biomech. Eng., 119(2), pp. 166–174. [CrossRef]
Chanda, S. , Gupta, S. , and Pratihar, D. K. , 2015, “ A Genetic Algorithm-Based Multi-Objective Shape Optimization Scheme for Cementless Femoral Implant,” ASME J. Biomech. Eng., 137(3), p. 034502. [CrossRef]
Engh, C. A. , McGovern, T. , Bobyn, J. , and Harris, W. H. , 1992, “ A Quantitative Evaluation of Periprosthetic Bone-Remodeling After Cementless Total Hip Arthroplasty,” J. Bone Joint Surg. Am., 74(7), pp. 1009–1020. [CrossRef] [PubMed]
Fraldi, M. , Esposito, L. , Perrella, G. , Cutolo, A. , and Cowin, S. , 2010, “ Topological Optimization in Hip Prosthesis Design,” Biomech. Model. Mechanobiol., 9(4), pp. 389–402. [CrossRef] [PubMed]
Hoffman, O. , 1967, “ The Brittle Strength of Orthotropic Materials,” J. Compos. Mater., 1(2), pp. 200–206. [CrossRef]
Challis, V. , Roberts, A. , and Wilkins, A. , 2008, “ Design of Three Dimensional Isotropic Microstructures for Maximized Stiffness and Conductivity,” Int. J. Solids Struct., 45(14), pp. 4130–4146. [CrossRef]
Kang, H. , Lin, C.-Y. , and Hollister, S. J. , 2010, “ Topology Optimization of Three Dimensional Tissue Engineering Scaffold Architectures for Prescribed Bulk Modulus and Diffusivity,” Struct. Multidiscip. Optim., 42(4), pp. 633–644. [CrossRef]
Torquato, S. , and Donev, A. , 2004, “ Minimal Surfaces and Multifunctionality,” Proc. R. Soc. London, Ser. A, 460(2047), pp. 1849–1856. [CrossRef]
Head, W. C. , Bauk, D. J. , and Emerson, R. H., Jr. , 1995, “ Titanium as the Material of Choice for Cementless Femoral Components in Total Hip Arthroplasty,” Clin. Orthop. Relat. Res., 311, pp. 85–90.
Zhang, L. , Klemm, D. , Eckert, J. , Hao, Y. , and Sercombe, T. , 2011, “ Manufacture by Selective Laser Melting and Mechanical Behavior of a Biomedical Ti–24Nb–4Zr–8Sn Alloy,” Scr. Mater., 65(1), pp. 21–24. [CrossRef]
Challis, V. J. , Xu, X. , Zhang, L. C. , Roberts, A. P. , Grotowski, J. F. , and Sercombe, T. B. , 2014, “ High Specific Strength and Stiffness Structures Produced Using Selective Laser Melting,” Mater. Des., 63, pp. 783–788. [CrossRef]
Williams, M. , 1952, “ Stress Singularities Resulting From Various Boundary Conditions,” ASME J. Appl. Mech., 19(4), pp. 526–528.
Bendsøe, M. P. , and Kikuchi, N. , 1988, “ Generating Optimal Topologies in Structural Design Using a Homogenization Method,” Comput. Methods Appl. Mech. Eng., 71(2), pp. 197–224. [CrossRef]
Svanberg, K. , 1987, “ The Method of Moving Asymptotes—A New Method for Structural Optimization,” Int. J. Numer. Methods Eng., 24(2), pp. 359–373. [CrossRef]
Svanberg, K. , 2002, “ A Class of Globally Convergent Optimization Methods Based on Conservative Convex Seperable Approximations,” SIAM J. Optim., 12(2), pp. 555–573.
Camron, H. , Pilliar, R. , and Macnab, I. , 1976, “ The Rate of Bone Ingrowth Into Porous Metal,” J. Biomed. Mater. Res., 10(2), pp. 295–302. [CrossRef] [PubMed]
Boyan, B. D. , Hummert, T. W. , Dean, D. D. , and Schwartz, Z. , 1996, “ Role of Material Surfaces in Regulating Bone and Cartilage Cell Response,” Biomaterials, 17(2), pp. 137–146. [CrossRef] [PubMed]
Khanoki, S. A. , and Pasini, D. , 2013, “ Fatigue Design of a Mechanically Biocompatible Lattice for a Proof-of-Concept Femoral Stem,” J. Mech. Behav. Biomed. Mater., 22, pp. 65–83. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Visualization of the microstructure sets used. Base cells shown at volume fractions of 10%, 30%, 50%, 70%, and 90%. The front octant has been removed to show internal structure. The right most picture visualizes directional dependence of the Young's modulus (E) for the 50% volume fraction microstructure as a fraction of the base material Young's modulus (E0). The axes in the lower left corner indicate orientation of all subfigures. (a) microstructure set A, (b) microstructure set B, (c) microstructure set C, and (d) microstructure set D.

Grahic Jump Location
Fig. 2

Schematics of the domains used for the femoral implant problem. (a) Orthographic schematic of the idealized design domain, with proximal torsion and fixed distal boundary conditions shown, (b) Isometric sketch of ideal domain, (c) 2D slice of rasterized ideal domain (top right quarter only), (d) Isometric sketch of domain with collar added, and (e) 2D slice of rasterized domain with collar (top right quarter only).

Grahic Jump Location
Fig. 3

Prosthetic designs displaying volume fraction, ρ. Slices taken from the back of the prosthesis to the front with the torque forcing clockwise in the page (i.e., torque axis is out of page). All designs shown have a base material with Young's modulus 60 GPa. Black is bone, white is void (a) m = 1, microstructure set A, (b) m = 1, microstructure set C, and (c) m = 2, microstructure set C. Note the filament-like structures with microstructure set A.

Grahic Jump Location
Fig. 4

Stress distributions across the interface Π for: m = 1, microstructure set A (a) homogeneous and (b) inhomogeneous; m = 1, microstructure set C (c) homogeneous and (d) inhomogeneous; m = 2, microstructure set C (e) homogeneous and (f) inhomogeneous. An angle of 0 corresponds to the part of Π where y = 0, x > 0. Stress values have been truncated to the range [103, 106], however, some extreme values are outside this. Vertical striations are an artifact of the discretization. Note the reduction in shear stress concentrated at the proximal end (top) from the homogeneous ((b), (d), and (f)) to the inhomogeneous ((a), (c), and (e)) and the further reduction from the m = 1 case (c) to the m = 2 case (e).

Grahic Jump Location
Fig. 5

Shear stress averaged on circles around the interface surface Π. Microstructure set C with a 60 GPa base material.

Grahic Jump Location
Fig. 6

(a) Two renders of a potential design with microstructure cells 2.86 mm across, the left has one quadrant removed to show internal structure, the right is unmodified. (b) Two renders of a design similar to (a) with microstructure cells that are 800 μm across, giving pores at a scale relevant for bone in-growth [23,24]. (c) and (d) Photographs of the design in (a) manufactured using selective laser melting (SLM) [18].



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In