INGEOEXPERT

Curso de Introducción a la
Modelización con FLAC3D

Profesor: Angel Fco. Silvestre Ordaz

Módulo 1:
Introducción a FLAC3D
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

ÍNDICE

1 Presentación e introducción al curso ..................................... 1

2 Objetivos de este modulo ........................................................ 3
3 El proceso de modelización numérica en geotencia .............. 3
3.1 Generalidades ....................................................................................... 3
3.2 Las etapas del proceso de modelización numérica ............................... 5
3.2.1 Recopilación de la información disponible ............................................................ 5
3.2.2 Crear un modelo conceptual del problema real en cuestión ................................. 6
3.2.3 Definir los objetivos que se persiguen con el modelo ........................................... 6
3.2.4 Desarrollar modelos simplificados (si se considera necesario) ............................ 7
3.2.5 Combinar toda la información disponible .............................................................. 7
3.2.6 Crear un modelo matemático conceptual detallado .............................................. 8
3.2.7 Definir expectativas ............................................................................................... 9
3.2.8 Ejecutar el modelo ................................................................................................. 9
3.2.9 Interpretar, analizar y cuestionar los resultados ................................................. 10
3.2.10 Presentar los resultados ...................................................................................... 10

4 El programa de cálculo numérico FLAC3D ........................... 10

4.1 Generalidades ..................................................................................... 10
4.2 Características de FLAC3D ................................................................. 12
4.3 Conceptos fundamentales en FLAC3D ............................................... 14
5 Fundamentos del funcionamiento de FLAC3D ..................... 23
5.1 Formulación en diferencias finitas en FLAC3D ................................... 23
5.2 Comparación con el método de los elementos finitos. Ventajas y
desventajas ................................................................................................... 26
6 Evaluación ............................................................................... 28
7 Referencias bibliográficas ..................................................... 31
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

LISTADO DE FIGURAS
Figura 3.1 – Etapas del proceso de modelización numérica ........................................................ 5
Figura 3.2 – Modelo conceptual del problema real en cuestión ................................................... 6
Figura 3.3 – Modelo matemático conceptual detallado ................................................................ 8
Figura 4.1 – Ejemplo de malla, “zone” y “gridpoints” en FLAC3D - I .......................................... 13
Figura 4.2 – Ejemplo de malla, “zone” y “gridpoints” en FLAC3D - II ......................................... 13
Figura 4.3 – Ejemplo de malla, “zone” y “gridpoints” en FLAC3D - III ........................................ 14
Figura 4.4 – Sistema de ejes coordenados en FLAC3D ............................................................. 14
Figura 4.5 – Contornos de una malla de FLAC3D ...................................................................... 16
Figura 4.6 – Ejemplo de malla y submalla en FLAC3D .............................................................. 17
Figura 4.7 – Ejemplo de rangos en FLAC3D .............................................................................. 18
Figura 4.8 – Ejemplo de elementos de rangos en FLAC3D ....................................................... 18
Figura 4.9 – Ejemplo de group y slot en FLAC3D - I .................................................................. 19
Figura 4.10 – Ejemplo de group y slot en FLAC3D - II ............................................................... 19
Figura 4.11 – Ejemplos de primitivas en FLAC3D - II ................................................................. 22
Figura 5.1 – División de las zonas en elementos tetraédricos .................................................... 24
Figura 5.2 – Esquema de representación del ciclo de cálculo explícito básico en FLAC3D ...... 25

LISTADO DE TABLAS
Tabla 1.1 – Módulos del Curso de introducción a la modelización con FLAC3D ......................... 2
Tabla 3.1 – Interpretar, analizar y cuestionar resultados ............................................................ 10
Tabla 4.1 – Aplicaciones del software FLAC3D .......................................................................... 11
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

1 PRESENTACIÓN E INTRODUCCIÓN AL CURSO

Bienvenidos al Curso de Introducción a la Modelización con FLAC3D. A lo largo de los ocho
módulos en los que consiste este curso, el alumno se introducirá en el manejo del software
para la modelización de problemas propios del ámbito de la ingeniería geotécnica. Tal y como
figura en la presentación del curso en la página web de Ingeoexpert, es importante destacar
que el objetivo del curso no es formar a expertos en el uso de FLAC3D, sino iniciar a los
alumnos en el empleo de FLAC3D, proporcionándoles las herramientas y conocimientos
básicos y elementales necesarios para que, posteriormente, puedan adentrarse por su cuenta y
con éxito en los aspectos más sofisticados y complejos de FLAC3D.
Además de lo anterior, y tal y como descubrirá el alumno a medida que avance en el curso,
FLAC3D es un software tremendamente extenso, que ofrece infinidad de posibilidades y
opciones distintas (todas igualmente válidas, con sus ventajas y desventajas) para afrontar un
mismo problema geotécnico. Todas estas opciones son imposibles de abordar en su totalidad
en un curso de estas características, fundamentalmente por razones y requerimientos de
tiempo y espacio. Así pues, ante estos condicionantes, se decidió enfocar el curso hacia la
iniciación a FLAC3D, dejando la posibilidad a los alumnos de adentrarse (una vez finalizado el
curso o durante la realización del mismo) en los aspectos más complicados y avanzados por su
cuenta, ayudados por los conocimientos y herramientas adquiridos.
En relación con las múltiples opciones que ofrece FLAC3D para abordar un mismo problema
geotécnico, conviene destacar dos aspectos principales:
1) Dado que, como comprobará el alumno a lo largo del curso, FLAC3D es un software cuya
principal característica es la programación mediante comandos y a través del lenguaje de
programación interno FISH, existen numerosas (por no decir infinitas) formas de escribir
un código para la resolución y análisis de un mismo problema geotécnico. En este sentido,
el alumno debe tener en cuenta que, en este curso, el profesor no le está enseñando la
única forma de llevar a cabo una modelización con FLAC3D, sino tan solo le está
mostrando una de las posibilidades para llevarla a cabo y las herramientas y
conocimientos de los que debe disponer para ello. Así pues, es importante que, en este
sentido, el alumno tenga en cuenta que las estructuras de programación que se facilitan
con el curso responden, en cierta medida, a las preferencias, punto de vista y forma
particular de emplear el software del profesor, que éste ha desarrollado a lo largo de su
trayectoria profesional. Esto implica que el alumno, una vez adquiridos los conocimientos
necesarios, tiene total libertad para escribir los códigos y estructurar sus proyectos en
FLAC3D de la forma que le resulte más clara y conveniente, sin que ello implique que sea
una forma incorrecta o menos válida por el mero hecho de ser distinta a la forma en que lo
hace el profesor. De hecho, el profesor anima a los alumnos a que traten de escribir sus
códigos y organizar sus proyectos de la forma que a ellos les resulte más clara, intuitiva y
útil.
2) Tal y como se ha mencionado en el punto anterior, una de las principales características y
ventajas de FLAC3D, y posiblemente su pilar central y fundamental, es la programación
mediante comandos y el lenguaje de programación interno FISH. Este aspecto
proporciona al usuario gran versatilidad y flexibilidad, tal y como se verá a lo largo del
curso. Así pues, este curso se ha enfocado exclusivamente en mostrar a los alumnos
cómo llevar a cabo simulaciones numéricas de problemas geotécnicos a través de los
comandos de FLAC3D y el lenguaje FISH. En las últimas versiones de FLAC3D (versiones
5, 6 y 7), Itasca Consulting Group ha introducido nuevas posibilidades en el software como
la herramienta Building Blocks (para la generación de geometrías complejas de manera
visual e interactiva), la herramienta de Extrusion (generación de modelos 3D por medio de
una sección 2D) y la importación de geometrías desde otros softwares, entre otras
posibilidades. Sin embargo, ningún modelo numérico puede llevarse a cabo
completamente (definición de la geometría, asignación de propiedades, establecimiento de
condiciones iniciales y de contorno, fases de construcción etc.) a través de estas
herramientas sin valerse en absoluto de los comandos de FLAC3D y del lenguaje interno
FISH. Por tanto, puesto que la programación mediante comandos y el lenguaje de

1
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

programación FISH constituye la piedra angular de FLAC3D y permite llevar a cabo

modelizaciones completas (desde la generación de la geometría hasta la resolución del
modelo), el curso se ha enfocado exclusivamente y en su totalidad en enseñar al
alumno cómo llevar a cabo simulaciones numéricas por medio de los comandos y el
lenguaje FISH. Las herramientas como Building Blocks, Extrusion, etc. quedan fuera del
contenido fundamental del curso y se abordarán, según el caso, mediante
videoconferencias específicas. Ello responde también a los requerimientos y
condicionantes de espacio y de tiempo de un curso de estas características. Puesto que
es imposible abarcar todas las posibilidades de FLAC3D en único curso de iniciación como
éste, se ha considerado que resultaba más relevante y de mayor utilidad para el alumno
aprender y conocer la programación mediante comandos y el lenguaje FISH, que le
permitirá llevar a cabo modelizaciones completas en FLAC3D. Una vez finalizado el curso,
el alumno podrá investigar por su cuenta el uso de las herramientas mencionadas
anteriormente, sabiendo que necesitará los comandos y el lenguaje de programación FISH
para llevar a cabo el cálculo y que el curso le habrá proporcionado los conocimientos para
ello.
El Curso de Introducción a la Modelización con FLAC3D se estructura en los ocho módulos
siguientes:
Tabla 1.1 – Módulos del Curso de introducción a la modelización con FLAC3D

Módulo Título Descripción

Se presentarán los fundamentos científicos y teóricos del método de las
diferencias finitas en el que se basa FLAC3D y sus diferencias con el
1 Introducción a FLAC3D
método de los elementos finitos. Así mismo, se expondrá el esquema y
procedimiento de resolución que sigue el software.
La interfaz de usuario Se explicará el funcionamiento, la organización y la estructura de la interfaz
2 de FLAC3D y la de usuario de FLAC3D. Así mismo, se detallará la estructura de proyecto
estructura de proyectos que suele seguirse al modelizar con FLAC3D.
Los comandos de En este módulo, se tratarán tanto los comandos propios de FLAC3D como el
3 FLAC3D y el lenguaje lenguaje de programación interno FISH y la forma de combinar ambos para
FISH llevar a cabo cálculos.
Generación de Se estudiará cómo pueden emplearse los comandos de FLAC3D y el
4 geometrías con lenguaje FISH para generar geometrías en FLAC3D. Se detallará también la
FLAC3D creación de regiones y grupos dentro de una geometría dada.
Se describirá cómo asignar los modelos constitutivos y los valores de los
Estados iniciales,
distintos parámetros a las diferentes zonas del modelo. Así mismo, se
condiciones de
5 explicará la generación de tensiones y la asignación de desplazamientos y
contorno y modelos
velocidades a los nodos para la creación de las condiciones de contorno del
constitutivos
modelo.
En este módulo se explicará cómo considerar e implementar diferentes fases
6 Fases constructivas constructivas para llevar a cabo la simulación numérica de un problema
geotécnico real.
El alumno aprenderá el funcionamiento de los diferentes elementos
estructurales disponibles en FLAC3D, así como de los elementos interfase.
7 Elementos estructurales Se explicará cómo introducir estos elementos en la malla del modelo, su
activación y desactivación durante las fases constructivas y la asignación de
los valores de sus respectivos parámetros.
El agua es un elemento fundamental en los cálculos geotécnicos. Dada su
importancia y las distintas posibilidades que ofrece FLAC3D para su
8 El agua en FLAC3D consideración en el análisis, se destina un módulo completo del curso a este
tema. El alumno aprenderá a tener en cuenta el agua en las simulaciones
numéricas con FLAC3D.

El primer módulo del curso es fundamentalmente teórico, pero esencial para comprender bien
el funcionamiento de FLAC3D y ser conscientes en todo momento de lo que está haciendo con
el programa. De hecho, uno de los principales problemas ligados a la modelización numérica
(ya sea en ingeniería geotécnica, cálculo estructural u otro campo de la ingeniería) radica en
que se desconocen (o no se conocen lo suficiente) los principios y fundamentos teóricos y
matemáticos en los que se basan los programas de cálculo. Así pues, antes de lanzarse a
modelizar con FLAC3D, es preciso conocer su esquema y lógica de funcionamiento. El resto de

2
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

los módulos (del 2 al 8) son de carácter más práctico, y a través de ellos el alumno irá
aprendiendo gradualmente cómo llevar a cabo una simulación numérica con FLAC3D.
Para llevar a cabo este curso, el alumno dispone de la licencia web de FLAC3D que
Ingeoexpert incluye como parte fundamental del contenido del curso. Esta licencia tiene un
periodo de validez equivalente a la duración del curso.
Dicho todo lo anterior, empezamos el Curso de Introducción a la Modelización con FLAC3D.
Estoy a vuestra disposición para resolver todas las dudas y contestar a todas las preguntas que
tengáis.
Espero que el curso sea de vuestro agrado y que quedéis plenamente satisfechos. Sed de
nuevo muy bienvenidos y gracias por vuestra confianza. ¡Comienza el curso!

2 OBJETIVOS DE ESTE MODULO

El primer módulo del Curso de Introducción a la Modelización con FLAC3D persigue los
siguientes objetivos:
1) Dar a conocer (o recordar) al alumno el proceso de modelización numérica en ingeniería
geotecnia, con sus sucesivas etapas. Este procedimiento es válido para los análisis
numéricos ejecutados con FLAC3D y su implementación es necesaria a la hora de utilizar
el programa.
2) Presentar FLAC3D al alumno, resaltando sus principios básicos y fundamentales.
3) Explicar de forma conceptual, sin entrar en desarrollos matemáticos profundos que no son
objeto de este curso, el método de las diferencias finitas en el que se basa FLAC3D y
comprender el funcionamiento del programa. Entender el funcionamiento del software de
análisis numérico que se está utilizando para el estudio de un problema determinado es
imprescindible para no solo usar correctamente el software, sino también para comprender
los resultados obtenidos, identificar errores, introducir correctamente los parámetros, etc.
En muchas ocasiones, los errores que se cometen no se deben a una falta de
conocimientos por parte del ingeniero de la materia en cuestión, sino por su
desconocimiento del funcionamiento de la herramienta que está utilizando. Por tanto,
dedicar una parte del módulo a estudiar cómo funciona FLAC3D es fundamental para
asegurar y garantizar un uso eficiente y satisfactorio del programa.

3 EL PROCESO DE MODELIZACIÓN NUMÉRICA EN

GEOTENCIA
3.1 Generalidades

La modelización numérica, ya sea en geotecnia o en cualquier otra disciplina (economía,

cálculo estructural, hidráulica, transporte, mecánica de fluidos, finanzas, etc.), consiste
básicamente en, dado un problema real, “transformarlo” en un modelo matemático más o
menos simplificado (en ocasiones no es posible, ni deseable, tener en cuenta todas las
variables y condiciones del problema por motivos de tiempo, dinero, potencia de cálculo,
complejidad del modelo, funcionalidad, operatividad, eficiencia, etc.) por medio de diferentes
herramientas para poder analizarlo y comprender su funcionamiento con unos objetivos
determinados (diseño, concepción, dimensionamiento, comprobación, etc.).
Puede decirse que la modelización numérica es la forma más precisa y realista de abordar el
análisis de los problemas en ingeniería, si bien como contrapartida suele implicar un mayor
coste y tiempo (de este último cada vez menos gracias al avance de la capacidad y potencia de
cálculo de los ordenadores). Así pues, a la hora de estudiar un problema concreto, el ingeniero
debe preguntarse antes si la modelización numérica es el mejor camino o no para el análisis
del problema en cuestión, frente a otros métodos como los cálculos a mano o las soluciones
analíticas o empíricas. Estos dos métodos son menos precisos que el análisis numérico, pero
consumen menos recursos de tiempo y dinero, y, dependiendo del problema del que se trata,

3
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

pueden proporcionar soluciones casi tan exactas como los modelos numéricos (o
razonablemente exactas para el grado precisión que se persigue). Por tanto, es fundamental
que, antes de lanzarse directamente a llevar a cabo un complejo análisis numérico, el ingeniero
se pregunte si el problema en concreto podría analizarse con suficiente precisión y seguridad
por medio de cálculos hechos a mano o soluciones analíticas y empíricas, o empleando
herramientas informáticas menos complejas como hojas de cálculo (Excel) y programas de
cálculo matemático.
En relación con lo anterior, es fundamental tener en cuenta que la solución a un problema de
ingeniería no es más verídica ni más exacta por el mero hecho de que sea proporcionada o se
haya obtenido a través de un análisis numérico. Un modelo numérico, ya sea con FLAC3D o
con cualquier otro software, únicamente proporcionará una solución fiable y realista si, en
primer lugar, el ingeniero sabe utilizar el software en cuestión y comprende su funcionamiento
interno y, en segundo lugar, si ha utilizado la INFORMACIÓN y los valores de los
PARÁMETROS adecuados para alimentar el modelo. Una simulación numérica, por muy
compleja, sofisticada y espectacular que sea, no servirá de nada (más que para proporcionar
figuras muy aparentes con colores muy alegres) si la información y los valores de los
parámetros empleados no son los adecuados y/o la metodología que ha empleado el ingeniero
no es la apropiada (por ejemplo, considerar parámetros a corto plazo en una situación a largo
plazo, comportamiento drenado en suelos cohesivos, etc.).
La diferencia fundamental entre la modelización numérica en ingeniería geotécnica, frente a
otros campos de la ingeniería civil como el cálculo estructural, radica precisamente en la
cantidad y la fiabilidad de la información de la que se dispone para llevar a cabo el estudio
del problema geotécnico por medio de un modelo numérico. Generalmente, la información a la
que se tiene acceso a la hora de estudiar un problema o actuación geotécnica es limitada, está
sometida a un cierto grado de variación (el terreno que parece uniforme puede modificar su
comportamiento y características de manera repentina a escasos metros de donde se está
ejecutando la obra) y a un cierto nivel de incertidumbre. La presencia de un cierto nivel de
incertidumbre y variabilidad deriva del hecho fundamental e indiscutible de que es IMPOSIBLE
obtener una información 100% completa y fiable en actuaciones en suelo y roca (túneles,
cimentaciones, muros, etc.) por una cuestión de extensión (y, por tanto, de presupuesto) y de
que tanto suelo y roca son materiales naturales que el hombre no ha fabricado (como es el
caso del hormigón y el acero), y por lo tanto no se tiene ningún control sobre sus parámetros y
comportamiento.
El nivel de incertidumbre y la posibilidad de que exista un elevado grado de variación de las
propiedades del terreno que no se detecte son tanto mayores cuanto mayor es la complejidad
geológica/geotécnica del lugar donde se tiene que llevar a cabo la actuación y cuanto menor
sea el presupuesto destinado a investigación geotécnica. Un menor presupuesto implica menos
ensayos y sondeos, lo cual afecta de forma significativa a la cantidad y la fiabilidad de la
información. Por otro lado, la accesibilidad al lugar de la actuación (que determinará la
posibilidad de emplear determinada maquinaria y equipos), condicionará también de manera
notable la información disponible. Así mismo, la ingeniería geotécnica se “enfrenta” de manera
especial a un elemento en cierta medida incontrolable, que puede alterar de manera
considerable las características y comportamiento del terreno y que en ocasiones puede tener
un comportamiento impredecible como es el agua.
Todo lo anterior implica que, en el caso de la modelización numérica en geotecnia, el nivel y la
calidad de la información de la que se dispone determinará los objetivos del análisis numérico.
Si la información es escasa y con elevado grado de incertidumbre y/o variabilidad, entonces el
modelo numérico servirá más bien para comprender y estudiar los mecanismos de
comportamiento y funcionamiento del problema en cuestión, más que para hacer predicciones,
diseños y dimensionamientos precisos. En este sentido, el modelo tendría la función más bien
de “laboratorio numérico” para comprobar y evaluar diferentes ideas y posibilidades. Si la
información es escasa, el ingeniero debe preguntarse sobre el grado de complejidad del
modelo. Un modelo, por complejo que sea, será poco más útil que un modelo simple mientras
no pueda construirse con información precisa y fiable. Por el contrario, si la información de la
que se dispone es abundante y fiable, el modelo tendrá más bien fines predictivos, de diseño y
de dimensionamiento precisos. En esta situación, el ingeniero puede plantearse la posibilidad

4
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

de llevar a cabo modelos más complejos y sofisticados, puesto que dispone de la información
necesaria para alimentarlo. En este caso, se habla de modelos numéricos “predictivos”.
Así pues, observamos que el uso y los objetivos que persigamos con el empleo de un software
de modelización numérica en geotecnia como FLAC3D (“laboratorio numérico” o modelo
“predictivo”) no depende del software en sí sino de las condiciones de campo y del presupuesto
destinado a investigaciones geotécnicas. Por supuesto, entre los dos extremos (información
abundante y fiable e información escasa y poco precisa) existen casos intermedios.
Corresponde por tanto al ingeniero el determinar el uso que le dará el modelo numérico.
Así pues, a partir de todo lo anterior, podemos concluir que:
1) La modelización numérica es la forma más precisa y exacta de abordar el estudio de los
problemas geotécnicos, si bien es más costosa en cuanto a tiempo y presupuesto frente a
los cálculos a mano y las soluciones analíticas y empíricas u otros métodos (hoja de
cálculo, etc.).
2) El análisis numérico en geotecnia demostrará su valía en geotecnia frente a otros métodos
siempre y cuando el modelo esté alimentado con la información más fiable y precisa.
3) Es la cantidad y la fiabilidad de la información de la que se dispone la que determinará el
uso y los objetivos del análisis numérico (“laboratorio numérico” o modelo “predictivo”) y no
el software que se utilice.

3.2 Las etapas del proceso de modelización numérica

El proceso de modelización numérica en ingeniería geotécnica y apropiado para el uso de

FLAC3D se compone de las siguientes etapas:

Figura 3.1 – Etapas del proceso de modelización numérica

3.2.1 Recopilación de la información disponible

En esta primera etapa, y antes de empezar cualquier tipo de modelo numérico, el objetivo es
recopilar toda la información disponible sobre el problema en cuestión a analizar. Está
información incluye, entre otros, los planos de construcción y geometría, registros de sondeos,
informes geotécnicos y geológicos, etc. Es en esta primera etapa, en función de la cantidad y
de la fiabilidad de la información de la que se disponga, donde se determinará el uso y los
objetivos del modelo numérico, es decir, si se empleará más bien como “laboratorio numérico”

5
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

para estudiar los mecanismos de comportamiento del problema o bien como “modelo
predictivo” para realizar predicciones, diseños y dimensionamientos precisos.

3.2.2 Crear un modelo conceptual del problema real en cuestión

