MÉTODOS DE DOMINIO

Discretiza el interior del macizo en elementos geométricamente simples a los que se les asigna

las propiedades del macizo rocoso. Dentro de estos métodos están: métodos de elementos

finitos (MEF), métodos de las diferencias finitas (MDF), ambos se consideran el macizo rocoso

como un continuo, también pertenece el método de elementos distintos (MED).

Tenemos entre ellos:

1. METODO DE LOS ELEMENTOS FINITOS

El MEF es el más usado y popular método numérico en las ciencias de la ingeniería. Su

capacidad para resolver complejos problemas tenso-deformacionales ha hecho que sea

atractivo en la comunidad de ingenieros de diversas ramas de la ingeniería. La aplicación del

MEF en el diseño de proyectos de ingeniería civil y minería es principalmente un análisis de

tensiones que conduce a la estimación de tensiones, deformaciones y desplazamientos

inducidos por los procesos de excavación. Esta información proporciona importante guía hacia

la evaluación ingenieril de la estabilidad, seguridad y economía del proyecto.

Programas para el Método de Elementos Finitos:

RS2 y RS3 (www.rocscience.com)

PLAXIS 2D y PLAXIS 3D (www.plaxis.com)

ABAQUS www.3ds.com/es/productos-y-servicios/simulia/productos/abaqus

MIDAS GTS www.simulsoft-ingenieros.es/nuestras-soluciones/midas-gts-nx

Tipos de elementos:

Así, por ejemplo, un elemento línea con 2 ó 3 nodos representa una barra sin rigidez a flexión

en 1, 2 ó 3 dimensiones, con 1, 2 ó 3 grados de libertad por nodo respectivamente, o bien con

rigidez, en cuyo caso aparece un g.d.l. rotacional adicional por nodo si se trata de dos

dimensiones, o tres si se trata de 3. Análogamente, triángulos o cuadriláteros pueden

representar planos en 3D.

Matriz de
rigidez

La matriz de rigidez de los elementos, Ke., relaciona las fuerzas que actúan sobre los nodos con

los desplazamientos que producen, mediante la fórmula.

Consideraciones geotécnicas
El M.E.F. se puede utilizar en numerosas aplicaciones geotécnicas como método de cálculo en

la aplicación de distintos métodos de tratamiento de los problemas en que intervienen

tensiones y deformaciones.

 Método de las tensiones efectivas: Muchos de los análisis geotécnicos se basan en el

principio de las tensiones efectivas. Estas técnicas se aplican a procesos

independientes del tiempo.

 Cargas: Pueden ser fuerzas aplicadas en nodos, fuerzas en cuerpos o tracciones en

superficies. Las primeras se tratan directamente, y las dos últimas deben convertirse
 en fuerzas nodales equivalentes. Alternativamente los datos de partida pueden ser los

desplazamientos en vez de las cargas.

 Tensiones iniciales: Normalmente, ya existe en el terreno un estado tensional antes de

aplicar una carga o llevar a cabo una excavación.

Formulación del MEF

La deformación de la estructura subterránea se modeliza mediante una formulación

Lagrangiana de las ecuaciones de momentos y constitutivas de los elementos de la estructura.

Dicha formulación es dependiente del tiempo. Para mayor claridad, se expondrá el caso

bidimensional. Las ecuaciones de momentos de la estructura se pueden expresar según:

PROCESO DE SOLUCIÓN NUMÉRICA

A la hora de realizar el cálculo, las ecuaciones diferenciales (1) se aproximan por ecuaciones en

diferencias aplicadas a la malla definida en cada caso. Las magnitudes vectoriales, como

posición, velocidad y aceleración, se definen en los nodos de la red. Las magnitudes escalares,

como tensiones, permeabilidad, porosidad y presión en los poros, se definen en el centroide

de cada elemento. Un análisis tipo comienza normalmente con la aplicación de una carga

externa a un contorno del dominio de la solución. Siguen una serie de ciclos de deformación

hasta que el sistema alcanza el equilibrio estático. En esta etapa, se tratan todas las zonas
como elementos lagrangianas, es decir, la masa de cada zona permanece constante, mientras

que cambian el volumen, forma y posición.

Aplicaciones del Método de Elementos Finitos – MEF

Cuando es aplicado en la ingeniería de rocas, el MEF consiste en discretizar el macizo rocoso

que rodea una excavación en pequeños elementos conectados a través de nodos. En el

proceso se calculan los desplazamientos de los nodos y luego las tensiones en el interior de los

elementos.

Se aplica en Análisis tenso-deformacionales de automóviles, aeronaves, edificios, estructuras

de puentes, campos de flujos de calor, fluidos, magnéticos, filtraciones.

En Ingeniería Geotécnica ha sido aplicado con gran éxito en problemas de Mecánica de Suelos,

Mecánica de Rocas e Ingeniería de Rocas.

Implementación del Método de Elementos Finitos – MEF

•Se inicia con el análisis de los objetivos.

 Primero, es importante saber lo que se espera obtener como resultado de una

modelación numérica y cuáles son sus limitaciones.

 Segundo, se requiere definir y juntar los datos de ingreso necesarios para el uso de un

programa computacional.

 Tercero, la información obtenida debe ser analizada con propósitos y criterio

ingenieriles.

Visión general del MEF

En el análisis bidimensional, la estructura bajo consideración es divida en elementos que

pueden ser triangulares, rectangulares o regiones de cuatro lados, los que pueden ser

rectos o curvos. Se asume que los elementos están conectados en ciertos puntos llamados
 nodos. Los nodos pueden estar localizados en los lados o en el interior de los elementos. A

los elementos que tienen nodos en los lados, así como en sus vértices se les denomina

elementos parabólicos. Los elementos parabólicos se presentan como los más populares y

probablemente son los que menor tiempo computacional ocupan a la hora del proceso

computacional .Las ecuaciones que describen la forma de los elementos son conocidas

como funciones de forma La unión de varios elementos constituye la malla de elementos

finitos .La configuración de la malla se realiza de forma arbitraria en dependencia de las

necesidades del usuario, aunque en la práctica, se observa que mientras mayor es el

número de elementos, mayor podrán ser el nivel de exactitud

Estructura del Método de Elementos Finitos – MEF

 En el análisis bidimensional, la estructura es dividida en elementos que pueden ser

triangulares, rectangulares o regiones de cuatro lados, que pueden ser rectos o curvos.

 Se asume que los elementos están conectados en ciertos puntos llamados NODOS.

Aquí se calculan los desplazamientos.

 Los nodos pueden estar localizados en los lados o en el interior de los elementos. En el

interior de los elementos se calculan las tensiones.

 A los elementos que tienen nodos en los lados, así como en sus vértices se les

denominan ELEMENTOS PARABOLICOS, y ocupan menor tiempo computacional.


  La unión de varios elementos parabólicos constituye la MALLA DE ELEMENTOS

FINITOS.

 •Mientras mayor es el número de elementos parabólicos, mayor será el nivel de

exactitud.

 La malla no requiere ser uniforme, es decir, los elementos no requieren ser de la

misma forma o dimensión.

 La única regla por seguir es que donde se esperen concentraciones criticas de

tensiones, los elementos sean más pequeños, por ejemplo, la base de un talud o la

periferia de una excavación subterránea.

Figura
1 Malla de elementos finitos de un túnel a gran profundidad. a) Malla de elementos triangulares. b) Malla de
elementos de 4 lados ( Rocksciencie, 2002)..

Pasos para el análisis del método de los elementos finitos

1. Discretización o modelado de la estructura: La estructura es dividida en una cantidad finita

de elementos, con ayuda de un preprocesador. Este paso es uno de los más cruciales para

obtener una solución exacta del problema, de esta forma, determinar el tamaño o la cantidad

de elementos en cierta área o volumen del elemento a analizar representa una ventaja del
método, pero a la vez implica que el usuario debe estar muy conciente de esto para no generar

cálculos innecesarios o soluciones erróneas.

2. Definir las propiedades del elemento: En este paso el usuario debe definir las propiedades

del elemento.

3. Ensamblar las matrices de rigidez de los elementos: La matriz de rigidez de un elemento,

consiste de coeficientes los cuales pueden ser derivados del equilibrio, residuos ponderados o

métodos de energía. La matriz de rigidez del elemento se refiere a los desplazamientos

nodales al ser aplicadas fuerzas en los nodos. El ensamble de las matrices de rigidez, implica la

aplicación de equilibrio para toda la estructura.

4. Aplicación de las cargas: Fuerzas externas concentradas o fuerzas uniformes y momentos

son especificados en este paso.

5. Definir las condiciones de frontera: Las condiciones de apoyo deben ser dadas, por ejemplo,

si el desplazamiento de ciertos nodos es conocido. Usando los elementos de la frontera se

pueden determinar las reacciones en los mismos.

6. Solucionar el sistema de ecuaciones algebraicas lineales: La secuencial aplicación de los

pasos descritos, conduce a un sistema de ecuaciones algebraicas simultáneas, donde los

desplazamientos nodales son desconocidos.

7. Calcular los esfuerzos: El usuario puede entonces calcular los esfuerzos, reacciones,

deformaciones u otra información relevante. El post-procesador ayuda a visualizar la salida en

forma gráfica.

APLICACIÓN DEL METODO DE ELEMENTOS FINITOS EN LA

DETERMINACION DEL SOSTENIMIENTO EN UNA EXCAVACION
SUBTERRANEA
El modelo representa un túnel en forma de herradura, excavado en roca altamente diaclasada

y de mala calidad, descrita como fracturada en bloques y que requiere de sostenimiento para

evitar el colpaso. Se asume que el campo de tensiones es constante antes de la excavación y

que las características de la roca son:

Tabla 1 Características del material

OBJETIVO:

Evaluar la deformación de la excavación con y sin sostenimiento.

DATOS INGRESADOS: Datos del material, geometría de la excavación y condiciones de
contorno.

DATOS DE SALIDA: Campo de tensiones, campo de deformaciones y campo de

desplazamientos.

Figura 2 Mallas de elementos finitos. a) Modelo de excavaciónsin sostenimiento. b) Modelo de excavación con
sostemiento de hormigon lanzado, reforzado con fibras de acero y pernos de anclaje

RESULTADO APLICANDO PHASE 2

Figura 3 Contorno de vectores desplazamiento total. a) Excavación sin sostemiento. b) Excavación con sostemiento
de pernos de anclaje y hormigón lanzado reforzado con fibras de acero. (Rockscience, 2002)