Estudio de Casos

CAPÍTULO 8
Estudio de casos:
raíces de ecuaciones
La finalidad de este capítulo es utilizar los procedimientos numéricos analizados en
los capítulos 5, 6 y 7 para resolver problemas de ingeniería reales. Las técnicas
numéricas son importantes en aplicaciones prácticas, ya que con frecuencia los
ingenieros encuentran problemas que no es posible resolver usando técnicas
analíticas. Por ejemplo, modelos matemáticos simples que se pueden resolver
analíticamente quizá no sean aplicables cuando se trata de problemas reales. Debido
a esto, se deben utilizar modelos más complicados. En esta situación, es conveniente
implementar una solución numérica en una computadora. En otros casos, los
problemas de diseño en la ingeniería llegan a requerir soluciones de variables
implícitas en ecuaciones complicadas.
Las siguientes aplicaciones son típicas de aquellas que en forma rutinaria se
encuentran durante los últimos años de estudio y en estudios superiores. Más aún,
son problemas representativos de aquellos que se encontrarán en la vida
profesional. Los problemas provienen de las cuatro grandes ramas de la ingeniería:
química, civil, eléc- trica y mecánica. Dichas aplicaciones también sirven para ilustrar
las ventajas y desven- tajas de las diversas técnicas numéricas.
La primera aplicación, tomada de la ingeniería química, proporciona un excelente
ejemplo de cómo los métodos para determinar raíces permiten usar fórmulas realistas
en la ingeniería práctica; además, demuestra de qué manera la eficiencia del método
de Newton-Raphson se emplea cuando se requiere de un gran número de cálculos
como método para la localización de raíces.
Los siguientes problemas de diseño en ingeniería se toman de las ingenierías
civil, eléctrica y mecánica. En la sección 8.2 se usan tanto métodos cerrados como
abiertos para determinar la profundidad y velocidad del agua que fluye en un canal
abierto. En la sección 8.3 se explica cómo las raíces de ecuaciones trascendentes se
usan en el dise- ño de un circuito eléctrico. En las secciones 8.2 y 8.3 también se
muestra de qué forma los métodos gráficos ofrecen un conocimiento del proceso de
localización de raíces. Por último, la sección 8.4 usa la localización de raíces
polinominales para analizar las vibraciones de un automóvil.
8.1 LEYES DE LOS GASES IDEALES Y NO IDEALES

(INGENIERÍA QUÍMICA Y BIOQUÍMICA)
Antecedentes. La ley de los gases ideales está dada por

pV = nRT (8.1)
donde p es la presión absoluta, V es el volumen, n es el número de moles, R es la
constante universal de los gases y T es la temperatura absoluta. Aunque esta
ecuación se utiliza
20 ESTUDIO
8.1 LEYES
DE CASOS:
DE LOSRAÍCES
GASES DEIDEALES Y NO 20
ampliamente por los ingenieros y científicos, sólo es exacta en un rango limitado de

presión y temperatura. Además, la ecuación (8.1) es apropiada solamente para algunos
gases.
Una ecuación de estado alternativa para los gases está dada por:
⎛ a ⎞
p+ (v – b) = (8.2)
RT
⎝ v ⎠
2
conocida como la ecuación de van der Waals, donde v = V/n es el volumen molar, a y
b son constantes empíricas que dependen del gas que se analiza.
Un proyecto de diseño en ingeniería química requiere que se calcule exactamente
el volumen molar (v) del dióxido de carbono y del oxígeno para diferentes
combinaciones de temperatura y presión, de tal forma que los recipientes que
contengan dichos gases se puedan seleccionar apropiadamente. También es
importante examinar qué tan bien se apega cada gas a la ley de los gases ideales,
comparando el volumen molar calcu- lado con las ecuaciones (8.1) y (8.2). Se
proporcionan los siguientes datos:
R = 0.082054 L atm/(mol K)
a = 3.592
bióxido de carbono
b = 0.04267
a = 1.360
oxígeno
b = 0.03183
Las presiones de diseño de interés son de 1, 10 y 100 atmósferas para combinaciones
de temperatura de 300, 500 y 700 K.
Solución. Los volúmenes molares de ambos gases se calculan usando la ley de los
gases ideales, con n = 1. Por ejemplo, si p = 1 atm y T = 300 K,
V RT L am 300 K
v= = = 0.082054 = 24.6162 L/mol
n p mol K 1 atm
Estos cálculos se repiten para todas las combinaciones de presión y de temperatura
que se presentan en la tabla 8.1.
TABLA 8.1 Cálculos del volumen molar.
Volumen Volumen molar Volumen molar

molar (ley de los (van der Waals) (van der
Waals)
Temperatura, Presión, gases ideales), Dióxido de carbono, Oxígeno,
K atm L/mol L/mol L/mol
300 1 24.6162 24.5126 24.5928

10 2.4616 2.3545 2.4384
100 0.2462 0.0795 0.2264
500 1 41.0270 40.9821 41.0259
10 4.1027 4.0578 4.1016
100 0.4103 0.3663 0.4116
700 1 57.4378 57.4179 57.4460
10 5.7438 5.7242 5.7521
100 0.5744 0.5575 0.5842
20 ESTUDIO
8.1 LEYES
DE CASOS:
DE LOSRAÍCES
Los cálculos del volumen molar a partir de la ecuación de van der Waals se llevan
a cabo usando cualquiera de los métodos numéricos para la determinación de raíces de
ecuaciones analizados en los capítulos 5, 6 y 7, con
⎛ a ⎞
ƒ (v) = p + (v – b) – (8.3)
RT ⎝ v ⎠
2
En este caso, como la derivada de ƒ(v) se determina fácilmente, entonces es

conveniente y eficiente usar el método de Newton-Raphson. La derivada de ƒ(v)
respecto a v está dada por
a
ƒ′(v) = p – + (8.4)
2 ab
2 3
v v
El método de Newton-Raphson se describe mediante la ecuación (6.6):
ƒ (v i )
v i+1 = v i –
ƒ′(v i )
la cual se utiliza para estimar la raíz. Por ejemplo, usando como valor inicial 24.6162,
el volumen molar del bióxido de carbono a 300 K y 1 atmósfera es 24.5126 L/mol.
Este resultado se obtuvo después de sólo dos iteraciones y tiene un ea menor del 0.001
por ciento.
En la tabla 8.1 se muestran resultados similares para todas las combinaciones de
presión y de temperatura de ambos gases. Se observa que los resultados obtenidos con
la ecuación de los gases ideales difieren de aquellos obtenidos usando la ecuación de
van der Waals, para ambos gases, dependiendo de los valores específicos de p y T.
Ade- más, como algunos de dichos resultados son significativamente diferentes, el
diseño de los recipientes que contendrán a los gases podría ser muy diferente,
dependiendo de qué ecuación de estado se haya empleado.
En este problema, se examinó una complicada ecuación de estado con el método
de Newton-Raphson. En varios casos los resultados variaron de manera significativa
respecto a la ley de los gases ideales. Desde un punto de vista práctico, el método de
Newton-Raphson fue apropiado aquí, ya que ƒ′(v) resultó sencillo de calcular. De esta
manera, es factible explotar las propiedades de rápida convergencia del método de
Newton-Raphson.
Además de demostrar su poder en un solo cálculo, este problema de diseño
muestra cómo el método de Newton-Raphson es especialmente atractivo cuando se
requiere una gran cantidad de cálculos. Debido a la velocidad de las computadoras
digitales, la eficiencia de varios métodos numéricos en la solución para la mayoría de
las raíces de ecuaciones no se distingue en un cálculo único. Incluso una diferencia de
1 s entre el método de bisección y el eficiente método de Newton-Raphson no
significa pérdida de tiempo cuando se realiza sólo un cálculo. Sin embargo, suponga
que para resolver un problema se necesita calcular millones de raíces. En tal caso, la
eficiencia del método podría ser un factor decisivo al elegir una técnica.
Por ejemplo, suponga que se requiere diseñar un sistema de control
computarizado automático para un proceso de producción de sustancias químicas.
Dicho sistema requiere una estimación exacta de volúmenes molares sobre una base
20 esencialmente
ESTUDIO
8.1 LEYES DE continua,
DE CASOS:LOSRAÍCES
GASES para fabricar
DEIDEALES Y NOen forma conveniente el producto final.20Se
instalan medidores
20 ESTUDIO
8.1 LEYES
DE CASOS:
DE LOSRAÍCES
que proporcionan lecturas instantáneas de presión y temperatura. Se debe obtener

valores de v para diversos gases que se usan en el proceso.
Para una aplicación como ésta, los métodos cerrados, tales como el de bisección o
de la regla falsa, posiblemente consumirían mucho tiempo. Además, los dos valores
iniciales que se requieren en estos métodos generarían un retraso crítico en el procedi-
miento. Dicho inconveniente afecta de igual forma al método de la secante, que
también necesita dos valores iniciales.
En contraste, el método de Newton-Raphson requiere sólo de un valor inicial para
determinar la raíz. La ley de los gases ideales podría emplearse para obtener un valor
inicial del proceso. Después, suponiendo que el tiempo empleado sea lo bastante corto
como para que la presión y la temperatura no varíen mucho entre los cálculos, la solu-
ción de la raíz anterior se puede usar como un buen valor inicial para la siguiente
aplicación. De esta forma, se tendría de forma automática un valor inicial cercano a la
solución, que es requisito indispensable para la convergencia del método de Newton-
Raphson. Todas estas consideraciones favorecerán de buena manera la técnica de
New- ton-Raphson en estos problemas.
8.2 FLUJO EN UN CANAL ABIERTO

(INGENIERÍA CIVIL E INGENIERÍA AMBIENTAL)
Antecedentes. La ingeniería civil constituye una disciplina amplia que incluye

