LIMITING MODELS FOR CALCIFICATION IN FIBROUS TISSUES ADJACENT TO ORTHOPEDIC IMPLANTS - VARIATIONAL INDICATOR FUNCTIONS AND INFLUENCES OF IMPLANT STIFFNESS
Tp. Harrigan, LIMITING MODELS FOR CALCIFICATION IN FIBROUS TISSUES ADJACENT TO ORTHOPEDIC IMPLANTS - VARIATIONAL INDICATOR FUNCTIONS AND INFLUENCES OF IMPLANT STIFFNESS, Bulletin of mathematical biology, 60(4), 1998, pp. 615-646
Calcification and eventual integration of orthopedic implants into bon
e is important to many load-bearing devices, and the influence of load
and implant stiffness on this process are assessed in this mathematic
al modelling study. Three research questions are posed in this study,
First, can limiting material models provide useful information on the
overall behavior of the tissue adjacent to a loaded orthopedic implant
? Second, can the limiting models lead to optimization criteria? Third
, can an optimization approach be used to differentiate between the fo
ur prospective remodeling rate equations which are proposed? The answe
rs are yes, yes, and no, respectively. A two degree of freedom lumped
parameter model for axial loading of an intramedullary implant is cons
idered. Two limiting composite material models are used, and the strai
n energy density in the calcified and non-calcified phases are assesse
d as stimuli for calcification. The rate equations posed here assume t
hat the calcified material volume fraction decreases at high strain-en
ergy densities, and increases at small strain-energy densities. In all
four cases (both models, both phases) the steady states for these rat
e equations find equilibrium points of indicator functions which are a
weighted sum of total strain energy and the mass of calcified tissue
in the layer considered. The weights on strain-energy density and mass
differ in each case. This shows that for appropriate choices of param
eters, all four models can yield the same results, and it also shows t
hat an optimization approach does not uniquely determine the appropria
te rate equation in these cases. The rate equations showed complicated
dynamic behavior and a phase-plane analysis was used which led to upp
er bounds on load, which depended on implant stiffness and distal supp
ort. The predictions of the four cases studied are compared. (C) 1998
Society for Mathematical Biology.