La PizarraMatemática.

Matemáticas para ingeniería

Modelos de oscilaciones y vibraciones

Una masa unida a un resorte puede volver hacia su equilibrio, pasarse y regresar. Antes de resolver la ecuación, dibujaremos qué fuerza empuja y en qué sentido. Después veremos cómo posición y velocidad iniciales fijan la oscilación, qué cambia al amortiguar y por qué una fuerza periódica necesita distinguir respuesta permanente y arranque.

1. Elegimos el equilibrio y el sentido antes de poner signos a las fuerzas

Medimos x desde el equilibrio estático y elegimos sentido positivo. En el modelo lineal, el resorte ejerce −kx, el amortiguador viscoso −cx′ y la fuerza exterior es F(t).

Con m, k>0 y c≥0, Newton da mx″=F−cx′−kx, o mx″+cx′+kx=F.
Las fuerzas se escriben después de elegir coordenada y sentido
Medimos x desde el equilibrioCaso dibujado: x>0, x′>0mResorte: −kxAmortiguador: −cx′m x″ = F − cx′ − kxDos datos: x(0) y x′(0)El peso constante se absorbesi usamos el equilibrio estático.
La masa está conectada a resorte y amortiguador. El caso ilustrado tiene desplazamiento y velocidad positivos; ambas fuerzas internas apuntan en sentido negativo. La fuerza exterior se representa por F en la ecuación. El desplazamiento se mide desde el equilibrio estático, lo que permite absorber un peso constante si el modelo y la geometría lo justifican.

La ecuación necesita dos datos: x(0) y x′(0). Un peso constante puede desaparecer al desplazar la coordenada al equilibrio; no porque la gravedad haya dejado de actuar. Un resorte no lineal, fricción seca o una geometría diferente exigen revisar las fuerzas.

Calcular la aceleración con signos antes de buscar toda la trayectoria

m=0,5 kg, k=50 N/m, c=3 kg/s; x=0,1 m, x′=0,2 m/s y F=1 N.

  1. Fuerzas internas

    El resorte aporta −5 N y el amortiguador −0,6 N, porque desplazamiento y velocidad son positivos.

  2. Neto

    F total=1−5−0,6=−4,6 N.
  3. Aceleración

    x″=−4,6/0,5=−9,2 m/s² en ese instante.

    Es una cuenta local; no fija por sí sola toda x(t).

Resultado

El resorte se opone al desplazamiento y el amortiguador a la velocidad.

Lo comprobamos

Si la velocidad cambiara de signo, también cambiaría el signo de −cx′.

Explicar por qué el peso puede quedar absorbido

Resorte vertical con extensión y positiva hacia abajo: m y″+c y′+ky=mg+F(t).

  1. Equilibrio estático

    Sin movimiento ni fuerza dinámica, y equilibrio=mg/k.
    Para m=0,5 kg, k=50 N/m y g=9,81 m/s², vale 0,0981 m.
  2. Nueva variable

    Tomar x=y−mg/k conserva derivadas: y′=x′ e y″=x″.
  3. Sustituir

    El término k·mg/k cancela mg y queda mx″+cx′+kx=F(t). x=0 indica equilibrio estático, no resorte sin estirar.

Resultado

Cambiar la referencia explica el término que ya no aparece.

Lo comprobamos

La extensión y y la desviación x no tienen el mismo cero físico.

2. Frecuencia angular, frecuencia en ciclos y amortiguamiento tienen unidades diferentes

Los términos mx″, cx′ y kx tienen dimensión de fuerza.

Por eso [c]=M/T y [k]=M/T².

Definimos ω₀=√(k/m), frecuencia angular natural, y ζ=c/(2√(mk)), razón adimensional de amortiguamiento.

En el caso libre sin amortiguar, el periodo es 2π/ω₀ y la frecuencia de ciclos es f₀=ω₀/(2π). Los radianes se tratan como adimensionales, pero se conserva el nombre frecuencia angular para recordar el factor 2π. En el caso amortiguado la frecuencia de oscilación puede cambiar o dejar de existir.

Pasar de masa y rigidez a una vuelta de oscilación

