Propiedades Físicas y Mecánicas de los

Materiales Rocosos

Jorge Cortez Campaña 2017

 Introducción

Objetivos Generales

Habilitar al alumno en los métodos de análisis y diseño tanto de excavaciones

subterráneas, taludes rocosos y estructuras mineras.

Objetivos Específicos

- Determinar los esfuerzos inducidos y desplazamientos en excavaciones de

geometría simple tanto en rocas masivas como estratificadas, mediante técnicas
analíticas y numéricas.
- Diseñar pilares tanto en rocas masivas como estratificadas.
- Presentar criterios de soporte activos en rocas masivas y estratificadas.
- Presentar métodos de análisis de la estabilidad de taludes rocosos considerando
sus modos de fractura.
Introducción

Fuente: Apuntes Marek Cala, AGH.

Introducción
Capitulo 1: Conceptos Generales

Repaso de los conceptos vistos en mecánica de sólidos

• Esfuerzo y nomenclatura
• Tensor de esfuerzos
• Esfuerzos en 2D y 3D
• Ecuaciones de Equilibrio
• Deformaciones en 2D y 3D
• Roseta de Deformación
• Estados de esfuerzo plano y deformación plana
• Determinación de las constantes de elasticidad (E,λ,K,G,ν)
• Ecuaciones constitutivas de la teoría de la elasticidad
¿Qué fenómeno vamos a entender?
 Esfuerzos

∆𝑭
𝝈𝒓 = 𝐥𝐢𝐦
∆𝑨→𝟎 ∆𝑨

¿Qué implicaciones tiene que por definición

sea un límite?

El concepto de esfuerzo es usado para describir la intensidad de las fuerzas

internas en un cuerpo bajo la influencia de un conjunto de fuerzas en su
superficie.
Esfuerzos Normales y Tangenciales

∆𝑵 + ∆𝑻 = ∆𝑭

∆𝑵
𝝈𝒏 = 𝐥𝐢𝐦
∆𝑨→𝟎 ∆𝑨
∆𝑇
∆𝑵 ∆𝑻
𝝈𝒕 = 𝐥𝐢𝐦
∆𝑨→𝟎 ∆𝑨
 Esfuerzos en 2D

𝝈𝒕 𝝈𝒏
𝜎𝑥𝑥 + 𝜎𝑦𝑦 𝜎𝑥𝑥 − 𝜎𝑦𝑦 𝝈𝒙𝒚
𝜎𝑛 = + cos 2𝜃 + 𝜏𝑥𝑦 sin 2𝜃
2 2 𝝈𝒙𝒙 𝜃
𝜎𝑦𝑦 − 𝜎𝑥𝑥
𝜎𝑡 = sin 2𝜃 + 𝜏𝑥𝑦 cos 2𝜃
2

DEMOSTREMOSLO… 𝝈𝒙𝒚
En la doble notación de índices para los componentes de esfuerzo, el primer sub
𝝈𝒚𝒚
índice indica la dirección de la normal de la superficie y el segundo el sentido de
acción de la componente de esfuerzo actuando sobre esa superficie.
Convención de Signos
 Esfuerzos en 2D

Máximo de la función
de Esfuerzo Normal

𝝈𝒙𝒙 = 𝟏𝟎 𝑴𝑷𝒂

𝝈𝒙𝒚 = 𝟓 𝑴𝑷𝒂

La función
𝝈𝒙𝒚 = 𝟓 𝑴𝑷𝒂
de Esfuerzo
Tangencial toma el
valor de 0
𝝈𝒚𝒚 = 𝟖 𝑴𝑷𝒂
 Esfuerzos en 2D

𝜎𝑥 + 𝜎𝑦 𝜎𝑥 − 𝜎𝑦
𝜎𝑛 (𝜃) = + cos 2𝜃 + 𝜏𝑥𝑦 sin 2𝜃
2 2
𝜎𝑦 − 𝜎𝑥
𝜎𝑡 (𝜃) = sin 2𝜃 + 𝜏𝑥𝑦 cos 2𝜃
2
 Esfuerzos en 2D

Esfuerzo Normal
Mínimo de 4 MPa
Esfuerzo Normal
Máximo de 14 MPa

Esfuerzo Normal
Máximo de 14 MPa

Esfuerzo Normal
Mínimo de 4 MPa
Esfuerzos en 2D
Esfuerzos Principales

Valores Propios o Eigenvalores

𝜎𝑥𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧
𝜎 = 𝜎𝑦𝑥 𝜎𝑦𝑦 𝜎𝑦𝑧 𝐝𝐞𝐭 𝝈 − 𝝀𝑰 = 𝟎
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧𝑧
 Esfuerzos Principales

