MTODOS DE SOLUCIN DE ECUACIONES NO LINEALES.

ALGORTMOS GENERICOS

JORGE MARIO PEA CONSUEGRA

INGENIERO MECNICO

Pgina 1

JORGE MARIO PEA CONSUEGRA

INGENIERO MECNICO

Pgina 2

JORGE MARIO PEA CONSUEGRA

INGENIERO MECNICO

Pgina 3

JORGE MARIO PEA CONSUEGRA

INGENIERO MECNICO

Pgina 4

JORGE MARIO PEA CONSUEGRA

INGENIERO MECNICO

Pgina 5

BIBLIOGRAFIA
R. Burden, J.D. Faires. Analisis Numrico. 7 Edicin. Thomson Learning. Mexico. 2002.
S. C. Chapra, R. P. Canale. Mtodos Numricos para Ingenieros. 5 Edicin. McGraw-Hill. Mexico.
2006.

JORGE MARIO PEA CONSUEGRA

INGENIERO MECNICO

Pgina 6

PROBLEMAS

217

eficiencia de la conversin algunas veces se mejora recirculando

una porcin de la corriente del producto, de tal forma que regrese a la entrada para un paso adicional a travs del reactor (figura
P8.2). La razn de recirculando se define como
R=

volumen de fluido que regresa a la entrada

volumen que sale del sistema

c = cent(1 e0.04t) + c0e0.04t

Suponga que se est procesando una sustancia qumica A para

generar un producto B. Para el caso en que A forma a B de acuerdo
con una reaccin autocataltica (es decir, en la cual uno de los
productos acta como catalizador o estimulante en la reaccin),
es posible demostrar que una razn ptima de recirculacin
debe satisfacer
ln

1 + R(1 X A )
R(1 X A )

Reactor de flujo
tipo tapn

Producto

O2

Si se asume que sta es la nica reaccin que se lleva a cabo, la

fraccin molar x de H2O que se disocia se representa por
x
1 x

2 pt
2+x

donde la nomenclatura cn representa la concentracin del componente N. Suponga que se define una variable x que representa
el nmero de moles de C producido. La conservacin de la masa
se utiliza para reformular la relacin de equilibrio como
K=

(cc , 0 + x )
(ca, 0 2 x )2 (cb , 0 x )

donde el subndice 0 indica la concentracin inicial de cada componente. Si K = 0.016, ca, 0 = 42, cb, 0 = 28 y cc, 0 = 4, calcule x.
8.6 Las siguientes reacciones qumicas se llevan a cabo en un
sistema cerrado

En equilibrio, stas pueden caracterizarse por

8.3 En un proceso de ingeniera qumica el vapor de agua (H2O)

se calienta a temperaturas lo suficientemente altas para que una
porcin significativa del agua se disocie, o se rompa, para formar
oxgeno (O2) e hidrgeno (H2):

K=

cc
ca2 cb

 C
2A + B

 C
A+D

Figura P8.2
Representacin esquemtica de un reactor de flujo tipo
tapn con recirculacin.

1
2

C
2A + B

K=

Reciclaje

 H2 +
H 2O

Si la concentracin inicial es c0 = 5 y la concentracin de entrada es cent = 12, calcule el tiempo requerido para que c sea el 85%
de cent.
8.5 Una reaccin qumica reversible

se caracteriza por la relacin de equilibrio

R +1
R[1 + R(1 X A )]

donde XA es la fraccin del reactante A que se convierte en el

producto B. La razn ptima de recirculacin corresponde a
un reactor de tamao mnimo necesario para alcanzar el nivel
deseado de conversin. Utilice un mtodo numrico para determinar la razn de recirculacin necesaria, de manera que se
minimice el tamao del reactor para una conversin fraccional
de XA = 0.95.

Alimentacin

donde K = la constante de equilibrio de la reaccin y pt = la

presin total de la mezcla. Si pt = 3.5 atm y k = 0.04, determine
el valor de x que satisfaga la ecuacin (P8.3).
8.4 La siguiente ecuacin permite calcular la concentracin de
un qumico en un reactor donde se tiene una mezcla completa:

(P8.3)

K1 =

cc
ca2 cb

K2 =

cc
ca cd

