Problemas Resueltos PDF
Problemas Resueltos PDF
Problemas Resueltos PDF
CAÍDA DE PARACAIDISTAS
COEFICIENTE DE SEGUNDO ORDEN
(Chapra Canale, Cap.1 – Problema 1.5 - Pág. 23)
Aplicando la segunda Ley de Newton, determinar la velocidad final ( ) de caída libre de un cuerpo, partiendo de un estado de
reposo y cercano a la tierra, si se considera que la fuerza hacia arriba generada por la resistencia del aire es no lineal y se expresa
como:
Solución:
Utilizando la ecuación dada tenemos:
Por tanto:
Graficamos las y para comparar resultados, los cuales se muestran en las siguientes tabla y gráfica, donde puede
verse a los 30 seg. se alcanza la velocidad límite:
ANALÍTICA NUMÉRICA
t Ec.(a) Ec.(b)
seg m/s m/s
0 0.00 0.00
2 18.83 19.61
4 33.76 36.76
6 43.55 47.72
8 49.21 52.76
10 52.24 54.55
12 53.80 55.11
14 54.58 55.28
16 54.97 55.33
18 55.16 55.34
20 55.26 55.35
22 55.30 55.35
24 55.33 55.35
26 55.34 55.35
28 55.34 55.35
30 55.35 55.35
32 55.35 55.35
34 55.35 55.35
36 55.35 55.35
38 55.35 55.35
40 55.35 55.35
Determinar el valor del coeficiente de arrastre (c) para que un paracaidista de masa tenga una velocidad
después de una caída libre de Tomar la aceleración de la gravedad como
Solución:
Siendo la incógnita el valor de , debemos plantear una expresión que permita hallar este parámetro implícito:
Por lo tanto el valor de que hace que es la ráiz de la ecuación, es decir el valor que resuelve el problema.
Solución con EXCEL:
El valor de ceros se puede ver que es aproximadamente c=15 mt. Para calcular dicho valor aplicamos en la tabla la función BUSCAR
OBJETIVO obteniendo c=14.7940903914612.
Solución con EULER:
Se define la f(c):
>fc:= "(m*g/x)*(1-%e^-(x*t/m))-v";
S
Q, U
Sección
B
Transversal (Ac)
Si Q= 5 m3/seg; n=0.03; y s=0.02%, y B=Ancho del canal = 20.0 mt, se pide determinar H y U por cualquier método numérico
conocido y comprobar con dos funciones en EULER
Solución:
De la Ley de Conservación de la masa se deduce la Ecuación de Continuidad la cual establece que:
Donde es el área de la sección transversal del canal (m2) y representa la velocidad del flujo (m/s). Para el canal rectangular
mostrado tenemos:
Por otra parte, para un flujo uniforme (es decir aquel en que el flujo no varía ni con la distancia ni con el tiempo) Robert Manning
propuso una expresión semi-empírica para la velocidad:
En la cual:
Conocidos y la Ec.(d) no puede resolverse implícitamente para hallar , por lo cual debe resolverse en la forma:
Solución con EXCEL:
H f(H) - Q
0.00 -5.00
0.10 -4.7982
0.20 -4.3636
0.30 -3.7572
0.40 -3.0055
0.50 -2.1254
0.60 -1.1292
0.70 -0.0265
0.80 1.1748
0.90 2.4681
1.00 3.8477
El valor de ceros se puede ver que es aproximadamente 0.70 mt. Para calcular dicho valor aplicamos en la tabla la función BUSCAR
OBJETIVO obteniendo la altura de la lámina de agua, H=0.7023 mt.
Reemplazando dicho valor en la Ec.(a) obtenemos la velocidad del agua, U=0.3560 m/s.
Solución con EULER:
>fH:= "(s^(1/2)/n)*((B*x)^(5/3)/(B+2*x)^(2/3))-Q";
Se guarda en H (altura del canal) la solucion de la función:
>H:=solve(fH,1); print(H,unit=" mt.")
0.70 mt.
Se calcula el radio hidraulico (R) y la velocidad (U) y se grafica:
>R:=(B*H)/(B+2*H); print(R,unit=" mt.")
0.66 mt.
>U:= R^(2/3)*s^(1/2)/n; print(U,unit=" m/seg.")
0.36 m/seg.
>plot2d(fH,title="f(H) vs H", a=0,b=1,c=-5,d=5, color=red, grid=2; ...
xl="Nivel agua (H) mt.");
>plot2d(label("corte eje = H",H,0), add=1);
CANAL TRIANGULAR
Dadas la pendiente longitudinal (s) y el coeficiente de fricción (n) de un canal triangular y el caudal Q en un punto dado, determinar
la velocidad (U) y la profundidad (H) si el régimen es de flujo libre abierto, partiendo de la ecuación de continuidad derivada de la
Ley de Conservación de Masa.
Si Q= 7.5 m3/seg; n=0.03; y s=0.02%, , se pide determinar H y U por cualquier método numérico conocido y comprobar con dos
funciones en EULER
d H
s (pend. Long)
α α
Q,U
Área Transversal x
Solución:
Donde es el área de la sección transversal del canal (m2) y representa la velocidad del flujo (m/s). Para el canal triangular
mostrado tenemos:
Por otra parte, para un flujo uniforme (es decir aquel en que el flujo no varía ni con la distancia ni con el tiempo) Robert Manning
propuso una expresión semi-empírica para la velocidad:
En la cual:
Combinando las Ec.(a) y (b) obtenemos:
Conocidos y la Ec.(d) no puede resolverse implícitamente para hallar , por lo cual debe resolverse en la forma:
H f(H) - Q
0.00 -5.00
0.10 -4.7982
0.20 -4.3636
0.30 -3.7572
0.40 -3.0055
0.50 -2.1254
0.60 -1.1292
0.70 -0.0265
0.80 1.1748
0.90 2.4681
1.00 3.8477
CANAL TRAPEZOIDAL
En el canal trapezoidal con una pendiente longitudinal que se muestra en la figura, circula libremente un caudal de agua
m3/seg que a una profundidad (m) alcanza una velocidad (m/seg).
Dados además:
B = 20.0 m ; n = 0.03 ; s = 0.0002; y α = 45°
Se pide:
1. Determinar la altura del agua H (m) en la sección transversal y su velocidad U (m/s) correspondiente al caudal dado
2. Determinar la variación (en porcentaje) de H y U si el canal es rectangular (B x H).
s = pendiente
Q, U
Ac = Sección Transversal
d H
α α
B x
Solución:
Se parte de la ecuación de continuidad que se origina en la Ley de Conservación de Masa:
Siendo en ancho del canal (m) en la base del trapecio y el ángulo de la pared del canal con la horizontal.
Para calcular la velocidad del agua en flujo uniforme (es decir que no varía con la distancia ni con el tiempo, siendo por tanto
constante) Robert Manning propuso la siguiente fórmula empírica (conocida como ecuación de Manning):
Donde:
= Coeficiente de Manning de rugosidad del canal (tiene en cuenta la fricción del canal).
= Pendiente longitudinal del canal.
= Radio hidráulico (m):
Dado que la incógnita H está implícita la ec. (c) podemos reescribirla como función de H así:
Que es finalmente la ecuación que debemos utilizar para calcular la altura H que corresponde a un caudal Q dado.
Para el caso de:
Q = 5.0 m3/s
B = 20.0 m
n = 0.03
s = 0.0002
α = 45°
Utilizamos EULER para calcular la altura H (m), luego usamos este valor para calcular el radio hidráulico R (m) y la correspondiente
velocidad (m/s):
PROBLEMAS RESUELTOS
DETERMINAR RAÍCES REALES (1)
Chapra-Canale, Cap. 5 – Problema 5.1 - Pág.137 -
1. Analíticamente
2. Gráficamente
3. Con Excel (buscar objetivo)
4. Método de la bisección con 11 iteraciones (tome xl=-2, xu=0 para la primera raíz y xl=5, xu=10 para la segunda raíz).
Comparar con EULER (MAXIMA)
5. Método de Newton-Raphson, tomando -2 como valor inicial para la primera raíz y 7 para la segunda raíz. Comparar con
EULER (MAXIMA)
6. Método de Secante, comparar con EULER (MAXIMA)
1. Solución Analítica:
Aplicamos la cuadrática y obtenemos:
2. Solución Gráfica:
Utilizamos la utilidad en línea function evaluator and grapher v3.4 (www.zweigmedia.com) obtenemos una diferencia de
0.0021056 con el cero:
Vr Inicial= -2 Vr Inicial= 7
paso= 0.1 paso= 0.01
x f(x) x f(x)
-2.000 -1.3000 7.000 0.5000
-1.900 -0.9240 7.010 0.4660
-1.800 -0.5560 7.020 0.4318
-1.700 -0.1960 7.030 0.3976
-1.600 0.1560 7.040 0.3634
-1.500 0.5000 7.050 0.3290
7.060 0.2946
Con goal-seek se calcula la raíz (x1=): -1.644745296 7.070 0.2600
Con goal-seek se calcula la raíz (x2=): 7.144564250 7.080 0.2254
7.090 0.1908
7.100 0.1560
7.110 0.1212
7.120 0.0862
7.130 0.0512
7.140 0.0162
7.150 -0.0190
7.160 -0.0542
7.170 -0.0896
Como podemos ver con la segunda iteración se obtiene suficiente precisión, así como para la segunda raíz, cuya la tabla es:
Método de la Secante
Nótese que e incluye el valor inicial
>solve(f2,-2) //Se incluye el valor inicial (1ra raíz)
-1.64459895781
>solve(f2,8) //Se incluye el valor inicial (2da raíz)
7.14459895781
7. Comparación de resultados
Método Raíz 1 Raíz 2
1. Analítico -1.644599 7.144599
2. Gráfico -1.64 7.14
3. Excel (buscar objetivo) -1.644745 7.144564
4. Bisección -1.645019 7.144775
5. EULER (bisect) -1.644599 7.144599
6. Newton-Raphson -1.644599 7.144599
7. EULER (newton) -1.644599 7.144599
Método Raíz 1 Raíz 2
8. Secante -1.644599 7.144599
9. EULER (solve) -1.644599 7.144599
1. Analíticamente
2. Gráficamente
3. Con Excel (buscar objetivo)
4. Método de la bisección con 11 iteraciones (tome xl=-2, xu=0 para la primera raíz y xl=5, xu=10 para la segunda raíz).
Comparar con EULER (MAXIMA)
5. Método de Newton-Raphson, tomando -2 como valor inicial para la primera raíz y 7 para la segunda raíz. Comparar con
EULER (MAXIMA)
6. Método de Secante, comparar con EULER (MAXIMA)
1. Solución Analítica:
Factorizando la expresión dada con EULER (MAXIMA):
Aplicando la cuadrática:
4. Método de la bisección
Método de la bisección
Definimos y graficamos la función
>function f2:= "-2+7x-5x^2+6x^3";
>plot2d(f2,-0.5,0.5,title=f2,thickness=3,color=2)
[ -0.5 0.5 -7.5 1 ]
5. Método de Newton-Raphson
Método de Newton-Raphson
Definimos y graficamos la función
>function f2:= "-2+7x-5x^2+6x^3";
>plot2d(f2,-0.5,0.5,title=f2,thickness=3,color=2)
[ -0.5 0.5 -7.5 1 ]
Métodos de Newton-Raphson y de la Secante
Primero Derivamos la función numéricamente
>diff(f2,1)
15
Método de Newton-Raphson
>newton(f2,"diff(f2,1)",1)
0.333333333334
6. Método de la secante
Valores iniciales x(i-1) = 0.00
x(i) = 1.00
Iteración x(i-1) x(i) fx(i-1) fx(i) x(i+1) Error
0 0.00 1.0000 -2.0000 6.0000 0.250000 -75.0000%
1 1.00 0.2500 6.0000 -0.4688 0.304348 21.7391%
2 0.25 0.3043 -0.4688 -0.1636 0.333474 9.5699%
3 0.30 0.3335 -0.1636 0.0008 0.333333 -0.0422%
4 0.33 0.3333 0.0008 0.0000 0.333333 0.0002%
5 0.33 0.3333 0.0000 0.0000 0.333333 0.0000%
Método de la Secante
Definimos y graficamos la función
>function f2:= "-2+7x-5x^2+6x^3";
>plot2d(f2,-0.5,0.5,title=f2,thickness=3,color=2)
[ -0.5 0.5 -7.5 1 ]
Método de la Secante
>solve(f2,1)
0.333333333333
Con la función allroots podemos hallar todas las raíces de la ecuación y vemos que hay una raíz real y dos imaginarias:
>:: allroots(-2+7*x-5*x^2+6*x^3,x)
(%o39) [x = 0.33333333333333, x = 0.96824583655185 I + 0.25,
x = 0.25 - 0.96824583655185 I]
1. Gráficamente
2. Con Excel (buscar objetivo)
3. Método de la bisección (tome xl=0.50 para la primera raíz). Comparar con EULER (MAXIMA)
4. Método de Newton-Raphson, tomando 0.5 como valor inicial para la primera raíz. Comparar con EULER (MAXIMA)
5. Método de Secante, comparar con EULER (MAXIMA)
1. Solución Gráfica:
Utilizamos la utilidad en línea function evaluator and grapher v3.4 (www.zweigmedia.com).
Inicialmente tomamos un intervalo de -10 a +10 y vemos que la curva presenta una raíz real, reducimos el área del gráfico para
obtener mayor precisión: para x=0.57940 el valor de la función es de 0.001.
x f(x)
0.570 -0.1648
0.575 -0.0655
0.580 0.0326
0.585 0.1293
0.590 0.2249
0.595 0.3191
0.600 0.4121
Goal seek:
0.57834 0.0002
3. Método de la bisección
Método de la bisección:
Definimos y graficamos la función
>function f2:= "-26+82.3*x-88*x^2+45.4*x^3-9*x^4+0.65*x^5";
>plot2d(f2,0.5,0.6,title=f2,thickness=3,color=2)
[ 0.5 0.6 -1.7171875 0.390544 ]
4. Método de Newton-Raphson
Método de Newton-Raphson
Definimos y graficamos la función
>function f2:= "-26+82.3*x-88*x^2+45.4*x^3-9*x^4+0.65*x^5";
>plot2d(f2,0.5,0.6,title=f2,thickness=3,color=2)
[ 0.5 0.6 -1.7171875 0.390544 ]
5. Método de la Secante
Valores iniciales x(i-1) = 0.00
x(i) = 0.50
Iteración x(i-1) x(i) fx(i-1) fx(i) x(i+1) Error
0 0.00 0.5000 -26.00000 -1.71719 0.535358 7.0716%
1 0.50 0.5354 -1.71719 -0.90618 0.574866 7.3797%
2 0.54 0.5749 -0.90618 -0.08713 0.579069 0.7311%
3 0.57 0.5791 -0.08713 -0.00501 0.579325 0.0443%
4 0.58 0.5793 -0.00501 -0.00003 0.579327 0.0003%
Método de la Secante
Definimos y graficamos la función
>function f2:= "-26+82.3*x-88*x^2+45.4*x^3-9*x^4+0.65*x^5";
>plot2d(f2,0.5,0.6,title=f2,thickness=3,color=2)
[ 0.5 0.6 -1.7171875 0.390544 ]
Método de la Secante
>solve(f2,0.5)
0.579326592862
Con la función allroots podemos hallar todas las raíces de la ecuación y vemos que hay una real y cuatro imaginarias:
>:: allroots(-26+82.3*x-88*x^2+45.4*x^3-9*x^4+0.65*x^5,x)
(%o39) [x = 0.57932659286195, x = 0.86062151428295 I
+ 1.077269835444284, x = 1.077269835444284 - 0.86062151428295 I,
x = 2.333784494569321 I + 5.556143791201663,
x = 5.556143791201663 - 2.333784494569321 I]
6. Comparación de resultados
Método Raíz 1
1. Gráfico 0.57940
2. Excel (buscar objetivo) 0.57834
3. Bisección 0.57922
4. EULER (bisect) 0.579326592862
5. Newton-Raphson 0.579327
6. EULER (newton) 0.579326592862
7.Secante 0.579327
8.EULER (solve) 0.579326592862
Se desea diseñar un tanque esférico para atender las necesidades de agua de una pequeña población. Hallar una expresión para el
volumen de líquido en función de la altura del nivel del agua ( ).
Determinar para y .
Solución:
Consideramos la geometría del tanque:
El elemento diferencial del volumen lo podemos calcular como:
dV
V h
Reemplazando:
Integrando:
El desplazamiento de una estructura está definido por la siguiente ecuación para una vibración amortiguada:
Graficamos las ecuaciones y para determinar los rangos de corte con el eje , es decir las raíces de las ecuaciones. Para ello
utilizamos el software Euler:
Problema 8-15 – Desplazamiento Estructura
Problema 8-15 - Desplazamiento estructura con vibración amortiguada:
Se utiliza la ecuación:
Parámetros:
>k:=0.5; w:=3; t:=0:0.1:15;
Ecuación sin reducción de desplazamiento:
>ft:=8*%e^(-k*t)*cos(w*t);
>plot2d(t,ft, a=0,b=12,c=-10, d=8,thickness=1,color=green);
Ecuación con reducción de desplazamiento:
>ft1:=8*%e^(-k*t)*cos(w*t)-4;
>plot2d(t,ft1,thickness=1, xl="t", yl="despl.",color=red,>add):
Los tiempos (corte con el eje x) para cada uno de los casos son:
>bisect("8*%e^(-k*x)*cos(w*x)",0,0.6); print(%,4, unit=" seg.")
0.5236 seg.
>bisect("8*%e^(-k*x)*cos(w*x)-4",0,0.4); print(%,4, unit=" seg.")
0.3152 seg.
Esta diferencia se aprecia mejor si ampliamos las gráficas tal como se muestra en la siguiente figura entre 0 y 1 seg:
Un cable catenario es aquel que cuelga entre dos puntos que no están en la misma línea vertical, sujeto usualmente solo a su propio
peso, que actúa como carga uniformemente repartida por unidad de longitud del cable (pero no es uniforme para la distancia
horizontal :
Considerando el tramo AB del cable, y las
tensiones del mismo en sus extremos
entonces la ecuación que relaciona la altura
de la catenaria con la distancia horizontal es:
B
Siendo:
Para:
A y
Solución:
Para los valores dados podemos escribir:
a) Calcular d usando la Regla del Trapecio y las Reglas de Simpson 1/3 y 3/8.
b) Estime una ecuación polinómica que determinen la tendencia de los puntos dados e integre dicha ecuación para hallar d.
c) Compare resultados.
d) Determine el valor medio de la función mediante la expresión:
Solución:
a) Tabulamos los datos indicando los intervalos en donde pueden aplicarse los métodos pedidos:
Trapecio Simpson 1/3 Simpson 3/8
x f(x) Ec.3.8 Ec.3.10
1segmx2pt 2segmx3pt Ec.3.12 3segmx4pt
0.0 1.0
0.1 7.0 0.40 1.10
0.2 4.0 0.55 1.39
0.3 3.0 0.35 0.70
0.4 5.0 0.40 0.83
0.5 2.0 0.35 0.35
d= 2.05 2.15 2.22
b) Deducimos una fórmula polinómica para todo el intervalo utilizando la opción AGREGAR TENDENCIA de Excel, obteniendo
el polinomio:
Con
El polinomio anterior lo integramos con EULER entre los límites establecidos para obtener un valor de d=2.27:
Problema 21.14 Integración Numérica
Funciones de Integración de funciones y=f(x):
- integrate("funcion",a,b) -> metodo adaptativo
- romberg("funcion",a,b) -> frecuentemente el valor mas exacto
- gauss("funcion",a,b) -> solo 10 evaluaciones de la funcion
El método del polinomio representa en este caso una tendencia mayor pero considerando el obtenido puede reflejar el
mejor valor.
d) Para la determinación del valor medio se puede utilizar como numerador de la expresión dada el “d” obtenido de la
integración del polinomio:
Se sabe que
Solución:
e) Tabulamos los datos indicando los intervalos en donde pueden aplicarse los métodos pedidos:
Trapecio Simpson 1/3 Simpson 3/8
Volumen Presión Trabajo Solo Trapecio
3 Ec.3.8 Ec.3.10 Ec.3.12
m (kPa) Total Ec.3.14
1segmx2pt 2segmx3pt 3segmx4pt
0.5 420.0 0.0
591.00
2.0 368.0 591.0
3.0 333.0 675.33 350.5
4.0 326.0 329.5
6.0 326.0 652.0
1,861.50
8.0 312.0 638.0
10.0 242.0 554.0
224.50
11.0 207.0 224.5
Trabajo (kJ) = 815.50 675.33 1,861.50 3,352.33 3,339.5
En la gráfica siguiente se indica la fórmula aplicada en cada caso, según el número de segmentos iguales, cuya suma, es decir el
trabajo total, es 3,352.33 kJ.
Si calculamos solamente trapecios y sumamos los valores parciales, es decir aplicamos la Ec.3.14, obtenemos un valor de
3,339.5 kJ (última columna de la Tabla anterior), diferente al obtenido mediante la combinación de los métodos de Trapecio y
Simpson.
f) Estimamos una fórmula polinómica para todo el intervalo utilizando la opción AGREGAR TENDENCIA de Excel, obteniendo
el polinomio:
Con
El polinomio anterior lo integramos con EULER entre los límites establecidos para obtener un trabajo total de 3,345.0 kJ:
Ec.3.8
Ec.3.10
Ec.3.10
Ec.3.8
El método del polinomio representa en este caso una tendencia central y vemos que, con respecto a éste, el método
combinado se aparta mucho más que el método de solamente Trapecios. El valor del Trabajo puede estar representado,
para los datos usados, por un valor medio o si se prefiere por el valor del método del polinomio.
CONTROL DEPORTIVO
Chapra-Canale, Cap. 21 – Problema 21.27 - Pág.641 -
En un proceso de control deportivo se han tomado los siguientes datos a un corredor:
Velocidad (m/seg) 5.0 6.0 5.5 7.0 8.5 8.0 6.0 7.0 7.0 5.0
Tiempo (seg) 1.0 2.0 3.25 4.5 6.0 7.0 8.0 8.5 9.3 10.0
Se sabe que
Solución:
h) Tabulamos los datos indicando los intervalos en donde, según los datos dados, pueden aplicarse los métodos pedidos:
Como puede verse en esta Tabla según el número de segmentos iguales, la suma, es decir la distancia total, es 60.26 mt.
Si calculamos solamente trapecios y sumamos los valores parciales, es decir aplicamos la Ec.3.14, obtenemos un valor 60.43
mt (última columna de la Tabla anterior), ligeramente diferente al obtenido mediante la combinación de los métodos de
Trapecio y Simpson.
i) Estimamos dos fórmulas polinómicas (una de segundo grado y otra que se ajuste lo mejor posible a los datos) para todo el
intervalo considerado, utilizando la opción AGREGAR TENDENCIA de Excel, obteniendo los siguientes polinomios:
Como es de esperar, el polinomio de grado 5 se ajusta mejor a los datos del problema.
Los polinomios anteriores los integramos con EULER entre los límites establecidos para obtener la distancia total: 56.22 mt
para el polinomio de grado 5 y 60.56 mt para el de grado 2.
Problema 21.27 Integración Numérica
Problema 21.27 - Proceso deportivo
Definimos las funciones en EXCEL con la opción "agregar tendencia"
Integramos y graficamos
j) Comparación y conclusiones:
Método Distancia
Combinado Trapecio + Simpson 1/3 60.26 mt
Solo Trapecio 60.42 mt
Polinomio grado 2 60.56 mt
Polinomio grado 5 56.23 mt
La diferencia en los resultados de los dos polinomios es significativa (7%). Se puede pensar que ya que el polinomio grado 5
se ajusta más a los datos (ver gráfica) su resultado debe acercarse más al resultado de la combinación de los métodos
trapecio y Simpson. Sin embargo los resultados dicen lo contrario, es el de grado menor el que más se acerca al resultado
de los métodos numéricos utilizados, debido ello a un mejor ajuste estadístico de la curva del polinomio de grado 2.
Calcular las reacciones ( ) y la fuerza de tensión (en un plano) en el cable o templete de un mástil vertical (que permanece
vertical con respecto a la horizontal como en el caso de un bote) fijado en su base sin producir momento, sometido a una fuerza de
viento que varía con su altura según la siguiente expresión:
Es decir:
Graficamos la Ec.(b) anterior para diferentes valores de (entre z=0 y la altura H=30 ft) con un tamaño de paso=1:
i z f(z) [lb/ft] Apromedio
0.0 0.00
1 1.0 31.18 15.59
2 2.0 50.01 40.60
3 3.0 61.40 55.71
4 4.0 68.08 64.74
5 5.0 71.65 69.87
6 6.0 73.13 72.39
7 7.0 73.16 73.14
8 8.0 72.20 72.68
9 9.0 70.56 71.38
10 10.0 68.46 69.51
11 11.0 66.04 67.25
12 12.0 63.43 64.74
13 13.0 60.72 62.08
14 14.0 57.95 59.33
15 15.0 55.18 56.57
16 16.0 52.44 53.81
17 17.0 49.76 51.10
18 18.0 47.14 48.45
19 19.0 44.61 45.88
20 20.0 42.18 43.39
21 21.0 39.83 41.01
22 22.0 37.59 38.71
23 23.0 35.46 36.52
24 24.0 33.42 34.44
25 25.0 31.48 32.45
26 26.0 29.64 30.56
27 27.0 27.89 28.77
28 28.0 26.24 27.07
29 29.0 24.68 25.46
30 30.0 23.20 23.94
En la tabla anterior se muestra el cálculo de las áreas de cada segmento para un tamaño de paso de 1 (columna Apromedio) cuya suma,
es decir la integral, es . Repitiendo este cálculo para un tamaño de paso de 0.5 obtenemos .
Por Método de Simpson 1/3:
Aplicamos la Ec.(3.11) para la aplicación múltiple este método
EC.(3.11)
Repitiendo este procedimiento para tamaño de paso de 0.5 obtenemos que podemos tomar como valor final, es decir
Paso 2: Calcular :
i z z*f(z) Apromedio
0.0 0.00
1 1.0 31.18 15.59
2 2.0 100.02 65.60
3 3.0 184.21 142.12
4 4.0 272.33 228.27
5 5.0 358.27 315.30
6 6.0 438.75 398.51
7 7.0 512.12 475.44
8 8.0 577.62 544.87
9 9.0 635.05 606.34
10 10.0 684.56 659.80
11 11.0 726.46 705.51
12 12.0 761.22 743.84
13 13.0 789.32 775.27
14 14.0 811.32 800.32
15 15.0 827.73 819.52
16 16.0 839.08 833.40
17 17.0 845.87 842.48
18 18.0 848.58 847.23
19 19.0 847.66 848.12
20 20.0 843.51 845.58
21 21.0 836.53 840.02
22 22.0 827.08 831.81
23 23.0 815.47 821.28
24 24.0 802.02 808.74
25 25.0 786.98 794.50
26 26.0 770.62 778.80
27 27.0 753.14 761.88
28 28.0 734.77 743.95
29 29.0 715.67 725.22
30 30.0 696.01 705.84
Por Método del Trapecio Simple:
En la tabla anterior se muestra el cálculo de las áreas de cada segmento para un paso de 1 (columna A promedio) cuya suma, es decir la
integral, es . Repitiendo este cálculo para un tamaño de paso de 0.5
obtenemos .
Por Método de Simpson 1/3 aplicación múltiple:
Similarmente de la Tabla anterior obtenemos los valores para la Ec.(3.11):
Repitiendo este procedimiento para tamaño de paso de 0.5 obtenemos que podemos tomar como valor
.
De tal manera que la Ec.(c) será:
En la figura inicial del problema se muestra el diagrama de cuerpo libre con base en el cual se plantean las ecuaciones de equilibrio
en el siguiente paso.
Reemplazando en Ec.2:
Resolvemos con las funciones matriciales de Excel (MINVERSE y MMULT), obteniendo resultados iguales a los anteriores:
0.099504 1 0 1,480.56
[A] = -0.99504 0 1 {b} = 0.00
0 0 3 19,326.92
Punto de aplicación de F:
>h:= num/F
13.05
En el Ejemplo 24.2 anterior utilizar los métodos del Trapecio de aplicación múltiple, de Simpson 1/3 de aplicación simple, de
Simpson 3/8 y de Boole para calcular las integrales y determinar las incógnitas. Considerar en cada caso paso de 1 y de 0.5.
Solución:
por los diferentes métodos pedidos y aplicar sus resultados para hallar los valores de las incógnitas de acuerdo al procedimiento del
ejemplo 24.2 anterior.
Como puede verse en la tabla anterior, a excepción del método de Boole para 6 segmentos y seis puntos (Ec.3.13b), los demás
resultados tienen valores muy similares.
PROBLEMA 24.16 – CAP.24 – PÁG. 690
Chapra Canale – Cap.24 – Problema 24.14- Pág. 690
Determinar el trabajo efectuado al desplazar horizontalmente un bloque aplicándole una fuerza variable con un ángulo también
variable durante un trayecto determinado. Los datos de la fuerza aplicada y del ángulo de la fuerza leídos durante el recorrido
se muestran en la siguiente tabla:
Si la función de fuerza es fácilmente integrable, entonces el trabajo se encuentra integrando entre los límites del
recorrido.
Para el caso que nos ocupa la función de fuerza y su ángulo de aplicación son variables, por lo tanto la ecuación para hallar el trabajo
sería del tipo:
La función de fuerza anterior no puede integrarse directamente ya que no conocemos la ecuación que la represente y por lo tanto
debe recurrirse a una integración numérica. Para ello a partir de los datos elaboramos la siguiente tabla y gráfica que la representa:
x (ft) F(x).cosθ(x)
0 0.0000
5 1.5297
10 9.5120
15 8.7025
20 2.8087
25 1.0881
30 0.3537
La gráfica anterior muestra la función de fuerza F(x).cosθ(x). Con la tabla mostrada calculamos el área bajo la curva o trabajo por los
métodos del trapecio y de Simpson 1/3 para un tamaño de paso de 5:
Trapecio:
x (ft) F(x).cosθ(x) Apromedio
0 0.0000
5 1.5297 3.824261
10 9.5120 27.60415
15 8.7025 45.53624
20 2.8087 28.77819
25 1.0881 9.741993
30 0.3537 3.604364
Suma Pares 12.320693 Área= 119.0892
Suma Impares 11.3203033
Sin embargo un tamaño de paso grande no dá confianza en los resultados, que entre sí tienen una diferencia de cerca del 2%. Por
esta razón es conveniente hallar el trabajo con un tamaño de paso menor, digamos 1, para lo cual es necesario obtener las
coordenadas intermedias. Como no disponemos de la ecuación de la función de trabajo, recurrimos a leer de las curvas dadas los
valores cada paso y a partir de allí proceder a calcular . En la tabla siguiente se muestran los valores leídos.
x (ft) F(x) lb θ rad
x (ft) F(x) lb θ rad 16 13.50 0.96
0 0.00 0.50 17 12.60 1.05
1 1.90 0.70 18 11.60 1.15
2 3.80 0.95 19 10.80 1.24
3 5.50 1.15 20 10.50 1.30
4 7.40 1.32 21 10.50 1.35
5 9.00 1.40 22 11.10 1.40
6 10.00 1.35 23 11.70 1.44
7 11.00 1.20 24 12.10 1.46
8 11.90 1.02 25 12.00 1.48
9 12.50 0.85 26 11.10 1.50
10 13.00 0.75 27 10.00 1.50
11 13.40 0.73 28 8.30 1.50
12 13.90 0.75 29 6.50 1.50
13 14.05 0.78 30 5.00 1.50
14 14.20 0.84
15 14.00 0.90
Comparamos gráficamente estos valores leídos con los datos del problema y vemos que se ajustan bastante bien tal como se aprecia
en las figuras siguientes.
Hecha la anterior comprobación procedemos a calcular el área bajo la curva para calcular el trabajo realizado por los mismos
métodos empleados anteriormente (trapecio y Simpson 1/3) cuyos resultados se muestran en la siguiente tabla:
x (ft) F(x)*cosθ(x) Apromedio
x (ft) F(x)*cosθ(x) Apromedio 16 7.7425 8.22253
0 0.0000 17 6.2694 7.005958
1 1.4532 0.7266 18 4.7385 5.503925
2 2.2104 1.831798 19 3.5078 4.123127
3 2.2467 2.228538 20 2.8087 3.158269
4 1.8365 2.04159 21 2.2996 2.554154
5 1.5297 1.683101 22 1.8866 2.093103
6 2.1901 1.859886 23 1.5260 1.706296
7 3.9859 3.088001 24 1.3379 1.431926
8 6.2281 5.106995 25 1.0881 1.212977
9 8.2498 7.238922 26 0.7852 0.936621
10 9.5120 8.880872 27 0.7074 0.746277
11 9.9853 9.748646 28 0.5871 0.647245
12 10.1705 10.07791 29 0.4598 0.523455
13 9.9883 10.07941 30 0.3537 0.406739
14 9.4780 9.733154
15 8.7025 9.090256
De tal manera que por el método del trapecio obtenemos un valor de trabajo (suma de las áreas promedio -columna Apromedio-):
.
a= 0.0
b= 30.0
fa= 0.0000
fb= 0.3537
n= 30
Tenemos entonces:
En la gráfica siguiente se muestra la curva de la función de fuerza (tabla anterior), que integramos para obtener el trabajo, en la cual
podemos ver dos picos (en x=2.5 y x=12.5) no detectados en el cálculo de tamaño de paso de 5. Por esta razón los valores obtenidos
para un tamaño de paso de 1 pueden tomarse como verdaderos.
Calcular la variación del caudal (o rapidez del flujo) con respecto al diámetro del tubo del agua (dados varios diámetros) que
circula por una tubería si se conoce la distribución de la velocidad en su sección transversal:
A
r0=2, 3, 4, 5, 6 cm
“tajada” infinitesimal
Distribución de la velocidad “v ”
Solución:
Partimos de la ecuación:
Donde:
Para el caso de una tubería de radio tenemos:
Por lo tanto:
En la Ec.1 anterior está dado por lo tanto podemos calcular el valor de correspondiente. De esta forma establecemos una
tabla de variación.
Para el primer dato tenemos que, reemplazando en la Ec.1:
Aplicamos EULER para resolver numéricamente esta ecuación (ver Código de Programa y Figura) obteniendo
>::vel:=ode2('diff(y,x)=(-1/6)*(1-x/R)^(-5/6),y,x)
R - x 1/6
(%o38) y = R (-----) + %c
R
>::ic1(vel,y=R,x=0), f(x):=rhs(%)
R - x 1/6
(%o39) y = R (-----)
R
R - x 1/6
(%o40) f(x) := R (-----)
R
Integramos entre 0 y R
>cm3:= integrate(Q,0,R)
19.8863229973
>lts:= cm3/1000
0.0198863229973
>integrate(Q,0,R)
19.8863229973
Este procedimiento lo repetimos para cada valor de y elaboramos la siguiente tabla que refleja la variación pedida:
4 2 0.020
6 3 0.067
8 4 0.159
10 5 0.311
12 6 0.537
Radio r (in) 0 1 2 3 4 5 6 7 8
Velocidad u (m/s) 0.92 0.89 0.85 0.79 0.72 0.54 0.43 0.20 0.00
Donde:
Calcule Q:
a) Ajustando una curva polinomial a los datos e integrando,
b) Utilizando un método numérico aplicable
c) Utilizando EULER o MATLAB.
Solución:
a) Ajustando curva a los datos de la velocidad (con Excel):
Radio r (m) 0 0.025 0.051 0.076 0.102 0.127 0.152 0.178 0.203
Velocidad u (m/s) 0.92 0.89 0.85 0.79 0.72 0.54 0.43 0.20 0.00
Graficamos en Excel la tabla anterior y obtenemos una ecuación del polinomio que refleja la curva de tendencia (ver Ejercicios
NYOP.xlsx/Probl. Cap 21 Integr Num):
La ecuación a integrar, entre será entonces:
c) Utilizando EULER:
Problema 24.45 Chapra-Canale
Definimos la función
>function f1 := "2*%pi*(-23.585x^3+0.3258x^2+0.9095x)"
Y graficamos:
>plot2d(f1,a,b,title=f1, thickness=2, color=3)
[ 0 0.203 0 0.459340641768 ]