m=2 kg y k=18 N/m.

  1. Angular

    ω₀=√(18/2)=3 rad/s.
  2. Periodo

    T=2π/3 s.
  3. Ciclos

    f₀=3/(2π) Hz.

    No se escriben tres hercios: tres era la frecuencia angular.

Resultado

Una vuelta necesita un cambio de fase 2π.

Lo comprobamos

Si se cuadruplicara m conservando k, ω₀ se dividiría por dos y el periodo se duplicaría.

La razón de amortiguamiento compara coeficientes con la misma dimensión

m=1, k=9 y c=2 en unidades coherentes.

  1. Natural

    ω₀=3.
  2. Razón

    ζ=2/(2√9)=1/3.

    Es un número sin unidades.

  3. Otros parámetros

    La tasa de la envolvente es c/(2m)=1 por unidad de tiempo y la frecuencia amortiguada es √(9−1)=√8 rad por unidad de tiempo.

Resultado

ζ no es una frecuencia ni el coeficiente c con otro nombre.

Lo comprobamos

c=6 daría ζ=1 para esta misma masa y rigidez.

3. La oscilación libre usa dos constantes que salen de los dos datos iniciales

Con c=0 y F=0, x″+ω₀²x=0 y x=Acos(ω₀t)+Bsin(ω₀t).
En t=0 se obtiene A=x(0); al derivar, x′(0)=Bω₀ y B=x′(0)/ω₀.

No se usa B igual a la velocidad sin ese divisor.

También se puede escribir Rcos(ω₀t+φ), con R=√(A²+B²), A=Rcosφ y B=−Rsinφ. La fase se elige con el cuadrante correcto. Es otra representación de los mismos dos datos, no una solución adicional.

Una velocidad inicial puede cambiar amplitud y fase

m=2, k=18; x(0)=1 y x′(0)=−3.

  1. Constantes

    ω₀=3, A=1 y B=−3/3=−1.
  2. Solución

    x(t)=cos3t−sin3t.
  3. Amplitud y fase

    R=√2 y φ=π/4, así x=√2cos(3t+π/4).

    La amplitud supera la posición inicial uno porque había también velocidad.

Resultado

Los dos datos participan en la oscilación.

Lo comprobamos

En t=0 la expresión da uno y su derivada da menos tres.

Partir del equilibrio no significa quedarse quieto

m=1, k=4; x(0)=0 y x′(0)=2.

  1. Constantes

    ω₀=2, A=0 y B=2/2=1.
  2. Solución

    x(t)=sin2t.
  3. Un punto

    En t=π/4, x=1 y x′=2cos(π/2)=0: la masa ha llegado a un extremo y cambia de sentido después.

Resultado

Se puede iniciar la oscilación solo con velocidad.

Lo comprobamos

La solución identicamente cero habría incumplido el dato de velocidad dos.

4. Las raíces de la ecuación libre distinguen tres regímenes

Buscar x=exp(rt) da mr²+cr+k=0.

Si 0≤ζ<1, aparecen raíces −c/(2m)±iωd y una oscilación con envolvente exponencial, donde ωd=ω₀√(1−ζ²). Si ζ=1, la raíz es repetida; si ζ>1, hay dos raíces reales negativas.

Amortiguar no significa siempre mantener una oscilación
m=1, k=9; x(0)=1, x′(0)=0Misma masa y resorte; cambia c1001234Azul c=2: oscila y decaeVerde c=6: críticoRojo c=10: sobreamortiguadoTiempo t; desplazamiento x.
Se comparan tres respuestas libres con masa uno, rigidez nueve, posición inicial uno y velocidad inicial cero. Con c dos, la solución oscila con envolvente exp de menos t. Con c seis es crítica y vale uno más tres t por exp de menos tres t. Con c diez combina nueve exp de menos t menos exp de menos nueve t, dividido entre ocho, y no oscila.

En el caso crítico la forma es (A+Bt)exp(−ω₀t). En el sobreamortiguado se combinan dos exponenciales reales. No se llama frecuencia de oscilación real a una raíz cuadrada negativa: esos dos regímenes no tienen una oscilación libre sinusoidal.

Un caso subamortiguado conserva oscilación y pierde envolvente

m=1, k=9, c=2; x(0)=1, x′(0)=0.

  1. Raíces

    r²+2r+9=0 da r=−1±i√8.
  2. Constantes

    x=exp(−t)[Ccos(√8t)+Dsin(√8t)].
    C=1 y x′(0)=−C+√8D=0, así D=1/√8.
  3. Resultado

    x=exp(−t)[cos(√8t)+sin(√8t)/√8].

    La fase interior oscila con √8 y la envolvente decae con tasa uno.

Resultado

La frecuencia amortiguada es menor que la natural tres.

Lo comprobamos

La amplitud instantánea no se identifica simplemente con exp(−t) ignorando las constantes de seno y coseno.

Crítico y sobreamortiguado dan formas distintas con los mismos datos

Mantener m=1, k=9, x(0)=1 y x′(0)=0.

  1. Crítico c=6

    La raíz doble es −3.

    De A=1 y B−3A=0 resulta x=(1+3t)exp(−3t).
  2. Sobreamortiguado c=10

    Las raíces son −1 y −9.

    Las ecuaciones A+B=1 y −A−9B=0 dan A=9/8 y B=−1/8.
  3. Comparar

    La segunda solución es [9exp(−t)−exp(−9t)]/8. Ambas regresan al equilibrio sin una sinusoidal libre; su evolución depende de estas raíces y datos, no solo de decir más amortiguamiento.

Resultado

La raíz repetida necesita el factor t; no se repite dos veces la misma exponencial.

Lo comprobamos

No se generaliza una comparación de rapidez entre estos datos a toda posible condición inicial o criterio de rendimiento.

5. El balance de energía separa lo que almacena el resorte, el movimiento y las pérdidas

La energía mecánica es E=mx′²/2+kx²/2.
Al derivar, E′=x′(mx″+kx)=F(t)x′−cx′².

La fuerza exterior aporta potencia F x′ y el amortiguador retira c x′² cuando c≥0.

Si F=0, E no aumenta.

Con c=0 es constante; con c>0 puede tener derivada cero en un instante donde x′=0 sin ser constante durante toda la trayectoria. El signo de la potencia exterior depende también de la velocidad, no solo del signo de F.

Posición y velocidad pueden repartirse la misma energía total

m=2, k=18, x=1 y x′=−3, sin fuerza ni amortiguamiento.

  1. Cinética

    mx′²/2=2·9/2=9.
  2. Potencial

    kx²/2=18·1/2=9.

    La energía total es 18.

  3. Extremos

    Cuando la velocidad es cero, |x|=√(2E/k)=√2.
    Cuando x=0, |x′|=√(2E/m)=3√2.

    La energía se intercambia sin cambiar su suma.

Resultado

La amplitud √2 coincide con la solución libre obtenida desde esos datos.

Lo comprobamos

El signo negativo inicial de velocidad no hace negativa su energía cinética.

La fuerza puede aportar o retirar energía aun con el mismo valor

F=5 N y c=2 kg/s.

  1. Velocidad positiva

    Con x′=2 m/s, el trabajo instantáneo es 10 W y la disipación 8 W.
    E′=2 W.
  2. Velocidad negativa

    Con x′=−2 m/s, la potencia de F es −10 W y la disipación sigue siendo 8 W.
    E′=−18 W.
  3. Instante de reposo

    Si x′=0, E′=0 en ese instante.

    El resorte puede acelerar la masa después, así ese dato aislado no demuestra ausencia de disipación futura.

Resultado

La tasa de energía requiere fuerza y velocidad en el mismo instante.

Lo comprobamos

Para calcular el valor de E también harían falta posición y coeficientes; aquí se ha calculado su derivada.

6. La respuesta forzada completa incluye un transitorio y una parte armónica

Con F=F₀cos(Ωt), F₀≥0 y c>0, una particular armónica es A cos(Ωt−δ).

A=F₀/√[(k−mΩ²)²+(cΩ)²] y δ=atan2(cΩ, k−mΩ²).

El cuadrante importa: para k−mΩ² negativo no basta la arctangente de un cociente sin contexto.

La solución completa añade una solución libre que ajusta los dos datos iniciales. En el modelo con c>0 y m, k>0 ese transitorio decae, quedando la particular. Ω es la frecuencia impuesta por la fuerza, no la natural ω₀ salvo que se indique esa igualdad.

Obtener amplitud y fase de una fuerza de frecuencia distinta de la natural

m=1, k=9, c=2 y F=5cos2t.

  1. Denominador

    k−mΩ²=5 y cΩ=4.
    Así A=5/√41 y δ=arctan(4/5), en el primer cuadrante.
  2. Forma por componentes

    x particular=(25cos2t+20sin2t)/41.

    Corresponde a A cos(2t−δ).

  3. Comprobar

    Sustituir en x″+2x′+9x da 5cos2t. En t=0 la particular vale 25/41 y su velocidad es 40/41: no cumple por sí sola reposo inicial.

Resultado

Una particular resuelve la ecuación, pero todavía no necesariamente sus datos iniciales.

Lo comprobamos

Si Ω=4, δ debe estar en el segundo cuadrante porque k−mΩ²=−7 y cΩ=8.

Ajustar el arranque sin cambiar la fuerza periódica

El mismo sistema, ahora con x(0)=x′(0)=0.

  1. Añadir libre

    Usar x=x particular+exp(−t)[Ccos(√8t)+Dsin(√8t)].
  2. Posición

    25/41+C=0, luego C=−25/41.
  3. Velocidad

    40/41−C+√8D=0 da D=−65/(41√8).

    Así quedan satisfechos ambos datos; la parte añadida desaparece al crecer t.

Resultado

La trayectoria de arranque y la respuesta permanente son dos partes de una misma solución.

Lo comprobamos

Eliminar el transitorio desde el principio habría incumplido las condiciones iniciales.

7. Un máximo de amplitud y una resonancia ideal no son el mismo caso

Para desplazamiento forzado con c>0, si 0≤ζ<1/√2 el máximo a frecuencia positiva ocurre en Ω pico=ω₀√(1−2ζ²). Para ζ≥1/√2 la amplitud no presenta ese máximo interior positivo y su mayor valor está en Ω=0. Esto se refiere a la respuesta de desplazamiento, no a cualquier variable de salida.

Sin amortiguamiento y con Ω=ω₀, el denominador armónico se anula y no existe esa particular acotada. La solución incluye un término proporcional a t sin(ω₀t). Su envolvente crece dentro del modelo ideal; no se usa una fórmula estacionaria que ya ha perdido sus hipótesis.

El máximo amortiguado no tiene que estar en la frecuencia natural

m=1, k=9, c=2 y amplitud de fuerza cinco.

  1. Minimizar denominador

    El cuadrado del denominador es (9−Ω²)²+4Ω²=Ω⁴−14Ω²+81.

    Su derivada es 4Ω(Ω²−7).

  2. Pico

    El mínimo positivo está en Ω=√7, donde el cuadrado vale 32 y A pico=5/√32.

  3. Comparar

    La frecuencia natural es tres, donde A=5/6.

    Esa amplitud es menor que 5/√32; natural y pico no son idénticos aquí.

Resultado

La localización del pico sale de la función de amplitud del sistema.

Lo comprobamos

ζ=1/3 cumple el requisito para ese máximo positivo.

En resonancia no amortiguada aparece crecimiento temporal

x″+9x=3cos3t, con x(0)=x′(0)=0.

  1. Coincidencia

    La fuerza tiene Ω=3 igual a ω₀, y c=0.
  2. Probar otra forma

    x=(t/2)sin3t tiene x″=3cos3t−(9t/2)sin3t.

    Al sumar 9x queda exactamente la fuerza.

  3. Datos y tamaño

    La posición y velocidad iniciales son cero. La envolvente crece como t/2, de modo que no hay una amplitud permanente acotada de la forma anterior.

Resultado

La fórmula armónica singular se sustituye por una solución que sí satisface la ecuación.

Lo comprobamos

El crecimiento se interpreta dentro de las hipótesis de oscilador lineal sin disipación.

8. La comprobación vuelve a la ecuación y a los datos; varios grados de libertad añaden modos

Tras resolver, calculamos x′ y x″, sustituimos en la ecuación original y comprobamos los dos datos. El balance de energía aporta otra verificación. Una fórmula que satisface la ecuación puede corresponder a otra trayectoria si falla un dato inicial.

Con varias masas aparecen matrices de masa, amortiguamiento y rigidez. En problemas sin amortiguamiento, masa SPD y rigidez simétrica, los modos se obtienen de Kφ=ω²Mφ y permiten coordenadas adecuadas. Un amortiguamiento general no queda necesariamente desacoplado por esos mismos modos; se revisa esa hipótesis.

Verificar una solución libre en las tres condiciones necesarias

m=2, k=18 y x=cos3t−sin3t.

  1. Derivadas

    x′=−3sin3t−3cos3t y x″=−9cos3t+9sin3t.
  2. Ecuación

    2x″+18x=0 en todo tiempo.
  3. Datos

    x(0)=1 y x′(0)=−3.

    La fórmula resuelve exactamente el problema inicial indicado; otros datos exigirían otras constantes.

Resultado

La ecuación y sus datos se comprueban por separado.

Lo comprobamos

La energía inicial 18 coincide con el balance conservativo calculado antes.

Dos modos explican dos masas acopladas sin amortiguamiento

M=I, K=[[2,−1],[−1,2]], posición inicial (1,0) y velocidad inicial cero.

  1. Modos

    K(1,1)ᵀ=(1,1)ᵀ y K(1,−1)ᵀ=3(1,−1)ᵀ.

    Las frecuencias son uno y √3.

  2. Descomponer

    (1,0)=(1/2)(1,1)+(1/2)(1,−1).

    Cada modo recibe velocidad inicial cero.

  3. Trayectoria

    x₁=(cos t+cos(√3t))/2 y x₂=(cos t−cos(√3t))/2.

    Al sustituir se cumplen ambas ecuaciones; la energía inicial total es uno.

Resultado

Las formas de vibrar colectivamente pueden simplificar un sistema acoplado.

Lo comprobamos

No se trata cada masa como un oscilador independiente ignorando los términos fuera de diagonal de K.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Una oscilación libre

m=1,k=4,c=0,F=0, x(0)=3,x′(0)=0. Halla x(t).

  1. Ecuación

    Con m=1, k=4, c=0 y F=0 queda x″+4x=0.
  2. Frecuencia

    ω₀=√4=2.
  3. Forma general

    x=Acos2t+Bsin2t.
  4. Posición

    x(0)=A=3.
  5. Velocidad

    x′(0)=2B=0, luego B=0 y x(t)=3cos2t.
  6. Control

    x″=−12cos2t y 4x=12cos2t cancelan.

    También se cumplen posición tres y velocidad cero al inicio.

Resultado

x(t)=3cos2t.

Lo comprobamos

x″+4x=−12cos2t+12cos2t=0.

EJEMPLO RESUELTO 2

Energía inicial

m=2,k=8, x=1,x′=3. Calcula E.

  1. Datos

    m=2, k=8, posición x=1 y velocidad x′=3, en unidades coherentes.
  2. Fórmula

    E=mx′²/2+kx²/2.
  3. Movimiento

    La parte cinética es 2·3²/2=9.
  4. Resorte

    La parte potencial es 8·1²/2=4.
  5. Total

    E=9+4=13 unidades de energía.
  6. Control

    Ambas contribuciones son no negativas. Sin especificar el sistema de unidades no se inventa una etiqueta física adicional; con SI coherente serían julios.

Resultado

E=13 unidades de energía.

Lo comprobamos

Ambos términos son no negativos con m,k positivos.

Los símbolos que hemos usado

m>0
Masa: coeficiente de la aceleración.
k>0
Rigidez del resorte lineal.
c≥0
Coeficiente de amortiguamiento viscoso pasivo.
ω₀=√(k/m)
Frecuencia angular natural sin amortiguamiento.
ζ=c/(2√(mk))
Razón de amortiguamiento, adimensional.
E
Suma de energía cinética y energía potencial elástica.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar