Despite recent need-based advances in orthopedic scaffold design, current implants are unsuitable as “total” scaffold replacements. Both mechanical requirements of stiffness/strength and biological stipulations dictating cellular behavior (attachment, differentiation) should be included. The amount of mechanical stimulation in the form of stresses, strains, and energies most suitable toward implant design is presently unknown. Additionally unknown is if whole-bone optimization goals such as uniform and non-uniform driving forces are applicable to a scaffold-bone interface. At the very least, scaffolds ready for implantation should exhibit mechanical distributions (dependent on loading type) on the surface within the typical mechanical usage window. Scaffold micro-architectures can be strategically shifted into that window. The overall goal of this study was to produce microarchitectures tailored to a more uniform mechanical distribution, while maintaining the morphological properties necessary to sustain its mechanical integrity. The mechanical adjustment stimuli investigated were von Mises stress, strain energy density, maximum principle strain, and volumetric strain. Scaffold models of a similar volume fraction were generated of three initial architectures (Rhombitruncated Cuboctahedron, hollow sphere, and trabecular-like bone cube) using high resolution voxel mapping. The resulting voxels were translated into finite element meshes and solved, with a specially written iterative solver created in Fortran90, under confined displacement boundary conditions. The result was verified against a commercial software. Once the mechanical distributions were identified one of two methods was chosen to alter the configuration of material in Cartesian space. The success of the alteration was judged through a diagnostic based on the histogram of mechanical values present on the surface of the micro-architecture. The first method used a compliant approach and, for the case of stress, reinforced locations on the surface with large stresses with extra material (strategically taken from the least stressed portions). The second method used a simulated annealing approach to randomly mutate the initial state in a “temperature” dependent manner. Results indicate that the mechanical distributions of the initial scaffold designs vary significantly. Additionally, the end state of the adjustment demonstrated anisotropy shifts toward the direction of loading. Moreover, the adjustment methods were found to be sensitive both to the mechanical parameter used for adjustment and the portion of the surface adjusted at each increment. In conclusion, scaffolds may be adjusted using a mechanical surface-based objective, as the surface of the scaffold is crucial toward its in vivo acceptance. This technique provides some mathematical specificity toward the whole of computer-aided tissue engineering.

1.
Mullender
M. G.
and
Huiskes
R.
,
Proposal for the regulatory mechanism of Wolff’s law
.
J Orthop Res
,
1995
.
13
(
4
): p.
503
12
.
2.
Hart, R., Bone Modeling and Remodeling: Theories and Computation, in Bone Biomechanics Handbook, S.C. Cowin, Editor. 2001, CRC Press LLC: Boca Raton. p. (21)1–(31)42.
3.
Adachi
T.
,
Tanaka
M.
, and
Tomita
Y.
,
Uniform stress state in bone structure with residual stress
.
J Biomech Eng
,
1998
.
120
(
3
): p.
342
7
.
4.
Huiskes
R.
, et al.,
Effects of mechanical forces on maintenance and adaptation of form in trabecular bone
.
Nature
,
2000
.
405
(
6787
): p.
704
6
.
5.
Cowin
S. C.
and
Weinbaum
S.
,
Strain amplification in the bone mechanosensory system
.
Am J Med Sci
,
1998
.
316
(
3
): p.
184
8
.
6.
Hung
C. T.
, et al.,
Real-time calcium response of cultured bone cells to fluid flow
.
Clin Orthop
,
1995
(
313
): p.
256
69
.
7.
Mattheck
C.
,
Design and growth rules for biological structures and their application to engineering
.
Fatigue & Fracture of Engineering Materials and Structure
,
1990
.
13
(
5
): p.
535
550
.
8.
Mattheck, C., Design in Nature, ed. W. Linnard. 1998: Springer-Verlag Berlin and Heidelberg GmbH & Co.
9.
Hollister
S. J.
,
Maddox
R. D.
, and
Taboas
J. M.
,
Optimal design and fabrication of scaffolds to mimic tissue properties and satisfy biological constraints
.
Biomaterials
,
2002
.
23
(
20
): p.
4095
4103
.
10.
Lin
C. Y.
,
Kikuchi
N.
, and
Hollister
S. J.
,
A novel method for biomaterial scaffold internal architecture design to match bone elastic properties with desired porosity
.
Journal of Biomechanics
,
2004
.
37
(
5
): p.
623
636
.
11.
Lin
C. Y.
, et al.,
Interbody fusion cage design using integrated global layout and local microstructure topology optimization
.
Spine
,
2004
.
29
(
16
): p.
1747
1754
.
12.
Sun
W
and
Lal
P.
,
Recent development on computer aided tissue engineering–a review
.
Comput Methods Programs Biomed
,
2002
.
67
(
2
): p.
85
103
.
13.
Yang
S.
, et al.,
The design of scaffolds for use in tissue engineering. Part II. Rapid prototyping techniques
.
Tissue Eng
,
2002
.
8
(
1
): p.
1
11
.
14.
Mironov
V.
, et al.,
Organ printing: computer-aided jet-based 3D tissue engineering
.
Trends Biotechnol
,
2003
.
21
(
4
): p.
157
61
.
15.
Bendsoe, M.P. and O. Sigmund, Topology Optimization. Theory, methods and Applications. 2003, Berlin, Heidelberg, New York: Springer Verlag.
16.
Bendsoe
M. P.
and
Kikuchi
N.
,
Generating Optimal Topologies in Structural Design Using a Homogenization Method
.
Computer Methods in Applied Mechanics and Engineering
,
1988
.
71
(
2
): p.
197
224
.
17.
Sethian
J. A.
and
Wiegmann
A.
,
Structural boundary design via level set and immersed interface methods
.
Journal of Computational Physics
,
2000
.
163
(
2
): p.
489
528
.
18.
Belytschko
T.
,
Xiao
S. P.
, and
Parimi
C.
,
Topology optimization with implicit functions and regularization
.
International Journal for Numerical Methods in Engineering
,
2003
.
57
(
8
): p.
1177
1196
.
19.
Van Rietbergen
B.
, et al.,
Tissue stresses and strain in trabeculae of a canine proximal femur can be quantified from computer reconstructions
.
J Biomech
,
1999
.
32
(
4
): p.
443
51
.
20.
Hughes
T. J. R.
,
Levit
I.
, and
Winget
J.
,
An Elementby-Element Solution Algorithm for Problems of Structural and Solid Mechanics
.
Computer Methods in Applied Mechanics and Engineering
,
1983
.
36
(
2
): p.
241
254
.
21.
Kwon, Y.W. and H. Bang, The Finite Element Method using MATLAB. 1997, New York: CRC Press LLC.
22.
Ingber
L.
and
Rosen
B.
,
Genetic Algorithms and Very Fast Simulated Reannealing - a Comparison
.
Mathematical and Computer Modelling
,
1992
.
16
(
11
): p.
87
100
.
23.
Ingber
L.
,
Simulated Annealing - Practice Versus Theory
.
Mathematical and Computer Modelling
,
1993
.
18
(
11
): p.
29
57
.
24.
Kirkpatrick
S.
,
Gelatt
C. D.
, and
Vecchi
M. P.
,
Optimization by Simulated Annealing
.
Science
,
1983
.
220
(
4598
): p.
671
680
.
This content is only available via PDF.
You do not currently have access to this content.