MTODOS NUMRICOS

PROBLEMAS APLICADOS A LA INGENIERIA CIVIL

Semana 05

EJERCICIOS EXPLICATIVOS
Sesin 2

1. Se desea calcular a tirante en un escaln en un

canal por medio de la aplicacin de la ecuacin de 4. Calcular el tirante normal de un canal trapezoidal,
la energa especifica. para un gasto de Q=0.75 m/s, que tiene un ancho
Se tiene un flujo como se indica en la figura, del de plantilla b=1.70m, talud k=2, pendiente de S=
que se conoce el gasto (Q=0.5 m/s) y la tirante de 0.0001 y coeficiente de rugosidad n=0.10. (use
la seccin 1 (y=0.15m), aguas arriba del escaln, mtodo de newton raphson) . Considerar una
en la seccin 2. La altura del escaln es de 0.05m tolerancia de 0.0001
y el ancho del canal de seccin rectangular es de
1.5m. Considerar una tolerancia de 0.0001. (usar
m. de biseccin) 5. Con los datos del ejercicio anterior halle la tirante
critica del canal de seccin trapezoidal (use mtodo
de la secante) . Considerar una tolerancia de 0.0001

2. Use el mtodo de la regla falsa para el ejercicio

anterior.

3. En ingeniera estructural, la frmula de la secante

define la fuerza por unidad de rea P/A, que
ocasiona la tensin mxima en una columna que
tiene una razn de esbeltez L/K dada es:

P m

A ec P L
1 2 sec 0.5
k EA K

Donde:
ec/k= razn de excentricidad, E= mdulo de elasticidad
en una viga de acero, para el ejercicio considerar E=
200000 MPa, ec/k=0.4 y = 2. Calcule P/A para un
L/K= 50. Recuerde que secx= 1/cosx. Considerar un
error de 0.0055.

1 Mtodos Numricos
SOLUCIN La energa especifica de la seccin 1:

1. Usando mtodo de biseccin. - q2 0.3332

Principio de la conservacin de la energa E y
entonces
E 0.15 0.402m
2 gy 2 19.62 0.15 2
entre ambas secciones.

Para hallar el intervalo donde est la

E1 E2 Z .......(1) solucin a la ecuacin planteada se tiene
en cuenta el tirante critico que existe en
Energa especifica se define como: el canal.

v2 q 2 entonces 0.3332
E y .........(2) yc yc
3 3 0.224m
2g g 9.81
Gasto unitario:
Entonces el tirante buscado estar entre
los valores:
Q
q ............(3)
b 0.15 c 0.224
Energa especfica para un canal de
a 0.16
seccin rectangular: b 0.224

q2 Rpta. - Como se considera una tolerancia

E y ........(4) de 0.0001, en la tercera iteracin se obtiene
2 gy 2 la solucin y= 0.184
(ver tabla 1)

Reemplazando la ecuacin (4) en lugar de E2

en la ecuacin (1) 2. Usando mtodo de la regla falsa

q2 Este mtodo de la regla falsa o regula

E1 y 2 z falsi, es similar al de biseccin, aunque
2 gy 2 en lugar de utilizar un promedio
aritmtico de los lmites del intervalo de
Ecuacin a emplear bsqueda para evaluar la siguiente
aproximacin a la raz, se emplea una
relacin de tringulos semejantes.
q2
y (z E1 ) y
3
2 0
2
2
2g y 23 0.352 y 22 0.00566 0

y 23 0.352 y 22 0.00566 0

Hallamos el gasto unitario:

Q
q
b

0.5
q 0.333m 2 / s
1.5

2 Mtodos Numricos
 3. Mtodo de la secante. Ecuacin de continuidad

Q
V ...........(2)
P m A

A ec P L
1 ( 2
) sec 0.5 ( ) El radio hidrulico se define como
k EA K
P
P 250 R .........(3)


A 1 (0.4) sec 0.5 P / 200000 A(50) A
Donde A= rea hidrulica y P el permetro
P 250 mojado.


A 1 (0.4) 1 / cos(0.5 P / 10000000 A Para un canal de seccin trapezoidal se
P 250 tiene que :

A 1 (0.4 / 0.607628314( P / A)
A (b ky) y........(4.1)
P 0.507628314 P / A
P b 2 y 1 k 2 ..(4.2)
A 250(0.607628314 P / A 0.4
P Sustituyendo la ecuacin 3 en 2 y
(151.9070785 P / A 100 0.607628314 P ordenando:
A
P P P
f ( ) 151.9070785( ) 2 100.607628314 Qn A
5
3
A A A
1
2 3
2
s P
Formula de la secante
Es decir, la funcin a resolver es

5
3
A Qn
f ( xi )( xi xi 1 ) donde f ( y) 0.....(5)
xi 1 xi xi 1 1, xi 2
2
3
1 2

f ( xi 1 ) f ( xi )
P S
P P P
f ( ) 151.9070785( ) 2 100.607628314 En la cual A y P son funciones del tirante
A A A y. Para la aplicacin del mtodo de
Newton- Raphson se requiere obtener la
Rpta.- La fuerza por unidad de rea que soporta derivada de la funcin, que en este caso
la columna es de 38.60930523 MPa (ver tabla 3). es:

4. Mtodo de Newton Raphson 2 5 dA A dP

f ' ( y) A 3 .....(6)
De la ecuacin de Manning 3 dy P dy

1 2 3 12 Ordenando:
V R S .......(1)
n 2
Donde R es el radio hidrulico y V la A 3 5 dA A dP
velocidad media del flujo f ' ( y) 2 3 dy P dy ......(7)
P 3

3 Mtodos Numricos
 Por otra parte, las expresiones para el Esta es la ecuacin que debe resolverse para
rea hidrulica y permetro mojado para encontrar el tirante crtico.
una seccin transversal trapecial
resultan:
Rpta .- El tirante crtico y= 0.2445
dP
2 1 k 2 .......(8.1)
dy
dA
b 2ky..........(8.2)
dy

Se usar la ecuacin recursiva para este

mtodo ser:

f ( y)
yi 1 y1
f ' ( y)
Donde f(y) se calcula en la ecuacin 5 y
el f(y) con las ecuaciones 7, 8.1 y 8.2.

Rpta .- La tirante normal

y= 0.9513742

5. Mtodo de secante

El tirante critico ocurre cuando el nmero

de Froude es igual a la unidad, condicin
que puede escribirse como:

Q
1........(1)
A3
g
T
Donde T es el ancho de la superficie
libre, que para un canal de seccin
trapecial est determinada por la
ecuacin:

T b 2ky.........(2)

La ecuacin 1 se puede escribir en la

siguiente forma:

A3 Q 2
0
Tg

Sustituyendo las expresiones para A y

T:

f ( y)
(b ky) y Q 2
3
.......(3)
b 2ky g

4 Mtodos Numricos
 Tabla 1

Iteracion a b c f(a) f(b) f( c )

1 0.16 0.224 0.192 0.0007448 -0.00076253 -0.00023824

2 0.16 0.192 0.176 0.0007448 -0.00023824 0.000208224

3 0.176 0.192 0.184 0.00020822 -0.00023824 -2.780800E-05

Tabla 2

Iteracin a b c f(c )
1 0.16 0.224 0.1946 -0.00029
2 0.16 0.1916 0.1842 -0.000325

Tabla3

5 Mtodos Numricos