LIMITING MODELS FOR CALCIFICATION IN FIBROUS TISSUES ADJACENT TO ORTHOPEDIC IMPLANTS - VARIATIONAL INDICATOR FUNCTIONS AND INFLUENCES OF IMPLANT STIFFNESS

Authors
Citation
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
Citations number
26
Categorie Soggetti
Mathematics, Miscellaneous","Biology Miscellaneous","Mathematics, Miscellaneous
ISSN journal
00928240
Volume
60
Issue
4
Year of publication
1998
Pages
615 - 646
Database
ISI
SICI code
0092-8240(1998)60:4<615:LMFCIF>2.0.ZU;2-A
Abstract
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.