Una vez recopilada la información sobre el problema a analizar, la siguiente fase consiste en
combinar la información recopilada para definir un modelo conceptual del problema. Este
modelo conceptual debe mostrar sobre todo las características geométricas y las dimensiones
del problema a analizar. A continuación, se muestra un ejemplo de un modelo conceptual para
el análisis de la excavación y sostenimiento de un túnel. Este primer modelo conceptual
permitirá establecer de manera preliminar las dimensiones y geometría del modelo, los tipos de
materiales, las condiciones de contorno y el estado inicial. Así mismo, nos permitirá determinar,
junto con la información recopilada en la etapa anterior, si es verdaderamente necesario
abordar el problema con un modelo numérico 3D o si un modelo 2D sería suficiente para
obtener resultados lo suficientemente precisos y realistas. Por último, este modelo conceptual
inicial ayuda a plantearse los primeros interrogantes sobre el problema a analizar como, por
ejemplo: ¿puede ser el modelo inestable? ¿serán los movimientos grandes o pequeños en
comparación con las dimensiones de los elementos? ¿afectará el agua al comportamiento de
los materiales? ¿podemos tener en cuenta simetrías para simplificar el análisis?

Figura 3.2 – Modelo conceptual del problema real en cuestión

3.2.3 Definir los objetivos que se persiguen con el modelo

Tras definir el modelo conceptual, el siguiente paso es establecer los objetivos que se
persiguen con la realización del modelo numérico. Esta etapa está muy estrechamente
relacionada con la primera, pues dependiendo de la cantidad, fiabilidad y precisión de la
información disponible podremos determinar unos objetivos u otros. La definición correcta de
los objetivos es fundamental para posteriormente definir un modelo conceptual matemático
detallado. Por tanto, esta etapa es clave en el proceso de modelización numérica. Si no se
definen los objetivos correctamente, entonces la modelización numérica pierde su objetivo y se

6
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

está perdiendo el tiempo innecesariamente. Siguiendo con el ejemplo del túnel anterior, una
definición de los objetivos podría ser:
1) Análisis de las convergencias en la pared del túnel.
2) Estudio de la extrusión del frente.
3) Esfuerzos en los elementos de sostenimiento y revestimiento.
4) Estabilidad del frente.
Definiendo los objetivos que perseguimos con nuestro modelo numérico con antelación, a
modo de “lista de cosas que hacer” o “check-list”, nos aseguramos de que nuestro modelo
conceptual matemático detallado incluirá todos los elementos necesarios para alcanzar estos
objetivos y que ninguno quedará olvidado. La definición de los objetivos está también
íntimamente relacionada con el nivel de detalle y complejidad necesaria en el modelo
numérico. En ocasiones, puede resultar muy tentador incluir elementos o geometrías complejas
en el modelo por el mero hecho de que existen en la realidad. Sin embargo, estas
complejidades deben omitirse si su influencia en la respuesta y comportamiento del modelo va
a ser muy limitada o si resultan irrelevantes para el cumplimiento correcto de los objetivos que
se han determinado.

3.2.4 Desarrollar modelos simplificados (si se considera necesario)

En ocasiones, cuando se trata de modelos numéricos 3D de gran complejidad y sofisticación,

puede resultar muy conveniente llevar a cabo previamente modelos sencillos que nos ayuden a
comprender mejor el comportamiento del modelo completo complejo. Así mismo, estos
modelos sencillos nos pueden servir para identificar y resolver problemas antes de tener que
hacerlo sobre el modelo detallado, lo cual siempre será más difícil y llevará más tiempo (que se
sumará al tiempo que se ha perdido construyendo el modelo completo sin haber identificado y
resuelto el problema). En este sentido, los modelos sencillos pueden obligarnos a replantear
los objetivos y el modelo conceptual definidos en las etapas anteriores.
Por ejemplo, llevar a cabo modelos sencillos previos nos puede permitir determinar si los
parámetros del terreno escogidos son verdaderamente representativos de su comportamiento,
o bien si estos parámetros están infravalorados o sobrevalorados. Generalmente, los modelos
sencillos nos permitirán determinar un rango de parámetros razonable para los terrenos con el
fin de tener en cuenta la incertidumbre y la variabilidad. Así mismo, por medio de modelos
sencillos, podemos calibrar los parámetros que definen el comportamiento de los elementos
estructurales (por ejemplo, simulando el ensayo de arrancamiento de un anclaje o de carga de
hundimiento de un pilote) e identificar elementos complejos del mallado que pueden generar
inestabilidades. En todos los casos, los modelos numéricos sencillos nos permitirán resolver y
comprender mejor (y más rápido) estos problemas que el modelo complejo detallado.

3.2.5 Combinar toda la información disponible

Realizados los modelos sencillos (si estos se consideran necesarios), establecidos los
objetivos, definido el modelo conceptual y recopilada la información, en esta etapa el objetivo
es combinar toda la información para introducirla posteriormente en el modelo matemático
conceptual detallado, que constituye la base y soporte fundamental para la construcción del
modelo numérico. Entre la información relevante que es preciso tener en cuenta, destaca:
1) Detalles de la geometría y las dimensiones.
2) Comportamiento de los materiales (terreno, estructuras, discontinuidades, etc.)
3) Condiciones iniciales.
4) Condiciones de contorno.
5) Cargas externas.
6) Fases de construcción.

7
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

7) Resultados de los modelos sencillos.

En relación con el comportamiento de los materiales, en especial del terreno, es en esta etapa
en la que se debe analizar y estudiar el rango de parámetros geomecánicos que determinan el
comportamiento del terreno con el fin de tener en cuenta la incertidumbre y la variabilidad en
nuestros cálculos. En este sentido, si el análisis requiere asignar parámetros diferentes a una
misma zona del modelo para cuantificar la influencia que tiene este cambio de parámetros, el
modelo deberá construirse teniendo en cuenta esta posibilidad, permitiendo identificar y
“seleccionar” esta zona de forma sencilla para que la modificación de los parámetros pueda
hacerse de manera segura y rápida.

3.2.6 Crear un modelo matemático conceptual detallado

En esta fase, toda la información combinada en la etapa anterior debe introducirse en el

modelo conceptual para crear el modelo matemático conceptual detallado que servirá de base
para nuestro modelo numérico y permitirá alcanzar los objetivos establecidos. Este modelo
matemático conceptual detallado no tiene por qué ser muy elaborado. En muchos casos, basta
con un dibujo esquemático hecho a mano, tal y como muestra la figura que aparece a
continuación, siguiendo con el ejemplo del análisis del túnel.

Figura 3.3 – Modelo matemático conceptual detallado

Es muy importante no omitir esta etapa del proceso de modelización numérica. Crear un
modelo matemático conceptual detallado, aunque sea un dibujo esquemático hecho a mano,
permite evitar errores y, en caso de cometerlos, localizarlos antes (evitando así pérdidas de
tiempo considerables). Así mismo, nos permite recordar de manera fácil y rápida la estructura
del modelo en caso de que transcurra mucho tiempo desde que se utilizó. Por otro lado, en
caso de que una tercera persona deba hacer uso del modelo, disponer de un modelo
matemático conceptual del mismo le ayudará a comprenderlo y a emplearlo de manera sencilla
y efectiva, sin tener que destinar varias horas a “desentrañar” sus “misterios”.
Para la realización del modelo matemático conceptual detallado, es necesario tener en cuenta
las siguientes premias:
1) ¿Cuál es el tiempo necesario para el cálculo? Éste será función de la capacidad de cálculo
del ordenador/es empleado/s y del tamaño del modelo. A mayor tamaño del modelo,
mayor será el tiempo necesario para el cálculo. Por tanto, el tamaño del modelo debe
ajustarse tanto por el tiempo requerido para ejecutar el análisis como por la precisión de
resultados necesaria. En muchas ocasiones, un modelo de 100 Gb puede aportar
resultados respecto a un modelo de 50 Gb que solo se diferencian en el tercer decimal,
necesitando posiblemente el doble de tiempo para completarse. Así mismo, contra más
fases contemple el modelo, mayor será el tiempo de cálculo necesario para completar el
análisis, por lo que es necesario asegurar que las fases consideradas cumplen
estrictamente con los objetivos definidos. Toda fase superflua que no contribuya a estos
objetivos debe eliminarse del modelo. Por otro lado, siempre que sea posible, resulta
conveniente disponer de varios ordenadores para poder llevar a cabo varios cálculos a la

8
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

vez y así ahorrar tiempo. Esto es particularmente útil cuando un mismo modelo debe
analizarse considerando distintos valores de los parámetros de los terrenos, elementos
estructurales, condiciones iniciales, casos de carga, etc.
2) Es importante que el estado del modelo pueda guardarse a lo largo de las fases
intermedias del mismo. De este modo, en caso de cometerse un error en alguna de las
fases, no es necesario empezar el cálculo desde el principio, sino que se puede reiniciar
desde la fase anterior al error. El guardar el modelo en fases intermedias del cálculo evita
también tener que empezar el análisis desde el principio en caso de que haya un problema
de desconexión o apagado repentino del ordenador. De producirse esta situación, el
modelo puede reiniciarse desde la fase intermedia anterior a la desconexión/apagado,
permitiendo ganar tiempo al no tener que empezar el cálculo desde el principio. Así
mismo, guardar el modelo a lo largo de las fases intermedias permite ir controlando y
analizando los resultados, detectando posibles inestabilidades, velocidades excesivas,
desplazamientos demasiado grandes o pequeños, etc. A la hora de guardar el modelo
durante las fases intermedias, es necesario comprobar que se dispone del espacio en el
disco suficiente. Dada esta limitación, es fundamental decidir cuáles son las fases
intermedias relevantes que es importante guardar, desechando aquellas que no sean
importantes o sean menos importantes.
3) A la hora de realizar una modelización numérica, y especialmente con FLAC3D, es
importante definir puntos de control en el modelo que nos permitan comprobar que se está
ejecutando correctamente. Estos puntos de control pueden ser desplazamientos y/o
velocidades en los nodos, tensiones en algunas zonas, etc. En el caso de FLAC3D, tal y
como se verá más adelante en el curso, resulta muy conveniente monitorizar y controlar
las “unbalanced forces” (fuerzas desequilibradas) y el ratio de convergencia. El
seguimiento de estas dos magnitudes permite confirmar la estabilidad del modelo y su
convergencia hacia un estado de equilibrio.

3.2.7 Definir expectativas

Definido el modelo matemático conceptual detallado, hay que preguntarse acerca de las
expectativas que tenemos sobre el mismo. Este hecho responde a la siguiente cuestión: si no
tenemos expectativas, ¿cómo sabemos que nuestro modelo es correcto? En modelización
numérica, nunca hay que dejarse convencer y distraer por las preciosas figuras llenas de
colores y escalas que nos permite obtener el software. Lo esencial es saber determinar si esas
figuras son correctas o no y si los resultados son realistas y tienen sentido desde el punto de
vista ingenieril. Para ello, es fundamental definir unas expectativas sobre los resultados del
modelo. Estas expectativas pueden materializarse de diferentes maneras, desde un escrito
donde se describe lo que se piensa que será el resultado, hasta la estimación de los resultados
por medio de soluciones analíticas o empíricas.
En este sentido, en el caso de la modelización del ensayo de carga de un pilote, sería
interesante comprobar que, en el hundimiento, toda la resistencia por fuste debe estar agotada.
Comprobar el valor total de la carga de resistencia por fuste que proporciona el modelo con el
valor esperado (calculado a mano o con hoja de cálculo) sería un ejemplo de expectativa.
También sería un ejemplo de expectativa comparar la convergencia de la pared de un túnel
que proporciona el modelo con los resultados de soluciones analíticas, para comprobar que
está dentro del orden de magnitud esperado. De esta forma, definiendo las expectativas con
antelación, podremos determinar posteriormente si los resultados del modelo son coherentes o
no.

