Métodos Numéricos

UNIVERSIDAD NACIONAL DE HUANCAVELICA
FACULTAD DE CIENCIAS DE LA INGENIERÍA

ESCUELA PROFESIONAL DE INGENIERÍA CIVIL-HUANCAVELICA
"Año de la lucha contra la corrupción y la impunidad"
DIMENSIONAMIENTO DE UNA PRESA DE GRAVEDAD

APLICANDO EL MÉTODO DE NEWTON RAPHSON Y
PUNTO FIJO
CÁTEDRA: MÉTODOS NUMÉRICOS

CATEDRÁTICO: Ing. Lincoln Condori Paytan
ESTUDIANTE: Angel Sullcaray Ichpas
CICLO: V
SECCIÓN: B
HUANCAVELICA 2019
DEDICATORIA
Dedico este trabajo al creador de todas las cosas, el que

me ha dado fortaleza para continuar cuando he estado
a punto de caer; por ello con toda humildad dedico en
primer lugar mi trabajo a Dios.
De igual forma, dedico este trabajo a mi madre

que ha sabido formarme con buenos sentimientos, hábi-
tos y valores, lo cual me ha ayudado a salir adelante en
los momentos más difı́ciles. Al hombre que me dio la vi-
da, el cual a pesar de todo siempre ha sido un excelente
padre para mı́, brindándome sustento en todo momento.
Y a mi familia en general, por brindarme su apo-

yo incondicional en los buenos y malos momentos.
MÉTODOS NUMÉRICOS 2 INGENIERÍA CIVIL

Índice general
0.1 DIMENSIONAMIENTO DE LA PRESA . . . . . . . . 7

0.2 PROCEDIMIENTO DE CÁLCULO . . . . . . . . . . 9

AGRADECIMIENTO
Agradezco al Ing. Anderson Lincol Condori Paitan por las recomendaciones

que me ha brindado para poder realizar este trabajo.
De la misma forma agrezco a los ingenieros Kennedy R. Gomez Tunque

y Carlos Gaspar Paco, quienes también aportaron con sus enseñanzas y re-
comendaciones.

INTRODUCCIÓN
Este trabajo esta orientado al lenguaje de programacion HP PPL, di-

señado para el predimensionamiento de una presa de gravedad aplicando
metodos numericos.
En este trabajo se aplicán los metodos abiertos como Newton-Raphson y

Punto Fijo para poder calcular la raiz de la ecuación no lineal generada por
las fuerzas que ejercen los agentes externos como, el agua(aguas arriba y
aguas abajo),el sismo, el hielo,los sedimentos y el suelo sobre la presa de
gravedad.

MARCO TEÓRICO
Las grandes fuerzas hidrostáticas que actúan sobre las presas tienden a
producir en las mismas:
1) un deslizamiento horizontal a lo largo de su base

2) el vuelco alrededor de la arista de aguas abajo (que se conoce como pie de
la presa) de la base.
Otro factor que puede afectar a la estabilidad de la presa es la fuerza hi-
drostática de levantamiento (o ascensional) que actúa sobre la base. de la
presa, producida por el agua filtrada bajo la misma.
La estabilidad de una presa se comprueba determinando:
•El coeficiente de seguridad contra el deslizamiento.

•El coeficiente de seguridad contra el vuelco.
•la presión sobre la base de la presa.
Figura 1: Represa Tablachaca(Quichuas-Huancavelica)

0.1. DIMENSIONAMIENTO DE LA PRESA

Sobre la presa mostrada actuan inmensas fuerzas ejercidas por el agua(aguas
abajo y aguas arriba), también se observa la fuerza ejercida por el hielo(Fhielo ).
Figura 2:
EL EFECTO DE LAS FUERZAS QUE ACTUAN SOBRE LA PRESA

Estas fuerzas pueden voltear la presa, por ellos debemos construir la mejor
presa para contrarrestar estas fuerzas, esto significa, calcular las dimensiones
apropiadas de la estructura, para esto emplearemos los métodos numéricos
ya que no es posible determinar por cálculos simples el valor de X .
Empleando el metodo de Newton Rapshon y Punto Fijo se estimará el valor

