Structure-function relationships for coralline hydroxyapatite bone substitute.

To improve the understanding of the functional requirements of trabecular bone substitutes, the structure-function relationships of coralline hydroxyapatite were determined and compared to those of trabecular bone from a variety of anatomic sites. Mechanical properties and permeability of cylindrical coralline hydroxyapatite specimens were measured and related to various morphological parameters that were obtained from analysis of high-resolution (20 microm) computer reconstructions of each specimen. Results indicated the average (+/-SD) Young's modulus (2900 +/- 1290 MPa, n = 20) and permeability (0.50 +/- 0.19 x 10(-9) m2, n = 21) of the coralline hydroxyapatite were within the range of values exhibited by high density trabecular bone; ultimate stress (5.87 +/- 1.92 MPa, n = 13), while in the range of mid-density trabecular bone, was low considering its high volume fraction (31.3 +/- 1.9%, n = 49); and ultimate strain (0.22 +/- 0.03%, n = 13) was much lower than that of trabecular bone from any anatomic site. The only correlation found between mechanical and morphological parameters was between Young's modulus and "fabric" (a scalar measure of architecture that combined the degree of microstructural anisotropy with orientation). These results provide insight into the in vivo performance of this implant, as well as the biomechanical requirements for successful trabecular bone substitutes in general.

Application of the Tsai-Wu quadratic multiaxial failure criterion to bovine trabecular bone.

As a first step toward development of a multiaxial failure criterion for human trabecular bone, the Tsai-Wu quadratic failure criterion was modified as a function of apparent density and applied to bovine tibial trabecular bone. Previous data from uniaxial compressive, tensile, and torsion tests (n = 139 total) were combined with those from new triaxial tests (n = 17) to calibrate and then verify the criterion. Combinations of axial compression and radial pressure were used to produce the triaxial compressive stress states. All tests were performed with minimal end artifacts in the principal material coordinate system of the trabecular network. Results indicated that the stress interaction term F12 exhibited a strong nonlinear dependence on apparent density (r2 > 0.99), ranging from -0.126 MPa-2 at low densities (0.29 g/cm3) to 0.005 MPa-2 at high densities (0.63 g/cm3). After calibration and when used to predict behavior of new-specimens without any curve-fitting, the Tsai-Wu criterion had a mean (+/- SD) error of -32.6 +/- 10.6 percent. Except for the highest density triaxial specimens, most (15/17 specimens) failed at axial stresses close to their predicted uniaxial values, and some reinforcement for transverse loading was observed. We conclude that the Tsai-Wu quadratic criterion, as formulated here, is at best only a reasonable predictor of the multiaxial failure behavior of trabecular bone, and further work is required before it can be confidently applied to human bone.

Implementation of strain rate as a bone remodeling stimulus.

Strain rate is implemented as a stimulus for surface bone remodeling. Using idealized models for trabecular bone structures, the surface remodeling predictions using the strain rate as the stimulus are compared with the predictions using the peak strain magnitude as the stimulus. For a uniaxially loaded cruciform shape, the comparison shows that the two surface remodeling stimuli predict the same final shape under a periodic compressive load, but the two evolutionary paths to final shapes are different. Two biaxially loaded regular grid models of trabecular structure were considered, one a grid of square diamond shaped elements and the other a brick wall patterned grid. For both of these idealized trabecular structures, the comparison shows that the two surface remodeling stimuli predict the same final shape under a periodic compressive load, even from these distinctly different initial grid patterns, and the evolutionary paths to final shapes are quite different. In general the two stimuli do not predict the same remodeling and the conditions under which they do are derived. The models developed are also applied to the data from the animal experiments reported in Goldstein et al. (1991), and it is shown that the strain rate stimulus predicts bone remodeling similar to what was experimentally observed.

Chaos in the discrete-time algorithm for bone-density remodeling rate equations.

We compare the predictions of the differential equation form of a class of bone-density stress adaptation models with their associated discrete-time computational algorithms. Although our considerations apply to the class of adaptation models based on bulk or apparent bone-density remodeling, we focus attention on a particular model in this class, a model employed by Weinans et al. [Trans. Orthop. Res. Soc. 14, 310 (1989); Trans. First World Congress of Biomechanics, Vol. II, p. 75 (1990)]. We show that the discrete-time computational algorithm of that stress adaptation model has a well-known chaos mechanism for stress values of practical interest. Further, we obtain a condition on the discrete-time step that prevents the transition to chaos, and conditions that insure monotonic convergence. This chaos mechanism is only present in the discrete-time computational algorithm; we show that the corresponding differential equation form of the bone-density stress adaptation model is smooth, monotonic and nonchaotic.