10 5
𝜎=
5 8
𝝈𝒙𝒙 = 𝟏𝟎 𝑴𝑷𝒂
10 5 1 0
𝐝𝐞𝐭 −𝜆 =𝟎
5 8 0 1
𝝈𝒙𝒚 = 𝟓 𝑴𝑷𝒂

𝝈𝒙𝒚 = 𝟓 𝑴𝑷𝒂
𝜆1 = 3.9
𝜆2 = 14.1
𝝈𝒚𝒚 = 𝟖 𝑴𝑷𝒂
Repaso de Determinantes
Determinantes
 Direcciones Principales
Vectores Propios o Eigenvectores

𝝈𝒙𝒚
𝝈𝒙𝒙

𝝈𝒙𝒚

𝝈𝒚𝒚
 Direcciones Principales

𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦
a = det 𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1
𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦 1
det 𝑟1 = 𝑎 𝑖Ƹ + 𝑏 𝑗 Ƹ
𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1 𝑎2 + 𝑏2

𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦
b = −det 𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1

El mismo proceso puede ser llevado para determinar

la dirección principal del esfuerzo menor
Tensor de Esfuerzos en 3D

𝜎𝑥𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧

𝜎 = 𝜎𝑦𝑥 𝜎𝑦𝑦 𝜎𝑦𝑧
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧𝑧

𝐝𝐞𝐭 𝝈 − 𝝀𝑰 = 𝟎

Valores Propios o Eigenvalores

𝝈𝒊𝒋 el primer índice indica la dirección de la normal a la superficie

y el segundo índice indica el sentido de acción del esfuerzo
Rotación de Tensores
Donde 𝒍, 𝒎 y 𝒏 son vectores unitarios en
coordenadas cartesianas ortogonales entre si

𝑧 𝑚

La ortogonalidad se verifica
realizando el producto punto
𝑥 entre cada par de vectores

𝑦 𝑛
 Rotación de Tensores 𝑙𝑥 𝑙𝑦 𝑙𝑧
𝑅 = 𝑚𝑥 𝑚𝑦 𝑚𝑧
𝑛𝑥 𝑛𝑦 𝑛𝑧

𝜎𝑥𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧 𝜎𝑙𝑙 𝜎𝑙𝑚 𝜎𝑙𝑛

𝜎 𝑥𝑦𝑧 = 𝜎𝑦𝑥 𝜎𝑦𝑦 𝜎𝑦𝑧 𝜎 𝑙𝑚𝑛 = 𝜎𝑚𝑙 𝜎𝑚𝑚 𝜎𝑚𝑛
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧𝑧 𝜎𝑛𝑙 𝜎𝑛𝑚 𝜎𝑛𝑛

𝜎 𝑙𝑚𝑛 = 𝑅𝜎 𝑥𝑦𝑧 𝑅𝑇
 Direcciones Principales
Vectores Propios o Eigenvectores

𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦 𝜎𝑥𝑧

𝑎𝑖Ƹ + 𝑏𝑗Ƹ + 𝑐𝑘෠

a = det 𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1 𝜎𝑦𝑧
𝜎𝑥𝑧 𝜎𝑦𝑧 𝜎𝑧𝑧 − 𝜎1

𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦 𝜎𝑥𝑧

𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦 𝜎𝑥𝑧 𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1 𝜎𝑦𝑧
b = −det

𝑐2
det 𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1 𝜎𝑦𝑧 𝜎𝑥𝑧 𝜎𝑦𝑧 𝜎𝑧𝑧 − 𝜎1

+
𝜎𝑥𝑧 𝜎𝑦𝑧 𝜎𝑧𝑧 − 𝜎1

+ 𝑏2
1
𝜎𝑥𝑥 − 𝜎1 𝜎𝑥𝑦 𝜎𝑥𝑧

𝑎2
c = det 𝜎𝑥𝑦 𝜎𝑦𝑦 − 𝜎1 𝜎𝑦𝑧
𝜎𝑥𝑧 𝜎𝑦𝑧 𝜎𝑧𝑧 − 𝜎1

𝑟1 =
El mismo proceso puede ser llevado para determinar
la dirección principal del esfuerzo menor
 Deformaciones
Rosetas de Deformación

𝜀𝑥 + 𝜀𝑦 𝜀𝑥 − 𝜀𝑦 𝛾𝑥𝑦
𝜀𝑥′ = + cos 2𝜃 + sin(2𝜃)
2 2 2

𝜀𝑥 + 𝜀𝑦 𝜀𝑥 − 𝜀𝑦 𝛾𝑥𝑦
𝜀𝑦′ = − cos 2𝜃 + sin(2𝜃)
2 2 2

