DIPLOMADO EN GEOMECÁNICA SUBTERRÁNEA Y

SUPERFICIAL

INFORME TÉCNICO DE SUSTENTACIÓN DE CURSO

VIRTUAL OFFLINE (ITSCVO)

CURSO
SOFTWARE PHASE2

ESTUDIO

ANÁLISIS DE TENSIONES Y DEFORMACIONES POR EL

MÉTODO DE ELEMENTOS FINITOS (FEM) DEL MACIZO
ROCOSO EN EL TAJEO ESPERANZA Nv. 920

ALUMNA: KADIMA MEJIA ZANABRIA

DOCENTE: ING. GUILLERMO RODRIGUEZ CAYLLAHUA

AREQUIPA – PERÚ, MARZO 2019

 RESUMEN EJECUTIVO
El siguiente informe trata sobre el análisis de las tensiones-deformaciones in situ
e inducidos sobre el macizo rocoso y la evaluación de un Factor de seguridad
para el tajeo Esperanza Nv. 920. El proyecto se encuentra ubicado a 9 km al
sudeste de Tintaya en la Provincia de Espinar, región Cusco, entre 4,030 a 4,170
msnm.

Se va a evaluar el proyecto Tajeo Esperanza Nv. 920, con información

proporcionada por el área de Geología (litología) y data proporcionada por el
Laboratorio de Geotécnia de las propiedades elásticas y geomecánicas de cada
unidad de geotécnica (Pared colgante, pared pendiente, roca mineral y relleno
hidráulico) y por otra parte se cuenta con dimensiones de las excavaciones a
realizar, excavaciones existentes. Se requiere saber cómo se comportará el
macizo rocoso frente a etapas de excavación teniendo como base un Factor de
seguridad de 1.2 y estimar la altura mínima del pilar que cumpla con dicho factor
de seguridad del diseño.

El análisis se realizó con ayuda del software Phase2, basado en el Método de

Elementos Finitos que nos permite dividir un medio continuo en un número finito
de partes.

Palabras claves: Tensiones in situ, inducidos, factor de seguridad, Phase2,

Método de elementos finitos, pilar, medio continuo.

2
 INDICE

RESUMEN EJECUTIVO .................................................................................... 2

1. INTRODUCCION ......................................................................................... 8
2. METODOLOGÍA DE ESTUDIOS................................................................. 9
2.1 PLANTEAMIENTO DEL PROBLEMA .......................................................... 9
2.2 OBJETIVOS ................................................................................................. 9
2.2.1 OBJETIVO GENERAL ........................................................................... 9
2.2.2 OBJETIVO ESPECIFICO .................................................................. 9
2.3 ALCANCES ................................................................................................ 10
2.4 UBICACIÓN DEL PROYECTO .................................................................. 10
2.5. ASPECTOS GEOLOGICOS ..................................................................... 11
3. MARCO TEÓRICO ...................................................................................... 12
3.1 MÉTODO DE ELEMENTOS FINITOS .................................................... 12
3.2 ESTADO TENSIONAL DE MACIZO ROCOSO ...................................... 15
4. INVESTIGACIONES BÁSICAS .................................................................... 22
4.1 DATOS DE LA EXCAVACIÓN SUBTERRÁNEA .................................... 22
4.2 METODOLOGIA DE EXPLOTACIÓN DEL CUERPO ............................. 22
4.3 PERFIL TRANSVERSAL DE LA EXCAVACIÓN .................................... 22
4.4 PROPIEDADES ELÁSTICAS Y GEOMECÁNICAS ................................ 24
4.5 ESTADO TENSIONAL IN SITU .............................................................. 25
5. METODOLOGÍA DE ANÁLISIS ................................................................... 26
5.1. MÓDULO DE IMPORTACIÓN ............................................................... 26
5.2 CONFIGURACIÓN DEL PROYECTO .................................................... 27
5.3. ESTADO TENSIONAL ........................................................................... 31
5.4. PROPIEDADES DE LAS UNIDADES GEOTÉCNICAS ......................... 32
5.5. CREACIÓN DE LA MALLA EN ELEMENTOS FINITOS. ....................... 36
5.6. MODULO DE CÁLCULO ....................................................................... 38
5.8. MÓDULO DE INTERPRETACIÓN......................................................... 38
6. RESULTADOS ............................................................................................. 49
7. CONCLUSIONES Y RECOMENDACIONES ............................................... 50
BIBLIOGRAFÍA ................................................................................................ 51

3
ANEXOS .......................................................................................................... 52

INDICE DE FIGURAS
Figura 1. Mapa de ubicación del Proyecto Tajeo Esperanza…………………...10

Figura 2. Dominio y su Contorno – subdivisión en elementos…………………..14

Figura 3. Esquema del método de elementos finitos…………………………….14

Figura 4. Ejemplo de un elemento discretizado en un número finito de

elementos……………………………………………………………………………..15

Figura 5. Estado inicial de esfuerzos in situ en macizo rocoso no

perturbados…………………………………………………………………………...16

Figura 6. Estado posterior de los esfuerzos perturbados por excavaciones

abiertas………………………………………………………………………………..16

Figura 7. Esfuerzos principales inducidos en un elemento de roca cerca de un

túnel horizontal, sujeto a esfuerzos verticales y horizontales…………………...17

Figura 8. Mediciones de esfuerzos verticales en obras mineras – civiles……..18

Figura 9. Relación de la constante “k” y la profundidad………………………….20

Figura 10. Tensión Horizontal media y vertical en función de la profundidad…21

Figura 11. Perfil transversal de la excavación…………………………………….23

Figura 12. Estado Tensional In Situ………………………………………………..25

Figura 13. Opciones de Importación al software Phase2………………………..26

Figura 14. Archivo DXF importado al Phase2 con contornos der colores……..27

Figura 15. Configuración General del proyecto en Phase2……………………..28

Figura 16. Configuración de las etapas del proyecto en el software Phase2…29

Figura 17. Configuración del análisis de esfuerzos………………………………30

Figura 18. Resumen del proyecto………………………………………………….31

Figura 19. Introducción de los valores del campo de esfuerzos………………..31

4
Figura 20. Propiedades del material que presenta la pared colgante………….32