de X para que la presa no inicie su volteo.

la presencia de estas agua arriba y aguas abajo generan cargas en el

suelo(q1 y q2 ), estas se llaman fuerzas de subpresión.
los valores de q1 y q2 vienen dadas por:

q1 = Cγw h1
q2 = Cγw h2

donde:
C= coef. de presión = 0,5
0.2. PROCEDIMIENTO DE CÁLCULO

1) Se secciona el área en 3 secciones para facilitar el cálculo.
2) Se realiza el calculo de las fuerzas verticales y horizontales y su punto de
aplicacion(CG).
Datos:
γconcreto = 23T n/m3
Fhielo = 2T n/m que actua a una distancia de 20m con respecto al punto A.
γagua = 1T n/m3

realizando el cálculo
Figura 3: tabla1
1 8
A1 = (2)(4) = m2
3 3
8
W1 = γconcr. V1 = 23( )(1)
3
W1 = 61,33T n
3
x1 = (2)
4
x1 = 1,5m
A2 = 5(22) = 110m2
W2 = γconcr. V2 = 23(110)(1)
W2 = 2530T n
5
x2 = 2 + ( )
2
x2 = 4,5m
1
A3 = (x − 5)(22)(1) = 11(x − 5)m2
2
W3 = γconcr. V3 = 23(11)(x − 5)(1)
W3 = 253(x − 5)T n
1
x3 = 7 + (x − 5)( )
3

x + 16
x3 = m
3
Hallando la fuerzas verticales y horizontales ejercidas por aguas arriba
hallando la FH1
FH1 = γagua hcg Aproy

FH1 = 1(8)(16)
FH1 = 128T n−f
Hallando la ubicación de FH1

Ixx
Ycp = Ycg +
Ycg Aproy
(1)(16)3 /12
Ycp = 8 +
8(16)
Ycp = 10,67m
entones:
xH1 = 20 − 10,67
xH1 = 9,33m

hallando la FH2 y FV 1
FH2 = γagua hcg Aproy

FH2 = (1)(18)(4)
FH2 = 72T n − f
hallando la ubicación de FH2 respecto al punto A

hallando el area A.
dA = (4 − y)dx
Z Z 2
dA = (4 − x2 )dx
A 0
Z 2
A= (4 − x2 )dx
0
x3 2
A = [4x − ]
3 0
A = 5,33m2
hallando el punto de aplicación de FH2
γagua sen(θ)Ixx
ycp − ycg =
FH ycg Aproy
(1)sen(90◦ )(2(4)3 /12)
ycp − ycg =
72(2)(4)
ycp − ycg = 0,0185
xH2 = 1,98m
cálculo de la fuerza vertical
FV 1 = γagua V
FV 1 = γagua (A + Arect )(longitud)
FV 1 = (1)(5,33 + 16(2))(1)

FV 1 = 37,33T n − f
hallando la ubicación de FV 1 respecto al punto A

R2
0 xdA
xV 1 =
A
R2
0 x(4 − x2 )dx
xV 1 =
A
x4
[2x2 − ]20
xV 1 = 4
5,33
xV 1 = 0,75m
Hallando la fuerzas verticales y horizontales ejercidas por aguas abajo
hallando la FB
8
θ = arctan( )
n

22
θ = arctan( )
x−5
hcg = 4sen2 θ
entonces:
FB = γagua hcg Aproy
FB = (1)(4sen2 θ)(8senθ)
FB = 32sen3 θ
hallando la ubicacion de la fuerza FB conrespectoaO
Ixx
Ycp = Ycg +
Ycg Aproy
(8sen3 θ)/12
Ycp = 4senθ +
4sen2 θ(8senθ)
4
Ycp = 4senθ + senθ
3
entonces:
4
Ycp − Ycg = senθ
3
hallando FH3 y FV 2
FH3 = FB senθ
FH3 = 32sen4 θT n − f
FV 2 = FB cosθ
FV 2 = 32sen3 θcosθT n − f
hallando la ubicación de FH3 respecto a O.
8
xH3 = sen2 θ m
3
hallando la ubicación de FV 2 respecto a 0.
8
xV 2 = senθcosθ m
3