𝜀𝑥 − 𝜀𝑦 𝛾𝑥𝑦
𝛾𝑥′𝑦′ =− sin 2𝜃 + cos(2𝜃)
2 2
 Deformaciones

𝛾𝑥𝑦 𝛾𝑥𝑧
𝜀𝑥𝑥
𝜀𝑥𝑥 𝜀𝑥𝑦 𝜀𝑥𝑧 2 2
𝛾𝑥𝑦 𝛾𝑦𝑧
𝜺 = 𝜀𝑦𝑥 𝜀𝑦𝑦 𝜀𝑦𝑧 𝜺=
2
𝜀𝑦𝑦
2
𝜀𝑧𝑥 𝜀𝑧𝑦 𝜀𝑧𝑧 𝛾𝑥𝑧 𝛾𝑦𝑧
𝜀𝑧𝑧
2 2
 𝜎 𝜎
Constantes Elásticas 𝜎𝑝

∆𝑙

Esfuerzo Axial
El módulo de elasticidad (o modulo de 𝜎
∆𝑡
Young) 𝑬 y la razón de Poisson (𝝂) son las dos 𝜀
constantes elásticas mas usadas para
caracterizar el comportamiento de un sólido.

𝜀
𝜎 𝜀𝑡

Deformación Transversal
𝜀𝑡
𝐸= 𝜈= 𝜀 𝜎
𝜀 𝜀 ∆𝑙 ∆𝑡
𝜀= 𝜀𝑡 =
𝑙 𝑡

El módulo de Young relaciona el esfuerzo aplicado para

producir una deformación longitudinal, por lo tanto se
mide en unidades de esfuerzo [Pa]. La razón de Poisson 1 𝐸 Módulo de Cizalle
relaciona deformación longitudinal con deformación 𝜀𝑡 𝐺=
transversal, por lo tanto es adimensional. 2 (1 + 𝜈)
Constantes Elásticas
 Métodos Indirectos

Constantes Elásticas

Métodos Directos
Constantes Elásticas
Ecuaciones Constitutivas

𝝈𝒙𝒙 𝝈𝟏 𝒖
Ecuaciones Constitutivas
10 5
𝜎=
5 8
10 5
𝜎=
5 8
10 5
𝜎=
5 8
10 5
𝜎=
5 8
Ecuaciones Constitutivas 10 5
𝜎=
5 8
10 5
𝜎=
5 8
10 5
𝜎=
5 8
10 5
𝜎=
5 8
Ecuaciones Constitutivas

𝝈𝒙𝒙 𝝈𝟏 𝒖
Ecuaciones Constitutivas 10 5
𝜎=
5 8
10 5
𝜎=
5 8
10 5
𝜎=
5 8
10 5
𝜎=
5 8
Ecuaciones Constitutivas
𝜎𝑥𝑥 (𝑥, 𝑦) 𝜎𝑥𝑦 (𝑥, 𝑦)
𝜎=
𝜎𝑥𝑦 (𝑥, 𝑦) 𝜎𝑦𝑦 (𝑥, 𝑦)
 Ecuación de Continuidad
𝝏𝝈𝒙𝒚 𝝏𝝈𝒚𝒚
𝝈𝒙𝒚 + ∆𝒚 𝝈𝒚𝒚 + ∆𝒚
𝝏𝒚 𝝏𝒚

𝝏𝝈𝒙𝒚 El sistema está en

𝝈𝒙𝒚 𝝈𝒙𝒚 + ∆𝒙
𝝏𝒙 equilibrio, por lo tanto.

∆𝒚 𝝏𝝈𝒙𝒙
𝝈𝒙𝒙 𝝈𝒙𝒙 + ∆𝒙
𝝏𝒙
∑𝐹𝑥 = 0
𝝈𝒙𝒚 ∑𝐹𝑦 = 0
𝝈𝒚𝒚

∆𝒙
Ecuación de Continuidad

𝝏𝝈𝒙𝒙 𝝏𝝈𝒙𝒚
𝝈𝒙𝒙 ∆𝒚∆𝒛 + 𝝈𝒙𝒚 ∆𝒙∆𝒛 − 𝝈𝒙𝒙 + ∆𝒙 ∆𝒚∆𝒛 − 𝝈𝒙𝒚 + ∆𝒚 ∆𝒙∆𝒛 = 𝟎
𝝏𝒙 𝝏𝒚

𝝏𝝈𝒙𝒙 𝝏𝝈𝒙𝒚
+ =𝟎
𝝏𝒙 𝝏𝒚

