Human bones optimize their shape according to the mechanical load acting on them. That is, they grow/shrink in regions of high/low mechanical stress, which then again changes the stress distribution. We simulate the interplay between mechanical stress and local growth to improve understanding of this process and corresponding diseases, like osteoporosis. The geometrically complicated micro structure of the bone is described implicitly by a phase field model.