Figura 21. Ingresando parámetros del macizo rocoso de la pared colgante….32

Figura 22. Propiedades del material que presenta la pared pendiente………..33

Figura 23. Calculador de parámetros de macizo de la pared pendiente……....33

Figura 24. Propiedades del material que presenta el mineral…………………..34

Figura 25. Calculador de parámetros de macizo rocoso del mineral…………..34

Figura 26. Propiedades del material que presenta el relleno hidráulico……….35

Figura 27. Definición de las propiedades de la junta…………………………….35

Figura 28. Ingreso de datos para la generación de la malla…………………….36

Figura 29. Proyecto discretizado en elementos finitos…………………………..37

Figura 30. Computarización y cálculo de los elementos finitos…………………38

Figura 31. Valores de Sigma 1 corte 1, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….39

Figura 32. Factor De Seguridad Corte 1 en contornos de la excavación y

puente o pilar…………………………………………………………………………39

Figura 33. Valores de Sigma 1 corte 2, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….40

Figura 34. Factor De Seguridad Corte 2 en contornos de la excavación y

puente o pilar…………………………………………………………………………40

Figura 35. Valores de Sigma 1 corte 3, en contornos de las excavaciones…...41

Figura 36. Factor De Seguridad Corte 3 en contornos de la excavación y

puente o pilar…………………………………………………………………………41

Figura 37. Valores de Sigma 1 corte 4, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….42

Figura 38. Factor De Seguridad Corte 4 en contornos de la excavación y

puente o pilar…………………………………………………………………………42

Figura 39. Valores de Sigma 1 corte 5, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….43

5
Figura 40. Factor de Seguridad Corte 5 en contornos de la excavación y puente
o pilar…………………………………………………………………………………..43

Figura 41. Valores de Sigma 1 corte 6, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….44

Figura 42. Factor de Seguridad Corte 6 en contornos de la excavación y

puente………………………………………………………………………………....44

Figura 43. Valores de Sigma 1 corte 7, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….45

Figura 44. Factor de Seguridad Corte 7 en contornos de la excavación y puente

o pilar…………………………………………………………………………….........45

Figura 45. Valores de Sigma 1 corte 8, en contornos de las excavaciones y

Pilar…………………………………………………………………………………….46

Figura 46. Factor de Seguridad Corte 8 en contornos de la excavación y puente

o pilar…………………………………………………………………………………..46

Figura 47. Modelamiento de sigma 1 con relleno de hidráulico los pilares

explotados…………………………………………………………………………….47

Figura 48. Modelamiento del factor de seguridad con el relleno hidráulico en los
pilares explotados……………………………………………………………………47

Figura 49. Modelamiento de sigma 1 con relleno de hidráulico de todas las

excavaciones…………………………………………………………………………48

Figura 50. Modelamiento del factor de seguridad con el relleno hidráulico en

todas excavaciones del proyecto…………………………………………………..48

6
INDICE DE TABLAS
Tabla 1. Etapas de excavación, no incluye las etapas de relleno………………23

Tabla 2: Propiedades Geomecánicas y Elásticas de las paredes y roca

mineral………………………………………………………………………………...24

Tabla 3: Propiedades Geomecánicas y Elásticas del relleno hidráulico.………24

Tabla 4. Propiedades de las juntas asociadas al contacto entre la roca

encajonante y la zona mineralizada………………………………………………..24

Tabla 5. Resultados promedio del FS y Sigma 1 durante las diferentes

secuencias de minado y rellenos…………………………………………………..49

7
1. INTRODUCCION

Para garantizar la estabilidad de las excavaciones subterráneas la cual está

ligada con la seguridad del personal y de las maquinas que trabaja en estas
excavaciones, he ahí la importancia del estudio del macizo rocoso con el objeto
de conocer sus propiedades geomecánicas y asegurar su estabilidad; es por
ello la necesidad del uso y conocimiento del software Phase2.

Al momento de iniciar un diseño, es necesario tener suficientes antecedentes

geotécnicos y geo mecánicos para ser capaz de conocer, con la mayor precisión
posible, el comportamiento del macizo rocoso, además de elementos naturales
encontrados dentro y alrededor de una excavación planeada.

En el siguiente informe se presenta un análisis de las tensiones y deformaciones

que se desarrollan en el macizo rocoso, producidos por las labores de
explotación de corte y relleno en el Tajeo Esperanza Nv. 920. Con el fin, de
estimar una altura adecuada para el pilar, cumpliendo un factor de seguridad de
1.2; para lo cual se realizó el análisis por el método de elementos finitos.

El uso del software Phase2 por ende es una herramienta necesaria e importante,
ya que gracias a ella se tiene la facilidad de observar de manera casi exacta el
comportamiento de los esfuerzos presentes en el macizo y poder determinar el
Factor de seguridad, además como tener conocimiento de diversos escenarios
posibles relacionados a la estabilidad del macizo rocoso.

Este análisis tiene su base en el Método de Elementos Finitos en donde el

macizo rocoso es discretizado en una cantidad finita de elementos, integrando
las ecuaciones constitutivas relacionadas con la tensión-deformación.

8
2. METODOLOGÍA DE ESTUDIOS

2.1 PLANTEAMIENTO DEL PROBLEMA

En el proyecto Tajeo Esperanza Nv. 920, inicialmente en la parte superior del

cuerpo mineral se cuenta con una cámara previamente explotada de 35 m de
altura, lo cual al momento de realizar el tajeo de las diferentes etapas, el macizo
rocoso estará sometida a una redistribución de esfuerzos a medida que la altura
del pilar disminuya, esto nos lleva a determinar una altura mínima de pilar que
garantice la seguridad del personal y equipos dentro del proyecto.

¿Cuál es el Factor de Seguridad de la excavación subterránea en sus distintas

etapas?

¿Cuáles son los esfuerzos in-situ e inducidos de la excavación subterránea?

2.2 OBJETIVOS

2.2.1 OBJETIVO GENERAL

 Analizar la distribución de esfuerzos sigma 1 y determinar el Factor de

seguridad del pilar del Tajeo Esperanza Nv. 920 en las diversas etapas
de las excavaciones para cumplir con el Factor de Seguridad 1.2 de
diseño.

2.2.2 OBJETIVO ESPECIFICO

 Analizar en qué etapas la excavación se encuentra estable el pilar y en
