Tema 5

Modelos Físicos

5.1. Velocidad de escape

Por medio de la 2ª Ley de Newton aplicado a un cuerpo de masa m
sometido al campo gravitatorio

ma = −mg(r)
donde
dv d2 r
a= = 2
dt dt
y su teoría de gravitación
M
g(r) = G
r2
siendo

a Aceleración

v Velocidad

r Posición

R Radio del planeta (Tierra: R = 6371 km)

M Masa del planeta (Tierra: M = 5,972 · 1024 )

G Constante de gravitación universal G = 6,67 · 10−11 N 2 m2 /kg

1
5. Modelos Físicos 5.1. Velocidad de escape

g Gravedad (Gravedad en la superficie de la Tierra: para un radio

M
r = R → g(R) = G 2 = g ≈ 9,8 m/s)
R
mediante la regla de la cadena se obtiene la ecuación diferencial que modeliza
la velocidad de escape de un campo gravitatorio de un cuerpo a una posición
r > R respecto del centro del planeta

d2 r
m = −mg(r)
dt2
con la condición inicial de que la velocidad inicial del cuerpo desde la su-
perficie del planeta (a una posición R desde el centro del planeta) es vo :
v(R) = vo

dv dr dv M M R2 R2
= v = −mG 2 = −mG 2 2 = −mg 2
dr dt dr r R r r

dv R2
v = −g 2
dr r
Se trata de una ED de Variables Separables. Separando variables se ob-
tiene
dr
vdv = −gR2 2
r
e integrando a ambos lados de la ecuación se tiene la solución general
Z Z
dr v2 gR2
vdv = −gR2 → = +C
r2 2 r
y usando la condición inicial v(R) = vo

vo2 gR2 v2
= + C → C = o − gR
2 R 2
se obtiene la solución única que modeliza la velocidad

v2 gR2 vo2
= + − gR
2 r 2
s
2gR2 vo2
v(r) = ± + − gR
r 2

2
5. Modelos Físicos 5.2. Desintegración radiactiva

cuyo signo depende si el cuerpo se aleja del planeta (positivo) o se acerca

(negativo). En el caso de escape será la solución positiva
s
2gR2 vo2
v(r) = + − gR
r 2
Como los campos gravitatorios actúan en todo el espacio, salir del campo
será posible a una distancia infinita del planeta: r = ∞, luego fuera del campo
gravitatorio el cuerpo se encontrará en un estado en ausencia de fuerzas
externas, por tanto a distancia infinita permanecerá a velocidad constante
s
vo2
v(∞) = − gR = c
2
y la velocidad inicial vo mínima necesaria para escapar del campo será aquel
que a distancia infinita, la velocidad sea nula
s
vo2
v(∞) = − gR = 0
2
q
vo = 2gR
En el caso de la Tierra, la velocidad de escape es
q
vo = 2 · 9,8 m/s2 · 6371000 m ≈ 11,2 km/s

5.2. Desintegración radiactiva

Los átomos de ciertos elementos radiactivos se desintegran de forma es-
pontánea y se transforman en átomos de un nuevo elemento. Rutheford mos-
tró que la velocidad de desintegración dN/dt es proporcional (λ) a la cantidad
de átomos N , y esa observación se traduce en una ecuación diferencial que
relaciona una variable con su variación respecto al tiempo

Constante de
Variación
proporcionalidad
de la
entre la variable
variable
dN y su derivada Variable
= −λ · N
dt
3
5. Modelos Físicos 5.2. Desintegración radiactiva

donde λ es la constante de proporcionalidad.

Por tanto, si se conoce el dato de la cantidad de átomos No a un tiempo to ,
se tiene una condición inicial para el problema de desintegración radiactiva,
formando un problema de Cauchy de primer orden
dN
= −λN, N (to ) = No
dt
cuya solución única es
N (t) = No e−λ(t−to )
que modeliza la desintegración radiactiva de una sustancia en el tiempo,
teniendo a desintegrarse por completo para tiempo infinito.
Un dato que caracteriza la velocidad de desintegración es la media vida
(o vida media), que indica el tiempo transcurrido desde el inicia del proceso
de desintegración (to = 0) hasta que se desintegra la mitad de la sustancia
N (t = TM V ) = No /2

