UNIVERSIDAD TECNOLÓGICA NACIONAL

Facultad Regional Córdoba

Departamento Ingenierı́a Mecánica 2012 Cálculo Avanzado

Trabajo Práctico #11

Cálculo Numérico de Raı́ces de Ecuaciones
*Tipo de Trabajo: Grupal, con cinco integrantes por cada grupo.

*Objetivo: Cada grupo deberá resolver el problema fı́sico que le corresponda aplicando alguno de los métodos de
resolución visto en clases. Además el método debe ser programado en el software Octave o M atlab.

*A presentar: Planteo del método utilizado para resolver el problema, listado del código, resultados finales (solu-
ción lógica del problema, número de iteraciones (n), error obtenido (Er xn+1 ) y gráficas: n vs xn+1 y n vs Er xn+1 .

*Fecha de entrega: 24 de Octubre de 2012 (curso 3s2).

*Problemas: A continuación presentan los problemas fı́sicos a resolver. Su numeración corresponde al número de
cada grupo.

1. La ecuación de estado de Redlich-Kwong está dada por:

R.T a
p= − √
v − b v(v + b) T

donde R es la constante universal de los gases [0, 518KJ/(kg.K)], T temperatura absoluta [K], p presión absoluta
[kP a] y v el volumen de un kg de gas (m3 /kg). Los parámetros a y b se calculan mediante:

R2 Tc2,5
a = 0, 427
pc
Tc
b = 0, 0866R
pc
donde pc = 4580 kP a y Tc = 191 K. Determine v y la masa de metano que puede contener un depósito
si su volumen es de 3 m3 , temperatura de −50 o C y presión de 65000 kP a. El error admisible propuesto es
Eadm = 0,0001 %.
Método de resolución: Punto fijo o Falsa posición.

2. El volumen V de un lı́quido contenido en un tanque horizontal cilı́ndrico de radio r y longitud L está relacionado
con la profundidad del lı́quido h por:
    
−1 r − h
p
2 2
V = r cos − r − h 2rh − h L
r

Determine h para r = 2 m, L = 5 m y V = 5 m3 . El error admisible propuesto es Eadm = 0,0005 %.

Método de resolución: Bisección o Secante.

3. La ecuación de Ergun, sirve para describir el flujo de un lı́quido a través de un lecho.

∆pρ Dp ε3 (1 − ε)
= 150   + 1, 75
G2o L (1 − ε) Dp Go
µ

∆p es la caı́da de presión, ρ es la densidad del fluido, Go es la velocidad másica (coeficiente de flujo de masa
dividido entre el área de la sección transversal), Dp es el diámetro de las partı́culas dentro del lecho, mu es la

Práctico #11 1 Autor: Adrián D. Boccardo.

 UNIVERSIDAD TECNOLÓGICA NACIONAL
Facultad Regional Córdoba

viscosidad del fluido, L es la longitud del lecho y ε es la fracción vacı́a del lecho.
Dados los siguientes valores para los parámetros encuentre la fracción vacı́a ε del lecho.
Dp Go
= 1000
µ
∆pρDp
= 10
G2o L
El error admisible propuesto es Eadm = 0,0008 %.
Método de resolución: Newton-Raphson o Secante.

4. En una sección de tubo la caida de presión se calcula ası́:

LρV 2
∆p = f
2D
donde ∆p es la caı́da de presión [P a], f el factor de fricción, L longitud del tubo [m], ρ densidad [kg/m3 ], V
velocidad [m/s] y D diámetro del tubo [m]. Para un flujo turbulento la ecuación de Colebrook proporciona un
medio para calcular el factor de fricción.
 
1 ε 2, 51
√ = −2 log + √
f 3, 7D Re f

donde ε es la rugosidad [m] y Re número de Reynolds,

ρV D
Re =
µ

donde µ es la viscosidad dinámica [N.s/m2 ].

Determine ∆p para un tramo horizontal de tubo liso de 0, 2 m de longitud, dadas rho = 1, 23 kg/m3 , µ =
1, 79e−5 N.s/m2 , D = 0, 005 m, V = 40 m/s y ε = 0, 0015 mm. Los tubos lisos tiene un Re < 105 , por lo que un
valor inicial para comenzar a calcular f puede ser obtenido empleando la formula de Blasius: f = 0, 316/Re0,25 .
El error admisible propuesto es Eadm = 0,00001 %.
Método de resolución: Punto fijo o Secante.

5. El desplazamiento de una estructura está definido por la siguiente ecuación para una solución amortiguada:

y = 9e−k.t cos (ω.t)

donde k = 0, 7 y ω = 4. Determine el tiempo que se requiere para que el desplazamiento disminuya a 3, 5. El

error admisible propuesto es Eadm = 0,00008 %
Método de resolución: Punto fijo o Secante.

6. En ingenierı́a estructural, la fórmula de la secante define la fuerza por unidad de área P/A que ocaciona la
tensión máxima σm en una columna que tiene una razón de esbeltez L/k, es:

P σm
= p
A 1 + (ec/k 2 ) sec [0, 5(L/k) P/(E.A)]

donde ec/k 2 es la razón de excentricidad y E es el módulo elásticidad. Si para una viga de acero E = 200000 M P a,
ec/k 2 = 0, 4 y σm = 250 M P a, calcule P/A para L/k = 50. El error admisible propuesto es Eadm = 0,0001 %.
Método de resolución: Falsa posición o Secante.

Práctico #11 2 Autor: Adrián D. Boccardo.

 UNIVERSIDAD TECNOLÓGICA NACIONAL
Facultad Regional Córdoba

7. Una viga uniformemente sujeta a una carga distribuida uniformemente que crece en forma lineal con las vincu-
laciones mostradas en la figura tiene una ecuación de la curva elástica:
wo
y= (−x5 + 2L2 x3 − L4 x)
120E.I.L
Si L = 600 cm, E = 50000 kN/cm2 , I = 30000 cm4 y wo = 2, 5 kN/cm, determinar el punto para el cual
se tiene máxima deflexión (valor de x para el cual dy/dx = 0) y su valor. El error admisible propuesto es
Eadm = 0,0001 %.
Método de resolución: Newton-Raphson o Bisección.

Figura 1:

8. Los sistemas mecánicos reales involucran la deflexı́on de resortes no lineales. Una masa m se libera desde una
distancia h sobre un resorte no lineal. La fuerza de resistencia del resorte F del resorte está dada por la ecuación:

F = −(k1 d + k2 d3/2 )

Es posible usar la conservación de la energı́a para demostrar:

2k2 d5/2 1
0= + k1 d2 − m.g.d + m.g.h
5 2
Resuelva cuál serı́a el valor de d dados los valores siguientes de los parámetros: k1 = 50000 g/s2 , k2 =
40 g/(s2 m( 0, 5)), m = 90 g, g = 9, 81 m/s2 y h = 0, 45 m. El error admisible propuesto es Eadm = 0,0005 %.
Método de resolución: Secante o Bisección.

Figura 2:

Práctico #11 3 Autor: Adrián D. Boccardo.

 UNIVERSIDAD TECNOLÓGICA NACIONAL
Facultad Regional Córdoba

9. En ciertas ocasiones los ingenieros aeroespaciales deben calcular las trayectorias de proyectiles. Un problema
parecido tiene que ver con la trayectoria de una pelota que se lanza. Dicha trayectoria está definida por las
coordenadas (x, y) como se muestra en la imagen. La trayectoria se modela con la ecuación:
g
y = [Link](θo ) − x2 + 1, 8
2vo2 cos2 θo

Calcule el ángulo inicial θo apropiado si la velocidad inicial v0 = 20 m/s y la distancia x al catcher es de 35 m.

Obsérvese que la pelota sale de la mano del lanzador con una elevación yo = 2 m y el catcher la recibe a 1 m.
Exprese el resultado final en grados. Considerar g = 9, 81 m/s2 y un error admisible de Eadm = 0,0001 %.
Método de resolución: Newton-Raphson o Falsa posición.

Figura 3:

Práctico #11 4 Autor: Adrián D. Boccardo.