𝝏𝝈𝒙𝒚 𝝏𝝈𝒚𝒚
𝝈𝒚𝒚 ∆𝒙∆𝒛 + 𝝈𝒙𝒚 ∆𝒚∆𝒛 − 𝝈𝒙𝒚 + ∆𝒙 ∆𝒚∆𝒛 − 𝝈𝒚𝒚 + ∆𝒚 ∆𝒙∆𝒛 = 𝟎
𝝏𝒙 𝝏𝒚

𝝏𝝈𝒙𝒚 𝝏𝝈𝒚𝒚 Si se considera una fuerza

+ =𝟎 interna, como el peso del
𝝏𝒙 𝝏𝒚
elemento diferencial
Ecuación de Continuidad

𝝏𝝈𝒙𝒙 𝝏𝝈𝒙𝒚
𝝈𝒙𝒙 ∆𝒚∆𝒛 + 𝝈𝒙𝒚 ∆𝒙∆𝒛 − 𝝈𝒙𝒙 + ∆𝒙 ∆𝒚∆𝒛 − 𝝈𝒙𝒚 + ∆𝒚 ∆𝒙∆𝒛 = 𝟎
𝝏𝒙 𝝏𝒚

𝝏𝝈𝒙𝒙 𝝏𝝈𝒙𝒚
+ =𝟎
𝝏𝒙 𝝏𝒚

𝝏𝝈𝒙𝒚 𝝏𝝈𝒚𝒚
𝝈𝒚𝒚 ∆𝒙∆𝒛 + 𝝈𝒙𝒚 ∆𝒚∆𝒛 − 𝝈𝒙𝒚 + ∆𝒙 ∆𝒚∆𝒛 − 𝝈𝒚𝒚 + ∆𝒚 ∆𝒙∆𝒛 + 𝝆𝒈∆𝒙∆𝒚∆𝒛 = 𝟎
𝝏𝒙 𝝏𝒚

𝝏𝝈𝒙𝒚 𝝏𝝈𝒚𝒚
+ + 𝝆𝒈 = 𝟎
𝝏𝒙 𝝏𝒚
Generalizando las Ecuaciones
de Continuidad

𝝏𝝈𝒙𝒙 𝝏𝝈𝒙𝒚 𝝏𝝈𝒙𝒛

+ + + 𝜷𝒙 = 𝟎
𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏𝝈𝒙𝒚 𝝏𝝈𝒚𝒚 𝝏𝝈𝒚𝒛
+ + + 𝜷𝒚 = 𝟎
𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏𝝈𝒙𝒛 𝝏𝝈𝒚𝒛 𝝏𝝈𝒛𝒛
+ + + 𝜷𝒛 = 𝟎
𝝏𝒙 𝝏𝒚 𝝏𝒛
Ecuaciones de Strain

Sabemos que producto de la aplicación de

esfuerzos el sólido se está desplazando y
deformando
Ecuaciones de Strain
𝐴′ 𝐵′ − 𝐴𝐵
𝜀𝑥𝑥 =
𝐴𝐵

Demostración con cálculos

muy avanzados en
progreso…

𝜕𝑢
𝜀𝑥𝑥 =
𝜕𝑥
Ecuaciones de Strain
Versión Alternativa

𝜕𝑢 1 𝜕𝑢 𝜕𝑣 𝜕𝑢 𝜕𝑣
𝜀𝑥𝑥 = 𝜀𝑥𝑦 = + 𝛾𝑥𝑦 = +
𝜕𝑥 2 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑥
𝜕𝑣 1 𝜕𝑣 𝜕𝑤 𝜕𝑣 𝜕𝑤
𝜀𝑦𝑦 = 𝜀𝑦𝑧 = + 𝛾𝑦𝑧 = +
𝜕𝑦 2 𝜕𝑧 𝜕𝑦 𝜕𝑧 𝜕𝑦
𝜕𝑤 1 𝜕𝑢 𝜕𝑤 𝜕𝑢 𝜕𝑤
𝜀𝑧𝑧 = 𝜀𝑥𝑧 = + 𝛾𝑥𝑧 = +
𝜕𝑧 2 𝜕𝑧 𝜕𝑥 𝜕𝑧 𝜕𝑥
 Ecuaciones de
Compatibilidad de Strain

𝝏𝟐 𝜺𝒙𝒙 𝝏𝟐 𝜺𝒚𝒚 𝝏𝟐 𝜸𝒙𝒚 𝝏𝟐 𝜺𝒙𝒙 𝝏 𝝏𝜸𝒚𝒛 𝝏𝜸𝒛𝒙 𝝏𝜸𝒙𝒛

+ = 𝟐 = − + +
𝝏𝒚 𝟐 𝝏𝒙 𝟐 𝝏𝒙𝝏𝒚 𝝏𝒚𝝏𝒛 𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒛

𝝏𝟐 𝜺𝒙𝒙 𝝏𝟐 𝜺𝒛𝒛 𝝏𝟐 𝜸𝒙𝒛 𝝏𝟐 𝜺𝒚𝒚 𝝏 𝝏𝜸𝒚𝒛 𝝏𝜸𝒛𝒙 𝝏𝜸𝒙𝒛

𝟐
+ 𝟐
= 𝟐 = − +
𝝏𝒛 𝝏𝒙 𝝏𝒙𝝏𝒛 𝝏𝒙𝝏𝒛 𝝏𝒚 𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏𝟐 𝜺𝒚𝒚 𝝏𝟐 𝜺𝒛𝒛 𝝏𝟐 𝜸𝒚𝒛 𝝏𝟐 𝜺𝒛𝒛 𝝏 𝝏𝜸𝒚𝒛 𝝏𝜸𝒛𝒙 𝝏𝜸𝒙𝒛
𝟐
+ 𝟐
= 𝟐 = + −
𝝏𝒛 𝝏𝒚 𝝏𝒚𝝏𝒛 𝝏𝒙𝝏𝒚 𝝏𝒛 𝝏𝒙 𝝏𝒚 𝝏𝒛
Resolver Problemas en
Teoría de la Elasticidad
Mejor Opción

1. Resolver el 2. Resolver el 3. Resolver el

Problema mediante Problema mediante Problema mediante
esfuerzos deformaciones desplazamientos
Ecuaciones “Faltantes” De acuerdo al comportamiento de la roca
(elástica, frágil, dúctil) se definirá un conjunto
según el tipo de análisis de ecuaciones para relacionar los esfuerzos con
las deformaciones, según los 3 modelos de
análisis.

𝜎 𝜎 𝜎

𝜀 𝜀

Análisis Elástico Análisis Elastoplástico Análisis Plástico

 Notación Tensorial

La notación de subíndices es útil para resumir y simplicar ecuaciones. Con esta notación se pueden recordar
fácilmente relaciones mas complejas.

Se utiliza la convención de que si un subíndice aparece dos veces en el mismo termino, entonces se suma
 Notación Tensorial

El Delta de Kronecker es:

Ssubíndices separados por una “coma” indica derivada parcial con respecto a…

𝜕𝑎𝑖
𝑎𝑖,𝑗 =
𝜕𝑗
Ecuaciones Constitutivas

Ley de Hooke generalizada

(Para el caso de un análisis elástico)
1 𝜈
𝑒𝑖𝑗 = 𝜎𝑖𝑗 − 𝜎𝑘𝑘 𝛿𝑖𝑗
2𝐺 1+𝑣

𝜈
𝜎𝑖𝑗 = 2𝐺 𝑒𝑖𝑗 + 𝑒𝑘𝑘 𝛿𝑖𝑗
1 − 2𝑣
Modelos Matemáticos de Elasticidad Lineal

Ecuaciones de Equilibrio Ecuaciones Constitutivas Ecuaciones de Strain - Desplazamiento

𝜕𝜎𝑥𝑥 𝜕𝜎𝑥𝑦 𝜕𝜎𝑥𝑧 1 𝜈 𝜕𝑢𝑥

+ + + 𝛽𝑥 = 0 𝑒𝑥𝑥 = 𝜎 − (𝜎 + 𝜎𝑦𝑦 + 𝜎𝑧𝑧 ) 𝑒𝑥𝑥 =
𝜕𝑥 𝜕𝑦 𝜕𝑧 2𝐺 𝑥𝑥 1 + 𝑣 𝑥𝑥 𝜕𝑥
1 𝜈 𝜕𝑢𝑦
𝜕𝜎𝑥𝑦 𝜕𝜎𝑦𝑦 𝜕𝜎𝑦𝑧 𝑒𝑦𝑦 = 𝜎 − (𝜎 + 𝜎𝑦𝑦 + 𝜎𝑧𝑧 ) 𝑒𝑦𝑦 =
+ + + 𝛽𝑦 = 0 2𝐺 𝑦𝑦 1 + 𝜈 𝑥𝑥 𝜕𝑦
𝜕𝑥 𝜕𝑦 𝜕𝑧
1 𝜈 𝜕𝑢𝑧
𝜕𝜎𝑥𝑧 𝜕𝜎𝑦𝑧 𝜕𝜎𝑧𝑧 𝑒𝑧𝑧 = 𝜎𝑧𝑧 − (𝜎 + 𝜎𝑦𝑦 + 𝜎𝑧𝑧 ) 𝑒𝑧𝑧 =
+ + + 𝛽𝑧 = 0 2𝐺 1 + 𝜈 𝑥𝑥 𝜕𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧 1
𝑒𝑥𝑦 = 𝜎 1 𝜕𝑢𝑥 𝜕𝑢𝑦
2𝐺 𝑥𝑦 𝑒𝑥𝑦 = +
2 𝜕𝑦 𝜕𝑥
1
𝑒𝑥𝑧 = 𝜎 1 𝜕𝑢𝑦 𝜕𝑢𝑧
2𝐺 𝑥𝑧 𝑒𝑦𝑧 = +
1 2 𝜕𝑧 𝜕𝑦
𝑒𝑦𝑧 = 𝜎
2𝐺 𝑦𝑧 1 𝜕𝑢𝑥 𝜕𝑢𝑧
𝑒𝑥𝑧 = +
2 𝜕𝑧 𝜕𝑥
Para un problema en 3D tenemos 15 ecuaciones y 15 incógnitas.
Modelos Matemáticos de Elasticidad Lineal

Estado Plano de Esfuerzos Estado Plano de Deformaciones

Los esfuerzos sólo ocurren en el plano xy por lo Las deformaciones sólo ocurren en el plano xy por lo
que 𝜎𝑧 = 0 y 𝜎𝑖𝑗 (𝑥, 𝑦). Para utilizar este estado que 𝑢𝑧 = 0 y 𝒖(𝑥, 𝑦). Para utilizar este estado se
debe garantizar que la geometría no cambia a lo
se debe garantizar que el cuerpo es lo
largo de la dimensión z. Un ejemplo sería un túnel o
suficientemente delgado en la dimensión z. Se
un talud lo suficientemente largo.
ocupa para modelar los esfuerzos en placas
delgadas.
Pero los túneles o los taludes no son infinitos. Se ha
demostrado que la solución matemática es
adecuada cuando la geometría se mantiene
constante durante al menos 5 veces las dimensiones
en el plano xy.
 Problemas Elasticidad 1D

𝑑𝜎𝑥𝑥
+ 𝛽𝑥 = 0
𝑑𝑥
𝑑 𝑑𝑢𝑥
𝜎𝑥𝑥 = 𝐸𝑒𝑥𝑥 𝐸 + 𝛽𝑥 = 0 en Ω
𝑑𝑥 𝑑𝑥
𝑑𝑢𝑥
𝑒𝑥𝑥 =
𝑑𝑥

En cada uno de los contornos

del dominio debo elegir ya
Condición de Borde de Direchlet. 𝒖 = 𝒖𝟎 sea la condición de Newmann
en 𝜕Ω o Direchlet, nunca ambos en
𝜕𝒖
Condición de Borde de Newmann 𝐸 =𝒕 la frontera, tampoco sólo
𝜕𝑥 condiciones de Newmann.
 Funciones de Airy

Cuando las fuerzas de cuerpo son constantes en un sólido con estado de deformación o
esfuerzo plano, el problema elástico se simplifica considerablemente mediante el uso de
una función 𝜙 llamada función de Airy. Así, el problema del cálculo del tensor de esfuerzos
se reduce a la determinación de una función que sea biarmónica.

𝛻4𝜙 = 0
En coordenadas cartesianas
𝜕4𝜙 𝜕4𝜙 𝜕4𝜙
4
+2 2 2+ 4 =0
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
 Funciones de Airy

𝜕2𝜙 𝜕2𝜙 𝜕2𝜙

𝜎𝑥𝑥 = 𝜎𝑦𝑦 = 𝜎𝑥𝑦 =−
𝜕𝑦 2 𝜕𝑥 2 𝜕𝑥𝜕𝑦

1 𝜕𝜙 1 𝜕 2 𝜙 𝜕2𝜙 𝜕 1 𝜕𝜙
𝜎𝑟𝑟 = + 2 2 𝜎𝜃𝜃 = 2 𝜎𝜃𝜃 =−
𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃
 𝜕𝜎𝑟𝑟 1 𝜕𝜎𝑟𝜃 𝜎𝑟𝑟 − 𝜎𝜃𝜃
+ + =0
Excavación Circular 𝜕𝑟 𝑟 𝜕𝜃 𝑟
𝜕𝜎𝑟𝜃 1 𝜕𝜎𝜃𝜃 2𝜎𝑟𝜃
+ + =0
𝜕𝑟 𝑟 𝜕𝜃 𝑟
1 𝜕 2 𝜀𝑟𝑟 𝜕 2 𝜀𝜃𝜃 2 𝜕𝜀𝑟𝜃 1 𝜕𝜀𝑟𝑟 2 𝜕𝜀𝜃𝜃 2 𝜕𝜀𝑟𝜃
𝝈𝟎 𝑟 2 𝜕𝜃 2
+
𝜕𝑟 2
− −
𝑟 𝜕𝑟𝜕𝜃 𝑟 𝜕𝑟
+
𝑟 𝜕𝑟
− 2
𝑟 𝜕𝜃
=0

𝜕𝑢𝑟
𝜀𝑟𝑟 =
𝜕𝑟
1 𝜕𝑢𝜃 𝑢𝑟
𝜀𝜃𝜃 = +
𝑟 𝜕𝜃 𝑟
𝑟
1 1 𝜕𝑢𝑟 𝜕𝑢𝜃 𝑢𝜃
𝜀𝑟𝜃 = + −
2 𝑟 𝜕𝜃 𝜕𝑟 𝑟
𝜃
𝑎 1
𝒌𝝈𝟎 𝜀𝑟𝑟 = 𝜎 − 𝜈𝜎𝜃𝜃
𝐸 𝑟𝑟
1+𝜈
𝜀𝜃𝜃 = −𝜈𝜎𝑟𝑟 + (1 − 𝑣)𝜈𝜎𝜃𝜃
𝐸
1+𝜈
𝜀𝑟𝜃 = 𝜎𝑟𝜃
𝐸
Ecuaciones de Kirsch

1+𝑘 𝑎2 1−𝑘 3𝑎4

𝜎𝜃𝜃 = 𝜎0 1 + 2 + 𝜎0 1 + 4 cos 2𝜃
2 𝑟 2 𝑟

1+𝑘 𝑎2 1−𝑘 4𝑎2 3𝑎4

𝜎𝑟𝑟 = 𝜎0 1 − 2 − 𝜎0 1 − 2 + 4 cos 2𝜃
2 𝑟 2 𝑟 𝑟
1−𝑘 2𝑎 2 3𝑎4
𝜎𝑟𝜃 = 𝜎0 1 + 2 − 4 sin 2𝜃
2 𝑟 𝑟

𝜎0 𝑎2 𝑎2
𝑢𝑟 = − 1 + 𝑘 − 1 − 𝑘 4 1 − 𝜈 − 2 cos 2𝜃
4𝐺𝑟 𝑟
𝜎0 𝑎2 𝑎2
𝑢𝜃 = − 1 − 𝑘 2 1 − 2𝜈 + 2 𝑠𝑖𝑛 2𝜃
4𝐺𝑟 𝑟
 Ecuaciones de Kirsch

𝑘 = 1.5

𝑘 = 0.5
Ecuaciones de Kirsch

Sigma 1 y Sigma 3 pueden ser obtenidos

mediante la búsqueda a los autovalores
del tensor de esfuerzos en coordenadas
polares.

4 16 48 8 48
22.5 1 − 2 − 7.5 1 − 2 + 4 cos 2𝜃 7.5 1 + 2 − 4 sin 2𝜃
𝑟 𝑟 𝑟 𝑟 𝑟
𝑇 𝑟, 𝜃 =
8 48 4 48
7.5 1 + 2 − 4 sin 2𝜃 22.5 1 + 2 + 7.5 1 + 4 cos 2𝜃
𝑟 𝑟 𝑟 𝑟
 Ecuaciones de Kirsch
2 2
8 90 8 48 8 48 2
𝜎1 = 22.5 + 2 cos 2𝜃 + + 7.5 1 − 2 + 4 cos 2𝜃 + 56.25 1 + 2 − 4 sin 2𝜃
𝑟 𝑟2 𝑟 𝑟 𝑟 𝑟

2 2
8 90 8 48 8 48 2
𝜎3 = 22.5 + 2 cos 2𝜃 − + 7.5 1 − 2 + 4 cos 2𝜃 + 56.25 1 + 2 − 4 sin 2𝜃
𝑟 𝑟2 𝑟 𝑟 𝑟 𝑟

Solución Analítica
de Sigma 1 y 3
respectivamente
Ecuaciones de Kirsch

𝑘 = 1.5

𝑘 = 0.5
 Ecuaciones de Kirsch Para las condiciones de esfuerzo dado, las
ecuaciones de Kirsch en el contorno de la
excavación son:

𝜎𝜃𝜃 = 45 + 30 cos 2𝜃
Esfuerzos Máximos y Mínimos
𝜎𝑟𝑟 = 0
𝜎𝑟𝜃 = 0

Esfuerzo σ1
80

70

60

50

MPa
40

30

20

10

0
0 1 2 3 4 5 6
Angulo θ

Solución Analítica Solución Numérica

 Ecuaciones de Kirsch
𝜕𝜎𝜃𝜃 (𝑟, 𝜃)
=0 Sea 𝑝0 = (𝑟0 , 𝜃0 ) una
𝜕𝑟 raíz del sistema de
𝜕𝜎𝜃𝜃 (𝑟, 𝜃)
=0 ecuaciones anterior .
𝜕𝜃
Como vimos anteriormente, los esfuerzos normales y
cortantes son cero en el contorno de la excavación.
Por lo tanto, en esta región Sigma 3 es igual a cero y Matriz Hessiana
Sigma 1 toma el valor de Sigma Theta. ¿Pero son los
esfuerzos máximos en 𝑟 = 𝑎 para todo el dominio?
𝜕 2 𝜎𝜃𝜃 𝜕 2 𝜎𝜃𝜃
ℋ 𝑟, 𝜃 = 𝜕𝑟 2 𝜕𝑟𝜕𝜃
𝜕 2 𝜎𝜃𝜃 𝜕 2 𝜎𝜃𝜃
𝜕𝑟𝜕𝜃 𝜕𝜃 2
Teorema

Si los dos valores propios de ℋ𝑝0 son positivos, 𝑝0 es un mínimo relativo.

Si los dos valores propios de ℋ𝑝0 son negativos, 𝑝0 es un máximo relativo.
Si uno de los valores propios de ℋ𝑝0 es positivo y otro negativo, 𝑝0 es un punto silla.
Planos de Falla

𝐻
 𝑊+𝐻
𝑒0 =
𝑊−𝐻
Excavaciones Elípticas 4(𝑥12 + 𝑧12 )
𝑏=
𝑊 2 − 𝐻2
8(𝑥12 − 𝑧12 )
𝑑= −1
𝑥 𝝈𝟎 𝑊 2 − 𝐻2
𝑒0
𝑢=𝑏+ 𝑏2 − 𝑑
𝑒0
𝐻 𝑒0
𝑧 𝑒=𝑢+ 𝑢2 − 1
𝑒0
𝑊
𝑒 + 1 𝑧1
𝜓 = tan−1
𝑒 − 1 𝑥1
𝑥
𝛽 2
𝑒+1 𝑧1
𝜃= tan−1
𝒌𝝈𝟎 𝑒−1 𝑥1
𝑧1
𝐶 = 1 − 𝑒𝑒0
𝑥1
𝐽 = 1 + 𝑒 2 − 2𝑒 𝑐𝑜𝑠(2𝜓)
 Excavaciones Elípticas

𝑧1 𝜃

𝑚 (𝑥0 , 𝑧0 )
𝑥1

𝑙
𝜎0 𝑒0 − 𝑒 𝐶 𝐽
𝜎𝑙𝑙 = 1+𝑘 𝑒2 −1 + 1−𝑘 𝑒 − 𝑒0 + 𝐶𝑒 cos 2𝜓 + 2𝛽 − 𝐶 cos(2𝛽)
𝐽2 2𝑒0 2

𝜎0
𝜎𝑚𝑚 = 1 + 𝑘 𝑒 2 − 1 + 2 1 − 𝑘 𝑒0 𝑒 cos 2𝜓 + 2𝛽 − cos(2𝛽) − 𝜎𝑙𝑙
𝐽

𝜎0 𝑒0 − 𝑒 𝐶𝑒 𝑗
𝜎𝑙𝑚 = 1+𝑘 sin(2𝜓) + 1 − 𝑘 𝑒 𝑒0 + 𝑒 sin 2𝛽 + 𝑒𝑠𝑖𝑛 2𝜓 − 2𝛽 − 𝑒0 + 𝑒 + 𝑒 2 𝑒0 sin(2𝜓 + 2𝛽)
𝐽2 𝑒0 2
Excavaciones Elípticas

Las ecuaciones anteriores pueden ser reducidas al

sólo considerar los esfuerzos alrededor de la
excavación elíptica, por definición los esfuerzos
normales y cortantes (𝜎𝑙𝑙 y 𝜎𝑙𝑚 ) son cero.

𝜎0
𝜎= 1 + 𝑘 1 + 𝑞 2 + 1 − 𝑞 2 cos 2𝜓 − 2𝛽 − 1−𝑘 1+𝑞 2 cos 2𝜓 + 1 − 𝑞 2 cos 2𝛽
2𝑞

𝑊
𝑞=
𝐻