Hallando las fuerzas de subpresión generadas en la base de la presa por efecto

de aguas arriba y aguas abajo
los valores de q1 y q2 estan dadas por:

donde:
C= coef. de presión = 0,5
q1 = Cγw h1 (longitud)
q1 = 1(20)(0,5)(1)
Tn − f
q1 = 10
m
q2 = Cγw h2 (longitud)
q2 = 1(8)(0,5)(1)
Tn − f
q2 = 4
m

hallando la magnitud de las fuerzas F1 y F2

1
F1 = (x + 2)(q1 − q2 )
2
1
F1 = (x + 2)(6)T n − f
2
hallando el punto de aplicación de F1
1
d1 = (x + 2)m
3
F2 = (x + 2)(q2 )
F2 = 4(x + 2)T n − f
Se aplicarán momentos respecto al punto O como dice la teoria, donde el

momento tendrá que ser cero para que la presa no inicie su volteo.
para ello se hallará las nuevas distancias al punto de aplicación de las fuerzas
respecto al punto O, donde las distancias hacia las fuerzas horizontales se
mantienen.

x0V 1 = (x + 2) − xV 1 = (x + 1,25)m
a1 = (x + 2) − x1 = (x + 0,5)m
a2 = (x + 2) − x2 = (x − 2,5)m
2x − 10
a3 = (x + 2) − x3 = m
3
2
a4 = (x + 2) − d1 = (x + 2)m
3
(x + 2)
a5 = m
2