N (TM V ) = No e−λ(TM V −to )

9 9
No
= No e−λTM V
2
1 1
 
= e−λTM V → −λTM V = ln = 1 − ln 2 = − ln 2
ln
2 2
ln 2
TM V =
λ
Conociendo la vida media de elementos radiactivos, se puede determinar
la constante de proporcionalidad λ. Por ejemplo, la vida media del carbono-
14 es de unos TM V = 5,600 años, luego la constante del modelo es
ln 2 ln 2
λ= = = 1,2378 · 10−4
TM V 5600
Esta técnica se usa para datar materiales con contenido de carbono. Cono-
ciendo la cantidad de isótopos de carbono-14 a tiempo inicial: N (to = 0) = No
y a tiempo actual: N (t = tactual ) = pNo , donde

N (t = tactual ) Cantidad de isótopos de carbono-14 en el momento actual

p= =
N (to = 0) Cantidad de isótopos de carbono-14 en el momento inicial

4
5. Modelos Físicos 5.3. Dinámicas de poblaciones

entonces
N (t) = No e−λ(t−to )
9 9
pNo = No e−λt
ln p ln p
p = e−λt → −λt = ln p → t = − =−
λ 1,2378 · 10−4
Ejemplo 1. Para datar un fósil cuyas mediciones indicaron que se había
desintegrado un 80 % de su contenido en carbono, quedando por tanto un
p = 20 %, entonces la edad del fósil es de
ln p ln 0,2
t=− =− ≈ 13000 años
λ 1,2378 · 10−4

5.3. Dinámicas de poblaciones

5.3.1. Modelo de Malthus
Es el modelo más simple, que consiste en suponer que la variación de la
población es proporcional a la cantidad de población existente:
dp
= ap, p(to ) = po
dt
que consiste en un problema de Cauchy, siendo
p(t) cantidad de población en el instante t
a constante del modelo, con a > 0
to instante inicial
po población inicial
El modelo consiste en una ED de Vbles. Sep., cuya solución única es

p(t) = po ea(t−to )
Ejemplo 2. Sabiendo que la población humana existente en la Tierra en el
año 1965 era de 3340 millones de personas, aplicando el modelo de Malthus
con una constante del modelo a = 0,02, el modelo predice que para el año
1999 la población sería de
p(1999) = 3340 · e0,02(1999−1965) = 6592,8 millones

5
5. Modelos Físicos 5.3. Dinámicas de poblaciones

En este caso, el modelo predice con una buena aproximación el creci-

miento de la población mundial durante ese periodo, pues en el año 1999 se
superó la cantidad de 6000 millones de habitantes en el mundo. Sin embar-
go, el modelo predice un crecimiento exponencial sin límites a medida que
pasa el tiempo, luego deja de ser un modelo aceptable, necesitando alguna
modificación que limite el crecimiento.

5.3.2. Modelo de Verhulst (Ecuación Logística)

Es un modelo modificado del anterior modelo de Malthus, que consiste
en la ED
dp
= ap − bp2 , a, b > 0, p(to ) = po
dt
Para valores pequeños de p el modelo equivale al modelo de Malthus, pero
para valores grandes de p la derivada es negativa y la población decrece.
El modelo consiste en una ED de Vbles. Sep. y también de Bernoulli.
Resolviendo por Bernoulli, aplicando el cambio de vble u = p−1 la ED se
convierte en la EDL
du
+ au = b
dt
cuya solución es
b
u(t) = Ce−at +
a
y deshaciendo el cambio se tiene

1
p(t) =
b
Ce−at +
a
y usando la condición inicial p(to ) = po finalmente se tiene la solución única
del problema de Cauchy
1
po =
b
Ce−ato +
a
b
po Ce−ato = 1 − po
a
b
1 − po 1
!
b ato
C= a = − e
po e−ato po a

6
5. Modelos Físicos 5.3. Dinámicas de poblaciones

