Universidad Nacional de Trujillo

Facultad de Ingeniería
Departamento de Ingeniería de Minas

Criterios de falla y deformación plana

Resistencia del macizo rocoso

MIH. Eduardo Noriega

emnoriegav@[Link]
15/05/2024 MSc. Eduardo Noriega 1
 Agenda

1. Definiciones
2. Estado de plano de tensiones
3. Elementos finitos triangulares

15/05/2024 MSc. Eduardo Noriega 2

 Estado plano de tensiones
Tensiones y Deformaciones

Tensiones y Deformaciones

Figura 1: Tensor de tensiones.

15/05/2024 MSc. Eduardo Noriega 3

 Estado plano
Notación de Voigt

Notación de Voigt

15/05/2024 MSc. Eduardo Noriega 4

 Estado plano
Estado plano de tensiones

Estado Plano
1 Tensión Plana

2 Deformación Plana

Figura 2: S´olido en estado plano de solicitaciones[2].

15/05/2024 MSc. Eduardo Noriega 5

Estado plano
Tensión y deformación plana

Estructuras sujetas a solicitaciones planas

Figura 3: Ejemplo deformaci´on plana.

Figura 4: Ejemplo de tensi´on plana.
 Estado plano
Estado de Tensión Plana

Hip´otesis Fundamentales
Tensión Plana (Plane Stress)
1 La placa es plana y tiene un plano de simetría (plano medio).
2 Todas las cargas y condiciones de apoyo se ubican en el plano medio.
3 El espesor de la placa es más pequeña en comparación con las dos
dimensiones restantes.
4 Los desplazamientos, las deformaciones y tensiones son uniformes a través del
espesor.
5 Las tensiones transversales (σzz , σxz , σyz ) son despreciables.

15/05/2024 MSc. Eduardo Noriega 7

 Estado plano
Estado de Deformación Plana

Hip´otesis Fundamentales
Deformación Plana
1 Plano de simetría (plano medio).
2 Todas las cargas y condiciones de apoyo se ubican en el plano medio.
3 Una dimensión de la estructura es mayor que las dos dimensiones restantes.
4 Los desplazamientos, las deformaciones y tensiones son uniformes a través del
espesor.
5 Las deformaciones transversales (εzz , εxz , εyz ) son despreciables.

15/05/2024 MSc. Eduardo Noriega 8

 Estado plano
Estado de Deformación Plana

Estados de tensióny deformaciónplana

(1.4)

TensiónPlana Deformación Plana

σzz = σxz = σyz = 0 (1.5) σxz = σyz = 0 (1.7)

εxz = εyz = 0 (1.6) εzz = εxz = εyz = 0 (1.8)
Estado plano
Condiciones de Contorno

Condiciones de frontera en problemas planos

Figura 5: Equilibrio en frontera[2].

 Elasticidad Bidimensional
Campo de desplazamientos

Campo de desplazamientos

(1.9)

donde
u(x, y) = ux −→ desplazamiento en dirección X
v(x, y) = uy −→ desplazamiento en dirección Y

15/05/2024 MSc. Eduardo Noriega 11

Estado plano
Estados de tensióny deformaciónplana

∂ux
εxx =
∂x
∂uy
εyy =
∂y Convenios de signos
γxy = ∂ux + ∂uy εxx = εx = ε11
∂y ∂x εyy = εy = ε22
εzz = εz = ε33
γxy /2 = εxy = εyx = ε12
 Elasticidad Bidimensional
Relación tensión-deformación

Relacióntensión-deformación

σ = Dε
σ = D(ε − ε o ) + σo Condiciones iniciales

(1.10)

Tensiónplana Deformación Plana (1.11)

E E(1 − ν)
d11 = d 22 = (1.12) d11 = d22 = (1.15)
1 − ν2 (1 + ν)(1 − 2ν)
ν
d12 = d21 = νd11 (1.13) d12 = d21 = d11 (1.16)
1− ν
E
d33 = G = (1.14) E
2(1 + ν) d33 = G = (1.17)
2(1 + ν)
15/05/2024 MSc. Eduardo Noriega 13
Elasticidad Bidimensional
Principios de los Trabajos Virtuales (PTV)
Elemento Finito Triangular
Discretización del campo de los desplazamientos

