Se efectúa una amplia revisión de los estudios clínicos referentes a la evolución de la densidad mineral ósea (DMO) realizados sobre grupos de pacientes tratadas con distintos fármacos, dosis y plazos. Se define una curva de ajuste de tipo analítico para cada tendencia que, además de regularizar la misma, permite obtener valores numéricos suficientemente precisos en cualquier instante de la evolución analizada, de manera que en cada momento de la evolución se puede obtener la rigidez del hueso (módulo de Young) en función de su densidad. A partir de un modelo de daño mecánico y aplicando resultados conocidos de la Mecánica de la Fractura, se define la probabilidad de fractura. Se combinan todos los resultados anteriores para desarrollar e implementar un código específico que pueda emplearse en combinación con modelos de elementos finitos. Sobre dos modelos de elementos finitos, uno de cuello femoral y otro de columna lumbar, se realiza el correspondiente análisis con el modelo predictivo desarrollado, comparando, bajo diferentes condiciones, la efectividad de los distintos tratamientos farmacológicos aplicados.