diversas áreas como estructural, geotecnia, transporte, ambiental y abastecimiento
del agua. Las dos últimas especialidades tienen que ver con la contaminación y
suministro de agua y, por lo tanto, implican un uso extensivo de la ciencia de
mecánica de fluidos.
Un problema general se relaciona con el flujo de agua en canales abiertos, ríos y
canales. La velocidad de flujo, que se mide frecuentemente en la mayoría de los
ríos y arroyos, se define como el volumen de agua que pasa por un punto específico
de un canal por unidad de tiempo, Q (m3/s).
Aunque la velocidad de flujo es una cantidad útil, una cuestión adicional se
relaciona con lo que sucede cuando se tiene una velocidad de flujo específico en un
canal con pendiente (figura 8.l). De hecho, suceden dos cosas: el agua alcanzará una
profundidad específica H (m) y se moverá a una velocidad específica U (m/s). Los
ingenieros ambientales pueden estar interesados en conocer tales cantidades para
predecir el transporte y el destino de los contaminantes en un río. Así, la pregunta
general sería: si se tiene una velocidad de flujo para un canal, ¿cómo se calculan la
profundidad y la velocidad?
FIGURA 8.1
P
Ac
B
S
H Q, U
20 ESTUDIO
8.2 FLUJO
DE CASOS:
EN UNRAÍCES
CANAL DE 20
Solución. La relación fundamental entre flujo y profundidad es la ecuación de

continuidad
Q = UAc (8.5)
donde Ac = área de la sección transversal del canal (m2). Dependiendo de la forma del
canal, el área puede relacionarse con la profundidad por medio de varias expresiones
funcionales. Para el canal rectangular mostrado en la figura 8.1, Ac = BH. Al sustituir
esta expresión en la ecuación (8.5) se obtiene
Q = UBH (8.6)
donde B = ancho (m). Debe observarse que la ecuación de continuidad se obtiene de la

conservación de la masa (recuerde la tabla 1.1).
Ahora, aunque la ecuación (8.6) ciertamente relaciona los parámetros del canal,
no es suficiente para responder nuestra pregunta. Suponiendo que se conoce B, se
tiene una ecuación y dos incógnitas (U y H). Por lo tanto, se requiere una ecuación
adicional. Para flujo uniforme (significa que el flujo no varía con la distancia ni con el
tiempo), el ingeniero irlandés Robert Manning propuso la siguiente fórmula
semiempírica (llamada en forma apropiada ecuación de Manning)
1 2/3 1/2
U= R S (8.7)
n
donde n = coeficiente de rugosidad de Manning (un número adimensional que toma en

cuenta la fricción del canal), S = pendiente del canal (adimensional, metros de caída
por longitud en metros) y R = radio hidráulico (m), el cual se relaciona con los
parámetros fundamentales mediante
Ac
R= (8.8)
P
donde P = perímetro mojado (m). Como su nombre lo indica, el perímetro mojado es

la longitud de los lados y el fondo del canal que está bajo el agua. Por ejemplo, para un
canal rectangular, éste se define como
P = B + 2H (8.9)
Se debe observar que así como la ecuación de continuidad se obtiene de la

conservación de la masa, la ecuación de Manning es una expresión de la conservación
del momentum. En particular, indica cómo la velocidad depende de la rugosidad, una
manifestación de la fricción.
Aunque el sistema de ecuaciones no lineales (8.6 y 8.7) puede resolverse simultá-
neamente (por ejemplo, usando el método de Newton-Raphson multidimensional que
se describe en la sección 6.5.2), un método más simple sería la combinación de
ecuaciones. La ecuación (8.7) se sustituye en la ecuación (8.6) y se obtiene
BH 2/3 1/2
Q= R S (8.10)
n
20 ESTUDIO
8.2 FLUJO
DE CASOS:
EN UNRAÍCES
CANAL DE 20
Así, el radio hidráulico, ecuación (8.8), junto con las diferentes relaciones para
un canal rectangular, se sustituye:
1/2 5/3
S ( BH )
Q= 2/3 (8.11)
n (B + 2H)
De esta forma, la ecuación contiene ahora una sola incógnita H junto con el valor
dado de Q y los parámetros del canal (n, S y B).
Aunque se tiene una ecuación con una incógnita, es imposible resolverla en forma
explícita para encontrar H. Sin embargo, la profundidad se determina numéricamente,
al reformular la ecuación como un problema de raíces.
1/2
S ( BH )
5/3
ƒ(H) = 2/3 –Q=0 (8.12)
n (B + 2H)
La ecuación (8.12) se resuelve rápidamente con cualquiera de los métodos para

localizar raíces, descritos en los capítulos 5 y 6. Por ejemplo, si Q = 5 m3/s, B = 20 m,
n = 0.03 y S = 0.0002, la ecuación es
5/3
(20 H )
ƒ(H) = 0.471405 2/3 – 5 = 0 (8.13)
(20 + 2H)
Puede resolverse para H = 0.7023 m. El resultado se verifica sustituyéndolo en la

ecua- ción (8.13):
5/3
(20 × 0.7023) –5
ƒ(H) = 0.471405 2/3 – 5 = 7.8 ×10 (8.14)
(20 + 2 × 0.7023)
que se acerca bastante a cero.

La otra incógnita, la velocidad, ahora se determina por sustitución en la ecuación
(8.6),
Q 5
U= = = 0.356 m/s (8.15)
BH 20(0.7023)
Así, se tiene una solución satisfactoria para la profundidad y la velocidad.
Ahora se buscará analizar un poco más los aspectos numéricos de este problema.
Una pregunta pertinente sería: ¿Cómo hacer para obtener un buen valor inicial para el
método numérico? La respuesta depende del tipo de método.
Para los métodos cerrados, como el de bisección y el de la falsa posición, se deter-
minaría, si es posible, estimar valores iniciales inferiores y superiores que contengan
siempre una sola raíz. Un método conservador podría ser elegir cero como el límite in-
ferior. Y, si se conoce, la profundidad máxima posible que puede presentarse, este
valor serviría como valor inicial superior. Por ejemplo, todos los ríos, con excepción
de los más grandes del mundo, tienen menos de 10 metros de profundidad. Por lo
tanto, se toman 0 y 10 como límites del intervalo para H.
Si Q > 0 y H = 0, la ecuación (8.12) siempre será negativa para el valor inicial in-
ferior. Conforme H se incrementa, la ecuación (8.12) también se incrementará en
forma
20 ESTUDIO
8.2 FLUJO
DE CASOS:
EN UNRAÍCES
CANAL DE 20
monótona, y finalmente será positiva. Por lo tanto, los valores iniciales deberán
conte- ner una sola raíz en la mayoría de los casos que se estudian con ríos y arroyos
natu- rales.
Ahora, una técnica como la de bisección debería ser muy confiable en la
búsqueda de una raíz. ¿Pero qué precio se paga? Al usar tal ancho del intervalo y una
técnica como la de bisección, el número de iteraciones para obtener una precisión
deseada podría ser computacionalmente excesivo. Por ejemplo, si se elige una
tolerancia de 0.001 m, la ecuación (5.5) sirve para calcular
log(10 / 0.001)
n= log 2 =
13.3
Así, se requieren 14 iteraciones. Aunque esto ciertamente no sería costoso para un

solo cálculo, podría ser exorbitante si se efectuaran muchas de estas evaluaciones. Las
alter- nativas serían: estrechar el intervalo inicial (en base a un conocimiento
específico del sistema), usar un método cerrado más eficiente (como el de la falsa
posición) o confor- marse con una menor precisión.
Otra forma de tener una mejor eficiencia sería utilizar un método abierto como el
de Newton-Raphson o el de la secante. Por supuesto que en tales casos el problema de
los valores iniciales se complica al considerar la convergencia.
Se obtiene una mayor comprensión de este problema examinando al menos
eficiente de los métodos abiertos: iteración de punto fijo. Al analizar la ecuación
(8.11), se observa que hay dos modos sencillos para despejar H; esto es, se resuelve
tanto para H en el numerador,
3/5 2/5
(Qn) (B+ 2H)
H= 3/10 (8.16)
BS
como para H en el denominador,

3 5/2
1 ⎡ S ( BH ) ⎤
H= ⎢ 3/2
– B⎥ (8.17)
2⎣ (Qn) ⎦
Ahora, aquí es donde el razonamiento físico puede ayudar. En la mayoría de los

ríos y arroyos, el ancho es mucho mayor que la profundidad. Así, la cantidad B + 2H
no varía mucho. De hecho, debe ser aproximadamente igual a B. Por lo contrario, BH
es directamente proporcional a H. En consecuencia, la ecuación (8.16) deberá
converger más rápido a la raíz, lo cual se verifica al sustituir los límites del intervalo
H = 0 y 10 en ambas ecuaciones. Con la ecuación (8.16), los resultados son 0.6834 y
0.9012, que son cercanos a la raíz verdadera, 0.7023. En contraste, los resultados con
la ecuación (8.17) son –10 y 8 178, los cuales están alejados claramente de la raíz.
La superioridad de la ecuación (8.16) se manifiesta además al graficar sus
componentes (recuerde la figura 6.3). Como se observa en la figura 8.2, la
componente g(H) de la ecuación (8.16) es casi horizontal. Así, esta ecuación no
únicamente converge, sino que debe hacerlo con rapidez. En cambio, la componente
g(H) de la ecuación (8.17) es casi vertical, indicando así una fuerte y rápida
divergencia.
20 ESTUDIO
8.2 FLUJO
DE CASOS:
EN UNRAÍCES
CANAL DE 20
y y
4 4
y2 = g(H)
2 y1 = H 2 y1 = H
FIGURA 8.2 y2 = g(H)

Gráfica de los componentes
para dos casos de iteración
de punto fijo, uno que 0 0
0 1 2 H 0 1 2 H
converge [a), ecuación
(8.16)] y uno que diverge a) b)
[b), ecuación (8.17)].
Hay dos beneficios prácticos de este análisis:
1. En el caso de que se use un método abierto más detallado, la ecuación (8.16)

