Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Introducción

La  roca  en  profundidad  está  sujeta  a  tensiones  resultantes  del  peso  de  los  estratos  suprayacentes  y  
de  las  tensiones  de  origen  tectónico.  Cuando  se  excava  una  abertura  en  esta  roca,  el  campo  de  
tensiones  se  interrumpe  localmente  y  se  induce  un  nuevo  conjunto  de  tensiones  en  la  roca  que  rodea  
la  abertura.  El  conocimiento  de  las  magnitudes  y  direcciones  de  estos  esfuerzos  in  situ  e  inducidos  es  
un  componente  esencial  del  diseño  de  excavaciones  subterráneas  ya  que,  en  muchos  casos,  se  
excede  la  resistencia  de  la  roca  y  la  inestabilidad  resultante  puede  tener  graves  consecuencias  en  el  
comportamiento  de  las  excavaciones.

Este  capítulo  trata  de  la  cuestión  de  las  tensiones  in  situ  y  también  de  los  cambios  de  tensión  que  se  
inducen  cuando  se  excavan  túneles  o  cavernas  en  roca  sometida  a  tensión.  Los  problemas  asociados  
con  la  falla  de  la  roca  alrededor  de  las  aberturas  subterráneas  y  con  el  diseño  del  soporte  para  estas  
aberturas  se  tratarán  en  capítulos  posteriores.

La  presentación,  que  sigue,  pretende  cubrir  solo  aquellos  temas  que  son  esenciales  para  el  
conocimiento  del  lector  cuando  se  trata  del  análisis  de  la  inestabilidad  inducida  por  esfuerzos  y  el  
diseño  del  soporte  para  estabilizar  la  roca  en  estas  condiciones.

Tensiones  in  situ

Considere  un  elemento  de  roca  a  una  profundidad  de  1000  m  por  debajo  de  la  superficie.  El  peso  de  
la  columna  vertical  de  roca  que  descansa  sobre  este  elemento  es  el  producto  de  la  profundidad  y  el  
peso  unitario  del  macizo  rocoso  suprayacente  (típicamente  alrededor  de  2,7  toneladas/m3  o  0,027  MN/m3).
Por  lo  tanto,  la  tensión  vertical  sobre  el  elemento  es  de  2700  toneladas/m2  o  27  MPa.  Esta  tensión  se  
estima  a  partir  de  la  relación  simple:

sigma  subv  =  gamma
σv z (1)
sub­z
donde  sigma  sub  v  es  la  tensión  vertical
donde  σv  es  la  tensión  vertical  es  el  peso  

de  unitario  
la  roca  sduprayacente  
e  la  roca  suprayacente   y  gamma  ebs  
y  z  es  la  profundidad   el  pla  
ajo   eso  
unitario  z  
superficie.  
es  la  profundidad  debajo  de  la  superficie.

Las  mediciones  de  tensión  vertical  en  varios  sitios  de  minería  e  ingeniería  civil  en  todo  el  mundo  
confirman  que  esta  relación  es  válida  aunque,  como  se  ilustra  en  la  Figura  1,  hay  una  cantidad  
significativa  de  dispersión  en  las  mediciones.
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Figura  1:  Mediciones  de  tensión  vertical  de  proyectos  de  minería  e  ingeniería  civil  en  todo  el  mundo.  
(Después  de  Brown  y  Hoek  1978).

Los  esfuerzos  horizontales  que  actúan  sobre  un  elemento  de  roca  a  una  profundidad  z  por  debajo  de  
la  superficie  son  mucho  más  difíciles  de  estimar  que  los  esfuerzos  verticales.  Normalmente,  la  
relación  entre  el  esfuerzo  horizontal  promedio  y  el  esfuerzo  vertical  se  denota  con  la  letra  k  tal  que:

sub  h  =kk  
gamma     ks  
ub  
zs  
igma  
sσigma  
σ  hsub  v  =  k   (2)
z v