Campo de los desplazamientos

ux = N1u1 + N2u2 + N3u3 (2.1)
x. x x
1 2 3
uy = N1u + N2u + N3u (2.2)
y. y y

Figura 6: Elemento finito trinagular. N # significa

Nodo local 1.
 Elemento Finito Triangular
Funciones de forma

ai + bix + ciy
Ni = ∀i ∈ {1, 2, 3} (2.4)
2Ae
ai = xjyk − xkyj (2.5)
bi = yj − yk (2.6)
ci = xk − yj cíclico en ijk; es decir, 123, 231 y 312 (2.7)

15/05/2024 MSc. Eduardo Noriega 16

 Elemento Finito Triangular
Funciones de Forma

15/05/2024 Figura 7: MSc. Eduardo

Funciones de forma elemento Noriega
trinagular lineal 17
Elemento Finito Triangular
Discretización del campo de los desplazamientos

Campo de los desplazamientos

(2.8)

u = Nae (2.9)
 Elemento Finito Triangular
Discretización del campo de deformaciones

Campo de deformaciones

(2.10)

(2.11)

(2.12)

15/05/2024 MSc. Eduardo Noriega 19

 Elemento Finito Triangular
Discretización del campo de deformaciones

Campo de deformaciones

(2.13)

ai + bix + c iy
Ni = ∀i ∈ {1, 2,3} (2.14)
2Ae
Las deformaciones son constantes dentro del elemento finito

15/05/2024 MSc. Eduardo Noriega 20

 Elemento Finito Triangular
Discretización del campo de deformaciones

Campo de deformaciones

(2.15)

ai = xj yk − xk yj (2.16)
bi = yj − yk (2.17)
ci = xk − xj cíclico en ijk; es decir, 123, 231 y 312 (2.18)

u = Bae (2.19)

15/05/2024 MSc. Eduardo Noriega 21

Elemento Finito Triangular
Ecuación de equilibrio de la discretización

u = Nae (2.20)
PTV aplicado a un elemento finito

(2.21)

(2.22)

(2.23)
(2.24)
(2.25)
PTV aplicado a un elemento finito

(2.26)

(2.27)

(2.28)

(2.29)
 Elemento Finito Triangular
Ecuación de equilibrio de la discretización

(2.30)

(2.31)

Keij − > Matriz de 2× 2

15/05/2024 MSc. Eduardo Noriega 24

 Ensamblaje
Matriz de rigidez global y vector de fuerzas globales

• Sea K y F la matriz de rigidez global y vector global de fuerzas externas

• loop on element in model part

• Ke $F_eˆ{ext} => Compute the local
• stiffness matrix and external load vector
of the element
• Add elemental contribution to K and F (Ensamblaje
de l a matriz de rigidez)
Elemento
Formulación Básica
Finito Rectangular de 4 Nodos

Figura 8: Elemento cuadril´atero de 4 nodos.

Campo de los Desplazamientos

ux = N1u1 + N2u2 + N3u3 + N4u4
x. x x x
uy = 1 2 3
N1u + N2u + N3u + N4u4
y. y y y
Elemento
Formulación básica
Finito Rectangular de 4 Nodos
Campo de los Desplazamientos
Elemento Finito Rectangular de 4 Nodos
Formulación Básica

Funciones de Interpolaci´on
Elemento Finito Rectangular de 4 Nodos
Formulación básica

Campo de las Deformaciones

Elemento Finito Rectangular de 4 Nodos
Formulación básica

Matriz de Rigidez
Elemento Finito Rectangular Bilineal
Funciones de Forma

Figura 9: Funciones de forma elemento 4 nodos

Elemento Finito Rectangular Seréndipido Cuadrático de 8
Nodos
Funciones de Forma

Figura 10: Funciones de forma elemento 8 nodos

Elemento Finito Rectangular Langrangiano Cuadrático
de 9 Nodos
Funciones de Forma

Figura 11: Funciones de forma elemento 9 nodos

Integración Numérica en Dos Dimensiones Sobre Dominios
Rectangulares
Funciones de Forma

Figura 12: Cuadratura de Gauss-Legendre de 1x1, 2x2, 3x3 y 4x4.

Ejemplo. Integrar la siguiente función f (ξ, η) = ξ2η2 para un elemento

rectangular plano de dimensión 2a× 2b.
Elemento Finito Cuadrilátero
Discretización del campo de los desplazamientos

Campo de los desplazamientos.

Modos incompatibles
ux = N1 u1x + N2 u2x + N1 u3x + N4 u4x (3.1)

uy = N1 u1y + N2 u2y + N1 u3y + N4 u4y (3.2)

(3.3)

(3.4)

(3.5)
Figura 13: Elemento finito cuadril´atero. N#
significa Nodo local 1.
(3.6)

Muy bueno en cargas de compresión o tracción pura.

Poca precisi´on en problemas de flexión
Elemento Finito Cuadrilátero
Discretización del campo de los desplazamientos

Adición de modos incompatibles

Campo de los desplazamientos.

u6

u6

Mejoran sustancialmente el comportamiento a flexión Figura 14: Elemento finito cuadril´atero. N#

estos elementos significa Nodo local 1.
 Ejemplo
Tensión plaza

Origen

15/05/2024 MSc. Eduardo Noriega 37

 Ejemplo
Tensión plana

Matriz de rigidez global

15/05/2024 MSc. Eduardo Noriega 38

 REFERENCIAS
1. “Edelbro, C. (2004). Evaluation of rock mass criteria. Licentiate Thesis 72, LULEA
2. Hudson, J.A., and Harrison, J.P. (1997). Engineering Rock Mechanics.
3. Stewart, S.B.V., & Forsyth, W.W. (1995). The Mathews method for open stope design, CIM Bulletin, 88(992), 45- 53.
4. Palmstrom, A. (2005). Measurements of and correlations between block size and rock quality designation (RQD). Tunnels and Underground
Space Technology, 20, 362-377.
5. Palmstrom, A. (1974).Characterization of jointing density and the quality of rock masses (in Norwegian). Internal report, A.B. Berdal, Norway,
26 p.
6. Palmstrom, A. (1982).The volumetric joint count – A useful and simple measure of the degree of rock mass jointing. IAEG Congress, New
Delhi, 1982. pp. V.221-V.228.
7. Palmstrom, A. (1996).The weighted joint density method leads to improved characterization of jointing. Int. Conf. on Recent Advances in
Tunnelling Technology. New Delhi, India, 6 p.
8. Deere, D.U., Hendron, A.J., Patton, F.D., and Cording, E.J. (1967). Design of surface and near surface construction in rock. Proc. 8th U.S.
Symp. Rock Mech., AIME, New York, pp. 237-302.
9. Deere, D.U., and Deere, D.W. (1988). The RQD index in practice. Proc. Symp. Rock Classif. Eng. Purp., ASTM Special Technical Publication 984,
Philadelphia, pp. 91-101.
10. Bieniawski, Z.T. (1973). Engineering classification of jointed rock masses. Trans. S. Afr. Int. Civ. Engrs., 15, 335- 344.
11. Bieniawski, Z.T. (1974). Geomechanics classification of rock masses and its application in tunnelling. In advances in rock mechanics, 28A), 27-
32. Washington, D.C.: Nat. Acad. Sci.
12. Bieniawski, Z.T. (1976). Rock mass classification in rock engineering. In Exploration for rock engineering. Proc. Of the Symp. (ed. Z.T.
Bieniawski), 1, 97-106, Cape Town: Balkema.
13. Bieniawski, Z.T. (1978). Determining rock mass deformability – experiences from case histories. Int. J. Rock Mech. Min. Sci. & Geomech.
Abstr., 15, 237-247.
14. Bieniawski, Z.T. (1979). The geomechanics classification in rock engineering applications. Proc. 4th congr., Int. Soc. Rock Mech., Montreux, 2,
41-48.
15/05/2024 MSc. Eduardo Noriega 39
 Gracias

15/05/2024 MSc. Eduardo Noriega 40