donde la nomenclatura cn representa la concentracin del componente N. Si x1 y x2 son el nmero de moles de C que se producen debido a la primera y segunda reacciones, respectivamente,
emplee un mtodo similar al del problema 8.5 para reformular las
relaciones de equilibrio en trminos de las concentraciones iniciales de los componentes. Despus, use el mtodo de NewtonRaphson para resolver el par de ecuaciones simultneas no
lineales para x1 y x2 si K1 = 4 104, K2 = 3.7 102, ca,0 = 50,

ESTUDIO DE CASOS: RACES DE ECUACIONES

218

cb,0 = 20, cc,0 = 5 y cd,0 = 10. Utilice un mtodo grfico para proponer los valores iniciales.
8.7 La ecuacin de estado de Redlich-Kwong est dada por
p=

RT
a

v b v(v + b) T

donde R = la constante universal de los gases [= 0.518 kJ/(kg

K)], T = temperatura absoluta (K), p = presin absoluta (kPa) y
v = volumen de un kg de gas (m3/kg). Los parmetros a y b se
calculan mediante
a = 0.427

R2 Tc2.5
pc

b = 0.0866 R

Tc
pc

donde pc = 4 580 kPa y Tc = 191 K. Como ingeniero qumico, se

le pide determinar la cantidad de combustible metano que se
puede almacenar en un tanque de 3 m3 a una temperatura de
50C con una presin de 65 000 kPa. Emplee el mtodo de localizacin de races de su eleccin para calcular v y luego determine la masa de metano contenida en el tanque.
8.8 El volumen V de un lquido contenido en un tanque horizontal cilndrico de radio r y longitud L est relacionado con la
profundidad del lquido h por
r h

V = r 2 cos 1
(r h 2 rh h 2 L
r

Determine h para r = 2 m, L = 5 m y V = 8.5 m3. Observe que si

usted utiliza un lenguaje de programacin o herramienta de
software, el arco coseno se puede calcular como
cos 1 x =

tan 1

2
1 x2

8.9 El volumen V del lquido contenido en un tanque esfrico

de radio r est relacionado con la profundidad h del lquido por

k 2 (3r h)
3
Determine h para r = 1 m y V = 0.75 m3.
8.10 Para el tanque esfrico del problema 8.9, es posible desarrollar las siguientes frmulas para el mtodo de punto fijo:
V=

h=

h 3 + (3V / )
3r

y
V
h = 3 3 rh 2

Si r = 1 m y V = 0.75 m3, determine si cualquiera de las dos alturas es estable, y el rango de valores iniciales para los que s son
estables.
8.11 La ecuacin de Ergun, que se da abajo, sirve para describir
el flujo de un lquido a travs de un lecho empacado. P es la

cada de presin, r es la densidad del fluido, GO es la velocidad

msica (el cociente del flujo de masa dividido entre el rea de la
seccin transversal), Dp es el dimetro de las partculas dentro
del lecho, es la viscocidad del fluido, L es la longitud del lecho
y e es la fraccin vaca del lecho.
(1 )
p Dp 3
= 150
+ 1.75
Go2 L (1 )
Dp Go

Dados los siguientes valores para los parmetros encuentre la
fraccin vaca e del lecho.
Dp Go
= 1 000

PDp
= 10
Go2 L
8.12 En una seccin de tubo, la cada de presin se calcula as:
L V 2
2D
donde p = cada de presin (Pa), f = factor de friccin, L =
longitud del tubo [m], r = densidad (kg/m3), V = velocidad (m/s),
y D = dimetro (m). Para el flujo turbulento, la ecuacin de
Colebrook proporciona un medio para calcular el factor de friccin,

1
2.51
= 2.0 log
+
f
3
.
7
D
Re
f

p = f

donde e = rugosidad (m), y Re = nmero de Reynolds,

VD

donde m = viscosidad dinmica (N s/m2).

Re =

a) Determine p para un tramo horizontal de tubo liso de

0.2 m de longitud, dadas r = 1.23 kg/m3, m = 1.79 105
N s/m2, D = 0.005 m, V = 40 m/s, y e = 0.0015 mm. Utilice
un mtodo numrico para determinar el factor de friccin.
Obsrvese que los tubos lisos tienen Re < 105, un valor
inicial apropiado se obtiene con el uso de la frmula de
Blasius, f = 0.316/Re0.25.
b) Repita el clculo pero para un tubo de acero comercial ms
rugoso (e = 0.045 mm).
8.13 El pH del agua tiene gran importancia para los ingenieros
ambientales y qumicos. Se relaciona con procesos que van de
la corrosin de tubos de lluvia cida. El pH se relaciona con la
concentracin del ion de hidrgeno por medio de la ecuacin
siguiente:
pH = log10 [H+]