Terzaghi  y  Richart  (1952)  sugirieron  que,  para  un  macizo  rocoso  cargado  gravitacionalmente  en  el  
que  no  se  permitió  deformación  lateral  durante  la  formación  de  los  estratos  suprayacentes,  el  valor  
) , es  la  relación  de  Poisson  de  El  
de  k  es  independiente  de  la  profundidad  y  está  dado  por  k  (1  donde  
macizo  rocoso  Esta  relación  se  usó  ampliamente  en  los  primeros  días  de  la  mecánica  de  rocas  pero,  
como  se  analiza  a  continuación,  resultó  ser  inexacta  y  rara  vez  se  usa  en  la  actualidad.

Las  mediciones  de  tensiones  horizontales  en  obras  civiles  y  mineras  de  todo  el  mundo  muestran  que  
la  relación  k  tiende  a  ser  alta  a  poca  profundidad  y  que  disminuye  a  mayor  profundidad  (Brown  y  
Hoek,  1978,  Herget,  1988).  Para  comprender  la  razón  de  estas  variaciones  de  esfuerzos  horizontales,  
es  necesario  considerar  el  problema  en  una  escala  mucho  mayor  que  la  de  un  solo  sitio.

Sheorey  (1994)  desarrolló  un  modelo  de  estrés  térmico  elastostático  de  la  tierra.  Este  modelo  
considera  la  curvatura  de  la  corteza  y  la  variación  de  las  constantes  elásticas,  de  densidad  y  térmicas.

2
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

coeficientes  de  expansión  a  través  de  la  corteza  y  el  manto.  Una  discusión  detallada  sobre  el  modelo  
de  Sheorey  está  más  allá  del  alcance  de  este  capítulo,  pero  proporcionó  una  ecuación  simplificada  
que  se  puede  usar  para  estimar  la  relación  de  tensión  horizontal  a  vertical  k.  Esta  ecuación  es:

(3)

donde  z  (m)  es  la  profundidad  bajo  la  superficie  y  Eh  (GPa)  es  el  módulo  de  deformación  promedio  
de  la  parte  superior  de  la  corteza  terrestre  medido  en  dirección  horizontal.  Esta  dirección  de  medición  
es  importante  particularmente  en  rocas  sedimentarias  en  capas,  en  las  que  el  módulo  de  deformación  
puede  ser  significativamente  diferente  en  diferentes  direcciones.

En  la  Figura  2  se  muestra  un  gráfico  de  esta  ecuación  para  un  rango  de  módulos  de  deformación.  
Las  curvas  que  relacionan  k  con  la  profundidad  bajo  la  superficie  z  son  similares  a  las  publicadas  
por  Brown  y  Hoek  (1978),  Herget  (1988)  y  otros  para  esfuerzos  medidos  in  situ.  Por  tanto,  se  
considera  que  la  ecuación  3  proporciona  una  base  razonable  para  estimar  el  valor  de  k.

Figura  2:  Relación  de  tensión  horizontal  a  vertical  para  diferentes  módulos  de  deformación  basados  
en  la  ecuación  de  Sheorey.  (Después  de  Sheorey  1994).
Como  señaló  Sheorey,  su  trabajo  no  explica  la  ocurrencia  de  esfuerzos  verticales  medidos  que  son  
más  altos  que  la  presión  de  sobrecarga  calculada,  la  presencia  de  esfuerzos  horizontales  muy  altos  
en  algunos  lugares  o  por  qué  los  dos  esfuerzos  horizontales

3
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

rara  vez  son  iguales.  Estas  diferencias  probablemente  se  deban  a  características  topográficas  y  
geológicas  locales  que  no  pueden  tenerse  en  cuenta  en  un  modelo  a  gran  escala  como  el  propuesto  
por  Sheorey.

Cuando  los  estudios  de  sensibilidad  hayan  demostrado  que  es  probable  que  las  tensiones  in  situ  
tengan  una  influencia  significativa  en  el  comportamiento  de  las  aberturas  subterráneas,  se  
recomienda  medir  las  tensiones  in  situ.  Las  sugerencias  para  establecer  un  programa  de  medición  
de  estrés  se  analizan  más  adelante  en  este  capítulo.

El  mapa  mundial  de  estrés

El  proyecto  World  Stress  Map,  terminado  en  julio  de  1992,  involucró  a  más  de  30  científicos  de  18  
países  y  se  llevó  a  cabo  bajo  los  auspicios  del  International  Lithosphere  Project  (Zoback,  1992).  El  
objetivo  del  proyecto  era  compilar  una  base  de  datos  global  de  datos  de  estrés  tectónico  
contemporáneo.

El  World  Stress  Map  (WSM)  ahora  se  mantiene  y  ha  sido  ampliado  por  el  Instituto  Geofísico  de  la  
Universidad  de  Karlsruhe  como  un  proyecto  de  investigación  de  la  Academia  de  Ciencias  y  
Humanidades  de  Heidelberg.  La  versión  2005  del  mapa  contiene  aproximadamente  16.000  conjuntos  
de  datos  y  se  pueden  descargar  de  Internet  varias  versiones  del  mapa  para  el  mundo,  Europa,  
América,  África,  Asia  y  Australia.  El  WSM  es  una  base  de  datos  de  acceso  abierto  a  la  que  se  puede  
acceder  en  www.world­stress­map.org  (Reinecker  et  al,  2005)

El  Mapa  de  estrés  mundial  de  2005  se  reproduce  en  la  Figura  3,  mientras  que  un  mapa  de  estrés  
para  el  Mediterráneo  se  reproduce  en  la  Figura  4.

Los  mapas  de  tensión  muestran  las  orientaciones  de  la  máxima  tensión  de  compresión  horizontal.
La  longitud  de  los  símbolos  de  estrés  representa  la  calidad  de  los  datos,  siendo  A  la  mejor  calidad.  
Se  supone  que  los  datos  de  calidad  A  registran  la  orientación  del  esfuerzo  de  compresión  horizontal  
máximo  dentro  de  10°­15°,  los  datos  de  calidad  B  dentro  de  15°­20°  y  los  datos  de  calidad  C  dentro  
de  25°.  Se  considera  que  los  datos  de  calidad  D  dan  orientaciones  de  tensión  tectónica  cuestionables.

La  versión  de  1992  del  World  Stress  Map  se  derivó  principalmente  de  observaciones  geológicas  
sobre  mecanismos  focales  de  terremotos,  alineaciones  volcánicas  e  interpretaciones  de  
deslizamiento  de  fallas.  Menos  del  5%  de  los  datos  se  basaron  en  mediciones  de  fracturación  
hidráulica  o  de  excavación  del  tipo  comúnmente  utilizado  en  proyectos  de  minería  e  ingeniería  civil.  
Por  el  contrario,  la  versión  de  2005  del  mapa  incluye  un  número  significativamente  mayor  de  
observaciones  de  rupturas  de  pozos,  fracturación  hidráulica,  extracción  de  muestras  y  ranurado  de  
pozos.  Por  lo  tanto,  vale  la  pena  considerar  la  precisión  relativa  de  estas  medidas  en  comparación  
con  las  observaciones  geológicas  en  las  que  se  basó  el  mapa  original.

4
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Figura  3:  Mapa  mundial  de  tensiones  que  da  orientaciones  de  la  máxima  tensión  
de  compresión  horizontal.  De  www.world­stress­map.org.

5
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Figura  4:  Mapa  de  tensiones  del  Mediterráneo  dando  orientaciones  de  la  máxima  
tensión  de  compresión  horizontal.  De  www.world­stress­map.org.

6
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Al  analizar  las  mediciones  de  esfuerzos  de  fracturamiento  hidráulico  y  de  extracción  de  muestras,  Zoback  
(1992)  tiene  los  siguientes  comentarios:

'Las  pruebas  detalladas  de  fracturación  hidráulica  en  varios  pozos  que  comienzan  muy  cerca  de  la  
superficie  (10­20  m  de  profundidad)  han  revelado  cambios  marcados  en  las  orientaciones  de  los  
esfuerzos  y  las  magnitudes  relativas  con  la  profundidad  en  los  primeros  cientos  de  metros,  
posiblemente  relacionados  con  los  efectos  de  la  topografía  cercana  o  una  alto  grado  de  fractura  
cerca  de  la  superficie.

Incluidas  en  la  categoría  de  mediciones  de  tensión  de  'superposición  de  testigos'  hay  una  variedad  
de  técnicas  de  medición  de  alivio  de  tensión  o  tensión.  Estas  técnicas  involucran  una  medición  
tridimensional  del  alivio  de  tensión  en  un  cuerpo  de  roca  cuando  está  aislado  del  volumen  de  roca  
circundante;  el  tensor  de  tensión  tridimensional  se  puede  calcular  posteriormente  con  un  
conocimiento  del  tensor  de  cumplimiento  completo  de  la  roca.  Hay  dos  inconvenientes  principales  
con  esta  técnica  que  restringen  su  utilidad  como  indicador  de  tensión  tectónica:  las  mediciones  
deben  realizarse  cerca  de  una  superficie  libre  y  el  alivio  de  la  tensión  se  determina  en  áreas  muy  
pequeñas  (de  unos  pocos  milímetros  cuadrados  a  centímetros  cuadrados).  Además,  se  ha  
demostrado  que  las  mediciones  cercanas  a  la  superficie  (con  mucho,  las  más  comunes)  están  
sujetas  a  los  efectos  de  la  topografía  local,  la  anisotropía  de  las  rocas  y  la  fracturación  natural  
(Engelder  y  Sbar,  1984).  Además,  muchas  de  estas  mediciones  se  han  realizado  para  aplicaciones  
de  ingeniería  específicas  (p.  ej.,  evaluación  del  sitio  de  una  represa,  trabajos  de  minería),  lugares  
donde  la  topografía,  fracturamiento  o  excavaciones  cercanas  podrían  perturbar  fuertemente  el  
campo  de  tensión  regional.'

Obviamente,  desde  una  escala  global  o  incluso  regional,  el  tipo  de  mediciones  de  tensión  de  ingeniería  
realizadas  en  una  mina  o  en  un  sitio  de  ingeniería  civil  no  se  consideran  muy  confiables.  Por  el  contrario,  las  
versiones  del  World  Stress  Map  que  se  presentan  en  las  Figuras  3  y  4  solo  se  pueden  usar  para  dar  
estimaciones  de  primer  orden  de  las  direcciones  de  estrés  que  probablemente  se  encontrarán  en  un  sitio  
específico.  Dado  que  tanto  las  direcciones  como  las  magnitudes  de  las  tensiones  son  de  importancia  crítica  
en  el  diseño  de  excavaciones  subterráneas,  se  deduce  que  puede  ser  necesario  un  programa  de  medición  
de  tensiones  en  cualquier  proyecto  importante  de  minería  subterránea  o  de  ingeniería  civil.

Desarrollo  de  un  programa  de  medición  de  estrés

Considere  el  ejemplo  de  un  túnel  que  se  excavará  a  una  profundidad  de  1000  m  por  debajo  de  la  superficie  
en  un  entorno  de  roca  dura.  La  profundidad  del  túnel  es  tal  que  es  probable  que  las  tensiones  in  situ  e  
inducidas  sean  una  consideración  importante  en  el  diseño  de  la  excavación.
Los  pasos  típicos  que  podrían  seguirse  en  el  análisis  de  este  problema  son:

El  Mapa  Mundial  de  Tensiones  para  el  área  bajo  consideración  dará  una  buena  primera  indicación  de  la  
posible  complejidad  del  campo  de  tensiones  regional  y  posibles  direcciones  para  la  máxima  tensión  de  
compresión  horizontal.

1.  Durante  el  diseño  preliminar,  la  información  presentada  en  las  ecuaciones  1  y  3  se  puede  utilizar  para  obtener  una  primera  estimación  aproximada  de  la  tensión  vertical  y  
1.  Durante  
horizontal   el  
pd
media  en  las   iseño  pdreliminar,  
roximidades   el  túnel.  Para  ulna  
a  pirofundidad  
nformación   pmresentada  
de  1.000   en  dlan  
,  estas  ecuaciones   as   etensión  
la   cuaciones   1  ys  
vertical  sigma   3  vs  =e  
ub   pMuede  
  27   utilizar  
k  =  1p,3  
Pa,  la  relación   ara  
(para  Eh  =  75  GPa)  y,  por  tanto,  la  tensión  horizontal  media  sigma  sub  h=  35,1  MPa.  Un  análisis  preliminar  de  las  tensiones  inducidas  alrededor  del  túnel  propuesto  muestra  
obtener  una  primera  estimación  aproximada  de  la  vertical  y  horizontal  promedio
que  es  probable  que  estas  tensiones  inducidas  excedan  la  resistencia  de  la  roca  y  que  la  cuestión  de  la  medición  de  la  tensión  debe  considerarse  con  más  detalle.  Tenga  
en  cuenta  7  que  para  muchas  aberturas  en  roca  fuerte  a  poca  profundidad,  los  problemas  de  tensión  pueden  no  ser  significativos  y  el  análisis  no  necesita  continuar.
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

tensión  en  las  inmediaciones  del  túnel.  Para  una  profundidad  de  1.000  m,  estas  
ecuaciones  dan  la  tensión  
σ
v vertical  =  27  MPa,  la  relación  k  =  1,3  (para  Eh  =  75  GPa)  y,  
σ h  =  35,1  MPa.  Un  análisis  preliminar  de  las  
por  tanto,  la  tensión  horizontal  media  
tensiones  inducidas  alrededor  del  túnel  propuesto  muestra  que  es  probable  que  estas  
tensiones  inducidas  excedan  la  resistencia  de  la  roca  y  que  la  cuestión  de  la  medición  
de  la  tensión  debe  considerarse  con  más  detalle.  Tenga  en  cuenta  que  para  muchas  
aberturas  en  roca  fuerte  a  poca  profundidad,  los  problemas  de  tensión  pueden  no  ser  
significativos  y  el  análisis  no  necesita  continuar.

Para  este  caso  particular,  los  problemas  de  estrés  se  consideran  importantes.  Un  próximo  paso  
típico  sería  buscar  en  la  literatura  en  un  esfuerzo  por  determinar  si  los  resultados  de  los  
programas  de  medición  de  tensión  in  situ  están  disponibles  para  proyectos  de  ingeniería  civil  o  
minas  dentro  de  un  radio  de,  digamos,  50  km  del  sitio.  Con  suerte,  estarán  disponibles  algunos  
resultados  de  medición  de  tensión  para  la  región  en  la  que  se  encuentra  el  túnel  y  estos  
resultados  se  pueden  usar  para  refinar  el  análisis  discutido  anteriormente.

Suponiendo  que  los  resultados  del  análisis  de  las  tensiones  inducidas  en  la  roca  que  rodea  el  
túnel  propuesto  indican  que  es  probable  que  se  desarrollen  zonas  significativas  de  falla  de  la  
roca  y  que  los  costos  de  soporte  probablemente  sean  altos,  probablemente  sea  justificable  
establecer  un  proyecto  de  medición  de  tensiones.  en  el  sitio.  Estas  mediciones  se  pueden  
realizar  en  pozos  profundos  desde  la  superficie,  mediante  técnicas  de  fracturación  hidráulica,  o  
desde  accesos  subterráneos  mediante  métodos  de  overcoring.  La  elección  del  método  y  el  
número  de  mediciones  a  realizar  depende  de  la  urgencia  del  problema,  la  disponibilidad  de  
acceso  subterráneo  y  los  costos  involucrados  en  el  proyecto.  Tenga  en  cuenta  que  muy  pocas  
organizaciones  de  proyectos  tienen  acceso  al  equipo  necesario  para  llevar  a  cabo  un  proyecto  
de  medición  de  tensión  y,  en  lugar  de  comprar  este  equipo,  puede  valer  la  pena  traer  una  
organización  que  tenga  el  equipo  y  que  se  especialice  en  tales
mediciones.

2.  Cuando  es  probable  que  se  encuentren  características  tectónicas  regionales,  como  fallas  
importantes,  las  tensiones  in  situ  en  la  vecindad  de  la  característica  pueden  rotar  con  
respecto  al  campo  de  tensión  regional.  Las  tensiones  pueden  ser  significativamente  
diferentes  en  magnitud  de  los  valores  estimados  a  partir  de  las  tendencias  generales  
descritas  anteriormente.  Estas  diferencias  pueden  ser  muy  importantes  en  el  diseño  de  
las  aberturas  y  en  la  selección  del  soporte  y,  cuando  se  sospecha  que  es  probable  que  
sea  así,  las  mediciones  de  tensión  in  situ  se  convierten  en  un  componente  esencial  del  
proceso  de  diseño  general.

Análisis  de  tensiones  inducidas

Cuando  se  excava  una  abertura  subterránea  en  un  macizo  rocoso  tensionado,  las  tensiones  en  
la  vecindad  de  la  nueva  abertura  se  redistribuyen.  Considere  el  ejemplo  de  las  tensiones  
inducidas  en  la  roca  que  rodea  un  túnel  circular  horizontal  como  se  ilustra  en  la  Figura  5,  que  
muestra  un  corte  vertical  normal  al  eje  del  túnel.
Antes  de  que  se  excave  el  túnel ,  las  tensiones  in  situ  y  se  drebanada  
istribuyen  
, due  
niformemente  
roca   en  la  
bajo  consideración.  Después  de  retirar  la  roca  del  interior  del  túnel,  las  tensiones  en  las  
inmediaciones  del  túnel  se  modifican  y

8
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

, tensiones  
se  inducen  nuevas  tensiones.  En  la  Figura  5  se  muestran  sobre  
tres   principales  
un  elemento  típico  dqe  
ue  
actúan  
roca  
σ  1  
2  σ  σ3 .

La  convención  utilizada  en  la  ingeniería  de  rocas  es  que  los  esfuerzos  de  compresión  son  
siempre  positivos  y  los  tres  esfuerzos  principales  están  numerados  de  tal  manera  
σ1 que  es  el  
mayor  esfuerzo  de  compresión  
σ3 y  el  menor  esfuerzo  de  compresión  o  el  mayor  esfuerzo  de  
tracción  de  los  tres.

Figura  5:  Ilustración  de  las  tensiones  principales  inducidas  en  un  elemento  de  roca  cerca  de  un  
túnel  horizontal  sujeto  a  una  tensión  in  situ  vertical  σv  en  un   , una  tensión  in  situ  horizontal  
plano  normal  al  eje  del  túnel  y  un  eje  de  túnel  de  tensión  in  situ  horizontal. la  σh2σh1  paralela  a  

9
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Figura  6:  Direcciones  principales  de  tensión  en  la  roca  que  rodea  un  túnel  horizontal  sujeta  a  una  tensión  in  situ  horizontal  
sigma  sub  h1  igual  Figura  6:  Direcciones  principales  de  tensión  en  la  roca  que  rodea  un  túnel  horizontal  sujeta  a  
a  a  3  sigma  sub  v,  donde  sigma  sub  v  es  el  esfuerzo  
vertical  in  situ.  esfuerzo  horizontal   v situ  igual  a  3  d,onde  
σh1  σin   es  el  esfuerzo  vertical  in  situ.  
σv

Figura  7:  Contornos  de  las  magnitudes  de  tensión  principal  máxima  y  mínima  en  la  roca  que  rodea  Figura  7:  Contornos  de  las  magnitudes  de  tensión  principal  máxima  y  

σ v tensión  
mínima  en  la  roca  un  túnel  horizontal,  sujeto  a  una  tensión  vertical  in  situ  de  sigma  sub  v  y  una  horizontal  que  rodea  una  horizontal  
túnel,  isn  
ometido  
situ  vertical  
a  una  
de  y  una  

σ3v .
tensión  in  situ  horizontal  de     sigma  sub  tensión  v.  situ  de  3

10
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Los   tres  
Figura   6  eqsfuerzos  
inclinados   ue  
a  
mla  
uestra  principales  
dirección   ston  
las  ddirecciones  
e  la   mutuamente  
ensión  de  
aplicada   pin  
la  tensión   erpendiculares  
saitu.  
plicada  
Esto   in  
es   epvidente  
situ. ero  pueden  
en  la  estar  
Esto  
las   es  uen  
rodea  
h1  
tensiones  
igual   vidente  
at  únel  
tres  
en   ea  
n  
vhleces  la  lFa  
orizontal  
roca   igura  
qtue   6,  qvaue  
ensión  
sujeto  
rodea   mtuestra  
ertical  
  un  
na  
únel  
tensión   las  
in  shitu   tensiones  
orizontal  
shigma  
orizontal  
ssub   esL
ujeto  
ivn  n  
.     lua  
aitu  
a   droca  
tn   que  
ensión  
irecciones  
sigma  in  
sub  
situ  
de  
horizontal  de  las  barras  más  largas  es  igual  a  tres  veces  la  tensión  in  situ  vertical.  
principal  
Las  
las   dσirecciones  
h1  mσáxima  
v  más   dse  
largas  
igma  
la  tensión  
seub  
n  e1sta  
,  
pmrincipal  
ientras  
figura  rm epresentan  
qáxima  
ue  las  m barras  
ientras  
las  edn  
irecciones  
qe ue  
sta  
las  figura  
barras  
de  
representan  
la  
mtensión  
ás  
cortas  
cada  σ 1  dban  
las   las  
arras   mdirecciones  
ás  cortas   de  la  tensión  principal  mínima  sigma  sub  3  en  
dan  
en  cada   las  direcciones  
elemento   de  la  tcensión  
σ3  elemento   principal  
onsiderado.   En  m ínima  
este  caso  
particular,  sigma  sub  2  es  coaxial  con  el  esfuerzo  in  situ  sigma  es  coaxial  con  
el  esfuerzo  in  situ  σh2 ,

consideró.   , otros  
En  este  caso  particular,  pero  los  otros  σ2  sub  h2  pero  los  
, principales  sigma  sub  1  y  sigma  sub  3  están  inclinados  a  los  esfuerzos  
esfuerzos  
principales  
σv  sub  h1  ys  sigma  
ys  eub  
igma   stán  
v  ein  
nclinados  hacia  y  en  
las  inmediaciones   dlel  
as  inmediaciones  de  σ1  σ3  σh1  
túnel.  el  tunel.

Los   contornos  
tensión  
se  
principal  
figura  
da  
en  
m
principal  
luestra  
a   de  
mFínima  
igura  
qm lLas  
ue   mcsta  
ragnitudes  
7áxima  
os  
.  
la  
E ontornos  
edistribución  
yf  igura   dde  
la  tensión  
m e  lda  
uestra  
las   tm
σe  ensión  
1  la  
m prrincipal  
qagnitudes  
tue  
ínima  
ensión  
la   m
sedistribución  
igma  
principal  
de  lsa   áxima  
ub   3  sde  stigma  
e  
d an   slub  
ensiones  
en   1igura  
a  Fσ   3  
y  sla  
tensión  
e  7 .  Esta  
concentra  
tres   e n  
tensiones  
digamos  
el  
centro  tdres   l a  
r oca  
el  avgujero,  
eces   c ercana  
se  cinsignificante.
campo  
oncentra  
el  
la  
radio   a l  
perturbación  
de  desde  t únel  
etensión  
n  la  e
roca  
l  dy   
q
cin   ue,  
entro  
el   a
cscercana  
itu  
ampo     
u
deel  na  
s  ia d istancia  
nsignificante.  
dagujero,  
e  
l  ttúnel  
ensión  
yla   d e  
  qpiue,  
n   d igamos  
erturbación  
vseces  
itu  
a  uena  
s  
el  
dristancia  
adio  
del  desde  
de  

Kirsch  (1898)  publicó  una  solución  analítica  para  la  distribución  de  esfuerzos  en  una  
placa  elástica  estresada  que  contenía  un  agujero  circular  y  esto  formó  la  base  de  muchos  
estudios  tempranos  del  comportamiento  de  las  rocas  alrededor  de  túneles  y  pozos.  
Siguiendo  el  camino  iniciado  por  Kirsch,  investigadores  como  Love  (1927),  Muskhelishvili  
(1953)  y  Savin  (1961)  publicaron  soluciones  para  excavaciones  de  varias  formas  en  
placas  elásticas.  Brown  publicó  un  resumen  útil  de  estas  soluciones  y  su  aplicación  en  
la  mecánica  de  rocas  en  una  introducción  a  un  volumen  titulado  Analytical  and  
Computational  Methods  in  Engineering  Rock  Mechanics  (1987).

Las  soluciones  de  forma  cerrada  aún  poseen  un  gran  valor  para  la  comprensión  conceptual  del  
comportamiento  y  para  la  prueba  y  calibración  de  modelos  numéricos.  Sin  embargo,  para  
propósitos  de  diseño,  estos  modelos  están  restringidos  a  geometrías  y  modelos  de  materiales  muy  simples.
Tienen  un  valor  práctico  limitado.  Afortunadamente,  con  el  desarrollo  de  las  computadoras,  
muchos  programas  poderosos  que  brindan  soluciones  numéricas  a  estos  problemas  ahora  están  
disponibles.  A  continuación  se  presenta  una  breve  revisión  de  algunas  de  estas  soluciones  numéricas.

Métodos  numéricos  de  análisis  de  tensiones.

La  mayoría  de  las  excavaciones  subterráneas  son  de  forma  irregular  y  con  frecuencia  se  
agrupan  cerca  de  otras  excavaciones.  Estos  grupos  de  excavaciones  pueden  formar  un  
conjunto  de  formas  tridimensionales  complejas.  Además,  debido  a  la  presencia  de  
características  geológicas  como  fallas  y  diques,  las  propiedades  de  las  rocas  rara  vez  son  
uniformes  dentro  del  volumen  de  roca  de  interés.  En  consecuencia,  las  soluciones  de  forma  
cerrada  tienen  un  valor  limitado  en  el  cálculo  de  las  tensiones,  desplazamientos  y  fallas  del  
macizo  rocoso  que  rodea  las  excavaciones  subterráneas.  En  las  últimas  décadas  se  han  
desarrollado  varios  métodos  numéricos  basados  en  computadora  y  estos  métodos  
proporcionan  los  medios  para  obtener  soluciones  aproximadas  a  estos  problemas.

11
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Los  métodos  numéricos  para  el  análisis  de  problemas  impulsados  por  esfuerzos  en  la  mecánica  de  rocas  
se  pueden  dividir  en  dos  clases:

Métodos  de  discretización  de  límites,  en  los  que  solo  el  límite  de  la  excavación  se  divide  en  
elementos  y  el  interior  del  macizo  rocoso  se  representa  matemáticamente  como  un  continuo  
infinito.  Estos  métodos  normalmente  se  restringen  a  análisis  elásticos.

Métodos  de  discretización  de  dominios,  en  los  que  el  interior  del  macizo  rocoso  se  divide  en  
elementos  geométricamente  simples,  cada  uno  con  propiedades  asumidas.  El  comportamiento  
colectivo  y  la  interacción  de  estos  elementos  simplificados  modelan  el  comportamiento  general  
más  complejo  del  macizo  rocoso.  En  otras  palabras,  los  métodos  de  dominio  permiten  la  
consideración  de  modelos  de  materiales  más  complejos  que  los  métodos  de  límites.  Los  métodos  
de  elementos  finitos  y  diferencias  finitas  son  técnicas  de  dominio  que  tratan  el  macizo  rocoso  
como  un  continuo.  El  método  de  elementos  distintos  es  también  un  método  de  dominio  que  
modela  cada  bloque  individual  de  roca  como  un  elemento  único.

Estas  dos  clases  de  análisis  se  pueden  combinar  en  forma  de  modelos  híbridos  para  maximizar  las  ventajas  
y  minimizar  las  desventajas  de  cada  método.

Es  posible  hacer  algunas  observaciones  generales  sobre  los  dos  tipos  de  enfoques  discutidos  anteriormente.  
En  los  métodos  de  dominio,  se  requiere  una  gran  cantidad  de  esfuerzo  para  crear  la  malla  que  se  utiliza  
para  dividir  el  macizo  rocoso  en  elementos.  En  el  caso  de  modelos  complejos,  como  los  que  contienen  
varias  aberturas,  el  mallado  puede  volverse  extremadamente  difícil.  Por  el  contrario,  los  métodos  de  límite  
solo  requieren  que  se  discretice  el  límite  de  la  excavación  y  que  el  macizo  rocoso  circundante  se  trate  
como  un  continuo  infinito.  Dado  que  se  requieren  menos  elementos  en  el  método  de  límites,  se  reduce  la  
demanda  de  memoria  de  la  computadora  y  de  la  habilidad  y  experiencia  del  usuario.  La  disponibilidad  de  
generadores  de  malla  altamente  optimizados  en  muchos  modelos  de  dominio  ha  reducido  esta  diferencia  
hasta  el  punto  en  que  la  mayoría  de  los  usuarios  de  programas  de  dominio  desconocerían  los  problemas  
de  generación  de  malla  discutidos  anteriormente  y,  por  lo  tanto,  la  elección  de  modelos  puede  basarse  en  
otras  consideraciones.

En  el  caso  de  los  métodos  de  dominio,  los  límites  exteriores  del  modelo  deben  colocarse  lo  suficientemente  
lejos  de  las  excavaciones  para  que  los  errores  que  surjan  de  la  interacción  entre  estos  límites  exteriores  y  
las  excavaciones  se  reduzcan  a  un  mínimo  aceptable.  Por  otro  lado,  dado  que  los  métodos  de  límite  tratan  
el  macizo  rocoso  como  un  continuo  infinito,  las  condiciones  de  campo  lejano  solo  necesitan  especificarse  
como  esfuerzos  que  actúan  sobre  todo  el  macizo  rocoso  y  no  se  requieren  límites  exteriores.  La  principal  
fortaleza  de  los  métodos  de  límites  radica  en  la  simplicidad  lograda  al  representar  el  macizo  rocoso  como  
un  continuo  de  extensión  infinita.  Sin  embargo,  es  esta  representación  la  que  dificulta  la  incorporación  de  
propiedades  y  discontinuidades  variables  del  material,  como  juntas  y  fallas.

Si  bien  se  han  desarrollado  técnicas  para  permitir  el  modelado  de  elementos  de  contorno  de  propiedades  
variables  de  la  roca,  este  tipo  de  problemas  se  modelan  más  convenientemente  mediante  métodos  de  
dominio.

12
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Antes  de  seleccionar  la  técnica  de  modelado  adecuada  para  determinados  tipos  de  problemas,  es  necesario  
comprender  los  componentes  básicos  de  cada  técnica.

Método  de  elementos  de  contorno

El  método  de  elementos  de  contorno  deriva  su  nombre  del  hecho  de  que  solo  los  límites  de  la  geometría  del  problema  
se  dividen  en  elementos.  En  otras  palabras,  solo  las  superficies  de  excavación,  la  superficie  libre  para  problemas  poco  
profundos,  las  superficies  de  junta  donde  las  juntas  se  consideran  explícitamente  y  las  interfaces  de  materiales  para  
problemas  de  materiales  múltiples  se  dividen  en  elementos.  De  hecho,  varios  tipos  de  modelos  de  elementos  de  
contorno  se  conocen  colectivamente  como  "método  de  elementos  de  contorno" (Crouch  y  Starfield,  1983).  Estos  
modelos  se  pueden  agrupar  de  la  siguiente  manera:

Método  indirecto  (esfuerzo  ficticio),  llamado  así  porque  el  primer  paso  en  la  solución  es  encontrar  un  conjunto  de  
esfuerzos  ficticios  que  satisfagan  las  condiciones  de  contorno  prescritas.  Estas  tensiones  se  utilizan  luego  en  el  cálculo  
de  las  tensiones  y  desplazamientos  reales  en  la  roca.
masa.

Método  directo,  llamado  así  porque  los  desplazamientos  se  resuelven  directamente  para  las  condiciones  de  contorno  
especificadas.

Método  de  discontinuidad  de  desplazamiento,  llamado  así  porque  la  solución  se  basa  en  la  superposición  de  la  
solución  fundamental  de  una  rendija  alargada  en  un  continuo  elástico  y  desplazamientos  cortantes  y  normales  en  la  
dirección  de  la  rendija.

Las  diferencias  entre  los  dos  primeros  métodos  no  son  evidentes  para  el  usuario  del  programa.
El  método  directo  tiene  ciertas  ventajas  en  términos  de  desarrollo  de  programas,  como  se  discutirá  más  adelante  en  la  
sección  sobre  enfoques  híbridos.

El  hecho  de  que  un  modelo  de  elementos  de  contorno  se  extienda  'hasta  el  infinito'  también  puede  ser  una  desventaja.
Por  ejemplo,  una  masa  rocosa  heterogénea  consta  de  regiones  de  extensión  finita,  no  infinita.  Se  deben  utilizar  
técnicas  especiales  para  manejar  estas  situaciones.  Las  juntas  se  modelan  explícitamente  en  el  método  de  elementos  
de  contorno  utilizando  el  enfoque  de  discontinuidad  de  desplazamiento,  pero  esto  puede  resultar  en  un  aumento  
considerable  en  el  esfuerzo  computacional.
A  menudo  se  encuentra  que  la  convergencia  numérica  es  un  problema  para  los  modelos  que  incorporan  muchas  
uniones.  Por  estas  razones,  los  problemas  que  requieren  la  consideración  explícita  de  varias  articulaciones  y/o  un  
modelado  sofisticado  del  comportamiento  constitutivo  de  las  articulaciones,  a  menudo  se  manejan  mejor  con  uno  de  
los  métodos  de  dominio,  como  los  elementos  finitos.

Una  aplicación  ampliamente  utilizada  de  los  elementos  de  contorno  de  discontinuidad  de  desplazamiento  es  el  
modelado  de  cuerpos  de  mineral  tabulares.  Aquí,  toda  la  veta  de  mineral  se  representa  como  una  "discontinuidad"  que  
inicialmente  se  llena  con  mineral.  La  minería  se  simula  mediante  la  reducción  de  la  rigidez  del  mineral  a  cero  en  
aquellas  áreas  donde  se  ha  producido  la  minería,  y  se  puede  examinar  la  redistribución  de  tensiones  resultante  a  los  
pilares  circundantes  (Salamon,  1974,  von  Kimmelmann  et  al.,  1984).

13
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Métodos  de  elementos  finitos  y  diferencias  finitas

En  la  práctica,  el  método  de  los  elementos  finitos  no  suele  distinguirse  del  método  de  las  diferencias  
finitas;  por  lo  tanto,  serán  tratados  aquí  como  uno  y  el  mismo.  Para  el  método  de  elementos  de  
contorno,  se  vio  que  las  condiciones  en  un  límite  de  dominio  podrían  estar  relacionadas  con  el  estado  
en  todos  los  puntos  a  lo  largo  de  la  roca  restante,  incluso  hasta  el  infinito.  En  comparación,  el  método  
de  elementos  finitos  relaciona  las  condiciones  en  algunos  puntos  dentro  de  la  roca  (puntos  nodales)  
con  el  estado  dentro  de  una  región  cerrada  finita  formada  por  estos  puntos  (el  elemento).
En  el  método  de  elementos  finitos,  el  problema  físico  se  modela  numéricamente  dividiendo  toda  la  
región  del  problema  en  elementos.

El  método  de  elementos  finitos  es  muy  adecuado  para  resolver  problemas  que  involucran  propiedades  
de  materiales  heterogéneos  o  no  lineales,  ya  que  cada  elemento  modela  explícitamente  la  respuesta  
del  material  que  contiene.  Sin  embargo,  los  elementos  finitos  no  son  adecuados  para  modelar  límites  
infinitos,  como  ocurre  en  los  problemas  de  excavación  subterránea.  Una  técnica  para  manejar  límites  
infinitos  es  discretizar  más  allá  de  la  zona  de  influencia  de  la  excavación  y  aplicar  condiciones  de  
contorno  apropiadas  a  los  bordes  exteriores.  Otro  enfoque  ha  sido  desarrollar  elementos  para  los  que  
un  borde  se  extiende  hasta  el  infinito,  es  decir,  los  llamados  elementos  finitos  'infinitos'.  En  la  práctica,  
los  preprocesadores  y  posprocesadores  eficientes  permiten  al  usuario  realizar  análisis  paramétricos  y  
evaluar  la  influencia  de  las  condiciones  límite  aproximadas  del  campo  lejano.  El  tiempo  requerido  para  
este  proceso  es  insignificante  en  comparación  con  el  tiempo  total  de  análisis.

Las  uniones  se  pueden  representar  explícitamente  usando  'elementos  de  unión'  específicos.  Se  han  
propuesto  diferentes  técnicas  para  manejar  tales  elementos,  pero  ninguna  técnica  única  ha  encontrado  
el  favor  universal.  Las  interfaces  conjuntas  se  pueden  modelar  utilizando  relaciones  constitutivas  
bastante  generales,  aunque  posiblemente  con  un  mayor  gasto  computacional  dependiendo  de  la  
técnica  de  solución.

Una  vez  que  el  modelo  se  ha  dividido  en  elementos,  se  han  asignado  las  propiedades  de  los  materiales  
y  se  han  prescrito  las  cargas,  se  debe  utilizar  alguna  técnica  para  redistribuir  las  cargas  desequilibradas  
y  así  determinar  la  solución  al  nuevo  estado  de  equilibrio.
Las  técnicas  de  solución  disponibles  se  pueden  dividir  ampliamente  en  dos  clases:  implícitas  y  
explícitas.  Las  técnicas  implícitas  ensamblan  sistemas  de  ecuaciones  lineales  que  luego  se  resuelven  
utilizando  técnicas  estándar  de  reducción  de  matrices.  Cualquier  no  linealidad  del  material  se  tiene  en  
cuenta  modificando  los  coeficientes  de  rigidez  (enfoque  de  la  secante)  y/o  ajustando  las  variables  
prescritas  (enfoque  de  tensión  inicial  o  deformación  inicial).  Estos  cambios  se  realizan  de  manera  
iterativa  de  modo  que  todas  las  ecuaciones  constitutivas  y  de  equilibrio  se  cumplan  para  el  estado  de  
carga  dado.

La  respuesta  de  un  sistema  no  lineal  generalmente  depende  de  la  secuencia  de  carga.
Por  lo  tanto,  es  necesario  que  la  trayectoria  de  carga  modelada  sea  representativa  de  la  trayectoria  
de  carga  real  experimentada  por  el  cuerpo.  Esto  se  logra  dividiendo  la  carga  total  aplicada  en  
incrementos  de  carga,  siendo  cada  incremento  lo  suficientemente  pequeño,  de  modo  que  la  
convergencia  de  la  solución  para  el  incremento  se  logre  después  de  solo  unas  pocas  iteraciones.  Sin  
embargo,  a  medida  que  el  sistema  que  se  modela  se  vuelve  cada  vez  más  no  lineal  y  el  incremento  de  carga  representa  un
14
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

porción  más  pequeña  de  la  carga  total,  la  técnica  de  solución  incremental  se  vuelve  similar  a  modelar  el  
comportamiento  cuasi­dinámico  del  cuerpo,  ya  que  responde  a  la  aplicación  gradual  de  la  carga  total.

Para  superar  esto,  se  propuso  una  técnica  de  solución  de  "relajación  dinámica" (Otter  et  al.,  1966)  y  
Cundall  (1971)  la  aplicó  por  primera  vez  al  modelado  geomecánico.  En  esta  técnica  no  se  forman  
matrices.  Más  bien,  la  solución  procede  explícitamente:  las  fuerzas  desequilibradas,  que  actúan  en  un  
punto  de  integración  material,  dan  como  resultado  la  aceleración  de  la  masa  asociada  con  el  punto;  la  
aplicación  de  la  ley  de  movimiento  de  Newton  expresada  como  una  ecuación  en  diferencia  produce  
desplazamientos  incrementales,  la  aplicación  de  la  relación  constitutiva  apropiada  produce  el  nuevo  
conjunto  de  fuerzas,  y  así  sucesivamente,  marchando  en  el  tiempo,  para  cada  punto  de  integración  
material  en  el  modelo.  Esta  técnica  de  solución  tiene  la  ventaja  de  que  se  acomodan  las  no  linealidades  
tanto  geométricas  como  materiales,  con  un  esfuerzo  computacional  adicional  relativamente  pequeño  en  
comparación  con  un  análisis  lineal  correspondiente,  y  el  gasto  computacional  aumenta  solo  linealmente  
con  el  número  de  elementos  utilizados.
Otra  ventaja  práctica  reside  en  el  hecho  de  que  la  divergencia  numérica  normalmente  da  como  resultado  
que  el  modelo  prediga  un  comportamiento  físico  obviamente  anómalo.  Por  lo  tanto,  incluso  los  usuarios  
relativamente  inexpertos  pueden  reconocer  la  divergencia  numérica.

La  mayoría  de  los  paquetes  de  elementos  finitos  disponibles  en  el  mercado  utilizan  técnicas  de  solución  
implícitas  (es  decir,  matriciales).  Para  problemas  lineales  y  problemas  de  no  linealidad  moderada,  las  
técnicas  implícitas  tienden  a  funcionar  más  rápido  que  las  técnicas  de  solución  explícita.  Sin  embargo,  
a  medida  que  aumenta  el  grado  de  no  linealidad  del  sistema,  las  cargas  impuestas  deben  aplicarse  en  
incrementos  más  pequeños,  lo  que  implica  un  mayor  número  de  reformas  y  reducciones  de  la  matriz  y,  
por  lo  tanto,  un  mayor  gasto  de  cálculo.  Por  lo  tanto,  los  problemas  altamente  no  lineales  se  manejan  
mejor  con  paquetes  que  utilizan  una  técnica  de  solución  explícita.

Método  de  elementos  distintos

En  condiciones  de  terreno  convencionalmente  descritas  como  bloques  (es  decir,  donde  el  espaciamiento  
de  las  juntas  es  del  mismo  orden  de  magnitud  que  las  dimensiones  de  la  excavación),  las  juntas  que  se  
cruzan  forman  cuñas  de  roca  que  pueden  considerarse  como  cuerpos  rígidos.  Es  decir,  estas  piezas  
individuales  de  roca  pueden  rotar  y  trasladarse  libremente,  y  la  deformación  que  tiene  lugar  en  los  
contactos  del  bloque  puede  ser  significativamente  mayor  que  la  deformación  de  la  roca  intacta.  Por  lo  
tanto,  las  cuñas  individuales  pueden  considerarse  rígidas.  Para  tales  condiciones,  normalmente  es  
necesario  modelar  explícitamente  muchas  articulaciones.  Sin  embargo,  el  comportamiento  de  tales  
sistemas  es  tan  altamente  no  lineal  que  incluso  un  código  de  elementos  finitos  unidos,  que  emplea  una  
técnica  de  solución  explícita,  puede  funcionar  de  manera  relativamente  ineficiente.

Un  enfoque  de  modelado  alternativo  es  desarrollar  estructuras  de  datos  que  representen  la  naturaleza  
en  bloques  del  sistema  que  se  analiza.  Cada  bloque  se  considera  un  cuerpo  libre  único  que  puede  
interactuar  en  lugares  de  contacto  con  los  bloques  circundantes.  Los  contactos  pueden  representarse  
mediante  la  superposición  de  bloques  adyacentes,  evitando  así  la  necesidad  de  elementos  de  unión  
únicos.  Esto  tiene  la  ventaja  adicional  de  que  pueden  ocurrir  desplazamientos  relativos  arbitrariamente  
grandes  en  el  contacto,  una  situación  que  generalmente  no  es  tratable  en  códigos  de  elementos  finitos.

15
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Debido  al  alto  grado  de  no  linealidad  de  los  sistemas  que  se  modelan,  se  favorecen  las  técnicas  de  
solución  explícita  para  códigos  de  elementos  distintos.  Como  es  el  caso  de  los  códigos  de  elementos  
finitos  que  emplean  técnicas  de  solución  explícita,  esto  permite  un  modelado  constitutivo  muy  general  del  
comportamiento  conjunto  con  un  pequeño  aumento  en  el  esfuerzo  computacional  y  da  como  resultado  que  
el  tiempo  de  cálculo  dependa  solo  linealmente  del  número  de  elementos  utilizados.  El  uso  de  técnicas  de  
solución  explícita  exige  menos  habilidades  y  experiencia  que  el  uso  de  códigos  que  emplean  técnicas  de  
solución  implícita.

Aunque  el  método  de  elementos  distintivos  se  ha  utilizado  más  ampliamente  en  entornos  académicos  
hasta  la  fecha,  se  está  abriendo  camino  en  las  oficinas  de  consultores,  planificadores  y  diseñadores.  Se  
requiere  más  experiencia  en  la  aplicación  de  esta  poderosa  herramienta  de  modelado  a  situaciones  
prácticas  de  diseño  y  la  subsiguiente  documentación  de  estas  historias  de  casos,  de  modo  que  se  pueda  
desarrollar  una  comprensión  de  dónde,  cuándo  y  cómo  se  aplica  mejor  el  método  de  elementos  distintos.

Enfoques  híbridos

El  objetivo  de  un  método  híbrido  es  combinar  los  métodos  anteriores  para  eliminar  las  características  
indeseables  y  conservar  tantas  ventajas  como  sea  posible.  Por  ejemplo,  al  modelar  una  excavación  
subterránea,  la  mayor  parte  de  la  no  linealidad  ocurrirá  cerca  del  límite  de  la  excavación,  mientras  que  el  
macizo  rocoso  a  cierta  distancia  se  comportará  de  manera  elástica.  Por  lo  tanto,  el  macizo  rocoso  de  
campo  cercano  podría  modelarse  utilizando  un  método  de  elementos  distintos  o  de  elementos  finitos,  que  
luego  se  vincula  en  sus  límites  exteriores  a  un  modelo  de  elementos  de  contorno,  de  modo  que  las  
condiciones  de  contorno  de  campo  lejano  se  modelen  exactamente.  En  tal  enfoque,  se  favorece  la  técnica  
de  elemento  de  contorno  directo,  ya  que  da  como  resultado  una  mayor  eficiencia  de  programación  y  
solución.

Lorig  y  Brady  (1984)  utilizaron  un  modelo  híbrido  que  consiste  en  un  modelo  de  elementos  discretos  para  
el  campo  cercano  y  un  modelo  de  elementos  de  contorno  para  el  campo  lejano  en  una  masa  rocosa  que  
rodea  un  túnel  circular.

Modelos  bidimensionales  y  tridimensionales.

Se  puede  utilizar  un  modelo  bidimensional,  como  el  ilustrado  en  la  Figura  5,  para  el  análisis  de  tensiones  
y  desplazamientos  en  la  roca  que  rodea  un  túnel,  pozo  o  pozo,  donde  la  longitud  de  la  abertura  es  mucho  
mayor  que  las  dimensiones  de  su  sección  transversal. .  Los  esfuerzos  y  desplazamientos  en  un  plano,  
normal  al  eje  de  la  abertura,  no  están  influenciados  por  los  extremos  de  la  abertura,  siempre  que  estos  
extremos  estén  lo  suficientemente  alejados.

Por  otro  lado,  una  casa  de  máquinas  subterránea  o  cámara  trituradora  tiene  una  forma  mucho  más  
equidimensional  y  el  efecto  de  las  paredes  de  los  extremos  no  se  puede  despreciar.  En  este  caso,  es  
mucho  más  apropiado  realizar  un  análisis  tridimensional  de  las  tensiones  y  desplazamientos  en  el  macizo  
rocoso  circundante.  Desafortunadamente,  este  cambio  de  dos  a  tres  dimensiones  no  es  tan  simple  como  
parece  y  hay  relativamente  pocas  buenas  tres  dimensiones.

dieciséis
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

modelos  numéricos  dimensionales,  que  son  adecuados  para  el  trabajo  de  análisis  de  estrés  de  rutina  en  una  
oficina  de  diseño  de  ingeniería  típica.

EXAMINE3D  (www.rocscience.com)  es  un  programa  de  elemento  de  contorno  tridimensional  que  proporciona  
un  punto  de  partida  para  el  análisis  de  un  problema  en  el  que  la  geometría  tridimensional  de  las  aberturas  
es  importante.  Dichos  análisis  tridimensionales  proporcionan  indicaciones  claras  de  las  concentraciones  de  
tensión  y  de  la  influencia  de  la  geometría  tridimensional.  En  muchos  casos,  es  posible  simplificar  el  problema  
a  dos  dimensiones  considerando  las  tensiones  en  las  secciones  críticas  identificadas  en  el  modelo  
tridimensional.

Modelos  tridimensionales  de  elementos  finitos  más  sofisticados  como  FLAC3D  (www.itascacg.com)  están  
disponibles,  pero  la  definición  de  los  parámetros  de  entrada  y  la  interpretación  de  los  resultados  de  estos  
modelos  ampliarían  las  capacidades  de  todos  menos  los  modeladores  más  experimentados.  Probablemente  
lo  mejor  sea  dejar  este  tipo  de  modelado  en  manos  de  estos  especialistas.

Se  recomienda  que,  cuando  el  problema  que  se  está  considerando  sea  obviamente  tridimensional,  se  lleve  
a  cabo  un  análisis  elástico  preliminar  por  medio  de  uno  de  los  programas  de  elementos  de  contorno  
tridimensionales.  Luego,  los  resultados  se  pueden  usar  para  decidir  si  se  requieren  más  análisis  
tridimensionales  o  si  se  pueden  modelar  secciones  bidimensionales  apropiadas  usando  un  programa  como  
PHASE2  (www.rocscience.com),  un  programa  de  elementos  finitos  poderoso  pero  fácil  de  usar  que  
generalmente  cumple  con  las  necesidades  de  la  mayoría  de  los  proyectos  de  diseño  de  excavaciones  
subterráneas.

Ejemplos  de  análisis  de  tensión  bidimensional

Un  programa  de  elementos  de  contorno  llamado  EXAMINE2D  está  disponible  como  descarga  gratuita  desde  
www.rocscience.com.  Si  bien  este  programa  se  limita  a  los  análisis  elásticos,  puede  proporcionar  una  
introducción  muy  útil  para  aquellos  que  no  están  familiarizados  con  los  métodos  numéricos  de  análisis  de  
tensión  descritos  anteriormente.  Los  siguientes  ejemplos  demuestran  el  uso  de  este  programa  para  explorar  
algunos  problemas  comunes  en  la  construcción  de  túneles.

forma  de  túnel

A  la  mayoría  de  los  contratistas  les  gusta  una  forma  de  herradura  simple  para  los  túneles,  ya  que  proporciona  
un  piso  amplio  y  plano  para  el  equipo  utilizado  durante  la  construcción.  Para  túneles  relativamente  poco  
profundos  en  roca  de  buena  calidad,  esta  es  una  forma  de  túnel  apropiada  y  hay  muchos  cientos  de  kilómetros  
de  túneles  en  forma  de  herradura  en  todo  el  mundo.

En  macizos  rocosos  de  mala  calidad  o  en  túneles  de  gran  profundidad,  la  forma  de  herradura  simple  no  es  
una  buena  opción  debido  a  las  altas  concentraciones  de  esfuerzos  en  las  esquinas  donde  las  paredes  
laterales  se  encuentran  con  el  piso  o  se  invierten.  En  algunos  casos,  las  fallas  que  se  inician  en  estas  
esquinas  pueden  provocar  un  fuerte  levantamiento  del  piso  e  incluso  la  falla  de  todo  el  perímetro  del  túnel,  
como  se  muestra  en  la  Figura  8.

17
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Figura  8:  Falla  del  revestimiento  en  un  túnel  en  forma  de  herradura  en  un  macizo  rocoso  de  mala  calidad  y  
sometido  a  grandes  esfuerzos.  Esta  falla  se  inició  en  las  esquinas  donde  la  solera  se  encuentra  con  las  
paredes  laterales.

Figura  9:  Dimensiones  de  un  túnel  de  herradura  
modificado  de  10  m  de  luz  diseñado  para  
superar  algunos  de  los  problemas  ilustrados  en  
la  Figura  8.

La  distribución  de  tensiones  en  el  macizo  rocoso  que  rodea  el  túnel  se  puede  mejorar  modificando  la  forma  
de  herradura  como  se  muestra  en  la  Figura  9.  En  algunos  casos,  esto  puede  eliminar  o  minimizar  los  tipos  
de  falla  que  se  muestran  en  la  Figura  8,  mientras  que,  en  otros  casos,  puede  ser  necesario  utilizar  un  perfil  
de  túnel  circular.
18
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Tensiones  in  situ:

Esfuerzo  principal  mayor  sigma  sub  1  =  10  MPa  Esfuerzo  principal  menor  sigma  sub  1  =  10  MPa  

Esfuerzo  principal  sigma  sub  3  =  7  MPa  Esfuerzo  principal  menor  intermedio  σ3  =  7  MPa  

9  MPa
Esfuerzo  principal  sigma  sub  2  =  9  MPa  Esfuerzo  principal  intermedio  σ2  =

Inclinación  de  la  tensión  principal  principal  hacia  la  horizontal  Inclinación  de  la  tensión  

principal  
hacia  epl  rincipal  
eje  =  15eje  horizontal  =  15º

Propiedades  del  macizo  rocoso:

Ángulo  de  fricción  phi  =  35   Cohesión  c  =  1  MPa  Ángulo  de  fricción  Φ  =  35º

Resistencia  a  la  tracción  =  cero  Módulo  de  deformación  Cohesión  c  =  1  MPa

E  
=  4600  MPa  
Resistencia  
a  la  tracción  =  cero
Módulo  de  deformación  E  =  4600  MPa

Figura  10:  Comparación  de  tres  perfiles  de  
excavación  de  túneles  usando  EXAMINE2D.
Los  contornos  son  para  el  Factor  de  Resistencia  
definido  por  la  relación  entre  la  resistencia  del  
macizo  rocoso  y  la  tensión  inducida  en  cada  
punto.  El  perfil  límite  deformado  (exagerado)  
se  muestra  dentro  de  cada  excavación.

19
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

La  aplicación  del  programa  EXAMINE2D  para  comparar  tres  formas  de  túneles  se  ilustra  en  la  Figura  
10.  En  este  análisis  se  usaron  las  típicas  tensiones  in  situ  "promedio"  y  las  propiedades  del  macizo  
rocoso  y  las  tres  figuras  comparan  los  contornos  del  Factor  de  Resistencia  y  los  perfiles  de  excavación  
deformados  (exagerados)  para  el  tres  formas  de  túnel.

Está  claro  que  el  piso  plano  del  túnel  de  herradura  (figura  superior)  permite  el  desplazamiento  hacia  
arriba  o  el  levantamiento  del  piso.  Las  esquinas  afiladas  en  la  unión  entre  el  piso  y  las  paredes  laterales  
del  túnel  crean  altas  concentraciones  de  tensión  y  también  generan  grandes  momentos  de  flexión  en  
cualquier  revestimiento  instalado  en  el  túnel.  La  falla  del  piso  generalmente  se  inicia  en  estas  esquinas,  
como  se  ilustra  en  la  Figura  8.

El  levantamiento  del  piso  se  reduce  significativamente  por  la  curvatura  cóncava  del  piso  de  la  forma  de  
herradura  modificada  (figura  central).  En  casos  marginales,  estas  modificaciones  a  la  forma  de  
herradura  pueden  ser  suficientes  para  prevenir  o  al  menos  minimizar  el  tipo  de  daño  ilustrado  en  la  
Figura  8.  Sin  embargo,  en  casos  severos,  un  perfil  de  túnel  circular  es  invariablemente  la  mejor  opción,  
como  lo  muestra  el  factor  de  resistencia  uniforme.  contornos  y  la  forma  del  límite  del  túnel  deformado  
en  la  figura  inferior  de  la  Figura  10.

Grandes  cavernas  subterráneas

Un  complejo  subterráneo  típico  en  un  proyecto  hidroeléctrico  tiene  una  casa  de  máquinas  con  una  luz  
de  20  a  25  m  y  una  altura  de  40  a  50  m.  En  esta  caverna  se  alojan  de  cuatro  a  seis  conjuntos  de  
turbinas  y  generadores  y  en  la  Figura  11  se  muestra  un  boceto  en  corte  de  uno  de  estos  conjuntos.
Los  transformadores  se  alojan  con  frecuencia  en  una  cámara  o  galería  paralela  a  la  casa  de  máquinas.
Idealmente,  estas  dos  cavernas  deberían  estar  lo  más  cerca  posible  para  minimizar  la  longitud  de  las  
barras  colectoras  que  conectan  los  generadores  y  transformadores.  Esto  tiene  que  equilibrarse  con  el  
tamaño  y,  por  lo  tanto,  con  la  estabilidad  del  pilar  entre  las  cavernas.  La  ubicación  relativa  y  la  distancia  
entre  las  cavernas  se  exploran  en  la  serie  de  modelos  EXAMINE2D  que  se  muestran  en  la  Figura  12,  
utilizando  las  mismas  propiedades  del  macizo  rocoso  y  las  mismas  tensiones  in  situ  que  se  enumeran  
en  la  Figura  10.

Figura  11:  Croquis  en  corte  del  diseño  
de  una  caverna  subterránea  
de  una  central  eléctrica  y  una  galería  
paralela  de  transformadores.

20
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Tensiones  in  situ:

=  10  MPa
Tensión  principal  mayor  σ1  =  10  Tensión  principal  mayor  sigma  sub  1

Tensión  principal  menor  sigma  sub  3  =  7  MPa
MPa  Tensión  intermedia  sigma  sub  2  =  9  Tensión  principal  menor  σ3  =  7

MPa  Inclinación  de  la  tensión  principal  principal  
a  MPa  el  eje  horizontal  =  15
Esfuerzo  intermedio  σ2  =  9  MPa
Inclinación  de  la  tensión  principal  
mayor  al  eje  horizontal  =  15º

Propiedades  del  macizo  rocoso:

Ángulo  de  fricción  Phi  =  35   Cohesión  Ángulo  de  fricción  Φ  =  35º

=  1t  racción
MPa  
Cohesión  
a  
la   cR
  =esistencia  
  1  MPa  c
=  deformación
Resistencia  
cero  =  cero  a
M  lódulo  
a  tracción  
de  
Módulo  de  deformación  E  =  4600
E   =  4600  MPa
MPa

Figura  12:  Comparación  de  tres  
diseños  subterráneos  de  galerías  
de  transformadores  y  centrales  
eléctricas,  usando  EXAMINE2D.  
Los  contornos  son  para  el  Factor  
de  Resistencia  definido  por  la  
relación  entre  la  resistencia  del  
macizo  rocoso  y  la  tensión  
inducida  en  cada  punto.  El  perfil  
límite  deformado  (exagerado)  se  
muestra  dentro  de  cada  
excavación.

21
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Figura  13:  Vectores  de  desplazamiento  y  formas  de  excavación  deformadas  para  la  casa  de  máquinas  
subterránea  y  la  galería  de  transformadores.

Un  examen  más  detallado  de  las  deformaciones  inducidas  en  el  macizo  rocoso  por  la  excavación  de  la  
casa  de  máquinas  subterránea  y  la  galería  del  transformador,  en  la  Figura  13,  muestra  que  la  más  
pequeña  de  las  dos  excavaciones  es  atraída  hacia  la  caverna  más  grande  y  su  perfil  se  distorsiona  en  
este  proceso.  Esta  distorsión  se  puede  reducir  reubicando  la  galería  del  transformador  y  aumentando  
el  espacio  entre  las  galerías  como  se  ha  hecho  en  la  Figura  12.

Cuando  es  probable  que  la  combinación  de  la  resistencia  del  macizo  rocoso  y  los  esfuerzos  in  situ  
causen  esfuerzos  excesivos  alrededor  de  las  cavernas  y  en  el  pilar,  una  buena  regla  general  es  que  la  
distancia  entre  las  dos  cavernas  debe  ser  aproximadamente  igual  a  la  altura  de  la  caverna  más  grande.

Se  anima  al  lector  interesado  a  descargar  el  programa  EXAMINE2D  (gratis  de  www.rocscience.com)  y  
usarlo  para  explorar  el  problema,  como  los  que  se  ilustran  en  las  Figuras  10  y  12,  por  sí  mismos.

22
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

Referencias

Marrón,  ET  1987.  Introducción.  Métodos  analíticos  y  computacionales  en  ingeniería  mecánica  de  rocas,  
(ed.  ET  Brown),  1­31.  Londres:  Allen  y  Unwin.
Brown,  ET  y  Hoek,  E.  1978.  Tendencias  en  las  relaciones  entre  las  tensiones  in  situ  de  las  rocas  medidas  
y  la  profundidad.  En  t.  J.  roca  mecánica.  mín.  ciencia  y  Geomecánica.  Resumen  15,  págs.  211­215.
Crouch,  SL  y  Starfield,  AM  1983.  Métodos  de  elementos  de  contorno  en  mecánica  de  sólidos
.  Londres:  Allen  y  Unwin.

Cundall,  PA  1971.  Un  modelo  de  computadora  para  simular  movimientos  progresivos  a  gran  escala  en  
,
sistemas  de  rocas  en  bloques.  Proceso  de  fractura  en  roca.  sim.  ISRM,  
Nancy  1,  Documento  2­8.

Engelder,  T.  y  Sbar,  ML  1984.  Estrés  in  situ  cerca  de  la  superficie:  introducción.  J.  Geophys.
Res.  89,  págs.  9321­9322.  Princeton,  Nueva  Jersey:  Princeton  University  Press.
Herget,  G.  1988.  Tensiones  en  roca.  Róterdam:  Balkema.
Hoek,  E.,  Carranza  –  Torres,  C.  and  Corkum,  B.,  2002.  Hoek  ­  Criterio  de  falla  de  Brown  –  edición  2002.  
En  Actas  de  NARMS­TAC  2002,  Toronto  (eds.  Bawden,RW,  Curran,  J.,  Telesnicki,  M)  págs.  267­273.  
Descargar  de  www.rocscience.com  Kirsch,  G.,  1898.  Die  theorie  der  elastizitat  und  die  
bedurfnisse  der  festigkeitslehre.
Veit.  Deit.  En  g.  42  (28),  797­807.
Lorig,  LJ  y  Brady,  BHG  1984.  Un  esquema  computacional  híbrido  para  el  diseño  de  excavación  y  soporte  
en  medios  rocosos  articulados.  En  Diseño  y  ejecución  de  excavaciones  subterráneas,  (eds  ET  Brown  
y  JA  Hudson),  105­112.  Londres:  británico.  Geotecnología.  Soc.

Amor,  AEH  1927.  Tratado  sobre  la  teoría  matemática  de  la  elasticidad.  Nueva  York:
Dover.

Muskhelishvili,  NI  1953.  Algunos  problemas  básicos  de  la  teoría  matemática  de  la  elasticidad.  4ª  ed.,  
traducida  por  JRM  Radok.  Gronigen:  Noordhoff.
Otter,  JRH,  Cassell,  AC  y  Hobbs,  RE  1966.  Relajación  dinámica.  proc.  Instn  Civ.  Ings  35,  633­665.

Reinecker,  J.,  Heidbach,  O.,  Tingay,  M.,  Sperner,  B.  y  Müller,  B.  2005:  The  release  2005  of  the  World  
Stress  Map  (disponible  en  línea  en  www.world­stress­map.org).
Salamon,  MDG  1974.  Mecánica  de  rocas  de  excavaciones  subterráneas.  En  Avances  en  Proc.  3er  
mecánica  de  rocas  , Congreso  ISRM.,  Denver  1B,  951­1009.  Washington  DC:
Academia  Nacional  de  Ciencias  
Savin,  GN  1961.  Concentraciones  de  tensión  alrededor  de  agujeros.  Londres:  Pérgamo.
Sheory,  PR  1994.  Una  teoría  para  tensiones  in  situ  en  isótropos  y  transversalmente  isótropos
roca.  En  t.  J.  roca  mecánica.  mín.  ciencia  y  Geomecánica.  Resumen  31(1),  23­34.

Terzaghi,  K.  y  Richart,  FE  1952.  Tensiones  en  roca  sobre  cavidades.  geotécnica  3,
57­90.

23
Machine Translated by Google

Esfuerzos  in  situ  e  inducidos

von  Kimmelmann,  MR,  Hyde,  B.  y  Madgwick,  RJ  1984.  El  uso  de  aplicaciones  informáticas  en  
BCL  Limited  en  la  planificación  de  la  extracción  de  pilares  y  el  diseño  de  diseños  de  minas.  
En  Diseño  y  ejecución  de  excavaciones  subterráneas,  (eds  ET  Brown  y  JA  Hudson),  53­64.  
Londres:  británico.  Geotecnología.  Soc.
Zoback,  ML  1992.  Patrones  de  tensión  de  primer  y  segundo  orden  en  la  litosfera:  la
Proyecto  Mapa  Mundial  de  Estrés.  J.  Geophys.  Res.  97(B8),  11761­11782.

24