ofrece un medio para obtener un excelente valor inicial. Por ejemplo, si H se elige
como cero, la ecuación (8.12) toma la forma
3/5
H0 = (Qn/3/10
B)
S
donde H0 será el valor inicial utilizado en el método de Newton-Raphson o en el
de la secante.
2. Se ha demostrado que la iteración de punto fijo ofrece una opción viable para este
problema específico. Por ejemplo, usando como valor inicial H = 0, en la ecuación
(8.16) se obtienen seis dígitos de precisión en cuatro iteraciones para el caso que
se examina. La fórmula de iteración de punto fijo sería fácil de manipular en una
hoja de cálculo, ya que las hojas de cálculo son ideales para fórmulas iterativas
conver- gentes que dependen de una sola celda.
8.3 DISEÑO DE UN CIRCUITO ELÉCTRICO

(INGENIERÍA ELÉCTRICA)
Antecedentes. Los ingenieros eléctricos emplean las leyes de Kirchhoff para

estudiar el comportamiento de los circuitos eléctricos en estado estacionario (que no
varía con el tiempo). En la sección 12.3 se analiza el comportamiento de dichos
estados estacionarios. Otro problema importante tiene que ver con circuitos de
naturaleza transitoria, donde súbitamente ocurren cambios temporales. Esta
situación se presenta cuando se cierra el interruptor como en la figura 8.3. En tal caso,
existe un periodo de ajuste al cerrar el interruptor hasta que se alcance un nuevo
estado estacionario. La longitud de este pe-
20 ESTUDIO
8.3 DISEÑO
DE CASOS:
DE UN
RAÍCES
CIRCUITO
DE 20
Interruptor
i
– –
Batería V0 Capacitor Inductor
+ +
Resistor
FIGURA 8.3
Un circuito eléctrico. Cuando se cierra el interruptor, la corriente experimenta una serie de
oscilaciones hasta que se alcance un nuevo estado estacionario.
riodo de ajuste está íntimamente relacionada con las propiedades de almacenamiento

de energía, tanto del capacitor como del inductor. La energía almacenada puede
oscilar entre estos dos elementos durante un periodo transitorio. Sin embargo, la
resistencia en el circuito disipará la magnitud de las oscilaciones.
El flujo de corriente a través del resistor provoca una caída de voltaje (VR), dada
por
VR = iR
donde i = la corriente y R = la resistencia del resistor. Si las unidades de R e i son

ohms y amperes, respectivamente, entonces las unidades de VR son voltios.
De manera semejante, un inductor se opone a cambios de corriente tales que la
caída del voltaje a través del inductor V L es
di
VL = L
dt
donde L = la inductancia. Si las unidades de L e i son henrios y amperes,

respectivamente, entonces las de VL son voltios, y las de t son segundos.
La caída del voltaje a través del capacitor (VC) depende de la carga (q) sobre éste:
= q
V
C
C
donde C = la capacitancia. Si las unidades de carga se expresan en coulombios,

entonces la unidad de C es el faradio.
La segunda ley de Kirchhoff establece que la suma algebraica de las caídas de
voltaje alrededor de un circuito cerrado es cero. Así que, después de cerrar el
interruptor se tiene
di q
L + Ri + = 0
dt C
Sin embargo, como la corriente se relaciona con la carga de acuerdo con
dq
i = dt
20 ESTUDIO
8.3 DISEÑO
DE CASOS:
DE UN
RAÍCES
CIRCUITO
DE 20
Por lo tanto,
q(t)
q0
d + R dq + 1 q = 0
2
Ldt 2 (8.18)
dt C
Tiempo
Ésta es una ecuación diferencial ordinaria lineal de segundo orden que se resuelve
usando los métodos de cálculo (véase la sección 8.4). Esta solución está dada por
FIGURA 8.4 ⎡ 1 ⎛ R ⎞2
– Rt /(2 L)
La carga en un capacitor q(t) = q e cos ⎢ – t⎥ (8.19)
0 ⎤
LC ⎝ 2L ⎠
como función del tiempo
después de cerrar el
interruptor de la figura 8.3. si en t = 0, q = q0 = V0 C y V0 = el voltaje de la batería. La ecuación (8.19) describe la
variación de la carga en el capacitor. La solución q(t) se grafica en la figura 8.4.
Un problema de diseño típico en ingeniería eléctrica consistiría en la
determinación del resistor apropiado para disipar energía a una razón especificada,
con valores conocidos de L y C. En este problema, suponga que la carga se debe
disipar a 1% de su valor original (q/q 0 = 0.01) en t = 0.05 s, con L = 5 H y C = 10 –4F.
Solución. Es necesario despejar R de la ecuación (8.19) con valores conocidos

para q, q0, L y C. Sin embargo, debe emplear una técnica de aproximación numérica,
ya que R es una variable implícita en la ecuación (8.19). Se usará el método de
bisección para dicho propósito. Los otros métodos estudiados en los capítulos 5 y 6
también son apropiados; aunque el método de Newton-Raphson tiene el
inconveniente de que la derivada de la ecuación (8.19) es un poco complicada.
Reordenando la ecuación (8.19),
⎡ 1 q
2
– Rt /(2 L) ⎛ R ⎞
ƒ = ⎤
(R) e cos ⎢ – t⎥ –
⎢ LC ⎝ 2L ⎥ q0
⎣ ⎠ ⎦
Utilizando los valores numéricos dados,
cos [ 2 000 – 0.01R 2 (0.05)] – 0.01

–0.005 R
ƒ(R) = e (8.20)
Un examen de esta ecuación sugiere que un rango inicial razonable para R es 0 a 400
Ω (ya que 2 000 – 0.01R2 debe ser mayor que cero). La figura 8.5 es una gráfica
de la ecuación (8.20), que confirma lo anterior. Al hacer veintiún iteraciones con el
método de bisección se obtiene una raíz aproximada R = 328.1515 Ω, con un error
menor al
0.0001 por ciento.
De esta forma, se especifica un resistor con este valor para el circuito mostrado en
la figura 8.6 y se espera tener una disipación consistente con los requisitos del
problema. Este problema de diseño no se podría resolver eficientemente sin el uso de
los métodos numéricos vistos en los capítulos 5 y 6.
20 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 20
f (R)
Raíz 325
0.0
200 400 R
– 0.2
– 0.4
– 0.6
FIGURA 8.5
Gráfica de la ecuación (8.20) usada para obtener los valores iniciales
de R que contienen a la raíz.
8.4 ANÁLISIS DE VIBRACIONES (INGENIERÍA MECÁNICA

E INGENIERÍA AERONÁUTICA)
Antecedentes. Las ecuaciones diferenciales sirven para modelar la vibración de

sistemas en ingeniería. Algunos ejemplos (figura 8.6) son el péndulo simple, una
masa sujeta a un resorte y un circuito eléctrico con un inductor y un capacitor
(recuerde la sección 8.3). La vibración de estos sistemas puede amortiguarse por
medio de algún
FIGURA 8.6
Tres ejemplos de vibraciones armónicas simples. Las flechas dobles indican las vibraciones
en cada sistema.
Corriente
Péndulo Resorte/masa Circuito LC

21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
–x
Amortiguador
Resorte
+x
m
FIGURA 8.7
Un carro de masa m.
mecanismo que absorba la energía. Además, la vibración puede ser libre o sujeta a
algún disturbio periódico externo. En este último caso, se dice que el movimiento es
forzado. En esta sección se examinará la vibración libre y forzada del automóvil, que
se muestra en la figura 8.7. El tratamiento general es aplicable a muchos otros
problemas de inge- niería.
Como se observa en la figura 8.7, un carro de masa m se soporta por medio de re-
sortes y amortiguadores. Los amortiguadores presentan resistencia al movimiento, que
es proporcional a la velocidad vertical (movimiento ascendente-descendente). La
vibra- ción libre ocurre cuando el automóvil es perturbado de su condición de
equilibrio, como ocurre cuando se pasa por un bache (agujero en el camino). Un
instante después de pasar por el bache, las fuerzas netas que actúan sobre m son la
resistencia de los resortes y la fuerza de los amortiguadores. Tales fuerzas tienden a
regresar el carro al estado de equilibrio original. De acuerdo con la ley de Hooke, la
resistencia del resorte es proporcional a su constante k y a la distancia de la posición
de equilibrio x. Por lo tanto,
Fuerza del resorte = –kx
donde el signo negativo indica que la fuerza de restauración actúa regresando el auto-
móvil a su posición de equilibrio (es decir, la dirección x negativa). La fuerza para un
amortiguador está dada por
dx
Fuerza de amoriguación = –c
dt
donde c es el coeficiente de amortiguamiento y dx/dt es la velocidad vertical. El signo

negativo indica que la fuerza de amortiguamiento actúa en dirección opuesta a la
velocidad.
Las ecuaciones de movimiento para el sistema están dadas por la segunda ley de
Newton (F = ma), que en este problema se expresa como
d2x dx
m × = –c + (–kx)
dt2 dt
Masa × aceleración = fuerza de amortiguamiento + fuerza del resorte
21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
o bien
d x2 dx
m 2 +c + kx = 0
dt dt
Observe la similitud con la ecuación (8.18) que se desarrolló en la sección 8.3 para un
circuito eléctrico.
Si se supone que la solución toma la forma x(t) = ert, entonces se escribe la
ecuación
característica
mr2 + cr + k = 0 (8.21)
La incógnita r es la solución de la ecuación característica cuadrática que se puede

obtener, ya sea en forma analítica o numérica. En este problema de diseño, primero
se utiliza la solución analítica para ofrecer una idea general de la forma en que el
movimiento del sistema es afectado por los coeficientes del modelo: m, k y c.
También se usarán diferentes métodos numéricos para obtener las soluciones, y se
verificará la exactitud de los resultados con la solución analítica. Por último,
sentaremos las bases para problemas más complicados que se describirán más tarde
en el texto, donde los resultados analíticos son difíciles o imposibles de obtener.
La solución de la ecuación (8.21) para r está dada por la fórmula cuadrática
r1 –c ± c 2 – 4mk
= (8.22)
r2
2m
Note el significado de la magnitud de c al compararla con 2 km . Si c > 2 km ,

r1 y r 2 son números reales negativos, y la solución es de la forma
x(t) = Aer1t + Ber2t (8.23)
donde A y B son constantes que se deben determinar a partir de las condiciones