cuáles no.
 Estimar la mínima altura del pilar entre excavaciones que cumpla el factor
de seguridad de diseño de 1.2
 Determinar la cantidad de fases óptimas de minado para mantener el pilar
dentro del factor de seguridad requerido.
 Conocer la importancia y manejo adecuado del Software Phase2.

9
2.3 ALCANCES

En el presente informe se realizará el estudio y análisis de estabilidad en función

del factor de seguridad y esfuerzo principal mayor Sigma 1 para el Tajeo
Esperanza Nv. 920, para el desarrollo de éste análisis se empleará la información
brindada por el CGI.

2.4 UBICACIÓN DEL PROYECTO

El proyecto Tajo Esperanza Nv. 920 del CGI se encuentra ubicado a 9 km al

suroeste de Tintaya en la Provincia de Espinar, región Cusco, entre 4,030 a
4,170 msnm.

TAJEO ESPERANZA

PLANO DE UBICACIÓN :
TAJEO ESPERANZA
CGI

Figura 1. Mapa de ubicación del Proyecto Tajeo Esperanza

10
2.5. ASPECTOS GEOLOGICOS

El Proyecto Tajeo Esperanza de tipo Skarn-Pórfido se encuentra dentro de la

Franja Eocena-Oligocena del cinturón Andahuaylas-Yauri. Este cinturón alberga
varios proyectos importantes como Las Bambas, Antapaccay, Haquira, Quechua
y Constancia. Tajeo Esperanza es un depósito asociado a fases de intrusión con
mineralización de Cu-Au-Ag.

La mayor mineralización de sulfuros de cobre se encuentra distribuida en los

cuerpos de skarn de granate-piroxenomagnetita, con mineralización de
calcopirita y bornita, que se encuentran a 250 m. debajo de la superficie en forma
de un potente manto irregular, ubicado entre el hornfels (piso) y la diorita (techo)
inclinado hacia el SW intruído por un dique y stock de pórfido Monzonita que se
incrementa en profundidad.

La que trajo la mineralización es el pórfido monzonita hornbléndica, se presenta

generalmente con venillas de cuarzo tipo stock work débil con mineralización de
calcopirita de baja ley, la diorita también es mineralizada en contacto con el
pórfido Monzonita hornbléndica con algunas concentraciones de leyes de Cu, se
tiene otro pórfido de monzonita biotítica que es post mineral.

Todavía falta evaluar más este yacimiento para convertirlo de open pit a
subterráneo, para que sea viable su explotación a gran escala.

11
3. MARCO TEÓRICO

Para realizar el diseño de una excavación, a niveles superficiales o subterráneos,

es necesario contar con un estudio detallado desde el ámbito de la ingeniería
que permita realizar el análisis apropiado de las múltiples variables que
condicionan la estabilidad del suelo o del macizo rocoso sobre la cual se
pretende construir una nueva excavación, esto con el fin de diseñar
excavaciones o secuencias de minado en óptimas condiciones que reduzcan al
máximo los riesgos, garantizando siempre la seguridad no solo del personal y
maquinarias que labora en obra sino también de la eficacia y eficiencia en los
avances de la operación.

A través del tiempo se han definido varios métodos numéricos y ecuaciones

constitutivas que han permitido realizar el análisis de toda esta información
arrojando resultados que permiten sintetizar y minimizar la incertidumbre frente
a la estabilidad de una excavación permitiendo a los Ingenieros de las compañías
Mineras, Centrales Hidroeléctricas y Civiles, entre otros campos, diseñar labores
de múltiples formas, tamaños y geometrías enmarcados en los mejores
estándares de calidad; dentro de los métodos más usados y aplicados en la
actualidad se tiene el Método de Elementos Finitos que permite compilar gran
parte de la información mencionada anteriormente y nos lleva a evaluar y definir
la mejor forma de proceder frente al diseño de una excavación en términos de
su estabilidad global y local.

3.1 MÉTODO DE ELEMENTOS FINITOS

Este método constituye un método numérico destinado a resolver mediante

ecuaciones matriciales las ecuaciones diferenciales que se plantean en sistemas
discretos (estructuras) o continuos (campos). Actualmente, se considera al
método de las Diferencias Finitas como una subclase del método de los

12
Elementos Finitos y de hecho se puede demostrar [Silvester‐Chari] que el
método FEM se reduce al método DF cuando las mallas son regulares.

Las aplicaciones actuales del método son muy extensas e incluyen sistemas
lineales y no lineales, estáticos, dinámicos tales como Mecánica de Sólidos,
Teoría de la Elasticidad, Mecánica de Fluidos, Transmisión de Calor y
Electromagnetismo.

En el caso de sistemas continuos, el método consiste en discretizar el dominio

de interés en Elementos Finitos y resolver, mediante una función de prueba o de
aproximación, la ecuación que rige el sistema en cada EF para luego sumar
todas las soluciones.

Dado un recinto cerrado los pasos para la resolución son:

1. Dividir el recinto en Elementos Finitos: Triángulos (3 nodos), Tetraedros

(4 nodos), etc.
2. Deducir la ecuación que describe el potencial f dentro de un EF.
3. Plantear las ecuaciones que dan las condiciones de ajuste de las
soluciones en las fronteras de los EF.
4. Calcular los potenciales en los nodos de cada EF mediante algunos de
los métodos que luego de mencionarán.
5. Resolver las ecuaciones algebraicas planteadas.

En general se puede decir entonces que el método de los elementos finitos

constituye un procedimiento general de discretización de los problemas
continuos planteados por expresiones definidas matemáticamente. En este
sentido entonces dado un recinto cerrado los pasos para su resolución se
resumen a continuación:

13
  El continuo se divide en un número finito de partes (elementos), cuyo
comportamiento se especifica mediante un número finito de parámetros,
incluyendo sus conexiones nodales. (Triángulos, Tetraedros, nodos etc.)
 Se asume que los elementos están conectados entre sí mediante un
número discreto de puntos, que se llaman nodos, situados en sus
contornos. (Zienkiewicz 1982)
 Los desplazamientos de los nodos serán las incógnitas fundamentales del
problema. (Zienkiewicz 1982)
 Se toman un conjunto de funciones que definen el campo de
