| Sumario: | La remodelación ósea es un proceso realizado por los huesos para su optimización y reparación. Numerosas enfermedades se asocian a desequilibrios producidos durante el proceso, y para la mayoría de ellas, una detección precoz sumado a un correcto tratamiento puede suponer un punto de inflexión. Es aquí donde los modelos matemáticos y recursos que brinda la ingeniería juegan un papel fundamental, permitiendo el análisis, simulación y estudio de tratamientos cuya finalidad es evitar o reducir los efectos de dichas enfermedades. Una de dichas enfermedades es la osteoporosis. Esta produce una reducción de la densidad ósea, aumentando la porosidad de los huesos y el riesgo a sufrir fracturas. La osteoporosis afecta en mayor medida a las mujeres al cesar la actividad ovárica. Este tipo de osteoporosis se denomina primaria de tipo I u osteoporosis postmenopaúsica. Para combatir la enfermedad se hace uso de tratamientos farmacológicos. Uno de los fármacos empleados es el denosumab, que actúa retrasando el proceso de reabsorción. En el presente trabajo, se pretende mejorar un código que implementa un modelo matemático basado en el proceso de remodelación ósea interna. En el modelo se incluyen los efectos producidos por la osteoporosis postmenopaúsica y la acción del fármaco denosumab. El código de partida implementa las ecuaciones del modelo matemático mencionado en el programa Abaqus. Este programa aplica el método de los elementos finitos, que discretiza un dominio en elementos y resuelve las ecuaciones diferenciales planteadas en los puntos de integración presentes en los elementos. El programa Abaqus será utilizado para realizar las simulaciones y su posterior análisis. Para la simulación del proceso es necesario el uso de subrutinas en el código, ya que estas permiten modelar el comportamiento mecánico del hueso. Las mejoras a dicho código se realizan desde dos enfoques. Primero se realiza un estudio local y a continuación uno no local. El objetivo principal del estudio local consiste en la obtención de las distintas fases que tienen lugar en el proceso, aunque también se implementarán otras modificaciones secundarias. El estudio no local presenta como finalidad la activación-desactivación discreta del conjunto de células que desarrollan el proceso y la obtención de una frecuencia de activación de las mismas que se asemeje a la reflejada en la literatura. Este segundo enfoque no local se desarrolla estudiando el entorno cercano a cada punto de integración presente en el dominio. Además, se lleva a cabo un análisis de sensibilidad sobre la influencia de la variación de algunos parámetros que forman parte del modelo matemático. El análisis permitirá ajustar los valores para los modelos generados, obteniendo unas simulaciones que se asemejan en mayor medida al proceso real. Los resultados obtenidos para ambos modelos son satisfactorios. Aunque durante el desarrollo del modelo no local se encuentra un caso particular que permitirá nuevas líneas de desarrollo e investigación.
|
|---|