1
p(t) = !
1 b −a(t−to ) b
− e +
po a a
En este caso, el modelo es convergente a tiempo infinito
1 a
lı́m p(t) = ! =
t→∞ 1 b b b
− 0+
po a a
Ejemplo 3. Aplicando la ecuación logística ajustada con la evolución de la
población mundial , usando los datos de población en 1950 como condición
inicial, y en 1990 y 2020 para ajustar las constantes del modelo:
p(to = 1950) = po = 2,54 · 109 habitantes
p(t1 = 1990) = p1 = 5,33 · 109 habitantes
p(t2 = 2021) = p2 = 7,80 · 109 habitantes
con la ecuación logística se tiene
1
p(t) = !
1 b −a(t−to ) b
− e +
po a a
1

p(t1 ) = = p1 

! 
1 b −a(t1 −to ) b 

− e +



po a a





1 

p(t2 ) = = p2 

! 
1 b −a(t2 −to ) b



− e +



po a a 

1 b
 
−
−1 
p a
a=  1
ln 
!  
1 b −a(t1 −to ) b 1  t1 − to  1 b

− e + = 
 −
po a a p1 


 po a
→ k
1
!
b
  
1 b −a(t2 −to ) b 1 

− e + =

 − 
−1

po a a p2 p a
 
a= ln  2
t2 − to  1 b

−
po a

7
5. Modelos Físicos 5.4. Circuitos eléctricos

obteniendo una ecuación transcendente

1 b 1 b
   
−  − 
1 p

a 1 p

a
ln  1 = ln  2
t1 − to 1 b t2 − to 1 b

− −
  
po a po a
en el que se asume que la fracción de constantes a y b es la variable de la
ecuación. Resolviendo mediante una herramienta de álgebra computacional
devuelve la solución
b
≈ 8,039 · 10−11
a
Despejando las constantes se tiene
1 b
 
− 
−1 p

a
a= ln  2  ≈ 0,02681
t2 − to 1 b
−

po a

b = a · 8,358 · 10−11 ≈ 2,155 · 10−12

En base a este modelo, la población mundial convergerá a
a 0,02731
lı́m p(t) = = ≈ 12441 millones
t→∞ b 2,282 · 10−12
y en el año 2050 será de
1
p(2050) = ! ≈ 9,82 · 109 habitantes
1 b −a(2050−to ) b
− e +
po a a

que se asemeja a la predicción de la figura (??) de 9,74 · 109 habitantes.

5.4. Circuitos eléctricos

Los circuitos eléctricos se pueden esquematizar mediante un sistema inter-
conectado de diferentes componentes eléctricos. Los componentes más básicos
son
Resistencia Su capacidad R se mide en ohmios (Ω)

8
5. Modelos Físicos 5.4. Circuitos eléctricos

Bobina Su inductancia L se mide en henrios (H)

Condensador Su capacitancia C se mide en faradios (F )

Generador La fuerza electromotriz E se mide en voltios (V )

Estos componentes se diferencian por la manera en que interactúan con la
corriente eléctrica, que sigue ciertas reglas reflejadas en las siguientes leyes:
1. 1ª Ley de Kirchoff. La suma de las corrientes que entran en el nodo
es igual a la suma de las corrientes que salen del nodo (la suma de las
corrientes que pasan por el nodo es igual a cero)
X
Ik = 0
k

2. 2ª Ley de Kirchoff. En un circuito cerrado, la suma de todas las caídas

de tensión es igua a la tensión total suministrada
X
E= Vk
k

3. Ley de Ohm. La tensión entre los extremos de una resistencia es pro-

porcional a la corriente eléctrica que circula por la resistencia, siendo
la constante de proporcionalidad la capacidad de la resistencia R

VR = R · IR

4. Ley de Faraday. La tensión entre los extremos de una bobina es propor-

cional a la variación de la corriente eléctrica que circula por la bobina
con respecto al tiempo, cuya constante de proporcionalidad es la in-
ductancia de la bobina L
dIL
VL = L ·
dt

5. Ley de Coulomb. La tensión entre los extremos de un condensador es

proporcional a la carga del condensador Q, cuya constante de propor-
cionalidad es la inversa de la capacitancia C
1
VC = Q
C
9
5. Modelos Físicos 5.4. Circuitos eléctricos

y como la corriente eléctrica es, por definición, la cantidad de carga

eléctrica que atraviesa una sección de un conductor en cada unidad de
tiempo, aplicado al condensador, será
dQ
IC =
dt

5.4.1. Circuito RLC en serie

Figura 5.1: Circuito RLC en serie

Dado un circuito RLC en serie compuesto por una resistencia R, una bo-
bina de inductancia L, un condensador de capacitancia C, y alimentado por
un generador eléctrico con una fuerza electromotriz E , la corriente eléctrica
I(t) del circuito se puede modelizar aplicando la 2ª Ley de Kirchoff

VR + VL + VC = E

de forma que
dIL 1
L + RIR + Q = E
dt C
y derivando toda la ecuación respecto al tiempo se tiene
IC
 q 
d2 IL dIR 1 dQ
L 2
+R + =E
dt dt C dt
d2 IL dIR 1 dE
L 2
+R + IC =
dt dt C dt
10
5. Modelos Físicos 5.4. Circuitos eléctricos

pero por la 1ª Ley de Kirchoff, la corriente que circula por cada componente
es la misma
IR = IL = IC = I
de manera que se obtiene una ecuación diferencial lineal de segundo orden
que modeliza la corriente eléctrica del circuito

d2 I dI 1 dE
L 2
+R + I =
dt dt C dt

o también expresado en función de la carga eléctrica en el circuito

d2 Q dQ 1
L + R + Q=E
dt2 dt C

Ejemplo 4. Dado un circuito RLC en serie con una resistencia de R = 10 Ω

(ohmios), una bobina de L = 5 H (henrios) y un condensador de C = 0,2 F
(faradios), y una fuerza electromotriz de
125
E(t) = sin(2t)
2
Usando la segunda Ley de Kirchoff se obtiene la ecuación diferencial que
modeliza la corriente eléctrica

VR + VL + VC = E

donde
dVR dIR
VR = R · IR → = R
dt dt
dIL dVL d2 IR
VL = L · → = L 2
dt dt dt
1 dVC 1
VC = Q → = IC
C dt C
y por la primera Ley de Kirchoff la corriente que pasa por cada componente
es la misma, luego
dVR dVL dVC dE
+ + =
dt dt dt dt
d2 I dI 1 dE
L 2 +R + I =
dt dt C dt

11
5. Modelos Físicos 5.4. Circuitos eléctricos

d2 I dI
5 2
+ 10 + 5I = 125 cos(2t)
dt dt
Se trata de una EDL de segundo orden. El proceso de resolución se inicia
calculando la solución de la EDLH asociada, expresada en forma estándar

d2 I dI
2
+2 +I =0
dt dt
para luego obtener la solución particular de la EDL No Homogénea

d2 I dI
2
+ 2 + I = 25 cos(2t)
dt dt
mediante el método de los coeficientes indeterminados o del método de la
variación de las constantes.
La solución general de la EDL será la suma de la solución homogénea más
la solución particular, siendo en este caso

I(t) = c1 e−t + c2 te−t − 3 cos(2t) + 4 sin(2t)

5.4.1.1. Vibraciones libres sin amortiguamiento: R = 0, E = 0

d2 I 1 d2 I
L 2
+ I = 0 → 2
+ ωo2 I = 0
dt C dt
2
dQ 1 d2 Q
L 2 + Q = 0 → 2 + ωo2 Q = 0
dt C dt
2
siendo ωo = 1/(LC)
La solución general en ambas EDs abarcan un conjunto de soluciones
equivalentes
I(t) = c1 cos(ωo t) + c2 sin(ωo t)
Q(t) = c1 cos(ωo t) + c2 sin(ωo t)
que también se pueden expresar como
c1 c2 !
cos(ωo t) + sin(ωo t) =
q
I(t) = c21 + c22 q q
c21 + c22 c21 + c22
q  
= c21 + c22 cos α cos(ωo t) + sin α sin(ωo t)
= A cos(ωo t − α)

12
5. Modelos Físicos 5.4. Circuitos eléctricos

c1 c2 !
cos(ωo t) + q sin(ωo t) =
q
Q(t) = c21 + c22 q
c21 + c22 c21 + c22
q  
= c21 + c22 cos α cos(ωo t) + sin α sin(ωo t)
= A cos(ωo t − α)
en el que la solución representa un movimiento ondulatorio donde
A representa la amplitud de la onda

ωo representa la frecuencia angular de la onda

α representa el ángulo de fase de la onda

El periodo será T = 2π/ωo

5.4.1.2. Vibraciones libres amortiguadas: R > 0, E = 0

d2 I dI 1
L 2 +R + I =0
dt dt C
2
dQ dQ 1
L 2 +R + Q=0
dt dt C
Las raíces de la ecuación característica son
q
−R ± R2 − 4L/C
λ=
2L
Dependiendo del valor del discriminante R2 − 4L/C se tienen tres casos
Caso 1 Circuito sobreamortiguado (sobreamortiguamiento): R2 − 4L/C >
0. Las dos raíces son reales y distintas

I(t) = c1 eλ1 t + c2 eλ2 t

Q(t) = c1 eλ1 t + c2 eλ2 t

Caso 2 Circuito críticamente amortiguado (caso crítico): R2 − 4L/C = 0.

Las dos raíces son reales e iguales

I(t) = c1 eλ1 t + c2 teλ2 t

Q(t) = c1 eλ1 t + c2 teλ2 t

13
5. Modelos Físicos 5.4. Circuitos eléctricos

Caso 3 Circuito subamortiguado: R2 − 4L/C < 0. Las dos raíces son com-
plejas
R   R
I(t) = e− 2L t c1 cos(µt) + c2 sin(µt) = Ae− 2L t cos(µt − α)
q
R2 − 4L/C
µ=
2L

5.4.1.3. Vibraciones forzadas con amortiguamiento: R > 0, E 6= 0

Considerando una fuerza electromotriz E(t) = Eo cos(ωt)

d2 I dI 1
L 2
+ R + I = −ωEo sin(ωt)
dt dt C
2
dQ dQ 1
L 2 +R + Q = Eo cos(ωt)
dt dt C
Resolviendo la EDL No homogénea para la carga, la solución homogénea
equivale al obtenido en el apartado anterior. La solución particular será

Qp = M cos(ωt) + N sin(ωt) = A cos(ωt − α)

donde
1
 
− Lω 2 Eo
M =  C
2 2
1

− Lω 2 + R2 ω 2
C
RωEo
N =  2 2
1

− Lω 2 + R2 ω 2
C
y la amplitud de la respuesta en función de la frecuencia externa es
Eo
A(ω) = s 2
1
− Lω 2 + R2 ω 2
C

14
5. Modelos Físicos 5.4. Circuitos eléctricos

5.4.1.4. Vibración forzada sin amortiguamiento: R = 0, E 6= 0

Considerando una fuerza electromotriz E(t) = Eo cos(ωt)

d2 I 1 d2 I ωEo
L 2
+ I = −ωE o sin(ωt) → 2
+ ωo2 I = − sin(ωt)
dt C dt L
d2 Q 1 d2 Q Eo
L 2 + Q = Eo cos(ωt) → 2 + ωo2 Q = cos(ωt)
dt C dt L
donde ωo2 = 1/(LC)
Resolviendo la EDL No homogénea para la carga, la solución homogénea
equivale al obtenido en el apartado anterior. La solución particular dependerá
de dos casos

Caso 1 ω 6= ωo

Qp = M cos(ωt) + N sin(ωt) = A cos(ωt − α)

y la solución general será

Eo
Q(t) = Ao cos(ωo t − α) +   cos(ωt)
L ωo2 − ω 2

Caso 2 ω = ωo  
Qp = t M cos(ωt) + N sin(ωt)
y la solución general será
Eo
Q(t) = Ao cos(ωo t − α) + t sin(ωo t)
2ωo L
que reproduce el fenómeno de resonancia pura.

5.4.2. Circuito RLC en paralelo

Dado un circuito RLC en paralelo compuesto por una resistencia R, una
bobina de inductancia L, un condensador de capacitancia C, y alimenta-
do por un generador eléctrico con una fuerza electromotriz E , la corriente
eléctrica I(t) del circuito se puede modelizar aplicando la 1ª Ley de Kirchoff

IR + IL + IC = 0

15
5. Modelos Físicos 5.4. Circuitos eléctricos

Figura 5.2: Circuito RLC en paralelo

de forma que, por propiedades de las componentes

VR dIR 1 dVR
IR = → =
R dt R dt
dIL VL
=
dt L
1 dIC d2 V C
VC = Q → = C 2
C dt dt
y derivando respecto al tiempo en la primera Ley y sustituyendo se tiene

1 dVR VL d2 VC
+ +C 2 =0
R dt L dt
pero por la 2ª Ley de Kirchoff, la tensión entre los extremos de cada compo-
nente será el mismo
VR = VL = VC = V
por tanto, se obtiene la EDL de segundo orden

d2 V 1 dV V
C 2
+ + =0
dt R dt L

16
5. Modelos Físicos 5.4. Circuitos eléctricos

Por otro lado, de la primera Ley también se obtiene

VR dVC
+ IL + C =0
R dt
y por la 2ª Ley de Kirchoff las tensiones son equivalentes
dIL
VR = VC = VL = L
dt
obteniendo la EDL se segundo orden

d2 IL L dIL
LC + + IL = 0
dt2 R dt

Por último, otra forma de modelizar el circuito consiste en aplicar el

cambio de variable U = IL0 , entonces

dIL
=U
dt
dIL


 = U
dt




dU U IL




 = − −
dt RC LC
pero como
dIL V
=
dt L
se obtiene el sistema de EDL de primer orden

 dIL V

 =
dt L




dV V IL





 = − −
dt RC C

Ejemplo 5. Dado un sistema RLC en paralelo con una resistencia de R =

1 Ω (ohmios), una bobina de L = 1 H (henrios) y un condensador de C =
0,5 F (faradios), el sistema se modeliza mediante el sistema de EDLs

17
5. Modelos Físicos 5.5. Vibraciones mecánicas

dIL


 = V
dt




dV




 = −2V − 2IL
dt
cuyos valores propios del sistema son λ = −1 ± i, donde la solución es
! ! !
IL sin t cos t
= c1 e−t + c2 e−t
V cos t − sin t − cos t − sin t
De igual manera, resolviendo las EDL de segundo grado
d2 IL dIL
2
+2 + 2IL = 0
dt dt
d2 V dV
+ 2 + 2V = 0
dt2 dt
tienen las mismas soluciones generales
IL (t) = c1 e−t cos t + c2 e−t sin t
V (t) = C1 e−t cos t + C2 e−t sin t
donde coincide la solución de IL en la EDL de segundo orden con el sistema de
EDLs de primer orden. Sin embargo, la solución de V , aunque no coincidan,
se tratan de combinaciones lineales de soluciones linealmente independientes,
pero que pueden expresarse unas en función de las otras, puesto que
   
V (t) = c1 e−t cos t − sin t + c2 e−t − cos t − sin t =
   
= c1 − c2 e−t cos t + −c1 − c2 e−t sin t =
= C1 e−t cos t + C2 e−t sin t

5.5. Vibraciones mecánicas

Dado un sistema mecánico compuesto por una masa m, un resorte de
constante elástica k, y un amortiguador con constante c, la ecuación que
gobierna su movimiento está determinado por la 2ª Ley de Newton
d2 y X
m = F
dt2
donde las fuerzas que intervienen son

18
5. Modelos Físicos 5.5. Vibraciones mecánicas

Figura 5.3: Sistema de masa-resorte-amortiguador

Fuerza elástica Fk = −k · y(t)

dy
Fuerza de amortiguamientoFc = −c ·
dt
Fuerza externa Fext = F (t)

De este modo, la ecuación que gobierna el movimiento para un sistema de

masa - resorte - amortiguador será

d2 y dy
m 2
+ c + ky = F (t)
dt dt

5.5.1. Vibraciones libres sin amortiguamiento: c = 0, F =

0

d2 y d2 y
m + ky = 0 → + ωo2 y = 0
dt2 dt2
siendo ωo2 = k/m
La solución general será

y(t) = c1 cos(ωo t) + c2 sin(ωo t)

19
5. Modelos Físicos 5.5. Vibraciones mecánicas

que también se pueden expresar como

c1 c2 !
cos(ω t) + sin(ω t)
q
2 2 q o o
y(t) = c1 + c2 q
=
c21 + c22 c21 + c22
q  
= c21 + c22 cos α cos(ωo t) + sin α sin(ωo t)
= A cos(ωo t − α)
en el que la solución representa un movimiento ondulatorio donde
A representa la amplitud de la onda
ωo representa la frecuencia angular de la onda
α representa el ángulo de fase de la onda
El periodo será T = 2π/ωo

5.5.2. Vibraciones libres amortiguadas: c > 0, F = 0

d2 y dy
m+ c + ky = 0
dt2 dt
Las raíces de la ecuación característica son
√
−c ± c2 − 4mk
λ=
2m
Dependiendo del valor del discriminante c2 − 4mk se tienen tres casos
Caso 1 Sistema sobreamortiguado (sobreamortiguamiento): c2 − 4mk > 0.
Las dos raíces son reales y distintas
y(t) = c1 eλ1 t + c2 eλ2 t

Caso 2 Sistema críticamente amortiguado (caso crítico): c2 − 4mk = 0. Las

dos raíces son reales e iguales
y(t) = c1 eλ1 t + c2 teλ2 t

Caso 3 Sistema subamortiguado: c2 −4mk < 0. Las dos raíces son complejas
c   c
y(t) = e− 2m t c1 cos(µt) + c2 sin(µt) = Ae− 2m t cos(µt − α)
√
c2 − 4mk
µ=
2m
20
5. Modelos Físicos 5.5. Vibraciones mecánicas

5.5.3. Vibraciones forzadas con amortiguamiento: c >

0, F 6= 0
Considerando una fuerza externa F (t) = Fo cos(ωt)

d2 y dy
m 2
+ c + ky = Fo cos(ωt)
dt dt
Resolviendo la EDL No homogénea, la solución homogénea equivale al
obtenido en el apartado anterior. La solución particular será

yp = M cos(ωt) + N sin(ωt) = A cos(ωt − α)

donde  
k − mω 2 Fo
M =  2 2
2 2 2
k − mω +c ω
cωFo
N =  2 2
2 2 2
k − mω +c ω
y la amplitud de la respuesta en función de la frecuencia externa es
Fo
A(ω) = r 2
k − mω 2 + c2 ω 2

5.5.4. Vibración forzada sin amortiguamiento: c = 0, F 6=

0
Considerando una fuerza externa F (t) = Fo cos(ωt)

d2 y d2 y Fo
m 2
+ ky = F o cos(ωt) → 2
+ ωo2 y = cos(ωt)
dt dt m
donde ωo2 = k/m
Resolviendo la EDL No homogénea, la solución homogénea equivale al
obtenido en el apartado anterior. La solución particular dependerá de dos
casos

21
5. Modelos Físicos 5.5. Vibraciones mecánicas

Caso 1 ω 6= ωo

yp = M cos(ωt) + N sin(ωt) = A cos(ωt − α)

y la solución general será

Fo
y(t) = Ao cos(ωo t − α) +   cos(ωt)
m ωo2 − ω 2

Caso 2 ω = ωo  
yp = t M cos(ωt) + N sin(ωt)
y la solución general será
Fo
y(t) = Ao cos(ωo t − α) + t sin(ωo t)
2mωo
que reproduce el fenómeno de resonancia pura.

22