Computer simulations of bone remodeling in response to mechanical stresses can be used to understand normal growth and development of the skeleton or to predict the remodeling of bone in response to prosthetic devices. Using a previously derived bone maintenance theory, a technique for computing bone density distributions was applied to the proximal femur and tibia using two-dimensional, multiple-loading finite element models. The models initially represented solid, homogeneous structures. Using an iterative bone remodeling technique that relates bone apparent density to loading history, the internal distributions of apparent density and elastic modulus for the normal bones were predicted. The finite element models were then modified to represent bones in which porous-coated femoral surface replacements and tibial tray components had been implanted. The same iterative remodeling method was then applied to predict the distribution of bone around these components. The predicted bone density distributions for the natural femur and tibia agree with previously documented normal bone morphology. The predicted bone density distributions around various implanted prostheses were characteristic of the component under investigation and were consistent with clinical and experimental findings of other investigators. In the femoral head, stress shielding occurred underneath the metal surface replacement cup, resulting in lower densities in the femoral head. The addition of a central femoral cup fixation peg caused bone hypertrophy around the peg. In the tibia, the stress concentrations around the pegs also resulted in denser bone, with a concomitant decrease in bone density at more peripheral locations underneath the prosthetic tray. This remodeling technique has the potential to be an important tool in predicting the possible remodeling consequences of new implant design features.