Paleontología Virtual

Día 6: Análisis Biomecánico (FEA)

E. Miguel Díaz de León Muñoz

2026-04-24

Fundamentos del FEA

¿Qué es el Análisis de Elementos Finitos (FEA)?

Simulando la física de lo extinto

  • El FEA es un método computacional para predecir cómo un objeto responde a fuerzas físicas (Bright 2014).
  • Esfuerzo (Stress): Fuerza interna que las partículas de un material ejercen entre sí.
  • Deformación (Strain): Cambio en la forma del material debido al esfuerzo.
  • En Paleontología: Permite testear hipótesis sobre dieta, fuerza de mordida y resistencia mecánica de huesos fósiles.

Ejemplo de FEA en craneos de crocodilomorfos

El FEA nació en la ingeniería, llegó a los fósiles

De aviones a dinosaurios

  • 1940s–1950s: Desarrollado para análisis estructural en ingeniería aeroespacial y civil.
  • Turner, Clough, Martin & Topp (1956): Primera aplicación formal — análisis de esfuerzos en alas de avión.
  • Clough (1960): Acuñó el término “finite elements”.
  • Paleontología: Rayfield et al. (2001) — primer análisis en cráneo fósil (Allosaurus fragilis).
  • Hoy: Herramienta estándar para estimar fuerza de mordida, resistencia ósea y modos de locomoción.

Análisis de elementos finitos en estructura de ala de avión — Boeing (1950s)

Cómo funciona: la analogía de la hamaca

  • Resolver la física de un sólido continuo de forma exacta es matemáticamente inabordable.
  • Solución: Dividir el objeto en miles de pequeños elementos simples (discretización).
  • Cada elemento sigue las leyes de elasticidad de forma local.
  • La suma de todos los elementos reproduce el comportamiento global del sólido.
  • La malla es esa red: cuanto más fina, más fiel la simulación.

La malla del modelo FEA funciona como la red de una hamaca: la carga se distribuye por los nodos

Nodos, elementos y malla

Anatomía de un modelo FEA

  • Nodo: Punto de conexión entre elementos; almacena los resultados (desplazamiento, esfuerzo).
  • Elemento: Unidad mínima de cálculo — en modelos 3D, normalmente un tetraedro.
  • Malla: La totalidad de nodos + elementos que describe geométricamente el objeto.

Resolución vs. coste computacional

  • Alta densidad → resultados detallados, tiempo de cómputo elevado.
  • Baja densidad → rápido, pero pierde detalle en zonas de alto gradiente de esfuerzo.
  • Convergencia: Se aumenta la densidad hasta que el resultado no cambia significativamente.

Tres ingredientes de una simulación

1. Material

  • Módulo de Young (E): Rigidez del tejido ante la deformación.
  • Coeficiente de Poisson (ν): Deformación transversal al comprimir.
  • Hueso cortical: E ≈ 20 GPa, ν ≈ 0.3

2. Cargas (Loads)

  • Vectores de fuerza aplicados en nodos o superficies.
  • En paleontología: fuerzas musculares estimadas con MyoGenerator (Día 5).

3. Restricciones (Constraints)

  • Nodos fijos que simulan el punto de apoyo del hueso.
  • Ej.: cóndilos mandibulares como fulcros de la palanca ósea.

El resultado: mapas de Von Mises y Strains

  • Von Mises Stress: Criterio escalar que integra todos los componentes del esfuerzo — indica si el material está cerca de fallar mecánicamente.
  • Principal Strains: Dirección y magnitud de la deformación máxima y mínima.
  • Mapa de color: Zonas rojas = alto esfuerzo; azules = bajo esfuerzo.
  • Ley de Wolff: Los huesos tienden a ser más gruesos en las zonas de mayor esfuerzo histórico — el FEA permite verificarlo en fósiles.
  • Fuerza de mordida: En un análisis estático, la fuerza de mordida puede obtenerse como la fuerza de reacción en el punto de constricción si las fuerzas musculares son aplicadas.

Fuerza de mordida en Tyrannosaurus rex a distintas edades

Preguntas que el FEA puede responder en paleontología

  • ¿Qué tan fuerte mordía un animal extinto?
  • ¿Qué zonas del cráneo soportaban mayor estrés durante la captura de presas?
  • ¿Era una morfología resistente a torsión o a compresión pura?
  • ¿Cómo cambia la biomecánica entre un individuo juvenil y uno adulto?
  • ¿Qué dieta es mecánicamente compatible con la morfología del cráneo?

Herramientas de Simulación

El Software: Fossils

Un solver abierto y eficiente

  • Fossils (Chatar et al. 2023): Un protocolo de código abierto diseñado específicamente para simular cargas biomecánicas en huesos.
  • Ventajas:
    • Optimizado para mallas volumétricas de alta resolución (millones de tetraedros).
    • Manejo directo de fuerzas musculares distribuidas.
    • Basado en el motor Metafor, pero con un solver lineal estático rápido.
  • El reto: La preparación de datos para FEA suele ser la parte más lenta del proceso.

Fossils

Introducción a BFEX

El puente entre Blender y Fossils

  • BFEX (Blender Finite Element eXporter) (Díaz de León-Muñoz, Boman, y Ferreira 2025):
  • Un add-on para Blender diseñado para simplificar la creación de modelos de FEA.
  • Objetivo: Automatizar la limpieza de mallas y la selección de áreas para condiciones de contorno.
  • ¿Por qué Blender? Por su potente motor de manipulación de mallas y su accesibilidad.

Ejemplo de modelo construido con BFEX

Flujo de Trabajo con BFEX (I)

1. Preparación y Control de Calidad

  • La malla debe ser Manifold (estanca) y sin errores geométricos.
  • Uso del 3D Print Toolbox:
    • Check All y asegurar que Non-Manifold Edges y Intersecting Faces sean 0.
  • Escala: Aplicar transformaciones (Ctrl+A) para asegurar unidades correctas (usualmente mm).

Flujo de Trabajo con BFEX (II)

2. Configuración del Proyecto

  • Definir un directorio de trabajo y nombre del proyecto.
  • El add-on crea automáticamente la estructura de archivos necesaria.
  • Colección de BFEX: Organiza los elementos de la simulación de forma independiente al modelo original.

Flujo de Trabajo con BFEX (III)

3. Definición de Superficies y Cargas

  • Zonas de Interés (ROI): Selección directa en la malla de las áreas de inserción muscular y restricciones.
  • Boundary Conditions (BCs):
    • Supports: Nodos fijados (ej. cóndilos mandibulares).
    • Loads: Aplicación de fuerzas vectoriales (fuerza muscular calculada en el Día 5).
  • Sub-meshes: BFEX extrae estas selecciones como archivos separados para que Fossils las reconozca.

Flujo de Trabajo con BFEX (IV)

4. Propiedades y Parámetros

  • Materiales: Definición de propiedades elásticas.
    • Módulo de Young (E): Rigidez del hueso (ej. 20 GPa para hueso cortical).
    • Coeficiente de Poisson (\(\nu\)): Deformación transversal.
  • Script de Parámetros: BFEX genera un archivo .py que contiene toda la lógica de la simulación lista para Fossils.

Ejecución y Análisis

Del “Exportar” al “Resolver”

  1. Exportación: BFEX genera todos los archivos (mallas y scripts).
  2. Batch Mode: Se puede ejecutar Fossils directamente desde la interfaz de BFEX o vía terminal.
  3. Visualización:
    • Los resultados se pueden importar de vuelta a Blender o abrir en ParaView para análisis detallado de mapas de color (Von Mises Stress, Principal Strains).

BFEX

Práctica de hoy: Mi primer FEA

Objetivos:

  1. Instalar y configurar BFEX.
  2. Limpiar un modelo de mandíbula.
  3. Crear un proyecto de FEA.
  4. Seleccionar los apoyos (cóndilo) y las fuerzas (aductores).
  5. Configurar materiales y exportar el script para Fossils.

Enlaces y Recursos

Bibliografía

Bright, Jen A. 2014. «A Review of Paleontological Finite Element Models and Their Validity». Journal of Paleontology 88 (4): 760-69. https://doi.org/10.1666/13-090.
Chatar, Narimane, Romain Boman, Valentin Fallon Gaudichon, Jamie A. MacLaren, y Valentin Fischer. 2023. «“Fossils”: A New, Fast and Open-Source Protocol to Simulate Muscle-Driven Biomechanical Loading of Bone». Methods in Ecology and Evolution 14 (3): 848-59. https://doi.org/10.1111/2041-210X.14051.
Díaz de León-Muñoz, E. Miguel, Romain Boman, y Gabriel S. Ferreira. 2025. «BFEX: A Toolbox for Finite Element Analysis With Fossils and Blender». Ecology and Evolution 15 (3): e71093. https://doi.org/10.1002/ece3.71093.