PROBLEMAS

219

Ingeniera civil y ambiental

8.15 El desplazamiento de una estructura est definido por la
ecuacin siguiente para una oscilacin amortiguada:

Las cinco ecuaciones que siguen gobiernan las concentraciones

de una mezcla de dixido de carbono y agua para un sistema
cerrado.
K1 =

[H + ][HCO3 ]
[CO 2 ]

K2 =

[H + ][CO32 ]
[HCO3 ]

y = 9ekt cos wt
donde k = 0.7 y w = 4.
a) Utilice el mtodo grco para realizar una estimacin inicial del tiempo que se requiere para que el desplazamiento
disminuya a 3.5.
b) Emplee el mtodo de Newton-Raphson para determinar la
raz con es = 0.01%.
c) Use el mtodo de la secante para determinar la raz con es
= 0.01%.

Kw = [H + ][OH ]
cT = [CO 2 ] + [HCO3 ] + [CO32 ]
Alk = [HCO3 ] + 2[CO32 ] + [OH ] [H + ]
donde Alk = alcalinidad, cT = total de carbn inorgnico, y las K
son coeficientes de equilibrio. Las cinco incgnitas son [CO2] =
dixido de carbono, [HCO3 ] = bicarbonato, [CO2
3 ] = carbonato,
[H+] = ion hidrgeno, y [OH] = ion hidroxilo. Resuelva para las
cinco incgnitas dado que Alk = 2 103, cT = 3 103, K1 =
106.3, y K2 = 1010.3, y Kw = 1014. Asimismo, calcule el pH de
las soluciones.
8.14 La ecuacin que se presenta a continuacin, describe la
operacin de un reactor de flujo por inyeccin de densidad constante para la produccin de una sustancia por medio de una reaccin enzimtica, donde V es el volumen del reactor, F es la tasa
de flujo del reactivo C, Cent y Csal son las concentraciones del
reactivo que entra y sale del reactor, respectivamente, y K y kmx
son constantes. Para un reactor de 500 L, con una concentracin
en la toma de Cent = 0.5 M, tasa de entrada de flujo de 40 L/s,
kmx = 5 103s1, y K = 0.1 M, encuentre la concentracin de C
a la salida del reactor.
V
=
F

Csal

Cent

kmx C

1
kmx

8.16 En ingeniera estructural, la frmula de la secante define la

fuerza por unidad de rea, P/A, que ocasiona la tensin mxima
sm en una columna que tiene una razn de esbeltez L/k dada
es:
P
m
=
A 1 + (ec/k 2 ) sec [0.5 P/( EA) (L/k )]
donde ec/k2 = razn de excentricidad, y E = mdulo de elasticidad. Si para una viga de acero, E = 200 000 MPa, ec/k2 = 0.4 y
sm = 250 MPa, calcule P/A para L/k = 50. Recuerde que sec x =
1/cos x.
8.17 Un cable en forma catenaria es aquel que cuelga entre dos
puntos que no se encuentran sobre la misma lnea vertical. Como
se ilustra en la figura P8.17a, no est sujeta a ms carga que su
propio peso. As, su peso (N/m) acta como una carga uniforme
por unidad de longitud a lo largo del cable. En la figura P8.17b,
se ilustra un diagrama de cuerpo libre de una seccin AB, donde

dC

Figura P8.17
a) Fuerzas que actan sobre
una seccin AB de un cable
flexible que cuelga. La
carga es uniforme a lo largo
del cable (pero no uniforme
por la distancia horizontal
x). b) Diagrama de cuerpo
libre de la seccin AB.

TB

A
w

W = ws
TA

y0
x

a)

b)

ESTUDIO DE CASOS: RACES DE ECUACIONES

220

TA y TB son las fuerzas de tensin en el extremo. Con base en los

balances de fuerzas horizontal y vertical, se obtiene para el cable
el siguiente modelo de ecuacin diferencial:
d2y w
dy
=
1+
dx
dx 2 TA

Puede emplearse el clculo para resolver esta ecuacin para la

altura y del cable como funcin de la distancia x.
y=

w
TA
cosh
w
TA

T
x + y0 A
w

c = 10 20(e0.15x e0.5x)

donde el coseno hiperblico se calcula por medio de la ecuacin:

1
cosh x = (e x + e x )
2
Utilice un mtodo para calcular un valor para el parmetro TA
dados los valores de los parmetros w = 12 y y0 = 6, de modo
que el cable tenga una altura de y = 15 en x = 50.
8.18 En la figura P8.18a se muestra una viga uniforme sujeta a
una carga distribuida uniformemente que crece en forma lineal.
La ecuacin para la curva elstica resultante es la siguiente
(vase la figura P8.18b)
y=

Utilice el mtodo de la biseccin para determinar el punto de

mxima deflexin (es decir, el valor de x donde dy/dx = 0). Despus, sustituya este valor en la ecuacin (P8.18) a fin de determinar el valor de la deflexin mxima. En sus clculos, utilice
los valores siguientes para los parmetros: L = 600 cm, E = 50 000
kN/cm2, I = 30 000 cm4 y w0 = 2.5 kN/cm.
8.19 En la ingeniera ambiental (una especialidad de la ingeniera civil), la ecuacin siguiente se emplea para calcular el nivel
de oxgeno c (mg/L) en un ro aguas abajo de la descarga de un
drenaje:

w0
( x 5 + 2 L2 x 3 L4 x )
120 EIL

(P8.18)

donde x es la distancia aguas abajo en kilmetros.

a) Determine la distancia aguas abajo de la corriente, a la cual
el nivel de oxgeno cae hasta una lectura de 5 mg/L. (Recomendacin: est dentro de 2 km de la descarga.) Encuentre
la respuesta con un error de 1%. Obsrvese que los niveles
de oxgeno por debajo de 5 mg/L por lo general son dainos
para ciertas especies de pesca deportiva, como la trucha y
el salmn.
b) Calcule la distancia aguas abajo a la cual el oxgeno se
encuentra al mnimo. Cul es la concentracin en dicha
ubicacin?
8.20 La concentracin de bacterias contaminantes c en un lago
disminuye de acuerdo con la ecuacin
c = 75e1.5t + 20e0.075t
Determine el tiempo que se requiere para que la concentracin
de bacterias se reduzca a 15 con el uso de a) el mtodo grfico, y
b) el mtodo de Newton-Raphson, con un valor inicial de t = 6
y criterio de detencin de 0.5%. Compruebe los resultados que
obtenga.
8.21 En ingeniera oceanogrfica, la ecuacin de una ola estacionaria reflejada en un puerto est dada por l = 16, t = 12,
v = 48:

Figura P8.18

w0

2x
2tv x

h = h0 sen
cos
+e

a)

(x = L, y = 0)
(x = 0, y = 0)
x

b)

Resuelva para el valor positivo ms bajo de x, si h = 0.5 h0.

8.22 Suponga el lector que compra una pieza de equipo en
$25 000 como pago inicial y $5 500 por ao durante 6 aos. Qu
tasa de inters estara pagando? La frmula que relaciona el
valor presente P, los pagos anuales A, el nmero de aos n y la
tasa de inters i, es la que sigue:
A=P

i(1 + i )n
(1 + i )n 1

PROBLEMAS

221

20 kips/ft

150 kips-ft

15 kips

Figura P8.24

8.23 Muchos campos de la ingeniera requieren estimaciones

exactas de la poblacin. Por ejemplo, los ingenieros de transporte quizs encuentren necesario determinar por separado la tendencia del crecimiento de una ciudad y la de los suburbios. La
poblacin del rea urbana declina con el tiempo de acuerdo con
la ecuacin:
Pu(t) = Pu,mxekut + Pu,mn
en tanto que la poblacin suburbana crece segn:
ps (t ) =

M(x) = 10[x 02 x 52] + 15x 81

+ 150x 70 + 57x
Emplee un mtodo numrico para encontrar el (los) punto(s) en
los que el momento es igual a cero.
8.26 Con el uso de la viga con apoyo simple del problema 8.24,
la pendiente a lo largo de ella est dada por:
duy
dx

( x) =

Ps , mx
1 + [ Ps , mx / P0 1]e kst

donde Pu,mx, ku, Ps,mx, P0 y ks son parmetros que se obtienen en

