MÓDULO #5 DR. RICARDO LÓPEZ G.
3. Aplicaciones de las ecuaciones diferenciales de primer orden.
3.1 Trayectorias Ortogonales y Oblicuas.
Definición #1:
2 curvas f 1 y f 2 se dicen que son ortogonales en un punto, si y solo si, sus
tangentes T 1 y T 2 son perpendiculares en el punto de intersección.
T1
f2
f1 T2
90°
Método General:
y 1=f ( x , C )
y ' 1=f ' ( x , C )
y ' 1=h ( x )
−1
y '¿ =
h(x)
y 2=R( x , C)
Ejemplo:
1. Encuentre la trayectoria ortogonal de las siguientes familias de curvas:
c
a. y= 2
1+ x
b. 4 y + x2 +1+c e 2 y =0
2. Encuentre la constante para que las familias y 3=Cx y x 2+ a y 2=b 2, sean
ortogonales.
Solución:
c
a. 1.a. y=
1+ x 2
Primeramente, buscaremos la ecuación diferencial de esta familia de curvas.
Como solo hay un parámetro que es la c, derivamos una vez:
' −2 xc
y= 2
( 1+ x 2)
Eliminamos la c de la ecuación diferencial despejando c de la familia de curvas y la
reemplazamos en la ecuación diferencial,
2
c= y (1+ x )
−2 x ( y ( 1+ x ) )
2
'
y= 2
( 1+ x2 )
Simplificamos 1+ x 2, y nos queda:
−2 xy '
y=
1+ x 2
Que sería la pendiente de la familia de curvas, definimos la pendiente de la trayectoria
ortogonal, que sería el inverso negativo y tendríamos:
2
1+ x
y '¿ =
2 xy
Resolvemos la nueva ecuación diferencial, que nos sale de variables separables, donde
separamos las variables quedándonos:
1 1
2 ydy= dx+ xdx
2x 2
Integrando termino a término obtenemos la solución:
2 1 1 2
y = Lnx + x +C
2 4
b. 4 y + x2 +1+C e2 y =0
' ' 2y
4 y +2 x +2 y C e =0
' −2 x
y= 2y
4 +2C e
4 y + x 2 +1
Como C= , y al sustituirla en la ecuación diferencial
e2 y
' −2 x
y=
( )
2
2 y −4 y + x + 1
4 +e
e2 y
Reduciendo la Expresión,
' −2 x
y=
4−8 y −2 x 2−2
−x
¿
1−4 y−x 2
Definimos la pendiente de la trayectoria ortogonal
dy 4 y + x 2−1
=
dx x
Le damos la forma lineal
2
' −4 x −1
y= y=
x x
4
−∫ dx 1
F . I .=e x = 4
x
y x 2 +1 1
=∫ dx +c
x4 x x4
dx dx
¿ ∫ 3 −∫ 5 + C
x x
1 1
¿− ( x−2 )+ ( x−4 ) +C
2 4
2 1 4
y=−x + +C x
4
2, y 3=Cx y x 2+ a y 2=b 2
Para que dos familias de curvas sean ortogonales sus pendientes deben ser el
inverso negativo de la otra.
3
3 y
y =Cx C=
x
2 '
3 y y =C
' C
y= 2
3y
y3
x y
y'= 2 =
3y 3 x
2 2 2
x + a y =b
2x+2ayy’=0
' −x
y=
ay
−1
Y’=-
y'
y −1 ay
= =
3 x −x x
ay
1
Simplificando a =
3
Definición #2:
2 curvas f 1 y f 2 se dicen que son isogonales u oblicuas en un punto. Si y solo
si, Sus tangentes T 1 y T 2 no son perpendiculares en el punto de intersección.
f1
T1
≠ 90 °
T2
f2
Fórmula: Sí y’= f (x. y), entonces:
dy f ( x , y )−tanθ
=
dx TI 1+ f ( x , y ) tanθ
Ejemplo:
1. Encuentre la trayectoria isogonal u oblicua de las siguientes familias de curvas.
a. x 2=2 a ( y−x √ 3 ), θ=45 °
b. y 2=¿4 cx ¿ , θ=60 °
Solución:
a. x 2=2 a ( y−x √ 3 ), θ=45 °
Buscamos la pendiente de la familia de curvas, derivando y eliminando la
constante si aparece en la ecuación diferencial.
x2
2 x=2 a ( y ' −√ 3 ) , a=
2 ( y−x √3 )
Sustituimos a en la ecuación diferencial:
2
x
x=
2¿ ¿
Despejamos y’
2 y −x √ 3
y'=
x
Sustituimos en la fórmula de trayectoria oblicua:
2 y−x √3
−tan 45 °
dy x
=
dx TI 2 y−x √ 3
1+ tan 45°
x
Sustituyendo el valor de la tangente de 45° y reduciendo tanto el numerador como
el denominador a una sola fracción obtenemos:
dy
=
2 y−x ( √3− √22 )
dx x + y √6
dy 2 y−1.02 x
=
dx x +2.4 y
dy dv
Ecuación homogénea, y=vx =v + x
dx dx
dv 2 vx −1.02 x
v+ x =
dx x+ 2.4 vx
dv 2 v−1.02
x = −v
dx 1+ 2.4 v
1+2.4 v dx
dv =−2.4
2
v −0.41 v−0.42 x
Completamos el trinomio cuadrado perfecto y le damos la forma a v:
2.68 dv 2.4 ( v−0.7 ) dx
+ =2.4
2 2
( v−0.7 ) −0.92 ( v−0.7 ) −0.92 x
2.68
√ 0.92
tanh (
−1 v−0.7
√0.96 )
+1.2 ln [ ( v−0.7 ) −0.92=2.4 Lnx+C ]
2
Regresamos a la variable original sustituyendo v
( )
y
−0.7
[( ) ]
2
2.68 −1 x y
tanh + 1.2 ln −0.7 −0.92=2.4 Lnx+C
√ 0.92 √ 0.96 x
c. y 2=4 cx , θ=60 °
2
' y
2 y y =4 c c=
4x
y
y ' =2
x
y
2 −tan 60°
dy x
=
dx y
1+2 tan 60 °
x
dy 2 y−x √ 3
=
dx x+2 y √ 3
dy dv
y=vx =v + x
dx dx
dv 2 xv−x √3
v+ x =
dx x +2 xv √3
dv 2 v −√ 3
x = −v
dx 1+ 2 v √ 3
dv −2 v 2−v+ √ 3
x =
dx 2 v √ 3+1
2 ( v−6.92 ) 1.28 dx
+ =2 √3
2 2
( v−6.92 ) +0.47 ( v −6.92 ) + 0.47 x
1.28 v−6.92
ln [ ( v−6.92 )2 +0.47 ]+ tan−1 =2 √ 3 Lnx +C
√ 0.47 √ 0.47
3.2 Mecánica Elemental.
Definición:
La cantidad de movimiento (o ímpetu) de un cuerpo se define como el producto
mv de su masa m y su velocidad v.
Segunda ley de Newton:
La razón de cambio con respecto al tiempo de la cantidad de
movimiento de un cuerpo es proporcional a la fuerza resultante que actúa sobre el cuerpo
y está en la dirección de la fuerza resultante:
d
( mv ) =kF ó
dt
dv
∑ F=m dt
A. Problemas de caída libre.
Vamos a considerar algunos ejemplos de un cuerpo que cae a través del aire hacia
la tierra, o que se lanza verticalmente hacia arriba. En estas circunstancias el
cuerpo encuentra una resistencia del aire, donde la magnitud de la resistencia
depende de la velocidad del cuerpo, pero no se conoce una ley general que
exprese exactamente esa dependencia.
Caso #1: Caso #2
R W R W
+
{
Mov. x0
+
{
v0 Mov. x0
Superficie
v0
Superficie
dv dv
∑ F=m dt ; −R−W =m
dt
dv dv
∑ F=m dt ; W −R=m dt
Ejemplo:
1. Un objeto que pesa 32 libras se suelta desde el reposo a una altura de 50 pies
arriba de la superficie de un lago en calma. Antes de que el objeto llegue a la
superficie del lago, la resistencia del aire (en libras) está dado por 2v donde v es la
velocidad (en pies por segundos). Después de que el objeto empieza a hundirse en
el agua la resistencia (en libras) está dada por 6v. Además, sobre el objeto actúa
una fuerza de empuje hacia arriba de 8 libras. Calcule la velocidad 2 segundos
después que empezó a hundirse en el agua.
Análisis del problema:
Diagrama de cuerpo libre
Origen:
x 0=0
v 0=0
Caso #1 Mov
h=50pies +
W=32lbs R=2v
Caso #2
W=32lbs R=6v Fe=8lbs
Caso #1:
Analizaremos que pasa en el aire:
dv
∑ F=m dt
dv
W −R=m
dt
32 dv
32−2 v=
32 dt
dv
32−2 v=
dt
dv
=dt
32−2 v
−1
ln ( 32−2 v ) =t + c
2
−2 ( t +c )
32−2 v=e
32−e−2( t +c ) 32−e−2 c e−2 t
V= =
2 2
−2 c
32−e −2 c
Para v=0; t=0 0= , e =32
2
32−32 e−2t
v=
2
v=16 ( 1−e )
−2 t
Buscaremos la ecuación de desplazamiento ya que necesitaremos para el segundo
caso ese tiempo y la velocidad con que entre a la superficie de agua.
dx
v= =16 ( 1−e−2 t )
dt
( 1
x=16 t+ e−2 t +c
2 )
Como x=0 cuando t=0, obtenemos que c=-8
( 1
x=16 t+ e−2 t −8
2 )
Como la altura es de 50 pies, buscaremos el tiempo que tarde en recorrer esa
distancia.
X=50; t=?
( 1 −2t
50=16 t + e −8
2 )
Despejamos las variables de las constantes, para obtener
58 1 −2 t
=t+ e
16 2
Como no hay manera de despejar la t, debemos utilizar el método del tanteo para
encontrar t. Lo hacemos eliminando el término que tenga la exponencial y esa t
que nos sale empezamos a tantear.
58
t= =3.62, este valor lo sustituimos en la ecuación
16
58 −2 (3.62)
=3.62=3.62+e =3.62
16
Como nos dio el mismo valor no hubo necesidad de tantear y ese es el tiempo con
que llega a la superficie del agua. El cual lo sustituimos en la función velocidad y
encontramos la velocidad con que entra al agua.
v=16 ( 1−e )=15.99 P/s
−2 (3.62 )
Caso #2:
dv
∑ F=m dt
dv
32−8−6 v=
dt
dv
24−6 v=
dt
dv
=dt
24−6 v
−1
ln ( 24−6 v )=t+ c
6
Ln(24-6v)=-6(t+c)
24−6 v=e−6 (t +c )
24−e−6(t +c)
v=
6
Para t=0 y v=15.99 tenemos
24−e−6 c −6 c
15.99= , donde e =−69.96
6
24+69.96 e−6 t
v=
6
24+ 69.96 e−12
v ( 2 )= =4 p/seg
6
3
2. Una pelota que pesa lb se lanza verticalmente hacia arriba desde un punto que
4
se encuentra 6 pies arriba de la superficie terrestre y con una velocidad inicial de
20 p/seg. A medida que asciende actúa sobre ella la resistencia del aire que es
1
numéricamente igual a v , donde v es la velocidad (en pies por segundo) ¿a que
64
altura llegará la pelota?
3 1
W = lb R= v
4 64
Mov. Superficie
x 0=6 pies ; v 0=20 pies/ seg
dv
∑ F=m dt
3
−3 1 4 dv
− v=
4 64 32 dt
dv
=−dt
2
32+ v
3
3
2 ( 2
ln 32+ v =−(t+ c)
3 )
−2
( t+ c )
3
e −32
v=
−2 c
3
2
e −32
20= 3
2
{v=20
t =0
3
−2
c
3
e =45.33
−2
t
3
45.33 e −32
v=
2
3
−2
t
dx 45.33 e 3 −32
v= =
dt 2
3
−2
−3 t
( 45.33 ) e 3 −32t
2
x= +c
2
3
{
t=0 , entonces c=107.99
x=6
−2
t
x=101.99 e 3 −48 t +107.99
Para alcanzar la altura máxima en ese punto la velocidad se hace cero, o sea
v=0
t=?{
−2
t
45.33 e 3 −32
0=
2
3
−2
3
t32
e =
45.33
−2
3
t=ln
32
45.33 ( )
t=0.52seg.
−2
( .52)
h=101.99 e 3 −48 (.52 ) +107.99
h=155.14 pies
B. Fuerza de Rozamiento.
Supondremos que este movimiento se da sobre una superficie áspera y que
encontrará no solo la resistencia del aire sino también otra fuerza resistiva al
movimiento, debido a la rugosidad de la superficie. Esta fuerza adicional se
denominará fuerza de rozamiento, y está dada por F R =μN
μ=coeficiente de rozamiento o rugosidad de la superficie
N= Normal= Fuerza perpendicular que ejerce la superficie sobre el objeto
Plano Horizontal:
Movimiento
FE R
W N
Mov.
F1
+
∝
R
FR
W
θ
Diagrama de Cuerpo Libre: N F 1 sen ∝
F1
F∝1 cos ∝
R θ
FR Wcosθ
W Wsenθ
Haremos la sumatoria de fuerza que actúan perpendicularmente ya que no hay
movimiento, y la utilizaremos para definir la normal N.
∑ F=0
N−F 1 Sen ∝−Wcosθ=0
N=F 1 Sen ∝+ Wcosθ
Y con la normal definimos la fuerza de rozamiento F R =μN
Ahora haremos la sumatoria de fuerzas paralelas donde si hay movimiento y con esta se
define la ecuación diferencial a trabajar.
dv
∑ F=m dt
dv
Wsenθ + F1 cos ∝−F R −R=m
dt
Ejemplo:
1. Una caja de leche enlatada pesa 24 Lb y se suelta desde el reposo en la parte
superior de un plano metálico que mide 30 pies de longitud y está inclinado 45°
respecto de la horizontal. La resistencia del aire (en libras) es numéricamente igual
a un tercio de la velocidad (en pies por segundo) y el coeficiente de rozamiento es
0.4.
a. ¿Cuál es la velocidad con que se mueve la caja un segundo después de haberse
soltado?
b. ¿Cuál es la velocidad que lleva la caja en el momento de llagar a la parte
inferior del plano?
Solución:
Datos: Diagrama de cuerpo libre.
W=24 v
FR
x 0=0 N
v 0=0
1
L=30 pies R= v 45° 24 cos45°
3
θ=45 ° 45 ° 24sen45°
μ=0.4
∑ F=0
N=24 cos 45 °
F R =μN =( 0.4 )( 24 cos 45 ° )=6.788
dv
∑ F=m dt
dv
24 sen 45 °−F R −R=m
dt
1 24 dv
24 sen 45 °−6,788− v=
3 32 dt
dv
13.48−0.44 v =
dt
dv
=dt
13.58−0.44 v
−1
ln (13.58−0.44 v )=t +c
0.44
−0.44 ( t+ c )
13.58−e
v=
o .44
]
−0.44 c
t=0 → 0= 13.58−e
v =0 0.44
−0.44 c
e =13.58
−0.44 t
13.58−13.58 e
v=
0.44
v=30.86 ( 1−e )
−0.44 t
a. La velocidad de la caja un segundo después de haberse soltado es:
v=30.86 ( 1−e )=10.986 p/ s
−0.44
dx
=30.86 ( 1−e )
−0.44 t
b. v=
dt
x=30.86 t+ ( 1 −0.44 t
0.44
e +c)
]
t=0 → 0=30.86 0+ 1 +c
x=0 0.44 ( )
c=−70.145
x=30.86 t+ (
1 −0.44 t
0.44
e )
−70.145
‖ (
x=30 → 30=30.86 t+ 1 e−0.44 t
t=? 0.44 )−70.145
Reduciremos hasta donde podamos la expresión quedando
−0.44 t
3.25=t +2.27 e
Como podemos notar no hay forma de despejar t, debemos usar la regla de tanteo
que va a consistir en eliminar el término que contiene la exponencial y partir de
t=3.25, e ir tanteando si al sustituir en la expresión nos da muy alto entonces
empezamos a bajar o si nos da muy bajo empezamos a subir.
−0.44 ( 3.25 )
t=3.25 3.25=3.25+ 2.27 e =3.79
Como nos dio muy alto empezamos a bajar hasta encontrar un valor muy próximo
a 3.25
t=3 3.25=3+ 2.27 e−0.44 (3 )=3.6
−0.44 ( 2.25 )
t=2.5 3.25=2.5+ 2.27 e =3.255
−0.44 ( 2.49 )
t=2.49 3.25=2.49+ 2.27 e 3.2489
Entonces t=2.49 es el más próximo
Así que v (2.49) = 30.86 ( 1−e−0.44 (2.49 )) =¿20.54 p/s es la velocidad con que
llega a la parte inferior del plano.
2. Un muchacho se desliza con su trineo por una pendiente de 30°. El peso
combinado del muchacho y su trineo es 72 lb y la resistencia del aire (libras) es
numéricamente igual a 2 veces su velocidad (en pies por segundo). Si el
movimiento se inició desde el reposo y su velocidad después de 5 segundos es de
10 pies/seg. ¿Cuál es el coeficiente de rozamiento del trineo sobre la nieve?
N=72cos30°
N
F R =62.35 μ
FR
R=2v
30° 72cos30°
72 lb
72sen30°
30°
dv
∑ F=m dt
72 dv
72 sen 30 °−62.35 μ−2 v=
32 dt
dv
36−62.35 μ−24=2.25
dt
dv
16−27.71 μ−0.88 v =
dt
dv
=dt
16−27.71 μ−0.88 v
16−27,71 μ−e−0.88 ( t+ c )
v=
0.88
‖
−0.88 c
v =0 → 0= 16−27.71 μ−e
t=0 0.88
e−0.88 c =16−27.71 μ
(16−27.71 μ ) ( 1−t −0.88 t )
v=
0.88
t=5 →10= (16−27.71 μ) ( 1−e )
‖
−0.88 (5 )
v =10 0.88
μ=0.25