X
MA = 0
FV 1 (x0V 1 ) + W1 (a1) + W2 (a2) + W3 (a3) + FV 2 (xV 2 ) + FH3 (xH3 ) = Fhielo (h1 ) +
FH1 (xH1 ) + FH2 (xH2 ) + F1 (a4) + F2 (a5)
2x − 10
37,33(x + 1,25) + 61,33(x + 0,5) + 2530(x − 2,5) + 253(x − 5)( )+
3
8 8
32sen3 θcosθ( senθcosθ) + 32sen4 θ( sen2 θ) = 2(20) + 128(9,33) + 72(1,98) +
3 3
2 x+2
3(x + 2) (x + 2) + 4(x + 2)( )
3 2
La expresion es reducida a:
494 2 256
x − 925,993x + sen4 (θ) − 3423,806 = 0
3 3
donde
22
θ = (arctan( )
x−5
entonces tenemos:
494 2 256 22
x − 925,993x + sen4 (arctan( ) − 3423,806 = 0
3 3 x−5
para hallar el valor de x aplicaremos el Metodo de Newton Raphson
494 2 256 22
f (x) = x − 925,993x + sen4 (arctan( ) − 3423,806
3 3 x−5
∂f (x)
f 0 (x) =
∂x
aplicando el metodo de newthon Raphson
x1 = x0 − f (x)/f 0 (x)
como
22
θ = (arctan( )
x−5
entonces
xi = 6
Emax = 1e−10

PROCEDIMIENTO
ingreso de datos iniciales
RESULTADOS
Por lo tanto el valor de X es:

x = 8,122m

Ahora aplicaremos el metodo de punto fijo

494 2 256
f (x) = x − 925,993x + sen4 (θ) − 3423,806 + x
3 3
como
22
θ = (arctan( )
x−5
entonces
xi = 6
Emax = 1e−10
PROCEDIMIENTO
ingreso de datos iniciales
RESULTADOS
Por lo tanto el valor de X es:

x = 8,122m

Métodos Numéricos

Cargado por

Copyright:

Formatos disponibles

Métodos Numéricos

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Métodos Numéricos

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE HUANCAVELICA

FACULTAD DE CIENCIAS DE LA INGENIERÍA

"Año de la lucha contra la corrupción y la impunidad"

DIMENSIONAMIENTO DE UNA PRESA DE GRAVEDAD

CÁTEDRA: MÉTODOS NUMÉRICOS

Dedico este trabajo al creador de todas las cosas, el que

De igual forma, dedico este trabajo a mi madre

Y a mi familia en general, por brindarme su apo-

MÉTODOS NUMÉRICOS 2 INGENIERÍA CIVIL

0.1 DIMENSIONAMIENTO DE LA PRESA . . . . . . . . 7

MÉTODOS NUMÉRICOS 3 INGENIERÍA CIVIL

Agradezco al Ing. Anderson Lincol Condori Paitan por las recomendaciones

De la misma forma agrezco a los ingenieros Kennedy R. Gomez Tunque

MÉTODOS NUMÉRICOS 4 INGENIERÍA CIVIL

Este trabajo esta orientado al lenguaje de programacion HP PPL, di-

En este trabajo se aplicán los metodos abiertos como Newton-Raphson y

MÉTODOS NUMÉRICOS 5 INGENIERÍA CIVIL

1) un deslizamiento horizontal a lo largo de su base

•El coeficiente de seguridad contra el deslizamiento.

Figura 1: Represa Tablachaca(Quichuas-Huancavelica)

MÉTODOS NUMÉRICOS 6 INGENIERÍA CIVIL

0.1. DIMENSIONAMIENTO DE LA PRESA

EL EFECTO DE LAS FUERZAS QUE ACTUAN SOBRE LA PRESA

Empleando el metodo de Newton Rapshon y Punto Fijo se estimará el valor

MÉTODOS NUMÉRICOS 7 INGENIERÍA CIVIL

la presencia de estas agua arriba y aguas abajo generan cargas en el

los valores de q1 y q2 vienen dadas por:

MÉTODOS NUMÉRICOS 8 INGENIERÍA CIVIL

0.2. PROCEDIMIENTO DE CÁLCULO

MÉTODOS NUMÉRICOS 9 INGENIERÍA CIVIL

MÉTODOS NUMÉRICOS 10 INGENIERÍA CIVIL

FH1 = γagua hcg Aproy

FH1 = 128T n−f

Hallando la ubicación de FH1

MÉTODOS NUMÉRICOS 11 INGENIERÍA CIVIL

FH2 = γagua hcg Aproy

hallando la ubicación de FH2 respecto al punto A

MÉTODOS NUMÉRICOS 12 INGENIERÍA CIVIL

MÉTODOS NUMÉRICOS 13 INGENIERÍA CIVIL

hallando la ubicación de FV 1 respecto al punto A

Hallando la fuerzas verticales y horizontales ejercidas por aguas abajo

MÉTODOS NUMÉRICOS 14 INGENIERÍA CIVIL

hallando la ubicación de FH3 respecto a O.

MÉTODOS NUMÉRICOS 15 INGENIERÍA CIVIL

Hallando las fuerzas de subpresión generadas en la base de la presa por efecto

los valores de q1 y q2 estan dadas por:

MÉTODOS NUMÉRICOS 16 INGENIERÍA CIVIL

hallando la magnitud de las fuerzas F1 y F2

Se aplicarán momentos respecto al punto O como dice la teoria, donde el

MÉTODOS NUMÉRICOS 17 INGENIERÍA CIVIL

MÉTODOS NUMÉRICOS 18 INGENIERÍA CIVIL

MÉTODOS NUMÉRICOS 19 INGENIERÍA CIVIL

Por lo tanto el valor de X es:

MÉTODOS NUMÉRICOS 20 INGENIERÍA CIVIL

Ahora aplicaremos el metodo de punto fijo

Por lo tanto el valor de X es:

MÉTODOS NUMÉRICOS 21 INGENIERÍA CIVIL

MÉTODOS NUMÉRICOS 22 INGENIERÍA CIVIL