3.2.8 Ejecutar el modelo

Durante la ejecución del modelo, es importante controlar los resultados de los puntos de control
establecidos, así como el ratio de convergencia y las “unbalanced forces”. Ello permite
asegurar que el modelo se está ejecutando correctamente. Por otro lado, cada cierto tiempo, se
recomienda comprobar los estados intermedios para asegurar que todo está en orden. En caso

9
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

de identificarse un problema, el cálculo se detiene, se resuelve el problema y el cálculo se

reinicia desde la fase intermedia más avanzada no afectada por el problema solucionado.

3.2.9 Interpretar, analizar y cuestionar los resultados

Concluido el cálculo, y tal y como se ha dicho anteriormente, no hay que dejarse llevar por las
impresionantes figuras que produce el software y, simplemente por su aspecto y por tratarse de
un potente software de cálculo numérico, considerar que los resultados son correctos. Hay que
tener en cuenta que el modelo numérico siempre nos dará una solución al problema que
hayamos planteado atendiendo a la información, modelos constitutivos, propiedades,
características y fases constructivas que le hayamos proporcionado y las hipótesis que
hayamos considerado. Por tanto, esta solución no tiene por qué ser necesariamente correcta
si se ha cometido un error en la información proporcionada o en las hipótesis consideradas.
Finalizado el cálculo, es importante que el ingeniero interprete, analice y cuestione los
resultados del modelo atendiendo a la secuencia que se proporciona en la siguiente tabla.
Tabla 3.1 – Interpretar, analizar y cuestionar resultados

Interpretar, analizar y cuestionar los resultados

• Revisar el modelo en busca de errores.

• Comprobar que las expectativas son correctas.
• Volver al modelo matemático conceptual detallado y hacer los cambios
necesarios.
No
¡IMPORTANTE! Hacer un solo cambio a la vez en el modelo cuando se esté buscando
¿Los resultados
errores o la solución a los mismos. Si hacemos diversos cambios a la vez (por ejemplo,
del modelo son
modificar cohesión y módulo de Young a la vez), no sabremos determinar qué parte del
semejantes a las
cambio en los resultados es debida a la modificación de cada parámetro.
expectativas?
Volver a definir objetivos y realizar los cambios
¿Responde el modelo a los No oportunos al modelo conceptual matemático
Sí detallado. Repetir el proceso.
objetivos planteados?
Si Proceder a la presentación de resultados.

3.2.10 Presentar los resultados

En modelización numérica, la mejor manera de presentar los resultados es por medio de

figuras, tablas y gráficas que permiten apreciar, comprender y analizar la evolución y desarrollo
de los movimientos y las tensiones y/o esfuerzos. A la hora de presentar los resultados, el
ingeniero debe tener en cuenta que, como autor del modelo y de los cálculos, su nivel de
comprensión de los mismos es mucho más elevado que el de las personas a quienes debe
presentar los resultados. Por ello, el ingeniero debe esforzarse en presentar los resultados de
la manera más clara y sencilla posible, nunca dando nada por supuesto, dado que lo que
puede resultar evidente para él puede no resultar tan evidente para los demás.
Aunque pueda parecer laborioso, seguir este proceso a la hora de llevar a cabo modelizaciones
numéricas puede ahorrarnos mucho tiempo y evitar que cometamos errores. Así pues, se
recomienda a los alumnos, especialmente los menos experimentados en modelización
numérica, que empleen esta metodología en su día a día cuando se estén enfrentando a una
modelización numérica.

4 EL PROGRAMA DE CÁLCULO NUMÉRICO FLAC3D

4.1 Generalidades

FLAC3D (Fast Lagrangian Analysis of Continua) es un programa de cálculo numérico

tridimensional generalista especialmente concebido para analizar problemas geotécnicos,
geológicos y mineros. El programa se basa en el método explícito de las diferencias finitas,
también conocido como el método de los volúmenes finitos. El software permite tener en
cuenta cualquier geometría a través de elementos de forma arbitraria, con condiciones de
contorno más o menos generales (esto es, bloqueo de los desplazamientos en la dirección

10
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

normal a los planos de los contornos laterales y bloqueo de todos los desplazamientos en la
base del modelo).
El software permite analizar problemas estáticos y dinámicos significativamente no linéales en
pequeñas o grandes deformaciones, ya sea en suelo o en roca y teniendo en cuenta o no
(según el caso) el acoplamiento hidromecánico y la interacción suelo-estructura. Incorpora una
gran variedad de modelos constitutivos elásticos, plásticos y elasto-viscoplásticos, así como
elementos lineales y superficiales especialmente concebidos para modelizar vigas, pilotes,
anclajes, bulones y placas. Así mismo, permite considerar planos de deslizamiento y la
interacción suelo-estructura por medio de elementos de tipo “interfase”. El módulo principal del
software puede complementarse con tres módulos opcionales (dinámico, fluencia (creep) y
térmico).
Los diferentes materiales del terreno (suelos, rocas, etc.) se representan por medio de
elementos poliédricos en una malla tridimensional que el usuario ajusta para considerar la
forma y la geometría del problema en cuestión a analizar. Cada elemento de la malla se
comporta según una determinada relación tensión/deformación (ley constitutiva), de carácter
lineal o no lineal, en respuesta a las fuerzas aplicadas y a las condiciones del contorno
establecidas. Los materiales pueden plastificar y la malla deformarse (en el caso del modo
“large-strain” o “grandes desplazamientos”, que no es objeto de este curso) y moverse.
Tal y como ser verá más adelante, FLAC3D resuelve de forma total las ecuaciones dinámicas
del movimiento, aunque se trate de problemas cuasi estáticos. La principal ventaja de esta
forma de proceder radica en la posibilidad que ofrece FLAC3D al usuario de modelizar de
forma muy precisa problemas de inestabilidad física, como es el caso del colapso o de la
fluencia (deformación) plástica excesiva. Ello también es, como se detallará a continuación,
una de las principales ventajas del software frente a los programas de cálculo numérico
basados en el método de los elementos finitos, en los que es necesario una convergencia total
de la solución.
La tabla que se presenta a continuación muestra diversas aplicaciones que puede tener el
programa FLAC3D, clasificadas por campos o dominios:
Tabla 4.1 – Aplicaciones del software FLAC3D

Campo /
Tipo de Proyecto Problema a resolver con FLAC3D
Dominio
Diseño de sostenimiento y revestimiento.
Análisis del método constructivo.
Factor de seguridad del frente.
Túneles
Análisis de convergencias.
Análisis dinámico.
Redes de flujo hacia los túneles.
Análisis de la capacidad portante.
Estudio de asientos.
Cálculo dinámico.
Cimentaciones
Ingeniería Modelización de ensayos de carga de pilotes.
Civil Estudio de las deformaciones y movimientos inducidos por obras
subterráneas.
Verificación de los empujes y comprobación de esfuerzos.
Análisis de los movimientos del terreno y de la estructura.
Estructuras de
Cálculo del factor de seguridad.
contención de tierra
Análisis de estabilidad.
Redes de filtración y métodos de rebajamiento del freático.
Cálculo de tensiones y deformaciones.
Presas Análisis dinámico.
Redes de filtración y flujo.

11
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Campo /
Tipo de Proyecto Problema a resolver con FLAC3D
Dominio
Análisis de la estabilidad. Cálculo del factor de seguridad.
Excavaciones /
Análisis dinámico.
Taludes
Diseño de sistemas de estabilización (soil nailing, micropilotes, etc.)
Análisis de la estabilidad de los taludes.
Estudio del método constructivo.
Corta
Análisis dinámico.
Redes de filtración.
Minería Diseño de sostenimientos.
Análisis de convergencias.
Túneles Factor de seguridad.
Diseño de voladuras
Intersección de galerías
Fracturación hidráulica e inyección
Diseño de pozos
Oil & Gas Convencional
Diseño y dimensionamiento de prospecciones.
Compactación y análisis de subsidencias.
Cálculo de cimentaciones.
Nuclear
Análisis sísmico.
Cálculo de cimentaciones.
Aerogeneradores
Energía Análisis dinámico.
Cálculo de sostenimientos y revestimientos de túneles y cavernas.
Centrales
Diseño del método y las fases constructivas.
hidroeléctricas
Análisis de convergencias del terreno.

4.2 Características de FLAC3D

Un modelo de FLAC3D está compuesto de “zones” (“zonas”), desde unos cientos a unas miles
o millones de ellas. Por tanto, el modelo más pequeño que podría hacerse con FLAC3D estaría
constituido por una única “zone”. Generalmente, las zonas suelen ser de forma hexaédrica, es
decir, con ocho vértices y seis caras, aunque también pueden ser en forma de pirámide, cuña,
tetraedro o “d-brick”, todas ellas con menos vértices y caras que el hexaedro, tal y como
veremos más adelante.
Los vértices de las “zones” son los “gridpoints” (nodos). La combinación de las “zones” y los
“gridpoints” dan lugar a la malla de FLAC3D (“mesh” o “grid” en inglés). Durante el cálculo, las
ecuaciones mecánicas se resuelven para cada “zone” / “gridpoint”. La imagen que se presenta
a continuación muestra el detalle de una malla (“mesh” / “grid”) y las “zones” y “gridpoints” que
la componen.

12
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Figura 4.1 – Ejemplo de malla, “zone” y “gridpoints” en FLAC3D - I

En las figuras que se presentan a continuación se puede apreciar de nuevo los conceptos de
malla (mesh / grid), “zone” y “gridpoint” en un modelo de FLAC3D más complejo que el anterior.
El alumno puede apreciar que las “zones” en FLAC3D pueden ser de formas diversas y
tamaños distintos, no teniendo que ser todas iguales ni sus lados de la misma longitud.
Observará el alumno que las zonas, en todos los casos, son de forma hexaédrica, con seis
caras y ocho vértices.

Figura 4.2 – Ejemplo de malla, “zone” y “gridpoints” en FLAC3D - II

13
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Figura 4.3 – Ejemplo de malla, “zone” y “gridpoints” en FLAC3D - III

La malla de FLAC3D se construye en función del sistema global de referencia según los ejes
coordenados X, y, Z, tal y como se muestra en la figura que aparece a continuación. Todos los
“gridpoints” y los centroides de las “zone” está localizados y posicionados por su vector
posición (x, y, z). Así mismo, cada “gridpoint” y cada “zona” están identificados por un
identificador (ID) que permite hacer referencia a ellos mediante el lenguaje de programación
FISH, tal y como veremos más adelante. Así pues, podemos concluir que la generación de una
malla con FLAC3D implica combinar y ajustar las zonas para reproducir la geometría del
problema que queremos analizar. En el ejemplo del túnel mostrado en las figuras anteriores,
observamos como las zonas se distribuyen para tener en cuenta el macizo rocoso, el túnel
(núcleo) y el sostenimiento del mismo.

Figura 4.4 – Sistema de ejes coordenados en FLAC3D

4.3 Conceptos fundamentales en FLAC3D

Antes de empezar a utilizar el programa FLAC3D, conviene tener claros una serie de
conceptos fundamentales. Ello evitará confusiones y permitirá entender mejor todo lo que se
explicará en los siguientes módulos. Algunos de estos conceptos son generales y de aplicación
a la mayor parte de programas de cálculo numérico en ingeniería geotécnica, ya sean en
diferencias finitas o en elementos finitos. Otros son de aplicación exclusiva a FLAC3D.

14
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

• Modelo (model): Se denomina modelo a la combinación de una malla (“mesh” / “grid”)

junto con unas determinadas condiciones de contorno, condiciones iniciales, modelos
constitutivos de los elementos, etc. El modelo se diferencia de la malla en que este
permite abordar el análisis de un problema físico concreto, mientras que la malla es
simplemente una geometría. Es la combinación de la geometría (malla) con el resto de las
características (condiciones iniciales, de contorno, leyes constitutivas, parámetros, etc.) lo
que constituye el modelo.
• Zona (zone): Es el dominio geométrico más pequeño en el que se puede evaluar un
cambio en un estado determinado (por ejemplo tensión, deformación, etc.). Las zonas son
generalmente de tipo poliédrico (hexaedros) de seis caras y ocho vértices, aunque pueden
tener otras formas (tetraedro, “d-brick”, cuñas y pirámides). La combinación de un cierto
número de zonas (desde unos cientos hasta millones) permite construir la malla del
modelo. Tal y como se verá a continuación, cada zona contiene dos grupos de cinco
subzonas tetraédricas cada uno, es decir, cada zona está compuesta por diez tetraedros
(dos grupos de cinco). El usuario no suele ser consciente de estos tetraedros (ya que
FLAC3D hace la división de las zonas en tetraedros automáticamente), que forman parte
fundamental del esquema de resolución numérica que emplea FLAC3D pero no aparecen
en las gráficas ni figuras que proporciona el software, en las que sí aparecen las zonas (si
bien los tetraedros de las zonas son accesibles por medio del lenguaje de programación
interno FISH). El punto característico de una zona es su centroide, que está definido por
las coordenadas (x, y, z) según el sistema de referencia global del modelo. De hecho,
todas las propiedades de una zona (modelo constitutivo, parámetros, tensión,
deformación, etc.) están referidas a su centroide. Para referirnos a una zona, siempre lo
haremos a través de su centroide o de su identificador (ID).
• Nodo (gridpoint): Los gridpoints están asociados a los vértices (esquinas) de las zonas.
Su posición está definida por las coordenadas (x, y, z) según el sistema de referencia
global del modelo. Podemos acceder a la información de un nodo conociendo sus
coordenadas y/o su identificador (ID).
• Malla (mesh / grid): Es la combinación de zonas para dar lugar a una geometría en
concreto. Representa la zona “física” del modelo, y en ella se almacenan todas las
variables de estado del mismo. En lo que respecta a estas variables de estado, todas las
magnitudes vectoriales (fuerzas, velocidades y desplazamientos) se almacenan en los
nodos, mientras que las magnitudes escalares y tensoriales (propiedades materiales,
tensiones, etc.) se almacenan en las zonas.
• Contornos (límite) del modelo (model boundary): Los contornos del modelo son los
planos o superficies (regulares o irregulares) que constituyen la periferia de la malla.
Generalmente, suele hablarse del contorno superior (superficie libre) y el contorno inferior
(base), ambos generalmente perpendiculares al eje Z (la base normalmente será
perpendicular al eje Z, mientras que la superficie libre puede no serlo si es irregular), y los
contornos laterales, dos perpendiculares al eje X y dos perpendiculares al eje Y. Los
contornos interiores (por ejemplo, agujeros en el modelo que representan excavaciones)
también se denominan contornos del modelo. La figura que se presenta a continuación
muestra el detalle de los contornos de una malla de FLAC3D.

15
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Figura 4.5 – Contornos de una malla de FLAC3D

• Condiciones de contorno (boundary condition): Las condiciones de contorno son las

restricciones, requerimientos o condicionantes aplicados en los contornos del modelo.
Estas condiciones pueden ser de distinto tipo. Para los modelos mecánicos (los que se
estudiarán en el ámbito de este curso), las condiciones en los contornos suelen considerar
desplazamientos fijos en una o varias direcciones o fuerzas fijas en una o varias
direcciones. Así pues, por ejemplo, suele ser habitual bloquear los desplazamientos en
todas direcciones de los nodos de la base del modelo, y los desplazamientos en la
dirección perpendicular al plano del contorno en el caso de los contornos laterales. Esto
último constituye dos ejemplos de condiciones de contorno. En el caso de modelos de flujo
(fuera del ámbito de este curso), puede hablarse también de contornos permeables (el
agua puede atravesarlos) o impermeables (el agua no puede atravesarlos), lo cual
constituye otro ejemplo de condición de contorno.
• Condiciones iniciales (initial conditions): Son los valores de las variables del modelo
(tensiones, presiones intersticiales, etc.) en el estado inicial, es decir, antes de que se
produzca cualquier perturbación o fenómeno de carga y/o descarga. El estado o
condiciones iniciales suele coincidir con la situación de terreno virgen o en su estado
natural (condiciones de campo libre).
• Modelo constitutivo (constitutive model): Modelo matemático que determina la relación
tensión/deformación en las zonas del modelo. FLAC3D dispone de una amplia biblioteca
de modelos constitutivos que incluye modelos elásticos, plásticos y elasto-viscoplásticos.
El comportamiento elástico lineal y el modelo de Mohr-Coulomb son dos ejemplos de
modelos constitutivos. El modelo constitutivo puede asignarse de manera individual a cada
zona del modelo FLAC3D.
• Zona nula (null zone): Son zonas que tienen asignado el modelo constitutivo “nulo” (null).
Sirven para representar agujeros y excavaciones en la malla, es decir, regiones donde no
hay presencia de material.
• Submalla (subgrid): La malla de un modelo de FLAC3D puede estar compuesta de varias
submallas combinadas entre sí. La interacción entre submallas puede hacerse mediante
elementos “interfase” o mediante el comando “attach” (pegar / juntar). Por ejemplo, en la
figura que se muestra a continuación, la malla del modelo está compuesta por dos
submallas, una la del terreno, y otra la del encepado. La interacción entre ambas se
produce a través de la interfase colocada en los contornos del encepado (en verde en la
figura).

16
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Figura 4.6 – Ejemplo de malla y submalla en FLAC3D

• Interfases (interfaces): Las interfases son elementos bidimensionales (superficies) que

permite conectar submallas entre sí y que pueden experimentar separación (por
deslizamiento o por apertura) durante el cálculo. Las interfases son muy útiles a la hora de
representar discontinuidades físicas como las fallas en el terreno (en particular en los
macizos rocosos), planos de contacto entre formaciones distintas, o para tener en cuenta
la interacción entre una estructura y el terreno de su alrededor (interacción terreno-
estructura). Este último caso es el representado en la figura anterior.
• Rango (range): Un rango en FLAC3D es un “filtro” o “condicionante” impuesto por el
usuario que restringe las operaciones de los comandos de FLAC3D, es decir, los
comandos se aplican únicamente a los objetos (zonas, nodos, elementos estructurales,
interfases, etc.) que satisfacen y quedan dentro de la definición (o dominio) del rango. Por
ejemplo, si quisiéramos dividir una malla en dos estratos diferentes, siendo el superior una
arena y el inferior una roca, podríamos asignar los parámetros para cada tipo de terreno
aplicando rangos. Así pues, para definir la arena, asignaríamos las propiedades de esta a
las zonas que estuvieran dentro del rango z base,arena y ztop, donde zbase,arena es la cota de la
base del estrato de arena y ztop la cota de la superficie del modelo. Para definir la roca,
procederíamos de modo similar, pero asignando las propiedades de esta a las zonas que
estuvieran dentro del rango zbase,arena y zbot donde zbot sería la cota de la base del modelo.
La figura que se muestra a continuación muestra un detalle del ejemplo descrito. Un rango
puede o no estar asociado con objetos específicos del modelo. Si el rango está definido en
función de los identificadores (ID) de los elementos o por grupos (se explicarán a
continuación), dicho rango estará asociado con objetos específicos del modelo. Por el
contrario, en el caso del ejemplo anterior, si un rango está definido en un dominio espacial
(x, y, z), el rango está fijo en el espacio, pero los elementos dentro de él pueden variar
entre una fase y otra del análisis.

17
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Figura 4.7 – Ejemplo de rangos en FLAC3D

• Elementos de rango (range elements): Se denomina elemento de rango a cualquier
criterio que define un rango. Todo rango tendrá por lo menos un elemento de rango,
pudiéndose combinar con otros elementos de rango, con el objetivo de definir uniones e
intersecciones de los elementos de rango. Siguiendo con el ejemplo anterior de la malla
dividido en dos estratos, el superior arena y el inferior roca, la instrucción por la cual se
asignan las propiedades de la arena a las zonas entre z base,arena y ztop sería un elemento de
rango. El rango en cuestión estaría constituido por un solo elemento de rango (posición de
las zonas según el eje Z). Si además consideráramos que solo una parte a la izquierda del
modelo (según el eje X) comprendido entre zbase,arena y ztop fuese arena, deberíamos incluir
al rango la instrucción de sólo considerar las zonas entre xmin y xarena, donde xmin sería la
coordenada X del contorno de la malla y xarena la coordenada X del límite del estrato de
arena según este eje. En este caso, el rango en cuestión estaría constituido por dos
elementos de rango, uno determinando la posición de las zonas según el eje X, y el otro
determinando la posición según el eje Z. La figura que se presenta a continuación muestra
un detalle del ejemplo descrito.

Figura 4.8 – Ejemplo de elementos de rangos en FLAC3D

18
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

• Grupo y slot (group and slot): Un grupo en FLAC3D hace referencia a un nombre que
puede asignarse a uno o más objetos, ya sean estos zonas, nodos, interfases, elementos
estructurales, etc. Este nombre resulta muy útil cuando se quiere llamar a todos los
objetos (o una parte de ellos) que lo tienen asignado para dotarles de un nuevo valor en
alguna de sus propiedades, aplicarles una determinada carga, excavarlos, etc. Los grupos
pertenecen a una categoría o nivel superior denominado slot. Cada elemento del modelo
(nodo, zona, elemento estructural, etc.) solo puede pertenecer a un grupo dentro de cada
slot. Es decir, un mismo elemento puede estar asignado a varios grupos, pero cada uno de
estos grupos debe pertenecer a un slot diferente. En cierta forma, podríamos decir que
cada slot es un “armario” en el que introducimos todos los elementos del modelo. Cada
uno de estos “armarios” puede estar organizado en “cajones” de diferente manera, de
forma que un mismo objeto, en cada “armario”, tendrá asignado un “cajón” distinto. Cada
uno de estos cajones sería un grupo. Las diferencias (y relación) entre slot y grupo puede
apreciarse en las siguientes figuras:

Figura 4.9 – Ejemplo de group y slot en FLAC3D - I

Figura 4.10 – Ejemplo de group y slot en FLAC3D - II

En la primera de las figuras, vemos como en el slot “default” (slot por defecto), las zonas
del modelo se distribuyen en los grupos “encepado” y “terreno”. En la segunda figura, sin
embargo, vemos como en el slot “1” las zonas se distribuyen en los grupos “arena” y

19
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

“roca”. Así pues, vemos como una misma zona (o elemento del modelo) puede pertenecer
a varios grupos, pero cada grupo pertenecerá a un slot diferente. Es decir, un objeto solo
puede pertenecer a un solo grupo dentro de cada slot. A diferencia de los rangos, que
pueden hacer referencia a objetos en concreto o zonas del espacio, un grupo SIEMPRE
hará referencia a objetos específicos del modelo, no pudiendo asignarse a una
determinada zona del espacio como el rango. Los grupos permiten hacer referencia y
“llamar” fácilmente a determinadas zonas del modelo para operar en ellas a través de los
comandos o el lenguaje FISH. Esto se lleva a cabo indicando el nombre del grupo de
interés en el rango que va asociado al comando que se quiere ejecutar sobre los
elementos del grupo. La asignación de grupos, slots y rangos se verá en detalle en los
módulos siguientes.
• Número ID (ID number): Cada elemento individual de FLAC3D, sea del tipo que sea
(zona, nodo, elemento estructural, interfase, etc.) está identificado por un número ID,
denominado también “identificador” (identifier). Estos números ID son números únicos que
asigna el software a los elementos, permitiendo al usuario identificar y “llamar” a
elementos específicos del modelo. De modo similar a los números ID, FLAC3D también
asigna a los elementos individuales de los entes estructurales (nodos, elementos, enlaces
nodos/malla, etc.) que constituyen las vigas, pilotes, placas, geomallas, etc. los
denominados números id de componente (id-component). Así pues, un pilote de 20 metros
de longitud, formada por elementos de 1 metro de longitud (20 elementos en total), y por lo
tanto 21 nodos, tendrá un número id que identificará a la totalidad del pilote, mientras que
cada uno de los elementos y nodos tendrá su propio número id de componente. Es decir,
este pilote tendrá un número id y 20 + 21 = 41 números id de componente (uno para cada
elemento y nodo que lo constituye).
• Nombres (names): Para determinados elementos del modelo, como historias, tablas e
interfases (que se explicarán más adelante), resulta conveniente que el usuario les asigne
un nombre corto, claro y fácil de recordar para poder referirse a ellos a lo largo del cálculo.
Generalmente, esto se lleva a cabo asignando a estas entidades un nombre (name), que
suele asignarse cuando estos elementos son creados. Un nombre (name) es una variable
tipo “string”, normalmente corta, aunque puede tener cualquier longitud. Un nombre puede
estar formado por un número entero (integer, por ejemplo, 1000, 1500, etc.). En este caso,
el software convierte automáticamente el número entero en una variable tipo string.
• Elemento estructural (structural element): En primer lugar, debemos distinguir entre
elemento estructural y entidad estructural. Una entidad estructural está compuesta de una
combinación de elementos estructurales. Podemos decir que FLAC3D incorpora dos tipos
de elementos estructurales básicos. El primero de ellos son los elementos lineales, de dos
nodos (uno por extremo), que permiten modelizar vigas, cables y pilotes, que son
entidades lineales unidimensionales. Generalmente, una viga, pilote o cable estará
compuesto por varios elementos lineales. Por ejemplo, una viga de 10 metros puede
dividirse en elementos de 0,5 metros (un total de 20), lo que implicaría un total de 41
nodos. El segundo tipo son los elementos triangulares de 3 nodos, cuya combinación
permite modelizar placas, láminas, sostenimientos (de túnel, muro pantalla, etc.) y
geomallas, que son entidades superficiales bidimensionales. Los elementos estructurales
se utilizan para simular la interacción de los soportes estructurales (pilotes, anclajes,
pantallas, etc.) con el suelo o con el macizo rocoso. Todos los elementos estructurales
admiten un comportamiento elástico lineal, y con alguno de ellos puede tenerse en cuenta
un comportamiento plástico no lineal. Por tanto, cada entidad estructural en FLAC3D está
formada por tres componentes: nodos estructurales, elementos individuales y las
conexiones nodo estructural / malla. Las características de cada uno de estos
componentes, tal y como veremos, permite distinguir entre viga, cable, pilote, placa,
lámina, sostenimiento y geomalla.
• Paso de cálculo (step): Dado que FLAC3D es un código de cálculo explícito, la solución
a un problema requiere de un cierto número de pasos de cálculo computacionales (steps).
Durante el cálculo computacional, la información asociada con el fenómeno que se está
investigando (descarga, carga, excavación, etc.) se propaga por las zonas y nodos de la
malla. Por tanto, es necesario un cierto número de pasos de cálculo (steps) para llegar a

20
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

una situación de equilibrio en el caso de un problema estático. Al paso de cálculo (step) en

FLAC3D también se le suele denominar timestep (misma traducción en español, paso de
cálculo) y cycle (ciclo).Tal y como se verá en la sección 5 de este módulo, en el caso de
cálculos para situaciones estáticas con modelos elásticos o elastoplásticos, los pasos de
cálculo son cantidades “artificiales”, es decir, no representan el tiempo real, sino que
únicamente se emplean en la resolución de las ecuaciones del movimiento (ecuaciones
dinámicas dependientes del tiempo) para llegar a un estado de equilibrio estable.
• Solución estática (static solution): En FLAC3D, se considera que un modelo estático ha
alcanzado una solución estática (static solution) o de estado estable (steady-state) cuando
el ratio de cambio de la energía cinética del modelo se aproxima a un valor nulo, es decir,
los nodos del modelo tienden a una velocidad nula o constante (aceleración nula). Esto se
consigue mediante la amortiguación de las ecuaciones del movimiento. Al concluir un
cálculo estático, el modelo está o bien en un estado de equilibrio o bien en un estado de
fluencia plástica en el caso de que una parte (o la totalidad) del modelo sea inestable bajo
las acciones a las que ha sido sometido. Como ejemplo de lo anterior, podemos decir que
una situación de estado en equilibrio se consigue cuando, tras aplicar unas cargas, una
cimentación determinada asienta una cantidad limitada. Por el contrario, si se aplica una
carga igual o superior a la carga última de la cimentación, está asentará de forma
incontrolada (a velocidad constante), puesto que se habrá producido su hundimiento. Este
sería un caso de fluencia plástica. Lo mismo podría decirse en el caso de un problema de
estabilidad de taludes. Para hallar la rotura del talud, y por tanto su factor de seguridad, es
necesario “inducir” en el modelo un estado de fluencia plástica, en el que se produce la
inestabilidad de una zona del modelo (la rotura del talud), que llegado un cierto número de
pasos de cálculo se desplazará a velocidad constante. El cálculo de modelos estáticos es
el cálculo fundamental y por defecto de FLAC3D, y el que veremos en este curso, siendo
los cálculos dinámicos o de tipo “creep” más complejos y quedando fuera del alcance de
este curso.
• Fuerza no equilibrada (unbalanced force): La fuerza no equilibrada (unbalanced force)
indica al usuario cuando se alcanza un estado de equilibrio (o de fluencia plástica) en el
caso de un análisis estático. Se considera que el modelo está en equilibrio exacto cuando
la resultante de las fuerzas nodales en cada nodo (gridpoint) del modelo es nula. A la
resultante de las fuerzas nodales en los nodos (no nula) se le denomina fuerza no
equilibrada o fuerza desequilibrada (out of balance force o unbalanced force en inglés). Es
importante tener en cuenta que, para un cálculo estático en FLAC3D, la resultante de las
fuerzas nodales nunca será exactamente nula. Sin embargo, se considera que el modelo
está en equilibrio cuando esta resultante de las fuerzas desequilibradas es muy pequeña
comparada a las fuerzas totales aplicadas en el problema. En el caso de una situación de
equilibrio, la fuerza desequilibrada iría disminuyendo progresivamente a medida que
aumentamos los ciclos (steps) de cálculo (si bien no llegaría a una situación de cero
absoluto). Por el contrario, en una situación de fluencia plástica, la fuerza desequilibrada
alcanzaría un valor constante no nulo, independientemente del número de ciclos (steps)
de cálculo.
• Pequeñas y grandes deformaciones (small / large strains): Por defecto, FLACD
funciona en el modo de pequeñas deformaciones (small strain), que es el modo en el que
trabajaremos en este curso. En la modalidad de pequeñas deformaciones, las
coordenadas de los puntos no se actualizan con los desplazamientos, por muy grandes
que estos sean. Es decir, todos los nodos del modelo mantienen su posición espacial a lo
largo de todo el cálculo. La malla tiene siempre la misma forma. Por el contrario, en el
modo de grandes deformaciones (large strain), la posición de los nodos se actualiza en
cada paso de cálculo (step) según los desplazamientos que experimenta dicho nodo, de
forma que la posición espacial de los nodos varía continuamente a lo largo del cálculo, no
manteniendo la malla una geometría constante, sino variable con el cálculo al modificarse
continuamente la posición de los nodos.
• Proyecto (project): Por proyecto (project) en FLAC3D se hace referencia al conjunto de
ficheros de input (.dat, etc.) y de ouptut (.sav, etc.) que componen un modelo de FLAC3D.
Un proyecto en FLAC3D se almacena en un fichero de proyecto (project file) el cual, al ser

21
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

ejecutado, “llama” a todos los ficheros de input y output vinculados al proyecto para que
sean accesibles al usuario desde la interfaz gráfica de FLAC3D. En el segundo módulo de
este curso entraremos en el detalle de la estructura de proyecto de FLAC3D.
• Primitiva (primitive): Se denominan primitivas (primitives) en FLAC3D a un conjunto de
geometrías ya predefinidas en FLAC3D cuya combinación permite crear mallas para
construir los modelos. Las diferentes primitivas de FLAC3D se explicarán en detalle en el
cuarto módulo de este curso. La figura que se presenta a continuación muestra algunas de
estas primitivas, como el tetraedro (tetrahedron), la cuña (wedge), la pirámide (pyramid) y
el hexaedro (brick). Las primitivas de FLAC3D están formadas por zonas.

Figura 4.11 – Ejemplos de primitivas en FLAC3D - II

• Comando (command): Se denomina comando de FLAC3D a una determinada instrucción

dada por el usuario con el fin de introducir cambios en el modelo. Los comandos de
FLAC3D, que se verán con detalle en el módulo 3 del curso, sirven, entre otras funciones,
para crear geometrías, asignar parámetros, introducir condiciones de borde, inicializar
tensiones, etc.
• FISH (FISH): FISH (abreviación de FLAC-ISH) es el lenguaje interno de programación de
FLAC3D. Este lenguaje de programación es la clave de la versatilidad y la flexibilidad que
aporta el software al usuario, permitiéndole tener un control absoluto del modelo numérico
y de sus parámetros. Mediante FISH, el usuario puede definir funciones (también llamadas
rutinas) que, en combinación con los comandos, le permite llevar a cabo de forma
automática el lanzamiento de distintos ficheros de cálculo, salvar el estado del análisis,
asignar parámetros a determinadas zonas del modelo, definir secuencias de excavación y
construcción, entre otras posibilidades. Incluso, el lenguaje FISH permite al usuario
introducir nuevas variantes físicas y matemáticas al esquema de cálculo predefinido de
FLAC3D. Por ejemplo, mediante una rutina FISH debidamente definida, el usuario podría
hacer variar el módulo de deformación y/o la cohesión del terreno (al que se ha asignado
un modelo de comportamiento Mohr-Coulomb) con la deformación de corte. Aunque el
modelo Mohr-Coulomb “de serie” de FLAC3D no contemple esta posibilidad, el usuario
podría tenerla en cuenta por medio de una función FISH.

22
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

5 FUNDAMENTOS DEL FUNCIONAMIENTO DE FLAC3D

5.1 Formulación en diferencias finitas en FLAC3D

Tal y como se ha mencionado anteriormente, FLAC3D es un programa de cálculo numérico

generalista basado en el método de las diferencias finitas con un esquema de resolución
explícito. En este apartado, se describirá de modo conceptual, sin entrar en todo el rigor y
detalle matemático (que no es objeto de este curso), la formulación explicita en diferencias
finitas que emplea el software FLAC3D para el análisis y la resolución de problemas
geotécnicos, geológicos y mineros.
El método de las diferencias finitas es el método más antiguo para la resolución de un conjunto
de ecuaciones diferenciales en derivadas parciales, dados un conjunto da valores iniciales y
unas condiciones de contorno. En el método de las diferencias finitas, cada derivada del
conjunto de ecuaciones diferenciales se sustituye directamente por una expresión algebraica
que relaciona las variables de campo (desplazamientos, velocidades, etc.) en zonas
específicas y discretas de la malla (nodos y zonas). Por tanto, el método de las diferencias
finitas, como el de los elementos finitos, implica la discretización (o división) de un medio
continuo en “partes” discretas e individuales (nodos y zonas) en las que se resolverán las
ecuaciones del movimiento (nodos) y constitutivas (zonas) para obtener, por un lado,
desplazamientos y velocidades y por otro tensiones y deformaciones.
De manera simplificada, puede decirse que un método de análisis numérico de resolución de
ecuaciones diferenciales temporales (dependientes del tiempo) es explícito cuando este
permite calcular y determinar el estado (inicialmente desconocido) de un sistema, medio o
problema objeto de análisis en un tiempo futuro (t n+1), a partir del estado y las condiciones del
sistema, medio o problema en el tiempo presente o actual (tn), que resultan conocidos. Por el
contrario, el método implícito permite calcular y determinar el estado (inicialmente desconocido)
de un sistema, medio o problema objeto de análisis en un tiempo futuro (t n+1), a partir del
estado y las condiciones del sistema, medio o problema en el tiempo presente o actual (t n), que
son conocidas, y en el tiempo futuro (tn+1), que, como se ha dicho, son desconocidas.
A partir de lo anterior, podemos concluir que un método de análisis numérico de resolución de
ecuaciones diferenciales temporales es explícito cuando el estado del problema en un tiempo
futuro (tn+1) puede calcularse en términos de cantidades conocidas (el estado del sistema en el
tiempo tn). Cuando el estado del problema en un tiempo futuro (t n+1) se calcula a partir de
magnitudes conocidas (en tn) y desconocidas (tn+1), se dice que el esquema de resolución del
método de análisis numérico es implícito. En este último caso, es necesario emplear técnicas
matriciales o iterativas (caso de los elementos finitos, que usan fundamentalmente un esquema
de resolución implícito) para obtener la solución del problema.
Todo lo anterior podría expresarse matemáticamente de forma sencilla considerando la
siguiente ecuación diferencial:
𝑦̇ = 𝑓(𝑦, 𝑡)
Mediante un método de análisis numérico explícito, la anterior ecuación diferencial podría
expresarse como:
𝑦𝑛+1 = 𝑦𝑛 + ∆𝑡 ∙ 𝑓(𝑦𝑛 , 𝑡𝑛 )
Vemos como, efectivamente, el estado futuro del sistema (y n+1) depende exclusivamente del
estado presente o actual del mismo, que resulta conocido (y n y tn). La magnitud Δt es el paso de
tiempo entre el estado yn+1 y yn.
Por el contrario, para un método de análisis numérico implícito, la ecuación diferencial anterior
se expresaría como:
𝑦𝑛+1 = 𝑦𝑛 + ∆𝑡 ∙ 𝑓(𝑦𝑛+1 , 𝑡𝑛+1 )
Vemos pues que, para el método implícito, el estado futuro del sistema (yn+1) se calcula a partir
del estado presente o actual del mismo, que resulta conocido (y n), pero también a partir del
estado futuro del sistema (yn+1), que resulta desconocido y aparece en ambos lados de la

23
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

ecuación. Por tanto, para este último método es necesario utilizar técnicas de cálculo matricial
e iterativas.
Otra forma de expresar matemáticamente la diferencia entre un método explícito y un método
implícito sería la que figura a continuación, considerando que Y(t) es el estado actual del
sistema e Y(t+ Δt) es el estado futuro del sistema, siendo Δt es el paso de tiempo entre el
estado futuro y el actual. Teniendo en cuenta que, para ambos casos (explícito e implícito) es
necesario resolver para Y(t+ Δt), obtendríamos:
Método explícito: 𝑌(𝑡 + ∆𝑡) = 𝐹(𝑌(𝑡))
Método implícito: 𝐺(𝑌(𝑡), 𝑌(𝑡 + ∆𝑡)) = 0
Observamos de nuevo que, para el método explícito, la solución del estado futuro del sistema
depende exclusivamente del estado presente o actual que es conocido, mientras que para el
método implícito, la solución del estado futuro depende tanto del estado futuro (desconocido)
como del estado actual o presente (conocido).
En FLAC3D, es el usuario el que lleva a cabo la discretización del medio o problema objeto de
análisis en zonas hexaédricas. Posteriormente, FLAC3D, de manera automática, discretiza
estas zonas hexaédricas en grupos de tetraedro cuyos vértices son los nodos (gridpoints) de la
malla, tal y como se ha mencionado anteriormente. Generalmente, cada zona hexaédrica se
divide en dos grupos (overlays) de 5 tetraedros cada uno, tal y como muestra la siguiente
figura.

Figura 5.1 – División de las zonas en elementos tetraédricos

Así pues, las fuerzas nodales se evalúan a partir de la media entre los dos grupos de
tetraedros. Cada uno de estos tetraedros tiene una deformación y tensión constante (única y
homogénea) dentro del propio tetraedro, que irá variando a lo largo del cálculo. La deformación
volumétrica de la zona se computa como la media de los tetraedros que la componen, mientras
que la deformación desviadora se calcula para cada tetraedro independientemente. La división
de las zonas en tetraedros se debe a la versatilidad de este elemento geométrico a la hora de
aplicar las ecuaciones matemáticas del esquema de resolución numérica. Para el alcance de
este curso, es suficiente que el alumno conozca que las zonas del modelo se dividen
internamente en elementos tetraédricos a los cuales puede acceder mediante el lenguaje FISH.
No entraremos en el detalle de las ecuaciones aplicadas a estos elementos.
El método explícito de resolución en diferencias finitas de FLAC3D está basado en el principio
de propagación y disipación de la energía cinética en un cuerpo deformable en movimiento. La
solución numérica implementada en el software se fundamenta en el empleo de las ecuaciones
dinámicas del movimiento. Es decir, se trata de resolver un problema estático por medio de las
ecuaciones dinámicas del movimiento. El ciclo de cálculo explícito básico de resolución de
FLAC3D, referido a un paso de cálculo (step) Δt, se presenta de forma esquemática en la figura
que se muestra a continuación.

24
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Figura 5.2 – Esquema de representación del ciclo de cálculo explícito básico en FLAC3D
Tal y como puede apreciarse en la figura anterior, el ciclo de cálculo en FLAC3D se inicia a
partir del cálculo de las aceleraciones de los nodos del modelo, empleando para ellos las
ecuaciones dinámicas del movimiento, siendo conocidas las fuerzas que actúan en ellos
(primer recuadro de la figura anterior). Una primera integración de las aceleraciones en el paso
de cálculo Δt/2 permite obtener las velocidades de los nodos y una segunda integración en el
siguiente paso de cálculo Δt/2 permite obtener los desplazamientos de los nodos (segundo
recuadro de la figura anterior). A partir de las velocidades y los desplazamientos de los nodos,
se pueden calcular las velocidades de deformación y deformaciones de las zonas y, empleando
las ecuaciones constitutivas (relaciones tensión/deformación), se obtienen las tensiones en las
zonas del modelo (tercer recuadro de la figura anterior). Finalmente, a partir de las tensiones
calculadas en las zonas, se calculan las nuevas fuerzas en los nodos (cuarto recuadro de la
figura anterior), iniciándose así un nuevo ciclo de cálculo y repitiéndose el proceso descrito de
nuevo.
Observando el ciclo de cálculo de FLAC3D anterior, podría entenderse que el cálculo
continuará hasta que la resultante de las fuerzas actuando sobre los nodos del modelo sea
nula, y por tanto los nodos estén en equilibrio perfecto. Esto es así en el sentido estricto, formal
y teórico. Sin embargo, llegado un determinado momento del cálculo, las fuerzas
desequilibradas de los nodos serán muy pequeñas en comparación con las fuerzas que actúan
sobre ellos en ese momento. Estás pequeñas fuerzas desequilibradas producirán cambios muy
poco significativos en los desplazamientos y velocidades de los nodos y en las deformaciones y
tensiones de las zonas (cambios que únicamente se observarán a partir del cuarto o quinto
decimal). Por tanto, podría considerarse que el modelo ha alcanzado un estado de equilibrio,
aunque no en el sentido exacto del término, pues sigue existiendo una componente de fuerzas
desequilibrada en los nodos, aunque es despreciable. Así pues, en FLAC3D, es el usuario
quien define la magnitud de la fuerza desequilibrada en los nodos a partir de la cual se
considera que el modelo está en estado de equilibrio y que dicha fuerza desequilibrada ya no
tiene influencia significativa sobre la magnitud de los resultados. Generalmente, y tal y como se
explicará posteriormente en el curso, se considera que el modelo está en equilibrio cuando el
ratio entre las fuerzas desequilibradas y las fuerzas totales actuando sobre los nodos es igual o
inferior a 1 x 10-5. Dicho de otra manera, los nodos de un modelo de FLAC3D que se considera
en equilibrio no estarán totalmente “quietos” (velocidad nula), pero sí lo suficientemente
“quietos” para considerar que el modelo está en equilibrio.

25
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

Del esquema del ciclo de cálculo en FLAC3D es importante que el alumno retenga los
siguientes aspectos:
• Las aceleraciones, velocidades y desplazamientos se calculan y son atributos de los
NODOS del modelo. Las tensiones y deformaciones se calculan y son atributos de las
ZONAS del modelo.
• En cada uno de los recuadros de la figura anterior, el valor de cada una de las variables se
modifica a partir del valor de la variable en el ciclo de cálculo anterior, lo cual demuestra
que se trata de un esquema de cálculo explícito. De hecho, al iniciarse el cálculo (t 0), son
conocidas las fuerzas en los nodos, las tensiones en las zonas y las aceleraciones,
velocidades y desplazamientos de los nodos (todos nulos al tratarse del estado inicial). El
estado del modelo inmediatamente posterior (t1), se calcula por tanto a partir de las
magnitudes del estado anterior (t0), que son conocidas, lo que se corresponde, tal y como
se ha explicado anteriormente, con un método de resolución explícito.
• El paso de cálculo Δt es igual a la unidad (Δt=1). Ello es así puesto que un esquema de
resolución explicito precisa de pasos de cálculo pequeños para garantizar su estabilidad
(al contrario que el método implícito, que puede considerar pasos de cálculo mayores). Por
otro lado, empleando un paso de cálculo lo suficientemente pequeño, se considera
aceptable que, durante un mismo ciclo de cálculo, las magnitudes calculadas en los tres
primeros recuadros (aceleraciones, velocidades y desplazamientos en los nodos y
deformaciones en las zonas) permanezcan constantes a lo largo del ciclo aunque se
calculen nuevas tensiones y fuerzas nodales al final de dicho ciclo, con las que no serían
compatibles y estarían en desequilibrio. Por otro lado, puesto que el paso de cálculo es lo
suficientemente pequeño, se impide que la información de este “desequilibrio” dentro de la
zona en cuestión se propague a las zonas vecinas del modelo dentro del mismo ciclo de
cálculo, permitiendo así al modelo converger gradualmente. En otras palabras, podemos
decir que, para un ciclo de cálculo dado, las zonas están “aisladas” unas de otras, no
pudiéndose propagar los desequilibrios internos generados en las zonas (para el ciclo de
cálculo considerado) entre sus zonas vecinas/adyacentes en ese mismo ciclo de cálculo.
Estos desequilibrios solo se propagarán entre las zonas al final de cada ciclo de cálculo,
justo antes de aplicar de nuevo las ecuaciones del movimiento, que nos permitirá obtener
nuevas aceleraciones y, por tanto, velocidades, desplazamientos y deformaciones.
• Para los cálculos estáticos con modelos elastoplásticos, que constituyen el contenido de
este curso de introducción a la modelización con FLAC3D, los pasos de cálculo son una
magnitud “artificial”, es decir, no se corresponden con el tiempo real ni tienen sentido
físico. Es tan solo un artificio que emplea el software para resolver las ecuaciones del
movimiento, que son dependientes del tiempo. En cálculos en los que sí debe tenerse en
cuenta el tiempo como variable, como es el caso de los cálculos con fluencia (creep), que
están fuera del alcance de este curso, el paso de cálculo es una magnitud con sentido
físico que representa el tiempo real.

5.2 Comparación con el método de los elementos finitos. Ventajas y

desventajas

Las principales ventajas de FLAC3D con respecto a los softwares de cálculo numérico basados
en elementos finitos (generalmente con métodos de resolución implícitos), son las siguientes:
1) Resulta mucho más apropiado para modelizar problemas y sistemas donde se produzcan
grandes plastificaciones, importantes fenómenos de fluencia plástica, situaciones de
colapso (inestabilidad física) y grandes desplazamientos, deformaciones y rotaciones. Al
contrario que en los softwares con un esquema de resolución implícito (generalmente de
elementos finitos), donde es necesaria la convergencia estricta del cálculo, FLAC3D
siempre proporcionará una “solución” para estas situaciones, puesto que permite simular,
estudiar y observar los fenómenos de inestabilidad como si de una “película” se tratase,
dado que su esquema de resolución se basa en la resolución progresiva de las
ecuaciones del movimiento. En otras palabras, en FLAC3D los fenómenos de inestabilidad
física no llevan aparejado un problema de inestabilidad numérica, como sí ocurre en el

26
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

caso de los elementos finitos, donde es generalmente necesario llegar a una convergencia
estricta. En este sentido, FLAC3D resulta mucho más eficiente, sobre todo en lo que a
tiempos de cálculo se refiere.
2) FLAC3D resulta más adecuado para analizar problemas caracterizados por una no
linealidad física y/o material que la mayoría de softwares de elementos finitos. Es decir,
FLAC3D puede considerar relaciones tensión/deformación no lineales con el mismo
tiempo de cálculo que las relaciones lineales, mientras que los elementos finitos (con
métodos implícitos) pueden tardar considerablemente más para resolver problemas no
lineales.
3) FLAC3D puede considerar cualquier modelo constitutivo sin ningún cambio en el algoritmo
de resolución. Muchos softwares de elementos finitos precisan de algoritmos o técnicas de
resolución distintas para tener en cuenta diferentes modelos constitutivos.
4) Tal y como se ha explicado anteriormente, el método implícito (generalmente empleado en
elementos finitos), precisa de métodos iterativos y matriciales. Ello implica disponer de una
importante cantidad de memoria RAM para almacenar la matriz de rigidez del sistema, que
debe ser invertida para la resolución del sistema de ecuaciones. FLAC3D, basado en un
método explícito, no precisa de una importante cantidad de memoria RAM al no tener que
almacenar ni construir ninguna matriz de rigidez ni invertirla. El hecho de no tener que
almacenar ni construir una matriz de rigidez hace que FLAC3D sea más eficiente a la hora
de analizar problemas en grandes deformaciones. Dado que en este tipo de problemas la
posición de los nodos varía a lo largo del cálculo (la geometría de la malla se va
modificando), la matriz de rigidez debe ser actualizada de forma continua. Este paso no es
necesario en FLAC3D al no tener que emplear métodos matriciales.
Por otro lado, las principales desventajas de FLAC3D, frente a los elementos finitos basados en
métodos de resolución implícitos fundamentalmente:
1) Los tiempos de cálculo son mucho más largos en FLAC3D que en elementos finitos para
problemas lineales o en modelos de comportamiento estable (sin inestabilidad física ni
grandes plastificaciones). Ello se debe, entre otras razones, al pequeño paso de cálculo
que debe emplearse en FLAC3D para asegurar la estabilidad de la solución (método
explícito), en comparación con el paso de cálculo en elemento finitos (método implícito),
que puede ser mucho mayor.
2) FLAC3D también operará con tiempos de ejecución más largos en modelos donde haya
una gran diferencia de rigideces entre materiales adyacentes o de tamaño entre zonas
vecinas.

27
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

6 EVALUACIÓN
A continuación, se presenta una colección de 15 preguntas tipo test para evaluar los
conocimientos adquiridos por el alumno en este módulo. Cada pregunta tipo test tiene cuatro
opciones, siendo únicamente una de ellas la correcta. Cada pregunta correcta puntuará +0,67,
cada pregunta incorrecta restará -0,3 y las preguntas no contestadas no computarán.
Las respuestas a las preguntas deben subirse a la plataforma de Ingeoexpert en un fichero
Excel con el siguiente nombre:
Nombre alumno_Curso FLAC3D Ingeoexpert_Modulo 1_Test
En la columna A del fichero Excel, empezando por la fila 1 y hasta la fila 15, se debe indicar el
número de pregunta. En la columna B, empezando por la fila 1 y hasta la fila 15, se debe
indicar la respuesta a cada pregunta. ¡SUERTE!
1) La modelización numérica, comparada con los cálculos a mano y las soluciones
analíticas o empíricas es:
a. Más precisa, pero implica mayor tiempo y coste.
b. Menos precisa e implica mayor tiempo y coste.
c. Más precisa, con menor tiempo y coste.
d. Menos precisa, con menor tiempo y coste.
2) Una solución numérica a un problema de ingeniería:
a. Siempre será correcta.
b. Será correcta siempre y cuando el ingeniero comprenda el software que esté
utilizando y utilice la información y parámetros adecuados.
c. No suele ser nunca correcta.
d. Depende del tipo de software que se utilice.
3) La principal diferencia entre la modelización numérica en geotecnia y en otros
campos de la ingeniería radica en:
a. La dificultad matemática de la disciplina en cuestión.
b. Los softwares que se utilizan.
c. La potencia de los ordenadores necesaria.
d. La cantidad y la fiabilidad de la información de la que se dispone.
4) El modelo numérico servirá más bien como “laboratorio numérico” si:
a. La información es abundante y fiable.
b. El ordenador es de baja potencia de cálculo.
c. La información disponible es escasa y poco fiable.
d. Ninguna de las anteriores es correcta.
5) Antes de llevar a cabo una simulación numérica, es necesario:
a. Definir objetivos y expectativas, por este orden.
b. Definir expectativas y objetivos, por este orden.
c. No definir ni objetivos ni expectativas.
d. No es necesario llevar a cabo ningún procedimiento.

28
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

6) Los modelos numéricos simplificados pueden ser necesarios previamente cuando:

a. Se trata de un modelo complejo y muy sofisticado, para comprender mejor su
comportamiento y mecanismos.
b. Cuando no se disponga de la capacidad de cálculo suficiente.
c. Nunca son necesarios. Son una pérdida de tiempo.
d. Siempre son necesarios.
7) El modelo conceptual matemático detallado se diferencia del modelo conceptual en
que:
a. No se diferencian, son equivalentes.
b. Está hecho a mano y el modelo conceptual a ordenador.
c. Además de integrar la información recopilada y el modelo conceptual, considera
también los objetivos que se persiguen con el modelo y los resultados de los modelos
simplificados.
d. Es un paso previo al modelo conceptual.
8) Cuando el modelo numérico en cuestión no responde a las expectativas, los
cambios deben hacerse:
a. Todos a la vez, para ahorrar tiempo y costes.
b. De uno en uno, para comprender bien su influencia en los resultados.
c. No es necesario hacer cambios. Un modelo numérico siempre es correcto.
d. Siempre con el objetivo de obtener los resultados que me sean más favorables.
9) FLAC3D es un software de cálculo numérico:
a. En elementos finitos con método explícito.
b. En elementos finitos con método implícito.
c. En diferencias finitas con método implícito.
d. En diferencias finitas con método explícito.
10) Un modelo de FLAC3D está compuesto por:
a. Zonas (zones) generalmente hexaédricas, siendo sus vértices los nodos (gridpoints).
b. Únicamente de zonas (zones).
c. Únicamente de nodos (gridpoints).
d. Por elementos finitos.
11) La diferencia entre “malla” y “modelo” radica en:
a. Son lo mismo exactamente.
b. La malla integra, además de la geometría del problema a estudiar, las condiciones
iniciales, condiciones de contorno, modelos constitutivos, etc. y permite analizar un
problema físico completo.
c. El modelo integra, además de la geometría del problema a estudiar, las condiciones
iniciales, condiciones de contorno, modelos constitutivos, etc. y permite analizar un
problema físico completo.
d. Ninguna de las anteriores es correcta.

29
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

12) Los elementos estructurales en FLAC3D:

a. Son muy poco útiles.
b. Son de dos tipos, elementos lineales de 2 nodos y elementos triangulares de 3 nodos.
c. Sirven solo para modelizar elementos lineales como vigas, pilotes y anclajes.
d. Únicamente sirven para modelizar elementos superficiales como placas y geomallas.
13) La diferencia fundamental entre grupo y slot radica en:
a. Son conceptos equivalentes.
b. Grupo es un concepto superior a slot. Una misma zona puede pertenecer a slots
distintos en grupos diferentes.
c. Slot es un concepto superior a grupo. Una misma zona puede pertenecer a grupos
distintos en slots diferentes.
d. Grupo hacer referencia a un conjunto de zonas y slot a un dominio espacial.
14) FLAC3D calcula las tensiones y las deformaciones:
a. En los nodos.
b. En las zonas.
c. No los calcula.
d. En nodos y en zonas.
15) Una de las principales ventajas de FLAC3D frente a los elementos finitos es:
a. Más eficiente a la hora de calcular problemas con inestabilidad física y/o con grandes
desplazamientos.
b. No tienen ventajas. Son métodos equivalentes.
c. Más eficiente para calcular problemas lineales de comportamiento estable.
d. Su esquema de resolución implícito.

30
Curso de Introducción a la Modelización con FLAC3D
Modulo 1: Introducción a FLAC3D

7 REFERENCIAS BIBLIOGRÁFICAS
[1] FLAC3D 6.0 – FLAC3D MODELING (an excerpt from FLAC3D Help). Itasca Consulting
Group. © 2017.
[2] FLAC3D 6.0 – Theory and Background (an excerpt from FLAC3D Help). Itasca Consulting
Group. © 2017.
[3] Boidy, Eric. PhD Thesis (Tesis Doctoral). Modelisation numérique du comportement différé
des cavités souterraines. Fecha de lectura: 12 de abril 2002. Páginas 95 a 102.
[4] Debernardi, Daniele. PhD Thesis (Tesis Doctoral). Viscoplastic Behaviour and Design of
Tunnels. Mayo de 2008. Páginas 201 a 206.

31