EICAM-DIS. GEO. Y ESTR. DE TÚN.-VF-2024 - Clase
EICAM-DIS. GEO. Y ESTR. DE TÚN.-VF-2024 - Clase
EICAM-DIS. GEO. Y ESTR. DE TÚN.-VF-2024 - Clase
DISEÑO GEOMÉTRICO Y
ESTRUCTURAL DE TÚNELES
2. Diseño Geológico-Geotécnico
Examen
Diseño de túneles
Síntesis y especificación de
Nota: El alcance de un programa de alternativas
investigaciones geotécnicas para un túnel con
coberturas superiores a 1.0 km, es limitado.
Macizo rocoso
• Obtención de parámetros representativos de la roca intacta
mediante ensayos de laboratorio y/o ensayos in-situ.
Probablemente nunca se entienda bien las propiedades del macizo rocoso, sino
hasta el momento en que se construya la obra!!!
Δ
Secciones típicas
Características del túnel
Los túneles se caracterizan por su trazado y sección,
definidos por criterios geométricos de gálibo, pendiente,
radio de curvatura y otras consideraciones de proyecto.
El sostenimiento se refiere a los elementos estructurales de
sujeción del terreno, aplicados inmediatamente después de
la excavación, con el fin de asegurar su estabilidad durante
la construcción y después de ella, así como garantizar las
condiciones de seguridad en el trabajo.
El revestimiento se coloca con posterioridad al
sostenimiento y consiste en aplicar sobre el mismo una capa
de hormigón, u otros elementos estructurales, con el fin de
proporcionar resistencia a largo plazo y dar un acabado
regular, mejorando su funcionalidad (condiciones
aerodinámicas, impermeabilidad, luminosidad, albergar
instalaciones y propiciar la estética de la obra).
«cuando el terreno no se
investiga, el terreno es un
riesgo»
• La inversión adecuada en los estudios
geológico-geotécnicos depende de la
complejidad geológica, longitud del
túnel, espesor de recubrimientos, etc.
-2a
Diseño Geológico-Geotécnico (*)
El punto de partida del diseño práctico de excavaciones
subterráneas se basa en las clasificaciones geomecánicas.
Δ
Las características generales que se evalúan en los distintos
sistemas de clasificación del macizo rocoso son:
– Resistencia de la roca intacta
– Índice RQD
– Espaciamiento y persistencia de las discontinuidades
– Orientación de las discontinuidades
– Características generales de las discontinuidades:
rugosidad, relleno.
– Estructura geológica: fallas, zonas débiles
– Filtraciones y presencia de agua
– Estado tensional
Principales Clasificaciones geotécnicas
1. Clasificación de Terzaghi (1946). (*) Las más utilizados
en la industria minera
(ordenadas de más a
2. Clasificación de Stini y Lauffer (1958). menos frecuente)
3. Clasificación de Deere (1964).
4. Clasificación de Bieniawski, o CSIR (1973). (3)
5. Clasificación de Barton o NGI (1974). (2)
6. Clasificación de Hoek o GSI (1995). (4)
7. Clasificación de Pälmstrom o RMi (1995).
8. Clasificación de Laubscher (1975), Laubscher & Jakubec
(2001). Índices IRMR y MRMR. (1)
(σH>σV)
Tensiones importantes:
Laderas con valles profundos y pendientes > 25º
Laderas con valles profundos y rocas blandas
Campos de tensiones naturales
The stress maps display the orientations of the maximum horizontal compressive stress SH. The length of the stress
symbols represents the data quality, with A being the best quality. Quality A data are assumed to record the orientation
of SH to within 10°-15°, quality B data to within 15°-20°, and quality C data to within 25°. Quality D data are considered to
give questionable tectonic stress orientations. (Zoback, 1992; Sperner et al., 2003; Heidbach et al., 2007). The tectonic
regimes are NF for Normal Faulting, SS for Strike-Slip faulting, TF for Thrust Faulting, and U for an Unknown regime. FMS
data which are flagged as PBEs are not included in these maps. www.world-stress-map.org
The WSM compiles information of the contemporary crustal stress from a wide range of
stress indicators. Basis of the compilation is the WSM quality ranking scheme to guarantee
comparability of the various stress indicators that reflect the stress state of very different
rock volumes. The stress indicators are grouped into four categories:
1. Earthquake focal mechanisms
2. Well bore breakouts and drilling-induced fractures
3. In-situ stress measurements (overcoring, hydraulic fracturing, borehole slotter)
4. Young geologic data (from fault-slip analysis and volcanic vent alignments)
Campo Tensional Regional
Incluye topografía
Campo Tensional Local
1 / 3: 2,90
1 / 2: 2,30
1 : 9,4 MPa; horizontal
2 / 3: 1,30
2 : 4,1 MPa; vertical
v / h (a): 1,30
3 : 3,2 MPa horizontal
v / h (b): 0,40
Campo Tensional Inducido
Estimación de la Tensión vertical in situ
Se suele suponer que la tensión vertical es igual al peso del
material ubicado por encima: σv = γ z
σv = γ z Se considera una buena estimación bajo superficies
horizontales, pero no cuando hay pendientes pronunciadas. La
superficie del terreno siempre forma una trayectoria de tensiones
principales — una tensión principal es perpendicular a la vertiente
y vale cero, las otras dos son paralelas al plano de la vertiente.
Cerca del fondo del valle, las tensiones in situ pueden alcanzar qu
Tensiones verticales vs profundidad
σv = γ z
(Válida estrictamente
para macizo rocoso
continuo, homogéneo
y superficie horizontal)
σh = σx = σy = [ / (1 – ) ] σz
k = σh/ σz = / (1 – )
Donde es el coeficiente de Poisson.
Para valores de = 0. 25, la ecuación anterior muestra que
las tensiones horizontales son un tercio de las verticales.
Relación
k = σh / σz con la
profundidad
La erosión tiende a
aumentar el valor de
k, con lo cual la
tensión horizontal es
mayor que la vertical
a profundidades bajas
(Hoek & Brown, 1980)
Medición de las tensiones in situ
Se realizan in situ: en sondeos, afloramientos, paredes de
galerías subterráneas. Las tensiones se obtienen a partir de
medidas de desplazamientos.
En todos los métodos se procede a alterar la roca para crear
una respuesta que pueda ser medida.
La precisión de la mediciones es baja, aceptándose resultados
con errores menores a 0.3 MPa.
Interesa determinar los valores de las 6 componentes del
tensor de tensiones.
Métodos de medición
1. Métodos directos
• Métodos de sobreperforación (overcoring)
- Célula biaxial: CSIR Doorstopper , USBM
Deformation gage
- Célula Triaxial: CSIR triaxial, CSIRO Triaxial Hollow
Inclusion Cell, LNEC Triaxial Solid Inclusion Cell
• Métodos de presurización
- Gato plano (flat jack)
- Fractura hidráulica
2. Métodos indirectos
• Borehole breakout
• Emisiones acústicas
• etc.
1. Métodos directos
En los métodos de sobreperforación (overcoring):
• Las células biaxiales permiten estimar las tres componentes de tensión
sobre el plano perpendicular al eje del sondeo (definiendo el estado
bidimensional). Se requiere un mínimo de dos ensayos (en general se
realizan tres) para conocer las seis componentes del tensor de tensiones.
• Las células triaxiales pueden proporcionar las seis componentes del tensor
de tensiones en un solo ensayo.
En los métodos de presurización:
• El ensayo del gato plano proporciona la tensión actuante en una dirección.
• El ensayo de fracturación hidráulica aporta las tensiones principales en las
tres direcciones del espacio.
σ x τ xy σ x τ xy τ xz σ x σ1
σy σy τ yz σ2
σ z σ 3
Célula biaxial Célula triaxial Ensayo del Fracturación
CSIR Doorstopper, CSIRO triaxial, CSIRO Triaxial gato plano hidráulica
USBM Deformation gage Hollow Inclusion Cell, LNEC
Triaxial Solid Inclusion Cell
Métodos de sobreperforación (Overcoring)
doorstopper
2. Método del USBM Deformation gage (overcoring)
Procedimiento para una Célula triaxial
(a) (b)
• Overcoring alrededor del segundo agujero y del mismo
diámetro que el primero (liberación de cargas). (c)
(c) (d)
• Si es posible, el material extraído debe ser sometido a una
carga de compresión que devuelva la muestra a sus
dimensiones iniciales, siendo esta carga equivalente al estado
de tensiones inicial, previo al corte. (e)
(e) (f)
Triaxial Solid Inclusión Cell
Se determina: ; alternativamente:
Triaxial Solid Inclusión Cell
Triaxial Solid Inclusión Cell
Método del Gato Plano (Flat Jack)
El método del gato plano (flat jack test) es un método de presurización (o de
restitución de tensiones):
• Tincelin, France, 1952.
• El método consiste en realizar una entalladura (corte plano) en la pared de
una galería o de un talud en el que se introduce un gato hidráulico plano.
• La entalladura puede realizarse con una sierra de disco o bien con
perforaciones adyacentes.
• El estado tensional se obtiene en función del volumen de líquido (agua o
aceite) necesario para que la roca recupere el estado de tensiones previo a
la realización del corte.
• Especialmente diseñado para la medición de tensiones
paralelas a la superficie de una roca y en las proximidades de
la misma.
Entalladura
Métodos de Fracturación Hidráulica
Método estándar de medida de tensiones in situ a grandes
profundidades (mayores de 500 m, fuera del alcance de otros
métodos).
EFH EHFP
1. Ensayo de Fracturación Hidráulica (EFH)
Procedimiento: El ensayo se realiza en un sondeo siguiendo el
siguiente procedimiento:
Seleccionar el tramo donde se realizará el ensayo como
mínimo de 1 m de longitud.
Sellar el tramo con dos "packers“ (obturadores).
Inyectar agua a presión en el tramo sellado hasta que se
produzca la fractura de la roca.
Se mide la presión de fluido necesaria para generar,
propagar, mantener y reabrir la fractura.
Se pueden utilizar obturadores de impresión de la pared del
sondeo para conocer la orientación de la fractura, de manera
que con estos tres datos, situación de la fractura, Pf y Ps, se
puede conocer el estado de tensiones alrededor del sondeo.
• Se supone que la fractura se inicia cuando la presión del fluido
supera la tensión circunferencial en la pared del sondeo y luego
la resistencia a la tracción de la roca.
convergencia
extrusión
C. Análisis tenso-deformacionales.
Los casos B y C se dan durante la etapa de diseño, ya que se cuenta con tiempo y
recursos. Se utiliza uno u otro (o ambos) en función del problema y del interés en
que el proyecto cuente con un buen nivel de ingeniería.
5.1 ¿Cómo realizar el calculo estructural del túnel?
B. Análisis 2D o 3D
C. Niveles de incertidumbre
También puede haber una mayor dependencia de la rigidez del material, del
nivel de esfuerzos alcanzados por el medio, así como del nivel de deformaciones.
7.2 ¿Qué programa de cálculo utilizar en la modelación?
Por otro lado, conforme los modelos disponibles en el mercado ofrecen más
posibilidades, el número de incertidumbres crece y su uso exige campañas
geotécnicas que pocas veces se justifican en una obra subterránea convencional.
En los dos casos se combinan atributos del macizo rocoso (de tipo geológico,
geométrico y tensional) en un número único relacionado con la calidad global
de la roca.
A su vez, este número permite a través de la experiencia recogida en su
utilización en casos reales, la definición de un sostenimiento del túnel y la
estimación de otros parámetros o datos de interés (resistencia del macizo
rocoso, tiempo de estabilidad de una excavación no sostenida, etc.).
Relación entre claro activo y tiempo de sostén en función del tipo de roca
Ej: Sostenimiento
necesario para un túnel
de 4.27 m de luz
(Skinner, 1988)
Terzaghi
RQD
Trayectoria de la tensiones
principales mayores y menores
alrededor de un túnel circular,
obtenida para un estado de
tensiones inicial uniaxial (σv ≠0 ;
σh =0) mediante un modelo
elástico lineal, bajo hipótesis de
homogeneidad, isotropía y
deformación plana.
Deflexión de las líneas de corriente alrededor de una
obstrucción cilíndrica
Solución analítica :
Independiente de las constantes elásticas y del tamaño
de la excavación.
Resultado inaceptable desde un punto de vista
ingenieril.
Lo interesante del cálculo no es la solución en sí, sino lo
que se desprende de ella.
En primera aproximación, da una idea de que las
tensiones no están controladas por las características
del material sino por la geometría del túnel. (Muy
importante a la hora de proyectar un sostenimiento).
Variaciones de las tensiones de borde
tangencial en función de k = ph / pz
σθA=3pz
C
σθA=2pz
A σθC=2pz
k = 0.33
σθC=-pz
σθA = pz ( 3 – k )
k > 0.33 y k < 3, σθ = siempre compresión
k < 0.33, aparecen tracciones σθC = pz( 3 k - 1 )
Tensiones tangenciales en una
excavación elíptica
Solución analítica
(campo biaxial de tensiones)
ρA , ρC = radio de curvatura en
A y C.
k=0
Circular
σA = 3 pz
σA = 3 pz
σC = -pz
Tensiones máximas
(compresiones)
Tensiones máximas
(tracciones)
Disposición pésima para
las tensiones en hastiales.
Conclusiones sobre excavaciones simples
k = 0.5 k=1
Sección típica en carreteras y ferrocarril
k=1 k = 0.5
Geometría “ideal” en función de los estados de tensiones en techo y hastiales.
TECHO
HASTIAL
Si queremos conocer más soluciones analíticas podemos
consultar los trabajos de los siguientes autores:
Jaeger y Cook “Fundamentals of Rock Mechanics”. Chapman
Hall, 11
Muskhlishvili “Some basic problems of the Math. Theor. Of
Elast.”. Nordhoff
Savin (1961) “Stress concentrations around holes”. Pergamon
Comentarios finales:
Podemos adaptar la forma de la sección de excavación al
estado de tensiones pero en la práctica nadie diseña así los
túneles.
Quizá para un caso muy concreto podría llevarse a cabo, pero
carece de sentido ir modificando la sección en función de las
características de las litologías que vamos atravesando.
RESISTENCIA DE LA ROCA MATRIZ Y DE MACIZOS ROCOSOS
Interesa conocer las tensiones que pueden resistir tanto la roca matriz como el
macizo rocoso utilizando criterios de rotura.
También interesa conocer la influencia de los distintos factores que afectan la
estabilidad del túnel en función de la profundidad, ya que ello condicionará el
mecanismo de rotura.
Por lo tanto, desde el punto de vista de la rotura, resulta útil distinguir entre
túneles someros y túneles profundos en roca.
En un túnel somero, los factores a considerar dependerán de la estructura,
grado de alteración de las juntas y discontinuidades del macizo rocoso. Por
tanto, cuando hablemos de resistencia nos estaremos refiriendo
fundamentalmente a la resistencia del macizo rocoso.
Para el caso de túneles profundos en roca, el estado de tensiones será
relevante. Cuando hagamos referencia a la resistencia nos estaremos
refiriendo principalmente a la resistencia de la roca matriz.
INTERACCIÓN TÚNEL-SOSTENIMIENTO
Se considera el avance de la
excavación y cuatro secciones
significativas.
Lejos del frente, en la roca (sección
AA’), sobre el futuro contorno
teórico del túnel actúa la tensión p0.
Esta sección aún no se ha
deformado, de manera que el
desplazamiento radial, ui de los
puntos de la sección teórica del
túnel es nulo.
En la sección BB’, ya excavada y próxima
al frente, la tensión p0 ha desaparecido y
el contorno del túnel ha experimentado
un desplazamiento hacia el interior (ui).
Debido a la marcada tridimensionalidad
del problema no es posible en principio
efectuar un análisis bidimensional en
sección plana.
pi = k (ui – ud)
La cavidad esférica, aparte del interés que tiene en si misma para el análisis de
excavaciones subterráneas de formas diversas, es una aproximación
interesante al comportamiento en las proximidades del frente y proporciona
información útil para entender sus condiciones de estabilidad.
DETERMINACIÓN DE LA CURVA CARACTERÍSTICA
ELASTICIDAD. TÚNEL CIRCULAR EN DEFORMACIÓN PLANA
Ɛz=0
2. “La cavidad esférica, aparte del interés que tiene en si misma para el análisis de excavaciones
subterráneas de formas diversas, es una aproximación interesante al comportamiento en las
proximidades del frente y proporciona información útil para entender sus condiciones de
estabilidad.”
DETERMINACIÓN DE LA CURVA CARACTERÍSTICA
ELASTOPLASTICIDAD. TÚNEL CIRCULAR EN DEFORMACIÓN
PLANA. CRITERIO DE ROTURA DE MOHR-COULOMB
= tensión radial en el
contacto entre las zonas
elástica y plástica
En la frontera, r = re
Zona elastoplástica: (ri < r < re)
Si 1 y 3 son las tensiones principales mayor y menor, el criterio
de rotura de Mohr-Coulomb se escribe:
= ángulo de rozamiento interno
(1) c = cohesión
(1)
Por lo anterior, t = trayectoria de tensiones típica de los puntos del terreno próximos al contorno.
A partir de un estado isótropo de tensiones, σθ (identificable con σ1) aumenta, mientras que σr
(identificable con σ3) disminuye.
La condición límite se alcanza en la envolvente representada por la ecuación (1).
Se supondrá que tras alcanzar la superficie límite (fluencia) las deformaciones continúan a tensión
desviadora (σ1-σ3 ) constante (plasticidad perfecta).
Para calcular las deformaciones plásticas se supone una ley de plasticidad no asociada.
Se adopta un potencial plástico similar a (1) con un ángulo de dilatancia ψ:
(2) C = constante
resulta (3)
En el límite, cuando r re
las tensiones elásticas deben
encontrarse en rotura, es
decir, cumplir con el criterio
de Mohr-Coulomb:
j
Si los sostenimientos actúan en tiempos
o distancias al frente diferentes, su
composición debe tener en cuenta este
hecho.
Dovelas y juntas
CERCHAS METÁLICAS Eb : módulo elástico del material del
bloque de apoyo (acero o madera)
wb : anchura de cada bloque de apoyo
tb : espesor del bloque de apoyo
Es : módulo de elasticidad del acero
wb : anchura de cada bloque de apoyo
As : sección de la cercha (acero)
Is : momento de inercia de la cercha
S = espaciamiento entre cerchas en
dirección longitudinal
X : canto (alma del perfil) de las
cerchas
σys : resistencia a tracción del acero
(rigidez)
(Hoek-Brown)
Por ello, durante mucho tiempo los túneles han sido diseñados y
construidos sin tomar en cuenta los efectos sísmicos, y sólo desde los
60’s fueron incorporados, por primera vez, procedimientos de diseño
sísmicos en proyectos subterráneos (Wang, 1993).
ANÁLISIS DE TÚNELES EN ROCA
Sección transversal del metro de Kobe (m) en la estación Daikai, antes y luego
del sismo de Kobe (Liu, 2008).
ANÁLISIS DE TÚNELES EN ROCA
A la izquierda, el daño sufrido por la estación de metro Daikai luego del sismo de 1995, a
la derecha, la subsidencia del terreno sobre la estación de metro Daikai (Liu, 2008).
Distancia epicentral: 25 km
Medidas de refuerzo
Este enfoque resulta satisfactorio para túneles más flexibles que el terreno, pero
puede ser muy conservador en caso contrario.
Métodos analíticos.
Modelos numéricos.
Una vez establecidos los escenarios sísmicos de diseño ODE (Operating Design
Earthquake) y MDE (Maximum Design Earthquake), es necesario determinar las
aceleraciones, velocidades y desplazamientos del terreno asociados a cada uno
de ellos.
Este enfoque resulta satisfactorio para túneles más flexibles que el terreno, pero
puede ser muy conservador en caso contrario.
Análisis sísmico de túneles - Metodología
Modelos numéricos
Refinamiento de
la malla de
elementos finitos
ECUACIONES DEL MEF
Hipótesis:
Acción de la gravedad.
Se consideran 3 etapas de análisis:
• Etapa 1: excavación del túnel;
• Etapa 2: excavación de la zanja;
• Etapa 3: carga del edificio.
Secuencia de pasos:
1. MODEL: Construcción del modelo.
2. COMPUTE: Cálculo. Seguir el Tutorial
= 6143.7 MPa
Este método requiere construir un modelo 2D-deformación plana del túnel con
presión interna = 0 en el interior de la excavación y determinar:
La etapa final, con presión interna cero, se utiliza para determinar la cantidad de
deformación antes de la instalación del soporte (paso 1).
En el diseño del túnel, estas cargas se obtienen a partir del análisis numéricos
en los cuales los "elementos de viga“ son unidos al contorno del túnel y los
esfuerzos normales, momentos flectores y esfuerzos de corte inducidos en
estos elementos se calculan directamente.
Por otro lado, el uso de pernos de roca actúan como "refuerzo de túnel“
cambiando las propiedades mecánicas de la masa rocosa que lo rodea.
Sección de ancho b de un
revestimiento compuesto por
un conjunto de vigas de acero
espaciadas a una distancia s,
incrustadas en hormigón
proyectado.
Sección de ancho b de un
revestimiento compuesto por
un conjunto de vigas de acero
espaciadas a una distancia s,
incrustadas en hormigón
proyectado.
NOTA: La “división de carga” determina cómo se aplica la carga de los esfuerzos in situ
(y/o las fuerzas de masa, si corresponde) a los elementos sólidos de la malla de elementos
finitos. La división de carga NO determina (directamente) la carga que se aplica a los
elementos de soporte (es decir, bulones o revestimiento). Recuerde que la carga se genera
en los elementos de soporte debido a la deformación de la malla de elementos finitos. Los
esfuerzos in situ no se aplican directamente a los elementos de soporte. Por lo tanto, la
división de carga no afecta DIRECTAMENTE la carga del sostenimiento, sino
INDIRECTAMENTE, debido al efecto en la malla de elementos sólidos.
….. Ir al Tutorial RS2 – ADDING SUPPORT
FIN