forma emprica. Determine el tiempo y los valores correspondientes de Pu(t) y Ps(t) cuando los suburbios son 20% ms grandes
que la ciudad. Los valores de los parmetros son: Pu,mx = 75 000,
Ku = 0.045/ao, Pu,mn = 100 000 personas, Ps,mx = 300 000 personas, P0 = 10 000 personas, ks = 0.08/ao. Para obtener las soluciones utilice los mtodos a) grfico, b) de la falsa posicin, y
c) de la secante modificada.
8.24 En la figura P8.24 se muestra una viga apoyada en forma
sencilla que est cargada como se ilustra. Con el empleo de
funciones de singularidad, el esfuerzo cortante a lo largo de la
viga se expresa con la ecuacin:
V(x) = 20[x 01 x 51] 15x 80 57
Por definicin, la funcin de singularidad se expresa del modo
que sigue:
( x a) n
x a n =
0

cuando x > a

cuando x a

Utilice un mtodo numrico para encontrar el(los) punto(s) en

los que el esfuerzo cortante sea igual a cero.
8.25 Con el uso de la viga apoyada en forma simple del problema 8.24, el momento a lo largo de ella, M (x) est dada por:

10
15
[ x 0 3 x 5 3 ] + x 8 2
3
2
57
x 2 238.25
+ 150 x 71 +
2

Utilice un mtodo numrico para encontrar el(los) punto(s)

donde la pendiente es igual a cero.
8.27 Para la viga con apoyo simple del problema 8.24, el desplazamiento a lo largo de ella est dado por la ecuacin:
uy ( x ) =

5
15
[ x 0 4 x 5 4 ] + x 8 3
6
6
57 3
2
+ 75 x 7 + x 238.25 x
6

a) Calcule el (los) punto(s) donde el desplazamiento es igual

a cero.
b) Cmo se usara una tcnica de localizacin de races para
determinar la ubicacin del desplazamiento mnimo?
Ingeniera elctrica
8.28 Ejecute el mismo clculo que en la seccin 8.3, pero determine el valor de C que se requiere para que el circuito disipe 1%
de su valor original en t = 0.05 s, dado R = 280 , y L = 7.5 H.
Emplee a) un enfoque grfico, b) la biseccin, y c) software para
encontrar races, tales como Solver de Excel o la funcin fzero
de MATLAB.
8.29 La ecuacin i = 9et cos (2pt), describe una corriente oscilatoria en un circuito elctrico, donde t se expresa en segundos.
Determine todos los valores de t de modo que i = 3.

ESTUDIO DE CASOS: RACES DE ECUACIONES

222

8.32 En la figura P8.32 se muestra un circuito con una resistencia, un inductor y un capacitor en paralelo. Para expresar la
impedancia del sistema se emplean las leyes de Kirchhoff, as:

a
x
q

1
=
Z

Figura P8.31

8.30 La resistividad r de un lubricante de slice se basa en la

carga q en un electrn, la densidad del electrn n, y la movilidad
del electrn m. La densidad del electrn est dada en trminos de
la densidad del lubricante N, y la densidad intrnseca de acarreo
ni. La movilidad del electrn est descrita por la temperatura T,
la temperatura de referencia T0, y la movilidad de referencia 0.
Las ecuaciones que se requieren para calcular la resistividad son
las siguientes:

1
qn

donde
n=

1
N + N 2 + 4 ni2
2

T
= 0
T0

2.42

Determine N, dado que T0 = 300 K, T = 1 000 K, 0 = 1 350 cm2

(V s)1, q = 1.7 1019 C, ni = 6.21 109 cm3, y un valor deseable de r = 6.5 106 V s cm/C. Use los mtodos a) biseccin, y
b) la secante modificada.
8.31 Una carga total Q se encuentra distribuida en forma uniforme alrededor de un conductor en forma de anillo con radio a.
Una carga q se localiza a una distancia x del centro del anillo
(vase la figura P8.31). La fuerza que el anillo ejerce sobre la
carga est dada por la ecuacin
F=

1
qQx
4e0 ( x 2 + a 2 )3/ 2

donde e0 = 8.85 1012 C2/(N m2). Encuentre la distancia

x donde la fuerza es de 1.25 N, si q y Q son 2 105 C para un
anillo con un radio de 0.9 m.

1
1
+ C
R2
L

donde Z = impedancia () y w = frecuencia angular. Encuentre

la w que da como resultado una impedancia de 75 , con el uso
tanto del mtodo de la biseccin como el de la falsa posicin,
con valores iniciales de 1 y 1000 y los parmetros siguientes: R
= 225 , C = 0.6 106 F, y L = 0.5 H. Determine cuntas iteraciones son necesarias con cada tcnica a fin de encontrar la
respuesta con es = 0.1%. Utilice el enfoque grfico para explicar
cualesquiera dificultades que surjan.
Ingeniera mecnica y aeroespacial
8.33 Para la circulacin de fluidos en tubos, se describe a la
friccin por medio de un nmero adimensional, que es el factor
de friccin de Fanning f. El factor de friccin de Fanning depende de cierto nmero de parmetros relacionados con el tamao
del tubo y el fluido, que pueden representarse con otra cantidad
adimensional, el nmero de Reynolds Re. Una frmula que pronostica el valor de f dado Re es la ecuacin de von Karman.

1
= 4 log10 Re 0.4

Valores comunes del nmero de Reynolds para flujo turbulento

son 10 000 a 500 000, y del factor de friccin de Fanning son
0.001 a 0.01. Desarrolle una funcin que utilice el mtodo de
biseccin con objeto de resolver cul sera el factor de friccin
de Fanning f, dado un valor de Re proporcionado por el usuario
que est entre 2 500 y 1 000 000. Disee la funcin de modo que
se garantice que el error absoluto en el resultado sea de Ea,d <
0.000005.
8.34 Los sistemas mecnicos reales involucran la deflexin de
resortes no lineales. En la figura P8.34 se ilustra una masa m que
se libera por una distancia h sobre un resorte no lineal. La fuerza
de resistencia F del resorte est dada por la ecuacin

Figura P8.34

Figura P8.32
h
h+d
d

a)

b)

PROBLEMAS

223

F = (k1d + k2d3/2)
Es posible usar la conservacin de la energa para demostrar
que
0=

2 k2 d 5/ 2 1
+ k1d 2 mgd mgh
5
2

v = u ln

Resuelva cul sera el valor de d, dados los valores siguientes de

los parmetros: k1 = 50 000 g/s2, k2 = 40 g/(s2 m0.5), m = 90 g,
g = 9.81 m/s2, y h = 0.45 m.
8.35 Los ingenieros mecnicos, as como los de otras especialidades, utilizan mucho la termodinmica para realizar su trabajo.
El siguiente polinomio se emplea para relacionar el calor especfico a presin cero del aire seco, cp kJ/(kg K), a temperatura
(K):
cp = 0.99403 + 1.671 104 T + 9.7215 108T 2
9.5838 1011T 3 + 1.9520 1014T4
Determine la temperatura que corresponda a un calor especfico
de 1.1 kJ/(kg K).
8.36 En ciertas ocasiones, los ingenieros aerospaciales deben
calcular las trayectorias de proyectiles, como cohetes. 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 ilustra en la figura P8.36. La trayectoria se modela
con la ecuacin
y = ( tan 0 ) x

Para g, utilice un valor de 9.81 m/s2, y emplee el mtodo grfico

para elegir valores iniciales.
8.37 La velocidad vertical de un cohete se calcula con la frmula que sigue:

g
x 2 + 1.8
2v cos2 0
2
0

Calcule el ngulo inicial q0, apropiado si la velocidad inicial

v0 = 20 m/s y la distancia x al catcher es de 35 m. Obsrvese que
la pelota sale de la mano del lanzador con una elevacin y0 = 2 m,
y el catcher la recibe a 1 m. Exprese el resultado final en grados.

m0
gt
m0 qt

donde v = velocidad vertical, u = velocidad con la que se expele

el combustible, en relacin con el cohete, m0 = masa inicial del
cohete en el momento t = 0, q = tasa de consumo de combustible,
y g = aceleracin de la gravedad hacia abajo (se supone constante e igual a 9.81 m/s2). Si u = 2000 m/s, m0 = 150 000 kg, y q =
2 700 kg/s, calcule el momento en que v = a 750 m/s. (Sugerencia: El valor de t se encuentra entre 10 y 50 s.) Calcule el resultado de modo que est dentro de 1% del valor verdadero.
Compruebe su respuesta.
8.38 En la seccin 8.4, el ngulo de fase f entre la vibracin
forzada que ocasiona el camino rugoso y el movimiento del carro,
est dada por la ecuacin:
tan =

2(c/cc )( / p)
1 ( / p)2

Como ingeniero mecnico, le gustara saber si existen casos en

que f = w/3 1. Utilice los otros parmetros de la seccin con
objeto de plantear la ecuacin como un problema de clculo de
races, y resulvala para w.
8.39 Se mezclan dos fluidos con temperatura diferente de modo
que alcanzan la misma temperatura. La capacidad calorfica del
fluido A est dada por:
cp = 3.381 + 1.804 102T 4.300 106 T 2
y la capacidad calorfica del fluido B se obtiene con:
cp = 8.592 + 1.290 101T 4.078 105 T 2
donde cp se expresa en unidades de cal/mol K, y T est en unidades de K. Obsrvese que

Figura P8.36

H =

T2

T1

v0
0
x

c p dT

El fluido A entra al mezclador a 400C, y el B a 700C. Al entrar

al mezclador hay lo doble de fluido A que B. A qu temperatura salen los dos fluidos del mezclador?
8.40 Un compresor opera a una razn de compresin Rc de 3.0
(esto significa que la presin del gas en la salida es tres veces
mayor que en la entrada). Los requerimientos de energa del
compresor Hp se determinan por medio de la ecuacin que se da
a continuacin. Suponga que los requerimientos de energa del
compresor son exactamente iguales a zRT1/MW, y encuentre la
eficiencia politrpica n del compresor. El parmetro z es la compresibilidad del gas en las condiciones de operacin del compre-

ESTUDIO DE CASOS: RACES DE ECUACIONES

224

xx

xy
xz

xy xz
yy yz

yz zz

en la que los trminos en la diagonal principal representan esfuerzos a la tensin o a la compresin, y los trminos fuera de la
diagonal representan los esfuerzos cortantes. Un campo tensorial
(en MPa) est dado por la matriz que sigue:

T3

T0
T1

T2

Figura P8.41

10 14 25
14 7 15

25 15 16
Para resolver cules son los esfuerzos principales, es necesario
construir la matriz siguiente (de nuevo en MPa):
10
14

25

14
25
7
15
15 16

s1, s2 y s3 se obtienen con la ecuacin

sor, R es la constante de los gases, T1 es la temperatura del gas
en la entrada del compresor, y MW es el peso molecular del
gas.
zRT1 n
HP =
( Rc(n1)/ n 1)
MW n 1
8.41 En los envases trmicos que se ilustran en la figura P8.41,
el compartimiento interior est separado del medio por medio de
vaco. Hay una cubierta exterior alrededor de los envases. Esta
cubierta est separada de la capa media por una capa delgada de
aire. La superficie de afuera de la cubierta exterior est en contacto con el aire del ambiente. La transferencia de calor del
compartimiento interior a la capa siguiente q1 slo ocurre por
radiacin (ya que el espacio se encuentra vaco). La transferencia
de calor entre la capa media y la cubierta exterior q2 es por conveccin en un espacio pequeo. La transferencia de calor de la
cubierta exterior hacia el aire q3 sucede por conveccin natural.
El flujo de calor desde cada regin de los envases debe ser
igual, es decir, q1 = q2 = q3. Encuentre las temperaturas T1 y T2
en estado estable. T0 es de 450C y T3 = 25C.
q1 = 10 9 [(T0 + 273)4 (T1 + 273)4 ]
q2 = 4(T1 T2 )

3 I 2 + II III = 0
donde
I = xx + yy + zz
II = xx yy + xx zz + yy zz xy2 xz2 yz2
III = xx yy zz xx yz2 yy xz2 zz xy2 + 2 xy xz yz
I, II y III se conocen como las invariantes de esfuerzos. Encuentre
s1, s2 y s3 por medio de una tcnica de localizacin de races.
8.43 La figura P8.43 ilustra tres almacenamientos conectados
por medio de tubos circulares. Los tubos estn hechos de hierro

Figura P8.43
h1
A

h2
1
Q1

h3

Q2

q3 = 1.3(T2 T3 )4 / 3

3
Q3

8.42 La forma general para un campo tensorial de tres dimensiones es la siguiente:

PROBLEMAS

225

fundido recubierto con asfalto (e = 0.0012 m), y tienen las caractersticas siguientes:
Tubo
Longitud, m
Dimetro, m
Flujo, m3/s

1
1800
0.4
?

2
500
0.25
0.1

3
1400
0.2
?

Q1 = 1 m3/s y r = 1.23 kg/m3. Todos los tubos tienen D = 500

mm y f = 0.005. Las longitudes de los tubos son: L3 = L5 = L8 =
L9 = 2 m; L2 = L4 = L6 = 4 m; y L7 = 8 m.
8.45 Repita el problema 8.44, pero incorpore el hecho de que el
factor de friccin se calcula con la ecuacin de von Karman, que
es:
1
= 4 log10 (Re f ) 0.4
f

Si las elevaciones de la superficie del agua en los almacenamientos A y C son de 200 m y 172.5 m, respectivamente, determine
la elevacin que alcanza en el almacenamiento B y los flujos en
los tubos 1 y 3. Obsrvese que la viscosidad cinemtica del agua
es de 1 106 m2/s, y utilice la ecuacin de Colebrook para obtener el factor de friccin (consulte el problema 8.12).
8.44 Un fluido se bombea en la red de tubos que se muestra en
la figura P8.44. En estado estacionario, se cumplen los balances
de flujo siguientes:
Q1 = Q2 + Q3
Q3 = Q4 + Q5
Q5 = Q6 + Q7
donde Qi = flujo en el tubo i [m3/s]. Adems, la cada de presin
alrededor de los tres lazos en los que el flujo es hacia la derecha
debe ser igual a cero. La cada de presin en cada tramo de tubo
circular se calcula por medio de la ecuacin:
16 fL 2
Q
2 2 D5
donde P= cada de presin [Pa], f = factor de friccin [adimensional], L = longitud del tubo [m], r = densidad del fluido [kg/m3],
y D = dimetro del tubo [m]. Escriba un programa (o desarrolle
un algoritmo en algn paquete de software de matemticas) que
permita calcular el flujo en cada tramo de tubo, dado que
P =

Figura P8.44
Q1

Q3

Q2

Q10

Q5

Q4

Q9

Q6

Q8

Q7

donde Re = nmero de Reynolds

Re =

VD

donde V = velocidad del fluido en el tubo [m/s], y = viscosidad

dinmica (N s/m2). Obsrvese que para un tubo circular, V = 4Q/
pD2. Asimismo, suponga que el fluido tiene una viscosidad de
1.79 105 N s/m2.
8.46 Sobre el trasbordador espacial, al despegar de la plataforma,
actan cuatro fuerzas, las que se muestran en el diagrama de
cuerpo libre (vase la figura P8.46). El peso combinado de los
dos cohetes de combustible slido y del tanque exterior de este,
es de WB = 1.663 106 lb. El peso del orbitador con carga completa es de WS = 0.23 106 lb. El empuje combinado de los dos
cohetes de combustible slido es TB = 5.30 106 lb. El empuje
combinado de los tres motores de combustible lquido del orbitador es de TS = 1.125 106 lb.
Al despegar, el empuje del motor del orbitador se dirige con
un ngulo q para hacer que el momento resultante que acta sobre
el conjunto de la nave (tanque exterior, cohetes de combustible
slido y orbitador) sea igual a cero. Con el momento resultante
igual a cero, la nave no girara sobre su centro de gravedad G al
despegar. Con estas fuerzas, la nave experimentar una fuerza
resultante con componentes en direccin vertical y horizontal.
La componente vertical de la fuerza resultante, es la que permite
que la nave despegue de la plataforma y vuele verticalmente. La
componente horizontal de la fuerza resultante hace que la nave
vuele en forma horizontal. El momento resultante que acta sobre
la nave ser igual a cero cuando q se ajusta al valor apropiado.
Si este ngulo no se ajusta en forma adecuada y hubiera algn
momento que actuara sobre la nave, sta tendera a girar alrededor
de su centro de gravedad.
a) Resuelva el empuje del orbitador TS en las componentes
horizontal y vertical, y despus sume los momentos respecto
del punto G, centro de gravedad de la nave. Iguale a cero
la ecuacin del momento resultante. Ahora, sta puede
resolverse para el valor de q que se requiere durante el
despegue.
b) Obtenga una ecuacin para el momento resultante que
acta sobre la nave en trminos del ngulo q. Graque el

ESTUDIO DE CASOS: RACES DE ECUACIONES

226

28
4

Tanque externo

Cohete de
combustible
slido

Orbitador

38
WB

WS

TS
TB

Figura P8.46

momento resultante como funcin del ngulo q en el rango

de 5 radianes a +5 radianes.
c) Escriba un programa de computadora para resolver para el
ngulo q por medio del mtodo de Newton para encontrar la
raz de la ecuacin del momento resultante. Con el empleo
de la grca, elija un valor inicial para la raz de inters.
Interrumpa las iteraciones cuando el valor de q ya no mejore
con cinco cifras signicativas.
d) Repita el programa para el peso de la carga mnima del
orbitador, que es WS = 195 000 lb.