desplazamientos en los elementos según sea el desplazamiento de los
nodos. (Zienkiewicz 1982)
 Deducir la ecuación que describe el potencial f dentro de un EF. (CGI)
 Plantear las ecuaciones finales que dan las condiciones de ajuste de las
soluciones en las fronteras de los elementos finitos. (CGI)
 Resolver las ecuaciones algebraicas planteadas. (CGI)

Figura 2. Dominio y su Contorno – subdivisión en elementos (Fuente: Zienkiewicz 1982)

Figura 3. Esquema del método de elementos finitos. (Fuente: Sánchez 2016 a partir de Bobet
2010)

14
El Método de Elementos Finitos trata de un método general para la solución de
problemas de contorno gobernados por ecuaciones diferenciales ordinarias o
parciales. En esencia se trata de una técnica que sustituye el problema
diferencial por otro algebraico, aproximadamente equivalente, para el cual se
conocen técnicas generales de resolución. Para ello hace uso de la
"discretización" o subdivisión de una región sobre la cual están definidas las
ecuaciones en formas geométricas simples denominadas elementos finitos. Las
propiedades materiales y relaciones gobernantes en estos elementos se
expresan en función de los valores desconocidos en las "esquinas" de los
elementos o nodos.

PROCESO DE
DISCRETIZACIÓN

SISTEMA CONTINUO MODELO DISCRETO

Figura 4. Ejemplo de un elemento discretizado en un número finito de elementos.

3.2 ESTADO TENSIONAL DE MACIZO ROCOSO

Todas las rocas en profundidad o materiales del subsuelo se encuentran por

naturaleza sometidos a una serie de esfuerzos cuyo origen está ligado en
primera instancia a todo el peso del material o rocas supra yacentes, cargas
litostática o presiones de roca y en segunda instancia sin ser menos importante,
a esfuerzos causados por origen tectónico.

Cuando abrimos una excavación todo el campo de esfuerzos inicial alrededor de

la misma, es disturbado, generándose nuevos esfuerzos, pero inducidos en la

15
roca alrededor de la labor abierta, se pierde el equilibrio inicial que se tenía y los
esfuerzos en general cambian tanto en sus magnitudes como en sus direcciones,
estos últimos valores son desconocidos tanto desde sus valores in situ no
perturbados como los posteriores inducidos.

El conocimiento de estos esfuerzos es de gran importancia al momento de

pretender diseñar alguna excavación subterránea, dado que si estos presentan
valores demasiado elevados constituyen un problema en términos de la
estabilidad para las labores y pueden traer consecuencias graves.

En general estos se pueden denominar esfuerzos in situ y esfuerzos inducidos:

Figura 5. Estado inicial de esfuerzos in situ en macizo rocoso no perturbados. Fuente.

(Sociedad Nacional de Minería Petróleo Energía, 2004)

Figura 6. Estado posterior de los esfuerzos perturbados por excavaciones abiertas. Fuente.
(Sociedad Nacional de Minería Petróleo Energía, 2004).

16
El macizo rocoso de manera natural está sometido a esfuerzos “naturales” o “in
situ”, la importancia de estos esfuerzos y su estimación es porque influencian
especialmente sobre:

 La orientación de la excavación
 En el diseño de las obras
 En las voladuras
 En la forma de la excavación

Además, las tensiones naturales son definidas principalmente por ciertos

factores tales son los siguientes (Carnero, 2014):

 Gravedad
 Efecto de Poisson
 Topografía
 Tectónica de placas
 Estructuras geológicas
 Propiedades de las rocas
 Heterogeneidades
 Vulcanismo
 Subsidencia

Por otro lado, cuando se realiza una excavación se generan los esfuerzos
inducidos, ya que se alterará el equilibrio que existía en el macizo, provocándose
que los esfuerzos inducidos se sumen al esfuerzo natural dando lugar a un
esfuerzo luego de la excavación o un esfuerzo final.

Figura 7. Esfuerzos principales inducidos en un elemento de roca cerca de un túnel horizontal,

sujeto a esfuerzos verticales y horizontales. Fuente: (Hoek, 2007)

17
Tensión Vertical

El esfuerzo in-situ debido al peso de la columna vertical de roca, resulta de hallar

el producto de la profundidad por el peso unitario de la masa de roca sobre
yacente, como muestra la siguiente expresión:

σv= y.z

Dónde:

 σv : es esfuerzo vertical
 γ : es el peso unitario de la roca sobre yacente
 z: es la profundidad debajo de la superficie.

Las mediciones de los esfuerzos verticales en varias obras de ingeniería minera

y civil en el mundo confirman que esta relación es válida, tal como se ilustra en
la siguiente figura, aunque exista una significativa cantidad de dispersión en las
mediciones.

Figura 8. Mediciones de esfuerzos verticales en obras mineras – civiles (H&B 1978)

18
Determinación de la constante “k”

Los esfuerzos horizontales actuantes sobre un elemento de roca a una

profundidad z debajo de la superficie, son mucho más dificultosos de estimar que
los esfuerzos verticales. Normalmente, la relación del esfuerzo horizontal al
esfuerzo vertical promedios es denotada por la letra “k“:

Mediciones de los esfuerzos horizontales en obras civiles y mineras en el mundo,

muestran que la relación “k” tiende a ser alta a poca profundidad y que decrece
con la profundidad (Brown y Hoek, 1978, Herget, 1988).

Sheorey (1994) desarrolló un modelo de esfuerzo termal elasto-estático de la

tierra, mediante la cual proporcionó una ecuación simplificada que puede ser
utilizada para estimar la relación esfuerzo horizontal a esfuerzo vertical. De este
forma “k” puede ser estimado como:

Donde:

 z: es la profundidad expresada en metros debajo de la superficie

 Eh: es el módulo de deformación promedio de la parte alta de la corteza
terrestre medida en una dirección horizontal, y está expresada en GPa.

Un gráfico de esta ecuación está dado en la Figura B, para un rango del módulo
de deformación. Las curvas que relacionan a “k” con la profundidad debajo de la
superficie z son similares a las publicadas por Brown y Hoek (1978), Herget
(1988) y otros para esfuerzos in-situ medidos.

19
 Figura 9. Relación de la constante “k” y la profundidad (Sheorey 1994)

Tensión Horizontal

También existen esfuerzos horizontales que actúan sobre la roca a profundidad

aunque es un poco más complicado de estimar que el esfuerzo vertical, sin
embargo la relación tensión horizontal promedio respecto a la tensión vertical es
denotada por la letra K de tal forma que:

Sin embargo Terzaghi y Richart (1952) sugirieron que para macizos rocosos
gravitacionalmente cargados donde no se permitió ninguna tensión lateral
durante la formación se estratos superiores o suparyacentes el valor de K es
independiente de la profundidad y está dado por:

Donde:

 k: Relación que expresa el esfuerzo vertical y horizontal

 v : Modulo Poisson para la masa rocosa

20
La recopilación de varias lecturas del esfuerzo horizontal a través de varios
proyectos en el mundo muestra que el valor de K tiende a ser alto cerca de
superficie y que baja o decrece en profundidad (Brown and Hoek 1978, Herget
1988).

Figura 10. Tensión Horizontal media y vertical en función de la profundidad. (Fuente: CGI)

21
4. INVESTIGACIONES BÁSICAS

Para el estudio del proyecto Tajeo Esperanza Nv. 920 se cuenta con información
como la profundidad de las labores, las orientaciones que éstas tienen y el factor
de seguridad de diseño de 1.2

Para el análisis de los esfuerzos, deformación y factor de seguridad del macizo

rocoso con la excavación previa y luego de las etapas a realizar, se obtuvo de
estudios realizado por el Centro Geotécnico Internacional CGI.

4.1 DATOS DE LA EXCAVACIÓN SUBTERRÁNEA

 Nombre de la labor: Tajeo Esperanza Nv.920.

 Factor de Seguridad del diseño: 1.2

4.2 METODOLOGIA DE EXPLOTACIÓN DEL CUERPO

 Cuerpo Mineral subvertical con potencia o espesor de 15 m.

 Excavación de forma ascendente (Drift & Fill).
 En la parte superior del cuerpo existe una cámara previamente
explotada, sin relleno con altura de 35 m.
 Después de cada corte no se rellenará la cámara, solo se hará hasta el
final después del último tajeo.
 Se realizará relleno hidráulico. Solo al final.
 La altura de cada corte es variable, ver el diseño propuesto.
 No se aplicará sostenimiento como pernos o shotcrete.

4.3 PERFIL TRANSVERSAL DE LA EXCAVACIÓN

En el proyecto se tendrá un total de 8 fases de excavación o de cortes como se

detalla en la figura 7 cuya altura será variable para los primeros cortes. La altura
del pilar también varía según el número de etapas de excavación realizadas
(Tabla 1). El macizo rocoso encajonante del proyecto es una granodiorita
moderadamente alterada con un peso específico de 2.90 Ton/m3.

22
 Figura 11. Perfil transversal de la excavación (Fuente: CGI).

ETAPAS DE
W
EXCAVACIÓN
1 55
2 40
3 35
4 30
5 25
6 20
7 15
8 10
Tabla 1. Etapas de excavación, no incluye las etapas de relleno (fuente: CGI)

23
4.4 PROPIEDADES ELÁSTICAS Y GEOMECÁNICAS

Las siguientes tablas muestran las propiedades elásticas y geomecánicas de las

paredes y roca mineral (Tabla 2) y las propiedades elásticas y geomecánicas del
relleno hidráulico (Tabla 3).

Propiedades
Propiedades Elásticas
geomecánicas
Unidad Peso
Módulo
geotécnic Especifico Módulo
Criterio de de UCS
a (MN/m3) de Young GSI mi D
falla Poisson [MPa]
E [GPa] n
Pared
0.0258 40 0.25 150 50 20 0.2
colgante
Hoek-Brown
Pared
0.0258 Generalizado 45 0.25 150 48 18 0.2
pendiente
Roca
0.0294 28 0.27 100 45 12 0.2
mineral
Tabla 2: Propiedades Geomecánicas y Elásticas de las paredes y roca mineral (Fuente: CGI).

Unidad Peso Módulo Módulo Tensile Cohesio

Friction angle
geotécnic Especifico Criterio de falla de Young de Strengt n
a (MN/m3) E [GPa] Poisson [MPa] [MPa] [Grados]
Relleno
0.0196 Mohr Coulomb 0.3 0.3 0 32 0
Hidraulico
Tabla 3: Propiedades Geomecánicas y Elásticas del relleno hidráulico (Fuente: CGI).

Tabla 4. Propiedades de las juntas asociadas al contacto entre la roca encajonante y la zona
mineralizada (fuente: CGI)

24
4.5 ESTADO TENSIONAL IN SITU

Para el análisis vamos a considerar que el ángulo de giro respecto a la horizontal

es de 7°, por tanto al procesar los datos en el software se utilizaran de forma
directa.

Figura 12. Estado Tensional In Situ (fuente: CGI)

25
5. METODOLOGÍA DE ANÁLISIS

5.1. MÓDULO DE IMPORTACIÓN

Primero es necesario importar o crear los límites del proyecto; para éste caso se
cuenta con un diseño en AutoCAD proporcionado por el CGI como archivo DXF
lo cual se va importar al software Phase 2.

Cuando se diseñe o construye en el AutoCAD debemos generar capas que

tengan los nombres que se va a exportar al Phase 2, como por ejemplo:

 External: Contorno o limite externo.

 Excavations: Excavaciones del modelo
 Stages: Etapas con los cortes de cada fase.
 Joints: Juntas o contactos.

Figura 13. Opciones de Importación al software Phase2 (fuente: Elaboración Propia)

26
 Figura 14. Archivo DXF importado al Phase2 con contornos der colores (fuente: Elaboración
Propia)

5.2 CONFIGURACIÓN DEL PROYECTO

5.2.1. GENERAL

En este apartado se va asignar el tipo de análisis, tipo de solución y las unidades

tal como se detalla a continuación:

A. Tipo de Análisis: Plane Strain, es el más usado, un modelo de

deformación plana, el asimétrico o Axisymmetric normalmente se usa
para realizar análisis en 3D. Para este proyecto vamos a realizar el
modelado en 2D o Plane Strain.

27
 B. Tipo de Solución: En esta sección consta de varias opciones a través de
las cuales se quiera dar la solución a los cálculos, sin embargo el más
recomendado y menos complicado es el Gaussian Elimination
C. Unidades: En esta sección seleccionaremos las unidades con los cuales
vamos a realizar el proyecto. Trabajaremos en unidades métricas (Metric)
y los esfuerzos en Mega Pascales (Mpa).

Figura 15. Configuración General del proyecto en Phase2 (Fuente. Elaboración propia).

5.2.2. ETAPAS

En esta parte del proyecto se va a crear todas las etapas involucradas en el

proyecto. Para nuestro estudio se crearon 10 etapas (8 fases y 2 rellenos
hidráulicos) tal como se muestran en la Figura16.

28
 Figura 16. Configuración de las etapas del proyecto en el software Phase2. (Fuente:
Elaboración propia.)

5.2.3. ANALISIS DE ESFUERZOS

A. Número máximo de Iteraciones: El valor que nos muestra el software

phase2 por defecto es de 500 y para nuestro proyecto se analizará con
ese valor.
B. Tolerancia: El valor que nos muestra el software es de 0.001. Este valor
es recomendable cuando se trabaja con materiales de tipo elástico ya que
la solución es exacta; y también cuando no existen problemas de
convergencia.
C. Number of Load Steps: El valor que nos muestra el software por defecto
es “Auto” lo cual nos va generar un número de pasos de carga usados en
la etapa de cálculo para cada fase.

29
D. Tipo de convergencia: Se puede escoger entre dos tipos de
convergencia “Absolute Energy” (energía absoluta) y “Square Root
Energy”. Para nuestro proyecto se utilizó el valor configurado por defecto
“Absolute Energy”.
E. Tensile Failure Reduces Shear Strenght to Residual: La configuración
para esta sección se va trabajar por defecto.
F. Joint tension reduces joint stiffness by a factor of: Esta opción se
dejará activada con el valor por defecto que es de “0.01”.
G. Tensile failure reduces Hoek – Brown tensile strength to zero: Esta
opción se dejará activada debido a que trabajaremos con el criterio de
ruptura de Hoek & Brown generalizado.
H. Use Effective Stress Analysis: Esta casilla quedará desactivada; debido
a que no habrá cambios en la presión de poros y tampoco en la tensión
efectiva que afecten la deformación.

Figura 17. Configuración del análisis de esfuerzos (Fuente: Elaboración propia).

30
 5.2.4. RESUMÉN DEL PROYECTO

Figura 18. Resumen del proyecto. (Fuente: Elaboración propia).

5.3. ESTADO TENSIONAL

En este parte se cargan los esfuerzos los cuales fueron proporcionados por el
CGI. Estos datos pueden ser del orden Gravitacionales o constantes, para
nuestro caso particular tenemos información de esfuerzos constantes como se
muestra en la Figura 28.

No se va a considerar en este caso la información sobre cargas sísmicas, dado

que se usa principalmente para análisis más superficiales y no subterráneos.

Figura 19. Introducción de los valores del campo de esfuerzos (fuente: Elaboración propia).

31
5.4. PROPIEDADES DE LAS UNIDADES GEOTÉCNICAS
Para las siguiente 3 propiedades se va a trabajar con el criterio de ruptura de
Hoek & Brown generalizado.

5.4.1. PARED COLGANTE

Figura 20. Propiedades del material que presenta la pared colgante. (Fuente: Elaboración
propia).

Figura 21. Ingresando parámetros del macizo rocoso de la pared colgante. (Fuente:
Elaboración propia).

32
 5.4.2. PROPIEDADES DE LA PARED PENDIENTE

Figura 22. Propiedades del material que presenta la pared pendiente. (Fuente: Elaboración
propia).

Figura 23. Calculador de parámetros de macizo rocoso de la pared pendiente. (Fuente:

Elaboración propia).

33
 5.4.3. PROPIEDADES DE LA ROCA MINERALIZADA

Figura 24. Propiedades del material que presenta el mineral (Fuente: Elaboración propia).

Figura 25. Calculador de parámetros de macizo rocoso del mineral. (Fuente: Elaboración
propia).

34
 5.4.4. PROPIEDADES PARA EL RELLENO HIDRAÚLICO

Para el caso del relleno hidráulico se va hacer el análisis con el criterio de

ruptura de Mohr Coulomb.

Figura 26. Propiedades del material que presenta el relleno hidráulico. (Fuente: Elaboración
propia).

5.4.5. PROPIEDADES DE LA JUNTA

Las propiedades de la junta fueron proporcionadas por el CGI. Revisar Tabla 4.

Figura 27. Definición de las propiedades de la junta. (Fuente: Elaboración propia).

35
5.5. CREACIÓN DE LA MALLA EN ELEMENTOS FINITOS.

Una vez ingresado la información anterior procedemos con la generación de la

malla, para esto debemos tener presente que Phase2 puede crear diferentes
tipos de malla como:

 Graduada: El enmallado que se genera se hace más pequeño conforme

se acerca a la excavación, hacia los contornos de éstas y es la opción
más usada en Minería.
 Uniforme: Todos los elementos del enmallado son de tamaño uniforme.
 Radial: Normalmente usado en túneles circulares, la distribución de la
malla es radial.

En nuestro caso particular trabajaremos con: graduada, 3 nodos y triangular

Figura 28. Ingreso de datos para la generación de la malla. (Fuente: Elaboración propia).

36
Figura 29. Proyecto discretizado en elementos finitos. (Fuente: Elaboración propia).

37
5.6. MODULO DE CÁLCULO

Después de culminar con la introducción de toda la información al software

Phase2 procedemos a realizar los cálculos.

Figura 30. Computarización y cálculo de los elementos finitos. (Fuente. Elaboración propia).

5.8. MÓDULO DE INTERPRETACIÓN

Una vez realizado los cálculos por medio de los elementos finitos se procede a
la interpretación de los resultados, siendo el objetivo del proyecto: analizar el
factor de seguridad y de la sigma 1 de cada etapa.

38
FASE 1

Figura 31. Valores de Sigma 1 corte 1, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 32. Factor De Seguridad Corte 1 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

39
FASE 2

Figura 33. Valores de Sigma 1 corte 2, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 34. Factor De Seguridad Corte 2 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

40
FASE 3

Figura 35. Valores de Sigma 1 corte 3, en contornos de las excavaciones y (fuente: Elaboración
propia)

Figura 36. Factor De Seguridad Corte 3 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

41
FASE 4

Figura 37. Valores de Sigma 1 corte 4, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 38. Factor De Seguridad Corte 4 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

42
FASE 5

Figura 39. Valores de Sigma 1 corte 5, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 40. Factor de Seguridad Corte 5 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

43
FASE 6

Figura 41. Valores de Sigma 1 corte 6, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 42. Factor de Seguridad Corte 6 en contornos de la excavación y puente (fuente:

Elaboración propia)

44
FASE 7

Figura 43. Valores de Sigma 1 corte 7, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 44. Factor de Seguridad Corte 7 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

45
FASE 8

Figura 45. Valores de Sigma 1 corte 8, en contornos de las excavaciones y Pilar (fuente:
Elaboración propia)

Figura 46. Factor de Seguridad Corte 8 en contornos de la excavación y puente o pilar (fuente:
Elaboración propia)

46
FASE 9: RELLENO HIDRAULICO 1

Figura 47. Modelamiento de sigma 1 con relleno de hidráulico los pilares explotados. (Fuente:
Elaboración propia)

Figura 48. Modelamiento del factor de seguridad con el relleno hidráulico en los pilares
explotados. (Fuente: Elaboración propia).

47
FASE 10: RELLENO HIDRAULICO 2

Figura 49. Modelamiento de sigma 1 con relleno de hidráulico de todas las excavaciones.
(Fuente: Elaboración propia).

Figura 50. Modelamiento del factor de seguridad con el relleno hidráulico en todas
excavaciones del proyecto. (Fuente: Elaboración propia).

48
6. RESULTADOS
SIGMA(1)
F.S(Promedio
ETAPA DE EXCAVACIÓN W(m) Promedio en el CONDICIÓN
en el Pilar)
Pilar( Mpa)
FASE 1 55 43,4 2,89 Estable
FASE 2 40 47,3 1,78 Estable
FASE 3 35 46,72 1,57 Estable
FASE 4 30 49,79 1,08 Inestable
FASE 5 25 52,14 0,99 Inestable
FASE 6 20 55,39 0,88 Inestable
FASE 7 15 59,98 0,84 Inestable
FASE 8 10 66,94 0,67 Inestable
RELLENO HIDRAULICO 1 34,27 5,25 Estable
RELLENO HIDRAULICO 2 32,68 5,79 Estable
Tabla 5. Resultados promedio del FS y Sigma 1 durante las diferentes secuencias de minado y
rellenos.

De acuerdo a los resultados mostrados en la Tabla 5, podemos apreciar que:

 El Factor de seguridad presenta un descenso con forme se avanza la

secuencia de minado para las 8 fases.

 Podemos apreciar que hasta la fase 3 presenta una condición estable,

estando dentro del factor de diseño de 1.2.

 Por otra parte, la fase 4, 5, 6,7 y 8 presentan una condición inestable con
un factor de seguridad inferior al de diseño de1.2.

 Conforme avanza la secuencia de minado existe una redistribución de

tensiones.

 La aplicación del relleno hidráulico es un factor importante debido a que

eleva el factor de seguridad logrando así estabilizar el pilar.

 Se pudo apreciar que las magnitudes de la tensión mayor Sigma 1

aumentan conforme se avanza en la secuencia de minado, alcanzando
un valor de 66,94 Mpa en la fase 8.

 La importancia del relleno hidráulico radica en que va generar una

redistribución de tensiones como se puede apreciar en la tabla 10;
generando una disminución de la tensión principal y por ende el pilar se
va a encontrar con un factor de seguridad muy estable.

49
7. CONCLUSIONES Y RECOMENDACIONES

 El método de elementos finitos resulta muy útil, práctico y de rápido

procesamiento en 2D a la hora de hacer un análisis de una excavación
subterránea.
 Se concluye que el software phase2 es de gran utilidad a la hora de
determinar el comportamiento de las tensiones en cada etapa de minado
y el factor de seguridad; el cual es de suma ayuda para poder determinar
cómo se va comportar el macizo rocoso ya sea estable o inestable.
 Para una mayor exactitud en los resultados se recomienda discretizar la
malla en elementos cada vez más finos. Esto lo podemos configurar en la
opción “Number of Load Steps”. Debemos de tener presente que si
discretizamos la malla en elementos cada vez más pequeños el tiempo
de cálculo va a demorar mucho más.
 El software phase2 nos permite determinar el factor de seguridad en cada
etapa de minado lo cual nos va indicar si alguna zona requiere de una
técnica de sostenimiento para un minado seguro y evitar caída de rocas
o el colapso del pilar.
 Se recomienda tener un conocimiento previo sobre el método de
elementos finitos y terminologías tales como: Que es un nodo,
desplazamiento nodal, grados de libertad, interpolación de puntos de
Gauss, ecuaciones constitutivas, deformación, tensión, condiciones de
borde, matriz de rigidez, etc

50
BIBLIOGRAFÍA

G. Rodríguez. (2018). MODULO 4: Software Phase2, Centro Geotécnico

Internacional, Lima, Perú.

Gonzáles de Vallejo, Luis (2002) “Ingeniería geológica”. Madrid

Introducción al método de los elementos finitos: aplicación a la Mecánica de

sólidos.

Hoek, E., & Brown, E. (1980). Empirical strengh criterion for rock masses. ISRM
News Journal.

Zienkiewicz, O. C. (1982). Numerical and Analytical Methods in Geomechanics.

Barcelona.

Rocscience (2013), Excavation & Support Design Software, Phase2 8.0.

51
ANEXOS

Anexo 1. Información del análisis Tajo Esperanza Nv. 920

Phase2 Analysis Information

TAJEO ESPERANZA NV. 920

Project Summary

 File Name: TAJEO ESPERANZA NV.920.fez

 Last saved with Phase2 version: 8.005
 Project Title: TAJEO ESPERANZA NV. 920
 Analysis: TENSION-DEFORMACION
 Author: KADIMA MEJIA ZANABRIA
 Company: CENTRO GEOTECNICO INTERNACIONAL

General Settings

 Number of Stages: 10
 Analysis Type: Plane Strain
 Solver Type: Gaussian Elimination
 Units: Metric, stress as MPa

Analysis Options

 Maximum Number of Iterations: 500

 Tolerance: 0.001
 Number of Load Steps: Automatic
 Convergence Type: Absolute Energy
 Tensile Failure: Reduces Shear Strength
 Joint tension reduces joint stiffness by a factor of 0.01

Groundwater Analysis

 Method: Piezometric Lines

 Pore Fluid Unit Weight: 0.00981 MN/m3
 Probability: None

52
Field Stress

 Field stress: constant

 Sigma one: 35 MPa (compression positive)
 Sigma three: 27 MPa (compression positive)
 Sigma Z: 27 MPa (compression positive)
 Angle from the horizontal to sigma 1: 7 degrees (counter-clockwise)

Mesh

 Mesh type: graded

 Element type: 3 noded triangles
 Number of elements on FASE 1: 1875
 Number of nodes on FASE 1: 1037
 Number of elements on FASE 2: 1749
 Number of nodes on FASE 2: 983
 Number of elements on FASE 3: 1695
 Number of nodes on FASE 3: 959
 Number of elements on FASE 4: 1659
 Number of nodes on FASE 4: 944
 Number of elements on FASE 5: 1611
 Number of nodes on FASE 5: 923
 Number of elements on FASE 6: 1577
 Number of nodes on FASE 6: 909
 Number of elements on FASE 7: 1521
 Number of nodes on FASE 7: 884
 Number of elements on FASE 8: 1467
 Number of nodes on FASE 8: 860
 Number of elements on RELLENO HIDRAULICO 1: 1965
 Number of nodes on RELLENO HIDRAULICO 1: 1065
 Number of elements on RELLENO HIDRAULICO 2: 2165
 Number of nodes on RELLENO HIDRAULICO 2: 1136

Mesh Quality

 All elements are of good quality

Poor quality elements defined as:

 Side length ratio (maximum / minimum) > 30.00

 Minimum interior angle < 2.0 degrees
 Maximum interior angle > 175.0 degrees

Material Properties
53
Material: PARED COLGANTE

Color ___
Initial element loading field stress & body force
Unit weight 0.0258 MN/m3
Elastic type isotropic
Young's modulus 9250.1 MPa
Poisson's ratio 0.25
Failure Criterion Generalized Hoek-Brown
Material type Elastic
Compressive strength 150 MPa
mb parameter 2.75001
s parameter 0.0026
a parameter 0.505734
Piezo to use None
Ru value 0

Material: PARED PENDIENTE

Color ___
Initial element loading field stress & body force
Unit weight 0.0258 MN/m3
Elastic type isotropic
Young's modulus 9147.8 MPa
Poisson's ratio 0.25
Failure Criterion Generalized Hoek-Brown
Material type Elastic
Compressive strength 150 MPa
mb parameter 2.28617
s parameter 0.002049
a parameter 0.506582
Piezo to use None
Ru value 0

Material: ROCA MINERAL

Color ___
Initial element loading field stress only
Elastic type isotropic

54
 Young's modulus 4666.6 MPa
Poisson's ratio 0.27
Failure Criterion Generalized Hoek-Brown
Material type Elastic
Compressive strength 100 MPa
mb parameter 1.35306
s parameter 0.001434
a parameter 0.508086
Piezo to use None
Ru value 0

Material: RELLENO HIDRAULICO

Color ___
Initial element loading field stress & body force
Unit weight 0.0196 MN/m3
Elastic type isotropic
Young's modulus 300 MPa
Poisson's ratio 0.3
Failure criterion Mohr-Coulomb
Peak tensile strength 0 MPa
Residual tensile strength 0 MPa
Peak friction angle 32 degrees
Peak cohesion 0 MPa
Material type Elastic
Piezo to use None
Ru value 0

Joint Properties

Joint: JUNTA

Color ___
Normal stiffness 100000 MPa/m
Shear stiffness 10000 MPa/m
Initial joint deformation allowed
Pressure from Groundwater Analysis Not Included
Additional Pressure Inside Joint Not Included
Slip Criterion No slip allowed

55
List of All Coordinates

Excavation boundary

X Y
-28.0009 49.9858
-13.0009 49.9858
-10.8998 85.0028
-25.8998 85.0028

Excavation boundary

X Y
-32.2009 -20.0142
-17.2009 -20.0142
-16.3009 -5.01418
-15.4009 9.98582
-15.1009 14.9858
-14.8009 19.9858
-14.5009 24.9858
-14.2009 29.9858
-13.9009 34.9858
-13.6009 39.9858
-28.6009 39.9858
-28.9009 34.9858
-29.2009 29.9858
-29.5009 24.9858
-29.8009 19.9858
-30.1009 14.9858
-30.4009 9.98582
-31.3009 -5.01418

56
External boundary

X Y
40 -50
40 60.6884
40 100
32.4424 100
11.577 100
-10 100
-25 100
-30.017 100
-80 100
-80 -50
-39.9984 -50
-34 -50
-19.9726 -50
-19 -50
-3.98445 -50
16.2045 -50

Stage boundary

X Y
-28.9009 34.9858
-13.9009 34.9858

Stage boundary

X Y
-29.2009 29.9858
-14.2009 29.9858

Stage boundary

X Y
-29.5009 24.9858
-14.5009 24.9858

57
Stage boundary

X Y
-29.8009 19.9858
-14.8009 19.9858

Stage boundary

X Y
-30.1009 14.9858
-15.1009 14.9858

Stage boundary

X Y
-31.3009 -5.01418
-16.3009 -5.01418

Stage boundary

X Y
-30.4009 9.98582
-15.4009 9.98582

Joint boundary

X Y
-34 -50
-32.2009 -20.0142

Joint boundary

X Y
-25.8998 85.0028
-25 100

58
Joint boundary

X Y
-13.0009 49.9858
-13.6009 39.9858

Joint boundary

X Y
-28.6009 39.9858
-28.0009 49.9858

Joint boundary

X Y
-17.2009 -20.0142
-19 -50

Joint boundary

X Y
-10.8998 85.0028
-10 100

59