Modelos Matemáticos

Curso: Ecuaciones Diferenciales Ordinarias
Profesora; Aura Del C. Gómez
1
2
3
4
5
6
7
MÓDULO 3: APLICACIONES DE LAS ECUACIONES

DIFERENCIALES DE PRIMER ORDEN
3.1. TRAYECTORIAS ORTOGONALES Y OBLICUAS
Trayectorias Ortogonales: Se dice que una familia de curvas 𝑇(𝑥, 𝑦, 𝑘) = 0 con k

una constante arbitraria, es una trayectoria ortogonal para una familia de curvas
𝐹(𝑥, 𝑦, 𝐶) = 0 dada, si cualquier curva de la familia T corta a cada uno d ellos
miembros de la familia de curvas 𝐹(𝑥, 𝑦, 𝐶) = 0 bajo un ángulo recto.
Procedimiento:
1) Obtener la Ecuación Diferencial Ordinaria al haz de curvas 𝐹(𝑥, 𝑦, 𝐶) = 0 es
decir 𝐹(𝑥, 𝑦, 𝑦′) = 0
1
2) Debe sustituir y’ en el paso anterior por , y así tendremos la ecuación
−𝑦′
−1
diferencial asociada para encontrar las trayectorias ortogonales 𝐹 (𝑥, 𝑦, 𝑦′ ) = 0
3) Resuelva la ecuación diferencial obtenida, para obtener las trayectorias

ortogonales 𝑇(𝑥, 𝑦, 𝑘) = 0 .
Ejemplo
Encuentre las trayectorias ortogonales de la familia 𝑦 = 𝑐𝑥 3

Solución
Sea 𝑦 = 𝑐𝑥 3 , entonces; 𝑦´ = 3𝑐𝑥 2
3𝑦
𝑦´ =
𝑥
Consideremos
𝑑𝑦 −𝑥
=
𝑑𝑥 3𝑦
3𝑦𝑑𝑦 = −𝑥𝑑𝑥
3𝑦 2 = −𝑥 2 + 𝐶
3𝑦 2 + 𝑥 2 = 𝐶
8
Trayectorias Oblicuas: Sea 𝐹(𝑥, 𝑦, 𝐶) = 0 una familia uniparamétrica de curvas

que intersecan a las curvas de la familia 𝑇(𝑥, 𝑦, 𝑘) = 0 en un ángulo constante 𝛼 ≠
90°, se llama trayectoria oblicua de la familia dada.
La familia de trayectorias oblicuas se obtiene a través de la expresión
𝒅𝒚 𝑓(𝑥,𝑦)+𝑇𝑎𝑛𝛼
= 1−𝑓(𝑥,𝑦)𝑇𝑎𝑛𝛼
𝒅𝒙
𝑑𝑦
Primero encontraremos f(x,y) que es la Ecuación diferencial en términos de 𝑑𝑥
𝑑𝑦
asociada a la primitiva 𝐹(𝑥, 𝑦, 𝐶) = 0 , así; derivando tenemos = 𝑓(𝑥, 𝑦) la cual
𝑑𝑥
será sustituida en la expresión para las trayectorias oblicuas

𝑑𝑦 𝑓(𝑥,𝑦)+𝑇𝑎𝑛𝛼
= ; y luego se resuelve esta Ecuación Diferencial.
𝑑𝑥 1−𝑓(𝑥,𝑦)𝑇𝑎𝑛𝛼
Ejemplo
Encuentre las trayectorias oblicuas de y= Cx formando un ángulo de 30°
Solución
Sea y= Cx
𝑑𝑦
=𝐶
𝑑𝑥
𝑦 𝑑𝑦 𝑦
Como 𝐶=𝑥 entonces =
𝑑𝑥 𝑥
Luego reemplazando en la fórmula

𝑦
𝑑𝑦 + 𝑇𝑎𝑛30°
= 𝑥 𝑦
𝑑𝑥 1 − 𝑇𝑎𝑛30°
𝑥
𝑦 √3
𝑑𝑦
= 𝑥+ 3
𝑑𝑥 𝑦 √3
1−𝑥∙ 3
9
3𝑦 + √3𝑥
𝑑𝑦 3𝑥
=
𝑑𝑥 3𝑥 − √3𝑦
3𝑥
𝑑𝑦 3𝑦 + √3𝑥
=
𝑑𝑥 3𝑥 − √3𝑦
Como es una ecuación homogénea
Sea y=ux dy= udx+xdu
Sustitueyndo (3𝑥 − √3𝑦)𝑑𝑦 = (3𝑦 + √3𝑥)𝑑𝑥
(3𝑢 + √3)𝑑𝑥 − (3 − √3𝑢)(𝑢𝑑𝑥 + 𝑥𝑑𝑢) = 0
3𝑢𝑑𝑥 + √3𝑑𝑥 − 3𝑢𝑑𝑥 − 3𝑥𝑑𝑢 + √3𝑢2 𝑑𝑥 + √3𝑢𝑥𝑑𝑢 = 0
√3(1 + 𝑢2 )𝑑𝑥 − 𝑥(3 − √3𝑢)𝑑𝑢 = 0
3−√3𝑢 √3𝑑𝑥
Separando 𝑑𝑢 =
1+𝑢2 𝑥
3 √3𝑢 √3𝑑𝑥
𝑑𝑢 − 1+𝑢2 𝑑𝑢 =
1+𝑢2 𝑥
𝒚
𝟐√𝟑𝑻𝒂𝒏−𝟏 ( ) − 𝒍𝒏|𝒙𝟐 + 𝒚𝟐 | = 𝑪
𝒙
3.2 MECÁNICA ELEMENTAL
Se tendrá en cuenta el movimiento de un cuerpo a lo largo de una línea recta. Si el cuerpo
se mueve de forma que su velocidad “v” va variando (movimiento acelerado), su

𝑑𝑣
aceleración dada por , se debe a una o varias fuerzas que actúan en el sentido del
𝑑𝑡
movimiento o en sentido opuesto. La fuerza neta que actúa sobre la masa es la suma
(algebraica) de varias fuerzas.
Se utilizan leyes como la Segunda Ley de Newton para modelar los sistemas mecánicos.
10
Una de las aplicaciones de las Ecuaciones Diferenciales Ordinarias de primer orden es la
solución de problemas de temperatura, en los que un objeto absorbe calor del medio
circundante. Para dichos casos se puede establecer la ley de Calentamiento-Enfriamiento
de Newton que dice: “La temperatura de un cuerpo se modifica a una velocidad que es
proporcional a la diferencia de las temperaturas entre el cuerpo y el medio externo,

𝑑𝑇
siempre que el medio mantenga constante su temperatura” = 𝑘(𝑇 − 𝑇𝑚 )
𝑑𝑡
𝑑𝑖
Otra de las ecuaciones fundamentales en circuitos eléctricos es 𝐿 𝑑𝑡 + 𝑅𝑖 = 𝐸(𝑡), donde
L (henrios) se denomina la inductancia , R (ohmios) la resistencia, i (amperios) la
corriente y E (voltios) la fuerza electromotriz o FEM.
También veremos algunos modelados no lineales, para ciertos problemas como de
vaciado o drenado, mezclas y cuerdas colgantes.
3.2.1. Caída Libre
Al construir el modelo matemático de un cuerpo que se desplaza en un campo
de fuerzas, la Segunda Ley de Newton nos indica que cuando la fuerza neta que actúa
sobre un cuerpo no es igual a cero, entonces la fuerza neta es proporcional a su
aceleración “a”, o más exactamente F=ma, donde m representa la masa del cuerpo.
Ejemplo
Utilice un modelo de Ecuación Diferencial para resolver el siguiente problema
11
Una bola se lanza hacia arriba con una velocidad de 128 pies/seg.
a) ¿Cuál es la velocidad después de 2s?
𝐹 = 𝑚𝑎
𝑑𝑣
𝑚 = 𝑚𝑎
𝑑𝑡
𝑑𝑣
𝑚 = −𝑚𝑔
𝑑𝑡
𝑑𝑣
= −𝑔
𝑑𝑡
𝑑𝑣 = −𝑔𝑑𝑡
𝑣(𝑡) = −32𝑡 + 𝐶
𝐶𝑜𝑚𝑜 𝑣(0) = 0 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝐶 = 128
Así
𝑣(𝑡) = −32𝑡 + 128
𝑣(2) = −64 + 128 = 64𝑝/𝑠
b) ¿Cuándo regresará a su posición de partida?

𝑑𝑥
𝑣(𝑡) = = −32𝑡 + 128
𝑑𝑡
𝑥(𝑡) = −16𝑡 2 + 128𝑡 + 𝐶
Como x(0)=0 entonces C=0
Así
𝑥(𝑡) = −16𝑡 2 + 128𝑡
0 = −16𝑡(𝑡 − 8)
Entonces t=0 ó t=8
Luego regresará a su posición de partida a los 8 segundos.
c) ¿Cuál es la altura máxima que alcanza antes de regresar?

Al alcanzar la altura máxima cuando la velocidad es cero entonces,
𝑣(𝑡) = −32𝑡 + 128
0 = −32𝑡 + 128 entonces t=4seg.
Luego
𝑥(4) = −16(4)2 + 128(4)
x(4)= 256 pies
Alcanzará la altura máxima de 256 pies
12
3.2.2. Fuerza con Rozamiento
Cuando hay fuerzas de rozamiento, como las dadas por la resistencia del aire,
el modelo anterior cambia, al ser afectado por esta fuerza, de la siguiente
manera:
𝐹 = 𝑚𝑎
𝑑𝑣
𝑚 = 𝑚𝑔 − 𝑘𝑣
𝑑𝑡
En esta ecuación k es una constante positiva de proporcionalidad.
Ejemplo
La resistencia del aire del agua que actúa sobre un bote es proporcional a su
velocidad instantánea, y es tal que a 20p/seg. La resistencia del agua es 40 lbs. Si el
bote pesa 320 lbs. y el único pasajero pesa 160 lbs. y si el motor puede ejercer una
fuerza estable de 50 lbs. en la dirección del movimiento.
✓ Encuentre la máxima velocidad a la cual el bote puede viajar.

✓ Encuentre la distancia recorrida y la velocidad a cualquier tiempo, asumiendo
que el bote parte del reposo.
Solución
Consideremos el modelo representado por F = Fpropulsión – Resistencia
𝑑𝑣
F=ma = mg-R entonces 𝑚 𝑑𝑡 = 𝐹𝑝𝑟𝑜𝑝𝑢𝑙𝑠𝑖ó𝑛 – 𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎
Como 𝑅 = 𝑘𝑣 → 40 = 𝑘(20) → 𝑘=2

𝑤 (320+160)
𝑚= = = 15 𝑠𝑙𝑢𝑔
𝑔 32
13

𝑑𝑣
15 𝑑𝑡 = 50 − 2𝑣
𝑑𝑣
15 = 2(25 − 𝑣)
𝑑𝑡
𝑑𝑣 2𝑑𝑡
=
(25 − 𝑣) 15
2
−𝑙𝑛|25 − 𝑣| = 𝑡 + 𝐶1
15
Como v(0)=0 entonces
−𝑙𝑛|25| = 𝐶1
Así
2
−𝑙𝑛|25 − 𝑣| = 𝑡 − 𝑙𝑛|25|
15
2
𝑙𝑛|25| − 𝑙𝑛|25 − 𝑣| = 𝑡
15
|25| 2
𝑙𝑛 = 𝑡
|25 − 𝑣| 15
25 2
= 𝑒 15𝑡
25 − 𝑣
2
25 = (25 − 𝑣)𝑒 15𝑡
2 2
25 = 25𝑒 15𝑡 − 𝑣𝑒 15𝑡
2
25 − 25𝑒 15𝑡
𝑣= 2
𝑒 15𝑡
𝟐
𝒗 = 𝟐𝟓(𝟏 − 𝒆−𝟏𝟓𝒕 )
2
Ahora 𝑥(𝑡) = ∫ 25(1 − 𝑒 −15𝑡 )𝑑𝑡
2
375
𝑥(𝑡) = 25𝑡 + 𝑒 −15𝑡 + 𝐶2
2
375
Como x(0)=0 entonces 𝐶2 = − 2
𝟐
𝟑𝟕𝟓 𝟑𝟕𝟓
𝒙(𝒕) = 𝟐𝟓𝒕 + 𝒆−𝟏𝟓𝒕 −
𝟐 𝟐
Para la velocidad máxima 𝑡 → ∞ entonces 𝑣𝑚á𝑥 = 25𝑝/𝑠𝑒𝑔
14
3.2.3. Circuitos Eléctricos en Serie
Considere el circuito en serie simple, que contiene al inductor, resistor y capacitor,
como se muestra:
La corriente que circula en un circuito después que el interruptor se cierra se representa
mediante i(t); la carga sobre un capacitor en el tiempo t está señalado como q(t). Las letras
L , C, R representan inductancia, capacitancia y resistencia respectivamente, y por lo general
son constantes. Ahora, de acuerdo con la segunda Ley de Kirchhoff, el voltaje E(t) que se
genera en un lazo cerrado debe ser igual a la suma de las caídas de voltaje en el lazo.
Como la corriente i(t) está relacionada con la carga q(t), presente en el capacitor mediante
𝑑𝑞
𝑖= .
𝑑𝑡
Para utilizar la Ley en una ecuación de primer orden para determinar la corriente sería
𝑑𝑖 𝑑𝑞 1
𝐿 𝑑𝑡 + 𝑅𝑖 = 𝐸(𝑡), y para determinar la carga 𝑅 𝑑𝑡 + 𝐶 𝑞 = 𝐸(𝑡)
15
Ejemplo
Una fuerza electromotriz de 200 voltios se aplica a un circuito RC en serie con
resistencia de 1000 ohms y capacitancia de 5x10-6 faradios.
Encuentre la carga q(t) sobre el capacitor si i(0)=0.4.
Determine la carga y la corriente en t=0.005 seg.
Determine la carga conforme 𝑡 → ∞
Solución
𝑑𝑞 1
Según la Ley de Kirchoff 𝑅 𝑑𝑡 + 𝐶 𝑞 = 𝐸(𝑡)
Sean R=1000Ω, C=5x10-6 faradios, E(t)=200 voltios,
i(0)=0.4
𝑑𝑞 1
1000 𝑑𝑡 + 5𝑥10−6 𝑞 = 200
𝑑𝑞
+ 200𝑞 = 0.2 es una ecuación lineal de primer orden
𝑑𝑡
𝜇(𝑡) = 𝑒 200𝑡
𝑑𝑞
𝑒 200𝑡 + 200𝑒 200𝑡 𝑞 = 0.2𝑒 200𝑡
𝑑𝑡
𝑑[𝑒 200𝑡 𝑞]
= 0.2𝑒 200𝑡
𝑑𝑡
0.2
𝑒 200𝑡 𝑞 = 200 𝑒 200𝑡 + 𝐶
𝑞(𝑡) = 0.001 + 𝐶𝑒 −200𝑡 derivando 𝑖(𝑡) = −200𝐶𝑒 −200𝑡

16
Como i(0)=0.4 entonces C=-0.002

De allí que derivando 𝒊(𝒕) = 𝟎. 𝟒𝒆−𝟐𝟎𝟎𝒕
Ahora 𝑖(0.005) = 0.4𝑒 −200(0.005) = 0.147
𝒒(𝒕) = 𝟎. 𝟎𝟎𝟏 − 𝟎. 𝟎𝟎𝟐𝒆−𝟐𝟎𝟎𝒕
Para 𝑞(0.005) = 0.001 + (−0.002)𝑒 −200(0.005) = 0.003

Si 𝑡 → ∞ q=0.001
Otras aplicaciones
a) De acuerdo con la Ley empírica de enfriamiento de Newton (o de
calentamiento), la velocidad con la que la temperatura de un cuerpo
cambia es proporcional a la diferencia que hay entre la temperatura del
cuerpo y la del medio que lo rodea, la denominada temperatura ambiental.
Si 𝑇(𝑡) representa la temperatura de un cuerpo en el momento t, 𝑇𝑚 la
𝑑𝑇
temperatura del medio que lo rodea y la velocidad a la que cambia la
𝑑𝑡
temperatura del cuerpo, la Ley de Newton de enfriamiento y
calentamiento se traduce en el enunciado matemático:
𝑑𝑇 𝑑𝑇
∝ (𝑇 − 𝑇𝑚 ) ó = 𝑘(𝑇 − 𝑇𝑚 )
𝑑𝑡 𝑑𝑡
Donde k es una constante de proporcionalidad.
Ejemplo
Una pequeña barra metálica, cuya temperatura inicial era de 20°C, se deja caer en
un gran recipiente que contiene agua hirviendo.
17
Solución
El modelo de Enfriamiento-Calentamiento de Newton es:
𝑑𝑇
= 𝑘(𝑇𝑚 − 𝑇) ; 𝑘 > 0
𝑑𝑡
La cual es una ecuación diferencial lineal de primer orden, por lo tanto:
Sea 𝑇𝑚 = 100°𝐶 , 𝑇(0) = 20°𝐶 , 𝑇(1) = 22°𝐶
𝑑𝑇
+ 𝑘𝑇 = 100𝑘
𝑑𝑡
𝜇(𝑡) = 𝑒 𝑘𝑡
𝑑𝑇
𝑒 𝑘𝑡 + 𝑒 𝑘𝑡 𝑘𝑇 = 100𝑒 𝑘𝑡 𝑘
𝑑𝑡
𝑑[𝑇𝑒 𝑘𝑡 ]
= 100𝑒 𝑘𝑡 𝑘
𝑑𝑡
𝑇(𝑡) = 100 + 𝐶𝑒 −𝑘𝑡
Como T(0)=20 entonces C=-80
𝑇(𝑡) = 100 − 80𝑒 −𝑘𝑡
39
Como T(1)=22 entonces 𝑘 = 𝑙𝑛 (40)
𝟑𝟗
𝑻(𝒕) = 𝟏𝟎𝟎 − 𝟖𝟎𝒆−𝒍𝒏(𝟒𝟎)𝒕
Cuánto tiempo le llevará a la barra alcanzar los 90°C si se sabe que su temperatura
aumentó 2° en un segundo?
39
90°𝐶 = 100 − 80𝑒 −𝑙𝑛(40)𝑡 entonces t = 82.1 seg.
¿Cuánto le llevará alcanzar los 98°C?

39
98°𝐶 = 100 − 80𝑒 −𝑙𝑛(40)𝑡 entonces t = 145.7 seg.
b) La mezcla de dos soluciones salinas de distintas concentraciones da lugar
a una ecuación diferencial de primer orden para la cantidad de sal
contenida en la mezcla. Supongamos que un tanque grande de mezcla
inicialmente contiene cierta cantidad de galones de salmuera (es decir,
agua en la que se ha disuelto cierta cantidad de libras de sal). Otra
solución de salmuera se inyecta al tanque grande a una velocidad dada en
galones por minuto; en este flujo de entrada, la concentración de sal está

18
dada en libras por galón. Cuando la solución se mezcla bien en el tanque
se extrae a una velocidad que puede ser igual a la inyectada al entrar al
tanque o mayor o menor. Si A(t) indica la cantidad de sal (en libras) que
hay en el tanque, en el momento t, entonces la velocidad a que A(t) cambia
será:
𝑑𝐴
= 𝑉𝑒𝑙. 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑎𝑙𝑚𝑢𝑒𝑟𝑎 − 𝑉𝑒𝑙. 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑎𝑙𝑚𝑢𝑒𝑟𝑎
𝑑𝑡
= 𝑅𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑅𝑠𝑎𝑙𝑖𝑑𝑎
La velocidad de entrada 𝑅𝑒𝑛𝑡𝑟𝑎𝑑𝑎 a la que ingresa la sal al tanque es el producto de
la concentración de sal del flujo de entrada y la velocidad del flujo de entrada, esto
deberá ser medido en libras por minuto.
La velocidad de salida 𝑅𝑠𝑎𝑙𝑖𝑑𝑎 a la que sale la sal del tanque es el producto de la
𝐴(𝑡)
concentración de sal del flujo de salida dado por por la velocidad de
𝑉𝑜𝑙.𝑆𝑎𝑙𝑚𝑢𝑒𝑟𝑎
salida del flujo, esto deberá ser medido en libras por minuto. Note que el volumen
de salmuera en la velocidad de salida es la diferencia de las velocidades del flujo de
entrada y de salida multiplicadas por el tiempo t.
Esta ecuación diferencial se resuelve para determinar A(t).
Ejemplo
Un tanque grande se llena a toda su capacidad con 500 galones de agua pura. Hacia
el tanque se bombea salmuera, conteniendo 2 libras de sal, a velocidad de 5 galones
por minuto. Perfectamente mezclada, la solución se bombea hacia afuera a la misma
19
velocidad. Encuentre la cantidad A(t) de libras de sal presentes en el tanque en el
tiempo t.
Solución
Una ecuación diferencial de primer orden para la cantidad de sal contenida en una
mezcla (salmuera), está dada por:
𝑑𝐴
= 𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟. 𝑑𝑒 𝑠𝑎𝑙 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ∙ 𝑉𝑒𝑙. 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑎𝑙𝑚𝑢𝑒𝑟𝑎
𝑑𝑡
− 𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟. 𝑑𝑒 𝑠𝑎𝑙 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎
∙ 𝑉𝑒𝑙. 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑎𝑙𝑚𝑢𝑒𝑟𝑎
𝑑𝐴 5𝐴
= 5(2) − 500
𝑑𝑡
𝑑𝐴 𝐴
= 10 − 100 la cual es una ecuación lineal de primer orden
𝑑𝑡
𝑑𝐴 𝐴
+ 100 = 10
𝑑𝑡
1
𝜇(𝑡) = 𝑒 100𝑡
1
1 𝑡 1
𝑡 𝑑𝐴 𝐴𝑒 100
𝑒 100 + = 10𝑒 100𝑡
𝑑𝑡 100
1
𝑑 [𝐴𝑒 100𝑡 ] 1
= 10𝑒 100𝑡
𝑑𝑡
1 1
𝐴𝑒 100𝑡 = 1000𝑒 100𝑡 + 𝐶
−1
𝐴 = 1000 + 𝐶𝑒 100𝑡
Como A(0)=0 entonces C=-1000
−1
𝐴(𝑡) = 1000 − 1000𝑒 100𝑡
c) En hidrodinámica, la Ley de Torricelli establece que la velocidad “v” de

flujo de salida de agua a través de un orificio plano ubicado en la parte
inferior de un tanque lleno hasta una altura “h” será igual a la velocidad
que un cuerpo adquiriría en caída libre desde una altura “h”; es decir 𝑣 =
20
√2𝑔ℎ, donde g es la aceleración debida a la gravedad. Esta última

1
expresión proviene de la ecuación de energía cinética 2 𝑚𝑣 2 con la energía
potencial 𝑚𝑔ℎ y resolviendo para v. Suponga que un tanque lleno de

agua puede drenar mediante un orificio bajo la influencia de la gravedad.
Deseamos encontrar la altura “h” del agua restante en el tanque en el
tiempo t.
Si el área del orificio es 𝐴ℎ , y la velocidad del agua que sale del tanque
es 𝑣 = √2𝑔ℎ, , entonces el volumen de agua que abandona el tanque es
𝐴ℎ √2𝑔ℎ, de este modo si V(t) indica el volumen del agua que hay en el
𝑑𝑉
tanque en el tiempo t, entonces: = −𝐴ℎ √2𝑔ℎ donde el signo negativo
𝑑𝑡
indica que V disminuye.
Por otro lado el volumen del agua en el tiempo t, puede expresarse como
𝑉(𝑡) = 𝐴𝑤 ℎ ; donde 𝐴𝑤 es el área de la superficie del agua.
Si derivamos ambos miembros con respecto a t, tenemos :
𝑑𝑉 𝑑ℎ
= 𝐴𝑤 𝑑𝑡
𝑑𝑡
Igualando ambas ecuaciones tenemos

𝑑ℎ 𝐴
= −𝑐 𝐴 ℎ √2𝑔ℎ , siendo c la constante de fricción
𝑑𝑡 𝑤
Resolvemos esta ecuación diferencial para determinar la altura del agua en el

tiempo t.
Es importante observar que esta ecuación sigue siendo válida incluso cuando 𝐴𝑤
No es constante, por ejemplo si el tanque tiene forma cónica.
Ejemplo
Suponga que un tanque cónico está invertido, y el agua gotea por un orificio circular,
con radio de 2 pulgadas, localizado en el centro de la base circular. ¿El tiempo que
le lleva a este tanque lleno vaciarse será el mismo que para el tanque cuyo vértice está
hacia abajo? Asuma que el coeficiente de fricción es de c=0.6 y g=32p/seg.
21
Solución
El modelo matemático como ecuación diferencial de primer orden es:
𝑑ℎ 𝐴ℎ
= −𝑐 √2𝑔ℎ
𝑑𝑡 𝐴𝑤
𝜋
𝐴ℎ es el área del orificio 𝑟ℎ es el radio del orificio=36
𝐴𝑤 es el área de la superficie del agua 𝑟𝑤 es el radio de la superficie del

agua=?
𝑟𝑤 20 − ℎ 2(20 − ℎ)
𝑟𝑤 =? → = → 𝑟𝑤 =
8 20 5
2(20−ℎ 2 4𝜋(20−ℎ)2 1 2 𝜋
𝐴𝑤 = 𝜋 [ ] = ; 𝐴ℎ = 𝜋 (6) =
5 25 36
De allí que:
𝜋
𝑑ℎ (36)
= −0.6 √2(32)ℎ
𝑑𝑡 4𝜋(20 − ℎ)2
( )
25
Simplificando y separando
𝑑ℎ 5
(20 − ℎ)2 = − 𝑑𝑡
√ℎ 6
1 1 3 5
(400ℎ−2 − 40ℎ2 + ℎ2 ) 𝑑ℎ = − 𝑑𝑡
6
80 2 5
800√ℎ − √ℎ3 + 5 √ℎ5 = − 6 𝑡 + 𝐶
3
2560√5
Como h(0)=20 entonces 𝐶= 3
80 2 5 2560√5
Por lo tanto 800√ℎ − √ℎ3 + 5 √ℎ5 = − 6 𝑡 +
3 3
22
Además; para h=0 entonces 𝑡 = 1024√5seg. = 38.16 minutos
PROBLEMAS PROPUESTOS
1) Durante el verano, en el estacionamiento sin techo, de la Facultad de Eléctrica, la

temperatura dentro del auto de la profesora llega a 42°C, mientras que en el exterior
es constante e igual a 34°C. Cuando la profesora entra a la camioneta, enciende el
aire acondicionado con el termostato en 16°C. Si la constante de tiempo para la
camioneta es 1/K=2 horas y para la camioneta con el aire acondicionado es de
1/K1=1/3 hora. En qué momento llegará la temperatura dentro del carro a agradables
24°C.
𝑑𝑇
= 𝑘(𝑇 − 𝑇𝑚 ) + 𝑈(𝑡) ; 𝐾 < 0 , 𝐾𝑢 < 0
𝑑𝑡
Resp. 20 minutos aprox.
𝑥 2 +𝑦 2
2) Determine las trayectorias ortogonales de la familia 𝐶 = 𝑥+𝑦
𝑥 2 +𝑦 2
Resp. 𝐶 = 𝑦−𝑥
1 √3
3) Determine la curva que pasa por (2 , ) y corta a cada miembro de la familia 𝑥 2 +
2
𝑦 2 = 𝑐 2 , a un ángulo de 60°.
−1 (𝑦)−𝜋)
Resp. (𝑥 2 + 𝑦 2 )𝑒 (2√3𝑇𝑎𝑛 𝑥 3 =1
4) El aire del interior de un pequeño cuarto con dimensiones de 12x8x8 metros tiene 3%
de monóxido de carbono. Empezando en t=0, se sopla aire fresco que no contiene
23
monóxido de carbono, hacia el interior del cuarto a razón de 100m3/min. Si el aire

del cuarto sale al exterior a través de una abertura a la misma velocidad, ¿Cuándo
tendrá el aire del interior del cuarto 0.01% de monóxido de carbono?
Resp. t=43.805 minutos
5) Un electroimán industrial se puede modelar como un circuito RL, cuando se

energiza mediante una fuente de voltaje. Si la inductancia es de 10H yel
embobinado contiene 3Ω de resistencia. ¿Cuánto tiempo tarda un voltaje aplicado
en energizar el electroimán hasta 90% de su valor asintótico?
Resp. t=7.68 seg.
6) Una cochera sin calefacción ni aire acondicionado tiene una constante de tiempo de
2 horas. Si la temperatura exterior varía como una onda senoidal con un mínimo de
50°F a las 2 a.m. y un máximo de 80°F a las 2p.m., determine los instantes en que el
edificio alcanza sus temperaturas máximas y mínimas suponiendo que el término
exponencial se extingue.
Resp. 51°F a las 3.50 a.m. , 78.3°F a las 3.50 p.m.
7) Un objeto con masa de 60 Kg parte del reposo en la parte superior de un plano

inclinado a 45°. Suponga que el coeficiente de fricción cinética es de 0.05. Si la
fuerza debida a la resistencia del aire es proporcional a la velocidad del objeto,
digamos -3v, determine la ecuación de movimiento del objeto. Cuánto tiempo
tardará el objeto en llegar a la parte inferior del plano inclinado si la rampa mide 10
m de largo?
Resp. t= 1.768 seg.
8) La trayectoria de una señal eléctrica binaria entre compuertas en un circuito integrado

se puede modelar como un circuito RC. La fuente de voltaje modela la compuerta de
transmisión y el condensador modela la compuerta de recepción. Por lo general la
resistencia es de 100Ω y la capacitancia es muy pequeña, digamos un Picofaradio pF
( 10-12F). Si el condensador no tiene carga inicialmente y la compuerta de transmisión
24
cambia de manera instantánea de 0 a 5 V, cuánto tiempo tarda el voltaje en la

compuerta de recepción en alcanzar, digamos 3V?.
Resp. t= 9.2x10-11seg
9) Un Tanque contiene 200 litros de fluido en el cual se han disuelto 30 gramos de sal.
Si se bombea agua pura al tanque a velocidad de 4 litros por minuto; perfectamente
bien mezclada, la solución se bombea hacia afuera a la misma velocidad. Encuentre
el número A(t) de gramos de sal presentes en el tanque en el tiempo t.
1
Resp. 𝐴(𝑡) = 30𝑒 −50𝑡
10) Cuando la velocidad v de un objeto es muy grande, la magnitud de la fuerza debida a

la resistencia del aire es proporcional a 𝑣 2 y la fuerza actúa en dirección opuesta al
movimiento del objeto. Un proyectil con masa de 3 Kg se lanza desde el suelo, hacia
arriba y con una velocidad inicial de 500m/s. Si la magnitud de la fuerza debida a la
resistencia del aire es de 0.1𝑣 2 , ¿En qué momento alcanzará el proyectil su máxima
altura sobre el suelo? ¿Cuál es la máxima altura?
Resp. ℎ𝑚𝑎𝑥 = 101.2 𝑚𝑒𝑡𝑟𝑜𝑠
11) Una persona está en una esquina A de un embalse rectangular, tiene en la esquina
adyacente B una barca atada al extremo de una cuerda de 20 metros de longitud. La
persona se desplaza hacia C caminando por el borde del embalse y manteniendo tensa
la cuerda. Determine la situación de la persona y de la barca cuando ésta se encuentra
a 12 metros de AC
Sugerencia: Elija el sistema de coordenadas de modo que AC sea el eje de las abscisas y AB el
eje de las ordenadas. Considere P(x,y) la posición de la barca cuando la persona está en Q
Resp. 22 metros ; 6 metros
12) Un embudo de 10 cm de diámetro en la parte superior y 1cm de diámetro en la parte

inferior tiene una altura de 24cm. Si se llena de agua, hallar el tiempo que se tarda
en vaciar.
Resp. 13.7 seg.
25
Ecuaciones Diferenciales como Modelos Matemáticos

En términos generales, un modelo matemático es una descripción matemática de
algo. Dicha función puede ser algo tan simple como una función.
Expresar suposiciones en términos Formulación Matemática

Suposicione de ecuaciones diferenciales.
s
Si es necesario, modificar
las suposiciones o
Resolución de las
incrementar la resolución
ecuaciones diferenciales
del modelo
Verificar las
predicciones del Presentar predicciones del modelo, Obtención de soluciones
modelo con por ejemplo, de forma gráfica.
hechos conocidos
Consideremos los siguientes enunciados, y escriba un modelo matemático que los

represente
1) Escribir la ecuación diferencial de la curva definida, mediante las condiciones
dadas. En cada punto xy la pendiente de la tangente es igual al cuadrado de la
abscisa del punto:
𝒚´ = 𝒙𝟐
2) Expresar mediante ecuaciones diferenciales el principio físico “el radio se
desintegra a una velocidad proporcional a la cantidad Q del radio presente.
𝒅𝑸
= −𝒌𝑸
𝒅𝒕
3) “La diferencia de potencial е a través de un elemento de inductancia es igual al
producto de L por la velocidad de cambio de la corriente i en la inductancia”.
𝒅𝒊
𝑬 = 𝑳 𝒅𝒕
26

Modelos Matemáticos

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Modelos Matemáticos

Cargado por

Información del documento

Descripción 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

Modelos Matemáticos

Cargado por

Copyright:

Formatos disponibles

Curso: Ecuaciones Diferenciales Ordinarias

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

Profesora; Aura Del C. Gómez

MÓDULO 3: APLICACIONES DE LAS ECUACIONES

Trayectorias Ortogonales: Se dice que una familia de curvas 𝑇(𝑥, 𝑦, 𝑘) = 0 con k

3) Resuelva la ecuación diferencial obtenida, para obtener las trayectorias

Encuentre las trayectorias ortogonales de la familia 𝑦 = 𝑐𝑥 3

Profesora; Aura Del C. Gómez

Trayectorias Oblicuas: Sea 𝐹(𝑥, 𝑦, 𝐶) = 0 una familia uniparamétrica de curvas

será sustituida en la expresión para las trayectorias oblicuas

Encuentre las trayectorias oblicuas de y= Cx formando un ángulo de 30°

Luego reemplazando en la fórmula

Profesora; Aura Del C. Gómez

Sea y=ux dy= udx+xdu

Sustitueyndo (3𝑥 − √3𝑦)𝑑𝑦 = (3𝑦 + √3𝑥)𝑑𝑥

(3𝑢 + √3)𝑑𝑥 − (3 − √3𝑢)(𝑢𝑑𝑥 + 𝑥𝑑𝑢) = 0

3𝑢𝑑𝑥 + √3𝑑𝑥 − 3𝑢𝑑𝑥 − 3𝑥𝑑𝑢 + √3𝑢2 𝑑𝑥 + √3𝑢𝑥𝑑𝑢 = 0

√3(1 + 𝑢2 )𝑑𝑥 − 𝑥(3 − √3𝑢)𝑑𝑢 = 0

3.2 MECÁNICA ELEMENTAL

Se tendrá en cuenta el movimiento de un cuerpo a lo largo de una línea recta. Si el cuerpo

se mueve de forma que su velocidad “v” va variando (movimiento acelerado), su

(algebraica) de varias fuerzas.

Profesora; Aura Del C. Gómez

Una de las aplicaciones de las Ecuaciones Diferenciales Ordinarias de primer orden es la

circundante. Para dichos casos se puede establecer la ley de Calentamiento-Enfriamiento

proporcional a la diferencia de las temperaturas entre el cuerpo y el medio externo,

L (henrios) se denomina la inductancia , R (ohmios) la resistencia, i (amperios) la

corriente y E (voltios) la fuerza electromotriz o FEM.

También veremos algunos modelados no lineales, para ciertos problemas como de

vaciado o drenado, mezclas y cuerdas colgantes.

3.2.1. Caída Libre

Al construir el modelo matemático de un cuerpo que se desplaza en un campo

sobre un cuerpo no es igual a cero, entonces la fuerza neta es proporcional a su

Utilice un modelo de Ecuación Diferencial para resolver el siguiente problema

Profesora; Aura Del C. Gómez

𝑣(2) = −64 + 128 = 64𝑝/𝑠

b) ¿Cuándo regresará a su posición de partida?

c) ¿Cuál es la altura máxima que alcanza antes de regresar?

Profesora; Aura Del C. Gómez

3.2.2. Fuerza con Rozamiento

el modelo anterior cambia, al ser afectado por esta fuerza, de la siguiente

velocidad instantánea, y es tal que a 20p/seg. La resistencia del agua es 40 lbs. Si el

fuerza estable de 50 lbs. en la dirección del movimiento.

✓ Encuentre la máxima velocidad a la cual el bote puede viajar.

Consideremos el modelo representado por F = Fpropulsión – Resistencia

Como 𝑅 = 𝑘𝑣 → 40 = 𝑘(20) → 𝑘=2

Profesora; Aura Del C. Gómez

Para la velocidad máxima 𝑡 → ∞ entonces 𝑣𝑚á𝑥 = 25𝑝/𝑠𝑒𝑔

Profesora; Aura Del C. Gómez

3.2.3. Circuitos Eléctricos en Serie

Considere el circuito en serie simple, que contiene al inductor, resistor y capacitor,

La corriente que circula en un circuito después que el interruptor se cierra se representa