iniciales de x y dx/dt. Tales sistemas se denominan sobreamortiguados.
Si c < 2 km , las raíces son complejas,
r1
= λ ± µi
r2
donde
⏐c 2 – 4 mk⏐
µ=
2m
y la solución es de la forma
x(t) = e–lt (A cos µt + B sen µt) (8.24)
Tales sistemas se conocen como subamortiguados.

Por último, si c = 2 km , la ecuación característica tiene una raíz doble y la
solución es de la forma
x(t) = (A + Bt)e–lt (8.25)
21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
donde
c
λ=
2m
A tales sistemas se les llama críticamente amortiguados.
En los tres casos, x(t) se aproxima a cero cuando t tiende al infinito. Esto
significa que el automóvil siempre regresa a la posición de equilibrio después de pasar
por un bache (¡aunque esto parecería poco probable en algunas ciudades que hemos
visitado!). Estos casos se ilustran en la figura 8.8.
El coeficiente de amortiguamiento crítico cc es el valor de c que hace que el
radical de la ecuación (8.22) sea igual a cero,
cc = 2 km o cc = 2mp (8.26)
donde
k (8.27)
p= m
La relación c/cc se llama factor de amortiguamiento, y a p se le conoce como la

frecuencia natural de la vibración libre no amortiguada.
Ahora, consideremos el caso donde el automóvil está sujeto a una fuerza
periódica dada por
P = Pm sen wt o d = dm sen wt
donde dm = Pm /k = la deflexión estática del carro sujeto a una fuerza Pm. La

ecuación diferencial que rige este caso es
d x dx
m 22 + c + kx = Pm sen ω t
dt dt
La solución general de esta ecuación se obtiene al sumar una solución particular a

la solución por vibración libre, dada por las ecuaciones (8.23), (8.24) y (8.25).
Conside-
FIGURA 8.8
x(t)
Vibraciones a) sobreamor-
tiguadas, b) subamortigua- Amortiguamiento
das y c) amortiguadas crítico
críticamente. Sobreamortiguamiento
Subamortiguamiento
21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
remos el movimiento en estado estacionario del sistema forzado donde se ha

amortiguado el movimiento transitorio inicial. Si consideramos que esta solución
en estado estacionario tiene la forma
xss (t) = xm sen (wt – f)

se demuestra que
xm xm 1
= = (8.28)
2 2 2
La cantidad xm /dm llamada factor de amplificación de la amplitud depende tan sólo de

la razón del amortiguamiento real con el amortiguamiento crítico, y de la razón de la
frecuencia forzada con la frecuencia natural. Observe que cuando la frecuencia forzada
w se aproxima a cero, el factor de amplificación se aproxima a 1. Si, además, el
sistema es ligeramente amortiguado, es decir, si c/cc es pequeño, entonces el factor de
amplifica- ción se hace grande cuando w es cercano a p. Si el amortiguamiento es
cero, entonces el factor de amplificación tiende a infinito cuando w = p, y se dice que
la función de fuerza entra en resonancia con el sistema. Por último, conforme w/p se
vuelve muy grande, el factor de amplificación se aproxima a cero. La figura 8.9
muestra una gráfica del factor de amplificación como una función de w/p para
diversos factores de amortiguamiento.
Observe que el factor de amplificación se conserva pequeño al seleccionar un
factor de amortiguamiento grande, o manteniendo muy distantes las frecuencias
natural y forzada.
El diseño del sistema de suspensión del automóvil comprende una solución
interme- dia entre comodidad y estabilidad para todas las condiciones de manejo y
velocidad. Se pide determinar la estabilidad del carro para cierto diseño propuesto que
ofrezca comodidad sobre caminos irregulares. Si la masa del carro es m = 1.2 × 106
gramos y tiene un sistema de amortiguadores con un coeficiente de amortiguamiento c
= 1 × 107 g/s.
Suponga que la expectativa del público en cuanto a la comodidad se satisface si la
vibración libre del automóvil es subamortiguada y el primer cruce por la posición de
equilibrio tiene lugar en 0.05 s. Si en t = 0, el carro súbitamente se desplaza una
distan-
FIGURA 8.9
Gráfica del factor de 6 c/cc = 0
amplificación de la
amplitud xm/xd [ecuación 0.125
(8.28)] contra la frecuencia 4
w entre la frecuencia 0.25
xm
natural p para diversos 0.5

valores del coeficiente 2
1
de amortiguamiento c
entre el coeficiente de
amortiguamiento crítico cc. 0
0 1 2
/p
21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
cia x0, desde el equilibrio, y la velocidad es cero (dx/dt = 0), la solución de la ecuación
de movimiento está dada por la ecuación (8.24), con A = x0 y B = x0l/m. Por lo tanto,
– λt ⎛ λ ⎞
x(t) = x 0 e ⎜ cos µt + sen µt⎟
⎝ µ ⎠
Nuestras condiciones de diseño se satisfacen si

λ
x(t) = 0 = cos (0.05µ) + sen (0.05µ)
µ
o bien
⎛ c ⎞ ⎛ c ⎞
2 2
k c k
0 = cos ⎜ 0.05 – 2 ⎟ + sen ⎜ 0.05 – (8.29)
⎝ m 4m ⎠ 4km – c
2
⎝ m 4m 2 ⎟⎠
Dado que se conocen c y m, el problema de diseño consiste ahora en encontrar valores

apropiados de k que satisfagan la ecuación (8.29).
Solución. Se pueden utilizar los métodos de la bisección, de la falsa posición o de la

secante, ya que esos métodos no requieren la evaluación de la derivada de la ecuación
(8.29), la cual podría resultar algo difícil de calcular en este problema. La solución es
k = 1.397 × 109, con 12 iteraciones, utilizando el método de bisección con un intervalo
inicial que va de k = 1 × 109 a 2 × 109 (ea = 0.07305%).
Aunque este diseño satisface los requerimientos de vibración libre (después de
caer en un bache), también debe probarse bajo las condiciones de un camino
accidentado. La superficie del camino se puede aproximar como
2πx ⎞
d = d m sen ⎛⎝
D ⎠
donde d es la deflexión, dm es la máxima deflexión de 0.1 m y D es la distancia entre

los picos que es igual a 20 m. Si v es la velocidad horizontal del automóvil (m/s),
entonces la ecuación de movimiento del sistema se escribe como
d x dx 2πv ⎞
m 22 + c + kx = kdm sen ⎛ t
dt dt ⎝ D ⎠
donde w = 2pv/D es la frecuencia forzada.

La estabilidad del carro se considera satisfactoria si en estado estacionario la
máxi- ma distancia xm es inferior a 0.2 m para todas la velocidades de manejo. El
factor de amortiguamiento se calcula de acuerdo con la ecuación (8.26)
7
c 10 1 × 10
= = = 0.1221
cc 2 km 2 1.397 ×10 9 (1.2 ×10 6 )
21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
Ahora, se buscan valores w/p que satisfagan la ecuación (8.28),

1
2= 2 2 2 2 (8.30)
[1 – (ω/ p) ] + 4(0.1221) (ω/ p)
Si la ecuación (8.30) se expresa como un problema de raíces

ƒ(ω / p) = 2 [1 – (ω / p) 2 ]2 + 4(0.1221) 2 (ω / p) 2 – 1 = (8.31)
0
Vea que los valores w/p se determinan al encontrar las raíces de la ecuación (8.31).
Una gráfica de la ecuación (8.31) se presenta en la figura 8.10. En ésta se muestra
que la ecuación (8.31) tiene dos raíces positivas que se pueden determinar con el
méto- do de bisección, usando el software TOOLKIT. El valor más pequeño para w/p
es igual a 0.7300 en 18 iteraciones, con un error estimado de 0.000525% y con valores
iniciales superior e inferior de 0 y 1. El valor mayor que se encuentra para w/p es de
1.1864 en
17 iteraciones, con un error estimado de 0.00064% y con valores iniciales superior e
inferior de 1 y 2.
También es posible expresar la ecuación (8.30) como un polinomio:
⎛ ω ⎞4 ⎛ω⎞
2
–
⎜ ⎝ ⎟p ⎠ 1.9404 + 0.75 (8.32)
⎝⎜ p ⎟⎠
y usar MATLAB para determinar las raíces como sigue:
>> a=[l 0 -1.9404 0 .75];

>> roots (a)
ans =
1.1864
-1.1864
0.7300
-0.7300
Lo cual confirma el resultado obtenido con el método de bisección. Esto también

sugiere que, aunque la ecuación (8.32) es una ecuación de cuarto grado en w/p,
también es una ecuación cuadrática en (w/p)2.
El valor de la frecuencia natural p está dado por la ecuación (8.27),
9
1.397 × 10 –1
p= 6 = 34.12 s
1.2 ×10
Las frecuencias forzadas, para las que la máxima deflexión es 0.2 m, entonces se
calculan como
w = 0.7300(34.12) = 24.91 s–1

w = 1.1864(34.12) = 40.48 s–1
21 ESTUDIO
8.4 ANÁLISIS
DE CASOS:
DERAÍCES DE 21
f( 2
/p
0
1 2
FIGURA 8.10
Gráfica de la ecuación
(8.31) que indica dos raíces –2
positivas.
con lo cual se obtiene

ωD 24.91(20 ) m 3 600 s km
v= = = 79.29 × = 285 km/hr (= 177 mi / hr)
2π 2(3.14159) s hr 1 000 m
ωD 40.48(20) m 3 600 s km
v= = = 128.85 × = 464 km/hr (= 288 mi / hr)
2π 2(3.14159) s hr 1 000 m
Así, con los resultados anteriores y la figura 8.10, se determina que el diseño del
carro propuesto se comportará de forma aceptable para velocidades de manejo
aceptables. Es decir, el diseñador debe estar consciente de que el diseño podría no
cumplir los requerimientos cuando el automóvil viaje a velocidades extremadamente
altas (por ejemplo, en carreras).
Este problema de diseño ha presentado un ejemplo extremadamente simple, pero
que nos ha permitido obtener algunos resultados analíticos que se utilizaron para
evaluar la exactitud de nuestros métodos numéricos para encontrar raíces. Los casos
reales pueden volverse tan complicados que sólo se obtendrían las soluciones a éstos
emplean- do métodos numéricos.
PROBLEMAS
Ingeniería química/Ingeniería bioquímica quiera de los métodos numéricos analizados en los capítulos 5
8.1 Realice el mismo cálculo que en la sección 8.1, pero ahora y
con alcohol etílico (a = 12.02 y b = 0.08407) a una temperatura 6, y realice los cálculos. Justifique la elección de la técnica.
de 400 K y una presión P de 2.5 atm. Compare los resultados 8.2 En ingeniería química, los reactores de flujo tipo tapón (es
con la ley de los gases ideales. Si es posible, utilice el software decir, aquellos en que el fluido va de un extremo al otro con
de su computadora para determinar el volumen molar. Si no, una mezcla mínima a lo largo del eje longitudinal) se usan
use cual- para convertir reactantes en productos. Se ha determinado
que la
21 ESTUDIO
PROBLEM
DE CASOS: RAÍCES DE 21
eficiencia de la conversión algunas veces se mejora donde K = la constante de equilibrio de la reacción y pt = la
recirculando una porción de la corriente del producto, de tal presión total de la mezcla. Si pt = 3.5 atm y k = 0.04, determine
forma que regre- se a la entrada para un paso adicional a través el valor de x que satisfaga la ecuación (P8.3).
del reactor (figura P8.2). La razón de recirculando se define 8.4 La siguiente ecuación permite calcular la concentración de
como un químico en un reactor donde se tiene una mezcla completa:
volumen de fluido que regresa a la entrada c = cent(1 – e–0.04t) + c0e–0.04t

R= volumen que sale del sistema
Si la concentración inicial es c0 = 5 y la concentración de entra-
Suponga que se está procesando una sustancia química A para da es cent = 12, calcule el tiempo requerido para que c sea el
generar un producto B. Para el caso en que A forma a B de 85% de cent.
acuerdo con una reacción autocatalítica (es decir, en la cual uno 8.5 Una reacción química reversible
de los productos actúa como catalizador o estimulante en la
reacción), es posible demostrar que una razón óptima de 2A + B C
recirculación debe satisfacer
se caracteriza por la relación de equilibrio
1 + R(1 – X A ƒ) R +1
ln =
R(1 – X A ƒ R[1 + R(1 – X A ƒ cc
) )] K= 2
c
a b
donde XAƒ es la fracción del reactante A que se convierte en el
producto B. La razón óptima de recirculación corresponde a donde la nomenclatura cn representa la concentración del com-
un reactor de tamaño mínimo necesario para alcanzar el nivel ponente N. Suponga que se define una variable x que representa
deseado de conversión. Utilice un método numérico para deter- el número de moles de C producido. La conservación de la
minar la razón de recirculación necesaria, de manera que se masa se utiliza para reformular la relación de equilibrio como
minimice el tamaño del reactor para una conversión fraccional
de XAƒ = 0.95. ( cc , 0 + x )
K= 2
(c – 2x) (c – x)
a, 0
Alimentación Producto donde el subíndice 0 indica la concentración inicial de cada

Reactor de flujo tipo
tapón componente. Si K = 0.016, ca, 0 = 42, cb, 0 = 28 y cc, 0 = 4,
Reciclaje calcule x.
8.6 Las siguientes reacciones químicas se llevan a cabo en un
sistema cerrado
Figura P8.2
2A + B C
Representación esquemática de un reactor de flujo tipo
tapón con recirculación. A+D C
En equilibrio, éstas pueden caracterizarse por

cc
K1 = 2
8.3 En un proceso de ingeniería química el vapor de agua (H2O) c

a b
cc
se calienta a temperaturas lo suficientemente altas para que una K =
2 ca cd
porción significativa del agua se disocie, o se rompa, para
formar oxígeno (O2) e hidrógeno (H2):
donde la nomenclatura cn representa la concentración del com-
H2O H2 + 1
O ponente N. Si x1 y x2 son el número de moles de C que se
2 2
produ-
cen debido a la primera y segunda reacciones, respectivamente,
Si se asume que ésta es la única reacción que se lleva a cabo, la emplee un método similar al del problema 8.5 para reformular
fracción molar x de H2O que se disocia se representa por las relaciones de equilibrio en términos de las concentraciones
iniciales de los componentes. Después, use el método de
Newton-
21 x 2 pt ESTUDIO
PROBLEM
DE CASOS: RAÍCES Raphson para resolver el par de ecuaciones simultáneas21no
K= (P8.3)DE
1–x 2+x lineales para x1 y x2 si K1 = 4 × 10–4, K2 = 3.7 × 10–2, ca,0 =
50,
21 ESTUDIO
PROBLEM
cb,0 = 20, cc,0 = 5 y cd,0 = 10. Utilice un método gráfico para caída de presión, r es la densidad del fluido, GO es la
pro- poner los valores iniciales. velocidad másica (el cociente del flujo de masa dividido entre
8.7 La ecuación de estado de Redlich-Kwong está dada por el área de la sección transversal), Dp es el diámetro de las
RT a partículas dentro del lecho, µ es la viscocidad del fluido, L es la
p= – longitud del lecho y e es la fracción vacía del lecho.
v–b v(v + b) T
donde R = la constante universal de los gases [= 0.518 kJ/(kg ∆pρ Dp ε 3 (1 – ε )
2 = 150 + 1.75
K)], T = temperatura absoluta (K), p = presión absoluta (kPa) y (1 – ε ) ⎛ D G ⎞
Go L ⎜ p ⎟o
v = volumen de un kg de gas (m3/kg). Los parámetros a y b se
calculan mediante ⎝ µ
2 2.5 Dados los siguientes valores para los parámetros encuentre la
R Tc Tc
a = 0.427 b = 0.0866R fracción vacía e del lecho.
pc pc Dp Go
donde pc = 4 580 kPa y Tc = 191 K. Como ingeniero químico, = 1 000
µ
se ∆PρDp
le pide determinar la cantidad de combustible metano que se 2 = 10
G L
puede almacenar en un tanque de 3 m3 a una temperatura de o
–50°C con una presión de 65 000 kPa. Emplee el método de lo- 8.12 En una sección de tubo, la caída de presión se calcula así:
calización de raíces de su elección para calcular v y luego L ρV 2
determine la masa de metano contenida en el tanque. ∆p = f
2D
8.8 El volumen V de un líquido contenido en un tanque
horizontal cilíndrico de radio r y longitud L está relacionado donde ∆p = caída de presión (Pa), f = factor de fricción, L =
con la profundidad del líquido h por longitud del tubo [m], r = densidad (kg/m3), V = velocidad
(m/s),
⎡ ⎛r–h⎞ ⎤ y D = diámetro (m). Para el flujo turbulento, la ecuación de
V = r 2 cos –1 – (r – 2rh – h 2 L
h ⎢ ⎥
⎣ ⎝ r ⎠ ⎦ Colebrook proporciona un medio para calcular el factor de fric-
Determine h para r = 2 m, L = 5 m y V = 8.5 m3. Observe que ción,
si usted utiliza un lenguaje de programación o herramienta 1 ⎛ ε 2.51 ⎞
de 20
= − . log +
software, el arco coseno se puede calcular como f 3.7D Re f
π ⎛ x ⎞ donde e = rugosidad (m), y Re = número de Reynolds,
cos –1 x = – tan –1 2 ⎟
2 ⎝⎜ 1 – x ρVD
⎠ Re =
8.9 El volumen V del líquido contenido en un tanque esférico 8.11 La ecuación de Ergun, que se da abajo, sirve para
de radio r está relacionado con la profundidad h del líquido por describir el flujo de un líquido a través de un lecho empacado.
∆P es la
πk (2 3r – h)
V=
3
Determine h para r = 1 m y V = 0.75 m3.
8.10 Para el tanque esférico del problema 8.9, es posible desa-
rrollar las siguientes fórmulas para el método de punto fijo:
3
h + (3V /π
)
h=
3r
y
⎛ V⎞
h = 3 3⎝ rh 2 – ⎠
π
Si r = 1 m y V = 0.75 m3, determine si cualquiera de las dos al-

turas es estable, y el rango de valores iniciales para los que sí
son estables.
22 µ ESTUDIO
PROBLEM
donde m = viscosidad dinámica (N · s/m2).
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 × 10–5
N · s/m2, D = 0.005 m, V = 40 m/s, y e = 0.0015 mm.
Utilice un método numérico para determinar el factor de
fricción. Obsérvese que los tubos lisos tienen Re < 10 5, un
valor inicial apropiado se obtiene con el uso de la fórmula
de Blasius, f = 0.316/Re0.25.
b) Repita el cálculo pero para un tubo de acero comercial
más rugoso (e = 0.045 mm).
8.13 El pH del agua tiene gran importancia para los ingenieros

ambientales y químicos. Se relaciona con procesos que van de
la corrosión de tubos de lluvia ácida. El pH se relaciona con la
concentración del ion de hidrógeno por medio de la ecuación
siguiente:
pH = – log10 [H+]
22 ESTUDIO
PROBLEM
Las cinco ecuaciones que siguen gobiernan las concentraciones Ingeniería civil y ambiental
de una mezcla de dióxido de carbono y agua para un sistema 8.15 El desplazamiento de una estructura está definido por la
cerrado. ecuación siguiente para una oscilación amortiguada:
[H + ][HCO − ]
22 K1 = 3
ESTUDIO
PROBLEM
DE CASOS: RAÍCES DE y = 9e cos wt 22
[CO 2 ]
+ 2− donde k = 0.7 y w = 4.
K 2 = [H ][CO– 3]
[HCO ]3 a) Utilice el método gráfico para realizar una estimación ini-
Kw = [H +][OH ] cial del tiempo que se requiere para que el desplazamiento
disminuya a 3.5.
c = [CO ] + [HCO – ] + [CO 2−
]T 2 3 3
b) Emplee el método de Newton-Raphson para determinar la
Alk = [HCO–3 ] + 2[CO23− ] + [OH– ] − [H+ ] raíz con es = 0.01%.
c) Use el método de la secante para determinar la raíz con es
donde Alk = alcalinidad, cT = total de carbón inorgánico, y las
= 0.01%.
K
son coeficientes de equilibrio. Las cinco incógnitas son [CO2] =
dióxido de carbono, [HCO3 ] = bicarbonato, [CO3 ] = 8.16 En ingeniería estructural, la fórmula de la secante define
carbonato, la fuerza por unidad de área, P/A, que ocasiona la tensión
[H+] = ion hidrógeno, y [OH–] = ion hidroxilo. Resuelva para
las máxima
cinco incógnitas dado que Alk = 2 × 10–3, cT = 3 × 10–3, K1 = sm en una columna que tiene una razón de esbeltez L/k dada
10–6.3, y K2 = 10–10.3, y Kw = 10–14. Asimismo, calcule el pH es:
de
las soluciones. P σm
=
8.14 La ecuación que se presenta a continuación, describe la
2
operación de un reactor de flujo por inyección de densidad A 1+ (ec/k ) sec [0.5 P/(EA)(L/k)]
constante para la producción de una sustancia por medio de
una re-
acción enzimática, donde V es el volumen del reactor, F es la donde ec/k2 = razón de excentricidad, y E = módulo de elastici-
tasa de flujo del reactivo C, Cent y Csal son las concentraciones
dad. Si para una viga de acero, E = 200 000 MPa, ec/k2 = 0.4 y
del reactivo que entra y sale del reactor, respectivamente, y K y sm = 250 MPa, calcule P/A para L/k = 50. Recuerde que sec x =
kmáx son constantes. Para un reactor de 500 L, con una
1/cos x.
concentración en la toma de Cent = 0.5 M, tasa de entrada de 8.17 Un cable en forma catenaria es aquel que cuelga entre dos
flujo de 40 L/s, puntos que no se encuentran sobre la misma línea vertical.
kmáx = 5 × 10–3s–1, y K = 0.1 M, encuentre la concentración de
Como se ilustra en la figura P8.17a, no está sujeta a más carga
C
a la salida del reactor. que su propio peso. Así, su peso (N/m) actúa como una carga
uniforme
Csal
K 1 por unidad de longitud a lo largo del cable. En la figura P8.17b,
F
V= –
∫
Cent kmáx C
+
kmáx
dC
se ilustra un diagrama de cuerpo libre de una sección AB, donde
Figura P8.17
a) Fuerzas que actúan sobre y
una sección AB de un cable
flexible que cuelga. La
carga es uniforme a lo largo TB
B
del cable (pero no uniforme
por la distancia horizontal
x). b) Diagrama de cuerpo
libre de la sección AB.
A
W = ws
w TA
y0
x
a) b)
22 ESTUDIO
PROBLEM
TA y TB son las fuerzas de tensión en el extremo. Con base en Utilice el método de la bisección para determinar el punto de
los balances de fuerzas horizontal y vertical, se obtiene para el máxima deflexión (es decir, el valor de x donde dy/dx = 0).
cable el siguiente modelo de ecuación diferencial: Des- pués, sustituya este valor en la ecuación (P8.18) a fin de
determinar el valor de la deflexión máxima. En sus cálculos,
utilice
2
d2 y w ⎛ dy ⎞ los valores siguientes para los parámetros: L = 600 cm, E = 50
000
= 1+
dx 2 TA ⎝ dx kN/cm 2, I = 30 000 cm 4 y w0 = 2.5 kN/cm.
⎠ 8.19
En la ingeniería ambiental (una especialidad de la ingenie-
Puede emplearse el cálculo para resolver esta ecuación para la ría civil), la ecuación siguiente se emplea para calcular el nivel
altura y del cable como función de la distancia x. de oxígeno c (mg/L) en un río aguas abajo de la descarga de un
drenaje:
TA ⎛ w ⎞ T
y= cosh x +y − A
w ⎝ ⎜TA ⎠⎟ 0 c = 10 – 20(e–0.15x – e–0.5x)
w
donde el coseno hiperbólico se calcula por medio de la ecua- donde x es la distancia aguas abajo en kilómetros.
ción:
a) Determine la distancia aguas abajo de la corriente, a la
1 x –x cual el nivel de oxígeno cae hasta una lectura de 5 mg/L.
cosh x = (e + e )
2 (Reco- mendación: está dentro de 2 km de la descarga.)
Utilice un método para calcular un valor para el parámetro TA Encuentre la respuesta con un error de 1%. Obsérvese que
dados los valores de los parámetros w = 12 y y0 = 6, de modo los niveles de oxígeno por debajo de 5 mg/L por lo
que el cable tenga una altura de y = 15 en x = 50. general son dañinos para ciertas especies de pesca
8.18 En la figura P8.18a se muestra una viga uniforme sujeta a deportiva, como la trucha y el salmón.
una carga distribuida uniformemente que crece en forma lineal. b) Calcule la distancia aguas abajo a la cual el oxígeno se
La ecuación para la curva elástica resultante es la siguiente encuentra al mínimo. ¿Cuál es la concentración en dicha
(véase la figura P8.18b) ubicación?
w0 5 2 3 4 8.20 La concentración de bacterias contaminantes c en un lago
y= (− x + 2L x − L x)
disminuye de acuerdo con la ecuación
(P8.18)
120EIL
c = 75e–1.5t + 20e–0.075t
Determine el tiempo que se requiere para que la concentración

Figura P8.18 de bacterias se reduzca a 15 con el uso de a) el método gráfico,
y b) el método de Newton-Raphson, con un valor inicial de t =
6 y criterio de detención de 0.5%. Compruebe los resultados
que
w0 obtenga.
8.21 En ingeniería oceanográfica, la ecuación de una ola esta-
e por l = 16, t = 12,
cionaria reflejada en un puerto está dada
v = 48:
⎡ ⎛ 2πx ⎞ cos ⎛ 2πtv + −x ⎤

L h = ⎢sen ⎝ λ ⎠ ⎞⎝
λ ⎠ ⎥
0 ⎦
⎣
a)
Resuelva para el valor positivo más bajo de x, si h = 0.5 h0.
8.22 Suponga el lector que compra una pieza de equipo en
(x = L, y = 0) $25 000 como pago inicial y $5 500 por año durante 6 años.
(x = 0, y = 0) ¿Qué tasa de interés estaría pagando? La fórmula que relaciona
x el valor presente P, los pagos anuales A, el número de años n y
la tasa de interés i, es la que sigue:
22 ESTUDIO
PROBLEM
i(1 + i )n
b) A = P (1+ i)n −1
22 ESTUDIO
PROBLEM
20 kips/ft
150 kips-ft
15 kips
5’ 2’ 1’ 2’
Figura P8.24
8.23 Muchos campos de la ingeniería requieren estimaciones M(x) = –10[〈x – 0〉2 – 〈x – 5〉2] + 15〈x – 8〉1
exactas de la población. Por ejemplo, los ingenieros de
+ 150〈x – 7〉0 + 57x
transporte quizás encuentren necesario determinar por
separado la ten- dencia del crecimiento de una ciudad y la de
Emplee un método numérico para encontrar el (los) punto(s) en
los suburbios. La población del área urbana declina con el
los que el momento es igual a cero.
tiempo de acuerdo con la ecuación:
8.26 Con el uso de la viga con apoyo simple del problema
Pu(t) = Pu,máxe–kut + Pu,mín 8.24, la pendiente a lo largo de ella está dada por:
duy −10 3 3 15 2
en tanto que la población suburbana crece según: (x) = [〈x − 0〉 − 〈x − 5〉 ] + 〈x − 8〉
dx 3 2
Ps , máx 1 57 2
p (t) = −k t + 150 〈 x − 7〉 + x − 238.25
s 1 + [Ps, / P −1]e s 2
0
máx
donde Pu,máx, ku, Ps,máx, P0 y ks son parámetros que se obtienen Utilice un método numérico para encontrar el(los) punto(s)
en forma empírica. Determine el tiempo y los valores donde la pendiente es igual a cero.
correspondien- tes de Pu(t) y Ps(t) cuando los suburbios son 8.27 Para la viga con apoyo simple del problema 8.24, el des-
20% más grandes que la ciudad. Los valores de los parámetros plazamiento a lo largo de ella está dado por la ecuación:
son: Pu,máx = 75 000,
Ku = 0.045/año, Pu,mín = 100 000 personas, Ps,máx = 300 000 −5 15
[〈x − 0〉 4 − 〈x − 5〉 4 ] + 3
personas, P0 = 10 000 personas, ks = 0.08/año. Para obtener u y (x) = 〈x − 8〉
6 6
las soluciones utilice los métodos a) gráfico, b) de la falsa
posición, y
c) de la secante modificada. 2 57 3
+ 75〈x − 7〉 + x − 238.25x
8.24 En la figura P8.24 se muestra una viga apoyada en forma 6
sencilla que está cargada como se ilustra. Con el empleo de a) Calcule el (los) punto(s) donde el desplazamiento es igual
funciones de singularidad, el esfuerzo cortante a lo largo de la a cero.
viga se expresa con la ecuación: b) ¿Cómo se usaría una técnica de localización de raíces para
determinar la ubicación del desplazamiento mínimo?
V(x) = 20[〈x – 0〉1 – 〈x – 5〉1] – 15〈x – 8〉0 – 57
Por definición, la función de singularidad se expresa del modo Ingeniería eléctrica

que sigue: 8.28 Ejecute el mismo cálculo que en la sección 8.3, pero
determine el valor de C que se requiere para que el circuito
disipe 1%
⎧(x − a) n cuando x > a ⎫ de su valor original en t = 0.05 s, dado R = 280 Ω, y L = 7.5 H.
n
〈 x – a〉 = ⎨ 0 ⎬ Emplee a) un enfoque gráfico, b) la bisección, y c) software
⎩ cuando x ≤ a
⎭ para encontrar raíces, tales como Solver de Excel o la función
fzero
Utilice un método numérico para encontrar el(los) punto(s) en 8.25 Con el uso de la viga apoyada en forma simple del
los que el esfuerzo cortante sea igual a cero. problema 8.24, el momento a lo largo de ella, M (x) está dada
por:
de
22 ESTUDIO
PROBLEM
MATLAB.
8.29 La ecuación i = 9e–t cos (2pt), describe una corriente
osci- latoria en un circuito eléctrico, donde t se expresa en
segundos. Determine todos los valores de t de modo que i = 3.
22 ESTUDIO
PROBLEM
8.32 En la figura P8.32 se muestra un circuito con una resisten-

a cia, un inductor y un capacitor en paralelo. Para expresar la
impedancia del sistema se emplean las leyes de Kirchhoff, así:
x
q
1 1 ⎛ 1 ⎞2
= 2 + ωC –
Z R ⎝ ωL ⎠
Q
donde Z = impedancia (Ω) y w = frecuencia angular. Encuentre
Figura P8.31
la w que da como resultado una impedancia de 75 Ω, con el uso
tanto del método de la bisección como el de la falsa posición,
con valores iniciales de 1 y 1000 y los parámetros siguientes: R
8.30 La resistividad r de un lubricante de sílice se basa en la = 225 Ω, C = 0.6 × 10–6 F, y L = 0.5 H. Determine cuántas
carga q en un electrón, la densidad del electrón n, y la iteraciones son necesarias con cada técnica a fin de encontrar
movilidad del electrón m. La densidad del electrón está dada en la respuesta con es = 0.1%. Utilice el enfoque gráfico para
términos de la densidad del lubricante N, y la densidad explicar cualesquiera dificultades que surjan.
intrínseca de acarreo ni. La movilidad del electrón está descrita
por la temperatura T, la temperatura de referencia T0, y la Ingeniería mecánica y aeroespacial
movilidad de referencia µ0. Las ecuaciones que se requieren 8.33 Para la circulación de fluidos en tubos, se describe a la
para calcular la resistividad son las siguientes: fricción por medio de un número adimensional, que es el factor
de fricción de Fanning f. El factor de fricción de Fanning
1 depende de cierto número de parámetros relacionados con el
ρ=
qnµ tamaño del tubo y el fluido, que pueden representarse con otra
donde cantidad adimensional, el número de Reynolds Re. Una fórmula
que pro- nostica el valor de f dado Re es la ecuación de von
−2.42
Karman.
1 ⎛ T ⎞
n=
2
(N + 2
N + 4ni
2
) y µ = µ0
T
1
( )
ƒ= 4 log10 Re ƒ − 0.4
0
Determine N, dado que T0 = 300 K, T = 1 000 K, µ0 = 1 350

Valores comunes del número de Reynolds para flujo turbulento
cm2 (V s)–1, q = 1.7 × 10–19 C, ni = 6.21 × 109 cm–3, y un valor
son 10 000 a 500 000, y del factor de fricción de Fanning
desea- ble de r = 6.5 × 106 V s cm/C. Use los métodos a)
son
bisección, y b) la secante modificada.
0.001 a 0.01. Desarrolle una función que utilice el método de
8.31 Una carga total Q se encuentra distribuida en forma uni-
bisección con objeto de resolver cuál sería el factor de fricción
forme alrededor de un conductor en forma de anillo con radio
de Fanning f, dado un valor de Re proporcionado por el usuario
a. Una carga q se localiza a una distancia x del centro del anillo
que esté entre 2 500 y 1 000 000. Diseñe la función de modo
(véase la figura P8.31). La fuerza que el anillo ejerce sobre la
que se garantice que el error absoluto en el resultado sea de
carga está dada por la ecuación
Ea,d <
1 qQx 0.000005.
F=
8.34 Los sistemas mecánicos reales involucran la deflexión de
resortes no lineales. En la figura P8.34 se ilustra una masa m
que
2
4πe0 (x + a )
2 3/ 2
se libera por una distancia h sobre un resorte no lineal. La fuerza
de resistencia F del resorte está dada por la ecuación
donde e0 = 8.85 × 10–12 C2/(N m2). Encuentre la distancia
x donde la fuerza es de 1.25 N, si q y Q son 2 × 10–5 C para un
anillo con un radio de 0.9 m.
Figura P8.34
Figura P8.32
h
h+d
d
R L C
22 ESTUDIO
PROBLEM
DE CASOS: RAÍCES DE a) b) 22
22 ESTUDIO
PROBLEM
F = –(k1d + k2d3/2) Para g, utilice un valor de 9.81 m/s2, y emplee el método

gráfico para elegir valores iniciales.
Es posible usar la conservación de la energía para demostrar 8.37 La velocidad vertical de un cohete se calcula con la
que fórmu- la que sigue:
2 k 2d
5 /2
1 m0
0= +
2
k1 d − mgd – mgh v = u ln − gt
5 2 m0 − qt
Resuelva cuál sería el valor de d, dados los valores siguientes donde v = velocidad vertical, u = velocidad con la que se
de los parámetros: k1 = 50 000 g/s2, k2 = 40 g/(s2 m0.5), m = expele el combustible, en relación con el cohete, m0 = masa
90 g, g = 9.81 m/s2, y h = 0.45 m. inicial del cohete en el momento t = 0, q = tasa de consumo de
combustible, y g = aceleración de la gravedad hacia abajo (se
8.35 Los ingenieros mecánicos, así como los de otras
especialidades, utilizan mucho la termodinámica para realizar supone constante e igual a 9.81 m/s2). Si u = 2000 m/s, m0 =
su trabajo. El siguiente polinomio se emplea para relacionar el 150 000 kg, y q =
calor espe- cífico a presión cero del aire seco, cp kJ/(kg K), a 2 700 kg/s, calcule el momento en que v = a 750 m/s. (Sugeren-
cia: El valor de t se encuentra entre 10 y 50 s.) Calcule el resul-
temperatura (K):
tado de modo que esté dentro de 1% del valor verdadero.
cp = 0.99403 + 1.671 × 10–4 T + 9.7215 × 10–8T 2 Compruebe su respuesta.
–9.5838 × 10–11T 3 + 1.9520 × 10–14T4 8.38 En la sección 8.4, el ángulo de fase f entre la vibración
forzada que ocasiona el camino rugoso y el movimiento del
Determine la temperatura que corresponda a un calor carro, está dada por la ecuación:
específico 2(c/cc )(ω / p)
de 1.1 kJ/(kg K). tan φ = 2
1 – (ω/ p)
8.36 En ciertas ocasiones, los ingenieros aerospaciales deben
calcular las trayectorias de proyectiles, como cohetes. Un pro- Como ingeniero mecánico, le gustaría saber si existen casos en
blema parecido tiene que ver con la trayectoria de una pelota que f = w/3 – 1. Utilice los otros parámetros de la sección
que se lanza. Dicha trayectoria está definida por las con objeto de plantear la ecuación como un problema de
coordenadas (x, y), como se ilustra en la figura P8.36. La cálculo de raíces, y resuélvala para w.
trayectoria se modela con la ecuación 8.39 Se mezclan dos fluidos con temperatura diferente de
g modo que alcanzan la misma temperatura. La capacidad
2
y = (tan θ 0 )x − 2 2 x + 1.8 calorífica del fluido A está dada por:
2v cos θ
0 0
cp = 3.381 + 1.804 × 10–2T – 4.300 × 10–6 T 2
Calcule el ángulo inicial q0, apropiado si la velocidad inicial
v0 = 20 m/s y la distancia x al catcher es de 35 m. Obsérvese y la capacidad calorífica del fluido B se obtiene con:
que la pelota sale de la mano del lanzador con una elevación y0
= 2 m, y el catcher la recibe a 1 m. Exprese el resultado final cp = 8.592 + 1.290 × 10–1T – 4.078 × 10–5 T 2
en grados.
donde cp se expresa en unidades de cal/mol K, y T está en uni-
dades de K. Obsérvese que
T2
Figura P8.36 ∆H = ∫ c p dT
T
1
y El fluido A entra al mezclador a 400ºC, y el B a 700ºC. 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 razón de compresión Rc de 3.0
(esto significa que la presión del gas en la salida es tres veces
mayor que en la entrada). Los requerimientos de energía del
compresor Hp se determinan por medio de la ecuación que se da
v 0
a continuación. Suponga que los requerimientos de energía del
0 compresor son exactamente iguales a zRT1/MW, y encuentre la
x eficiencia politrópica n del compresor. El parámetro z es la
com- presibilidad del gas en las condiciones de operación del compre-
23 ESTUDIO
PROBLEM
⎡σ xx σ xy
σ xz ⎤
⎢ ⎥
⎢ σ xy σ σ yz
σ zz⎥
yy
⎢σ
⎣ xz σ yz
⎥
⎦
en la que los términos en la diagonal principal representan es-
fuerzos a la tensión o a la compresión, y los términos fuera de
la diagonal representan los esfuerzos cortantes. Un campo
tensorial (en MPa) está dado por la matriz que sigue:
⎡10 14 25⎤
⎢ ⎥
⎢14 7 15⎥
T3
T0 ⎣⎢25 15 16 ⎦⎥
T2
T1 Para resolver cuáles son los esfuerzos principales, es necesario
construir la matriz siguiente (de nuevo en MPa):
⎡10 − σ 14 25 ⎤
Figura P8.41 ⎢ ⎥
14 7−σ 15 ⎥
⎢
25 15 16 − σ
s1, s2 y s3 se obtienen con la ecuación
sor, R es la constante de los gases, T1 es la temperatura del gas 3 2
σ − Iσ + IIσ − III = 0
en la entrada del compresor, y MW es el peso molecular del
donde
gas.
I = σ xx + σ yy + σ zz
zRT1 n
HP =
MW n − 1
( Rc(n−1)/n − 1) II = σ xxσ yy + σ xx σ zz + σ yy σ zz − σ xy2 − σxz2 − σyz 2
2 2 2
8.41 En los envases térmicos que se ilustran en la figura P8.41, III = σxx σyy σzz − σxx σyz − σyy σxz − σ
zz σ
xy + 2σxy σxz yz
el compartimiento interior está separado del medio por medio

de vacío. Hay una cubierta exterior alrededor de los envases. I, II y III se conocen como las invariantes de esfuerzos.
Esta cubierta está separada de la capa media por una capa Encuentre
delgada de aire. La superficie de afuera de la cubierta exterior s1, s2 y s3 por medio de una técnica de localización de raíces.
está en con- tacto con el aire del ambiente. La transferencia de 8.43 La figura P8.43 ilustra tres almacenamientos conectados
calor del compartimiento interior a la capa siguiente q1 sólo por medio de tubos circulares. Los tubos están hechos de hierro
ocurre por radiación (ya que el espacio se encuentra vacío). La
transferencia de calor entre la capa media y la cubierta exterior
q2 es por con- vección en un espacio pequeño. La transferencia
de calor de la cubierta exterior hacia el aire q3 sucede por Figura P8.43
convección natural.
h1
igual, es decir, q1 = q2 = q3. Encuentre las temperaturas T1 y
T2
en estado estable. T0 es de 450ºC y T3 = 25ºC. A h2
1
q1 = 10 [(T0 + 273) − (T1 + 273) B
2
Q1 h3
]
q2 = 4(T1 − T2 ) Q2
4/3
q3 = 1.3(T2 − T3 3
) Q3
8.42 La forma general para un campo tensorial de tres dimen- C
siones es la siguiente:
23 ESTUDIO
PROBLEM
fundido recubierto con asfalto (e = 0.0012 m), y tienen las ca- Q1 = 1 m3/s y r = 1.23 kg/m3. Todos los tubos tienen D =
racterísticas siguientes: 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.
Tubo 1 2 3 8.45 Repita el problema 8.44, pero incorpore el hecho de que el
Longitud, m 1800 500 1400 factor de fricción se calcula con la ecuación de von Karman, que
Diámetro, m 0.4 0.25 0.2 es:
Flujo, m3/s ? 0.1 ?
1
= 4 log
10 (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, donde Re = número de Reynolds
respectivamente, determine la elevación que alcanza en el
ρVD
almacenamiento B y los flujos en los tubos 1 y 3. Obsérvese Re =
que la viscosidad cinemática del agua es de 1 × 10–6 m2/s, y µ
utilice la ecuación de Colebrook para obtener el factor de
donde V = velocidad del fluido en el tubo [m/s], y µ =
fricción (consulte el problema 8.12).
viscosidad dinámica (N ⋅ s/m2). Obsérvese que para un tubo
8.44 Un fluido se bombea en la red de tubos que se muestra en
circular, V = 4Q/ pD2. Asimismo, suponga que el fluido tiene
la figura P8.44. En estado estacionario, se cumplen los
una viscosidad de
balances de flujo siguientes:
1.79 × 10–5 N ⋅
Q1 = Q2 + Q3 s/m2.
8.46 Sobre el trasbordador espacial, al despegar de la
Q3 = Q4 + Q5 plataforma, actúan cuatro fuerzas, las que se muestran en el
Q5 = Q6 + Q7 diagrama de cuerpo libre (véase la figura P8.46). El peso
combinado de los dos cohetes de combustible sólido y del
tanque exterior de este,
es de WB = 1.663 × 106 lb. El peso del orbitador con carga com-
donde Qi = flujo en el tubo i [m3/s]. Además, la caída de
presión pleta es de WS = 0.23 × 106 lb. El empuje combinado de los dos
alrededor de los tres lazos en los que el flujo es hacia la derecha
debe ser igual a cero. La caída de presión en cada tramo de cohetes de combustible sólido es TB = 5.30 × 10 lb. El empuje
tubo
circular se calcula por medio de la ecuación: combinado de los tres motores de combustible líquido del orbi-
16 fL ρ 2 tador es de TS = 1.125 × 106 lb.
∆P = 2 Q Al despegar, el empuje del motor del orbitador se dirige con
π 2D 5 un ángulo q para hacer que el momento resultante que actúa
donde ∆P= caída de presión [Pa], f = factor de fricción sobre el conjunto de la nave (tanque exterior, cohetes de
[adimensional], L = longitud del tubo [m], r = densidad del combustible sólido y orbitador) sea igual a cero. Con el
fluido [kg/m3], y D = diámetro del tubo [m]. Escriba un momento resultante igual a cero, la nave no giraría sobre su
programa (o desarrolle un algoritmo en algún paquete de centro de gravedad G al despegar. Con estas fuerzas, la nave
software de matemáticas) que permita calcular el flujo en experimentará una fuerza resultante con componentes en
cada tramo de tubo, dado que dirección 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.
Figura P8.44 El momento resultante que actúa sobre la nave será igual a cero
cuando q se ajusta al valor apropiado. Si este ángulo no se
Q1 Q3 Q5 ajusta en forma adecuada y hubiera algún momento que actuara
sobre la nave, ésta tendería a girar alrededor de su centro de
gravedad.
Q10 Q9 Q8
Q2 Q4 Q6
23 a) Resuelva el empujeESTUDIO
delPROBLEM
orbitador TS en las componentes
horizontal y vertical, y después
sume los momentos respecto
Q7
del punto G, centro de gravedad
de la nave. Iguale a cero la
ecuación del momento
resultante. Ahora, ésta puede
resolverse para el valor de q que
se requiere durante el despegue.
b) Obtenga una ecuación para el
momento resultante que actúa
sobre la nave en términos del
ángulo q. Grafique el
23 ESTUDIO
PROBLEM
momento resultante como función del ángulo q en el

rango de –5 radianes a +5 radianes.
28’
c) Escriba un programa de computadora para resolver para el
’ ángulo q por medio del método de Newton para encontrar
la raíz de la ecuación del momento resultante. Con el
empleo de la gráfica, elija un valor inicial para la raíz de
Tanque interés. Interrumpa las iteraciones cuando el valor de q ya
externo no mejore con cinco cifras significativas.
d) Repita el programa para el peso de la carga mínima del
orbitador, que es WS = 195 000 lb.
Cohete de
combustibl
e sólido
Orbitado
r
38’
WB
S
W
TS
TB
Figura P8.46

Estudio de Casos

Cargado por

Copyright:

Formatos disponibles

Estudio de Casos

Cargado por

Información del documento

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

Estudio de Casos

Cargado por

Copyright:

Formatos disponibles

CAPÍTULO 8

8.1 LEYES DE LOS GASES IDEALES Y NO IDEALES

Antecedentes. La ley de los gases ideales está dada por

ampliamente por los ingenieros y científicos, sólo es exacta en un rango limitado de

TABLA 8.1 Cálculos del volumen molar.

Volumen Volumen molar Volumen molar

300 1 24.6162 24.5126 24.5928

En este caso, como la derivada de ƒ(v) se determina fácilmente, entonces es

que proporcionan lecturas instantáneas de presión y temperatura. Se debe obtener

8.2 FLUJO EN UN CANAL ABIERTO

Antecedentes. La ingeniería civil constituye una disciplina amplia que incluye

Solución. La relación fundamental entre flujo y profundidad es la ecuación de

donde B = ancho (m). Debe observarse que la ecuación de continuidad se obtiene de la

donde n = coeficiente de rugosidad de Manning (un número adimensional que toma en

donde P = perímetro mojado (m). Como su nombre lo indica, el perímetro mojado es

Se debe observar que así como la ecuación de continuidad se obtiene de la

La ecuación (8.12) se resuelve rápidamente con cualquiera de los métodos para

Puede resolverse para H = 0.7023 m. El resultado se verifica sustituyéndolo en la

que se acerca bastante a cero.

Así, se requieren 14 iteraciones. Aunque esto ciertamente no sería costoso para un

como para H en el denominador,

Ahora, aquí es donde el razonamiento físico puede ayudar. En la mayoría de los

FIGURA 8.2 y2 = g(H)

Hay dos beneficios prácticos de este análisis:

1. En el caso de que se use un método abierto más detallado, la ecuación (8.16)

8.3 DISEÑO DE UN CIRCUITO ELÉCTRICO

Antecedentes. Los ingenieros eléctricos emplean las leyes de Kirchhoff para

riodo de ajuste está íntimamente relacionada con las propiedades de almacenamiento

donde i = la corriente y R = la resistencia del resistor. Si las unidades de R e i son

donde L = la inductancia. Si las unidades de L e i son henrios y amperes,

donde C = la capacitancia. Si las unidades de carga se expresan en coulombios,

Sin embargo, como la corriente se relaciona con la carga de acuerdo con

Solución. Es necesario despejar R de la ecuación (8.19) con valores conocidos

Utilizando los valores numéricos dados,

cos [ 2 000 – 0.01R 2 (0.05)] – 0.01

8.4 ANÁLISIS DE VIBRACIONES (INGENIERÍA MECÁNICA

Antecedentes. Las ecuaciones diferenciales sirven para modelar la vibración de

Péndulo Resorte/masa Circuito LC

Fuerza del resorte = –kx

donde c es el coeficiente de amortiguamiento y dx/dt es la velocidad vertical. El signo

La incógnita r es la solución de la ecuación característica cuadrática que se puede

Note el significado de la magnitud de c al compararla con 2 km . Si c > 2 km ,

donde A y B son constantes que se deben determinar a partir de las condiciones

Tales sistemas se conocen como subamortiguados.

La relación c/cc se llama factor de amortiguamiento, y a p se le conoce como la

donde dm = Pm /k = la deflexión estática del carro sujeto a una fuerza Pm. La

La solución general de esta ecuación se obtiene al sumar una solución particular a

remos el movimiento en estado estacionario del sistema forzado donde se ha

xss (t) = xm sen (wt – f)

La cantidad xm /dm llamada factor de amplificación de la amplitud depende tan sólo de

natural p para diversos 0.5

Nuestras condiciones de diseño se satisfacen si

Dado que se conocen c y m, el problema de diseño consiste ahora en encontrar valores

Solución. Se pueden utilizar los métodos de la bisección, de la falsa posición o de la

donde d es la deflexión, dm es la máxima deflexión de 0.1 m y D es la distancia entre

donde w = 2pv/D es la frecuencia forzada.

Ahora, se buscan valores w/p que satisfagan la ecuación (8.28),

Si la ecuación (8.30) se expresa como un problema de raíces

y usar MATLAB para determinar las raíces como sigue: