Modelado numérico aplicado a minas subterráneas y a cielo abierto

(Parte 3 de 4)

Carlos Carranza-Torres, Ph.D., P.E.,

Professor of Geotechnical Engineering,
Department of Civil Engineering,
University of Minnesota, Duluth Campus.
Minnesota, USA
E-mail: [email protected]
Clasificación de los métodos numéricos rigurosos

1. Finite Element Method (FEM)

2. Finite Difference Method (FDM)
3. Boundary Element Method (BEM)
4. Discrete Element Method (DEM)

From, Brady B.H.G. and E.T. Brown, 2004, ‘Rock Mechanics for
Underground Mining’, 3rd Edition, Kluwer Academic Publishers.
Algunos softwares comerciales utilizados en mecánica de rocas

1. Finite Element Method (FEM)

→Phase2/RS2 (www.rocscience.com); PLAXIS (www.plaxis.com);
ABAQUS (Dassault Systems)
2. Finite Difference Method (FDM)
→FLAC/FLAC3D (www.itascacg.com)
3. Boundary Element Method (BEM)
→Examine2D/3D (www.rocscience.com); MAP3D (www.map3d.com)
4. Discrete Element Method (DEM)
→UDEC/3DEC/PFC/PFC3D (www.itascacg.com)
Clasificación de los métodos numéricos rigurosos

1. Finite Element Method (FEM)

2. Finite Difference Method (FDM)
3. Boundary Element Method (BEM)
4. Discrete Element Method (DEM)

From, Brady B.H.G. and E.T. Brown, 2004, ‘Rock Mechanics for
Underground Mining’, 3rd Edition, Kluwer Academic Publishers.
Finite Element Method (FEM)

From Zienkiewicz O.C. and R.L. Taylor, 2000, ‘The Finite Element Method’,
Volume I: The Basis. 5th Edition. Butterworth-Heinemann
Finite Element Method (FEM) - continuación
Algunos softwares comerciales utilizados en mecánica de rocas

1. Finite Element Method (FEM)

→Phase2/RS2/RS3 (www.rocscience.com); PLAXIS (www.plaxis.com);
ABAQUS (Dassault Systems)

Modelo Phase2 (RS2) de sección de túnel.

Modelo RS3 del avance de un túnel.

Finite Element Method (FEM). Características básicas

-El método requiere la discretización del medio en una ‘malla de elementos’ delimitado
por nodos.

-Las ecuaciones que gobiernan el problema mecánico (equilibrio, compatibilidad de

deformaciones y ley constitutiva del material) se expresan en forma matricial para
cada elemento (por ej. utilizando consideraciones de ‘trabajo mecánico virtual’) y se
ensamblan en una matriz global asociada a un sistema de ecuaciones linear global,
que cuando se resuelve, permite la computación de las incógnitas en los nodos de la
malla de elementos (las fuerzas y desplazamiento nodales). A partir de las cantidades
nodales mencionadas, se computan los componentes de los tensores de esfuerzo y
deformación para los distintos elementos (o nodos) de la malla.
Finite Element Method (FEM). Características básicas (continuación)

-El método permite incorporación de modelos de material no-lineares (p. ej. elasto-
plástico), a partir de un proceso iterativo que requiere resolver el sistema de
ecuaciones global repetidamente, hasta que la ley constitutiva del material se satisface
en todos los elementos/nodos de la malla.

-El método permite la incorporación de un número contado de planos de

discontinuidad (para representar diaclasamiento o fallas), pero (el método) no está
preparado para simular grandes desplazamientos a través de estos planos de
discontinuidad, o creación de nuevos contactos entre zonas limitadas por
discontinuidades, etc.

-El método permite realizar simulaciones estáticas y dinámicas (p. ej. para simular el
efecto de sismos en excavaciones).

-El método permite el cómputo de factores de seguridad (FS) para taludes,

implementando la técnica de reducción de la resistencia.
Técnica de reducción de la resistencia para el cómputo del factor de seguridad con FEM
Clasificación de los métodos numéricos rigurosos

1. Finite Element Method (FEM)

2. Finite Difference Method (FDM)
3. Boundary Element Method (BEM)
4. Discrete Element Method (DEM)

From, Brady B.H.G. and E.T. Brown, 2004, ‘Rock Mechanics for
Underground Mining’, 3rd Edition, Kluwer Academic Publishers.
Algunos softwares comerciales utilizados en mecánica de rocas

2. Finite Difference Method (FDM)

→FLAC/FLAC3D (www.itascacg.com)

Modelo FLAC para simular excavación Modelo FLAC3D para estudiar la estabilidad del
de sección de caverna en etapas, y verificar fondo de un rajo en la transición de minado de
sostenimiento. rajo a subterráneo.
Finite Difference Method (FDM). Características básicas

-También conocido como ‘Dynamic Relaxation Method’ (o ‘Método de Relajación

Dinámica’).

-Como en el caso del método de Elementos Finitos, el método de Diferencias Finitas

requiere discretizar el medio en una malla de elementos delimitados por nodos.

-Como en el caso del método de Elementos Discretos (as ser discutido

posteriormente), las ecuaciones de momento dinámico se aplican y se resuelven en
nodos y elementos utilizando un algoritmo de ‘marcha en el tiempo’ (integración tipo
Euler), con un esquema de amortiguamiento que tiene el propósito de que el sistema
converja a una solución pseudo-estática.

-El método no requiere ensamblar o resolver un sistema de ecuaciones, como en el

caso del método de los Elementos Finitos. En contraste, el método requiere la
integración de las ecuaciones en forma iterativa (como se mencionara en el punto
anterior).
Finite Difference Method (FDM). Características básicas (continuación)

-Modelos de material no-lineares (p. ej., modelos elasto-plásticos dilatantes) pueden

ser incorporados fácilmente con el método.

-Como en el caso de los Elementos Finitos, el método permite la incorporación de un

número contado de planos de discontinuidad (para representar diaclasamiento o
fallas), pero no está preparado para simular grandes desplazamientos a través de
estos planos de discontinuidad, o creación de nuevos contactos entre zonas limitadas
por discontinuidades, etc.

-El método permite realizar simulaciones estáticas y dinámicas (p. ej. para simular el
efecto de sismos en excavaciones).

-El método permite el cómputo de factores de seguridad (FS) para taludes,

implementando la técnica de reducción de la resistencia.
Análisis continuo (FLAC) de una sección de deformación plana de una caverna
superficial excavada en diversas etapas en roca competente (el material ha sido
representado con modelo constitutivo tipo Hoek-Brown).
Representación de contornos de magnitud de desplazamientos resultantes
en el modelo continuo (FLAC).
Modelo 3D complejo de una mina a rajo
abierto y su transición a subterránea. El
modelo se ha armado con un ‘mallador’
comercial de bajo costo, y ha sido
importado a FLAC3D (de Itasca) para el
cómputo. El objetivo del modelo ha sido
estudiar la concentración de esfuerzos en
el piso del rajo (y techo de los
caserones).
 Etapa de Excavación 1 de 5

Modelo 3D complejo de una mina a rajo

abierto y su transición a subterránea. El
modelo se ha armado con un ‘mallador’
comercial de bajo costo, y ha sido
importado a FLAC3D (de Itasca) para el
cómputo. El objetivo del modelo ha sido
estudiar la concentración de esfuerzos en
el piso del rajo (y techo de los
caserones).
 Etapa de Excavación 2 de 5

Modelo 3D complejo de una mina a rajo

abierto y su transición a subterránea. El
modelo se ha armado con un ‘mallador’
comercial de bajo costo, y ha sido
importado a FLAC3D (de Itasca) para el
cómputo. El objetivo del modelo ha sido
estudiar la concentración de esfuerzos en
el piso del rajo (y techo de los
caserones).
 Etapa de Excavación 3 de 5

Modelo 3D complejo de una mina a rajo

abierto y su transición a subterránea. El
modelo se ha armado con un ‘mallador’
comercial de bajo costo, y ha sido
importado a FLAC3D (de Itasca) para el
cómputo. El objetivo del modelo ha sido
estudiar la concentración de esfuerzos en
el piso del rajo (y techo de los
caserones).
 Etapa de Excavación 4 de 5

Modelo 3D complejo de una mina a rajo

abierto y su transición a subterránea. El
modelo se ha armado con un ‘mallador’
comercial de bajo costo, y ha sido
importado a FLAC3D (de Itasca) para el
cómputo. El objetivo del modelo ha sido
estudiar la concentración de esfuerzos en
el piso del rajo (y techo de los
caserones).
 Etapa de Excavación 5 de 5

Modelo 3D complejo de una mina a rajo

abierto y su transición a subterránea. El
modelo se ha armado con un ‘mallador’
comercial de bajo costo, y ha sido
importado a FLAC3D (de Itasca) para el
cómputo. El objetivo del modelo ha sido
estudiar la concentración de esfuerzos en
el piso del rajo (y techo de los
caserones).
Modelo 3D complejo de una presa de
enrocado con pantalla de concreto
generado con un mallador externo e
importado a FLAC3D (de Itasca). El
objetivo del modelo ha sido estudiar los
asentamientos del enrocado y el
comportamiento sísmico de la presa –el
modelo ha sido sujeto a un análisis
dinámico basado en un sismograma de
entrada.
Etapa de Construcción 1 de 4 Etapa de Construcción 2 de 4

Etapa de Construcción 3 de 4 Etapa de Construcción 4 de 4

Etapa de Construcción 1 de 4 Etapa de Construcción 2 de 4

Etapa de Construcción 3 de 4 Etapa de Construcción 4 de 4

Modelo 3D complejo de otra caverna de chancado
generado con un mallador externo e importado a
FLAC3D (de Itasca). El objetivo del modelo ha sido
estudiar el comportamiento del revestimiento de
concreto para la caverna.