Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 84 vistas

    PENDULO

    Cargado por

    ribert
    Este documento presenta un estudio de sistemas oscilatorios lineales mediante el método de Euler. Se analizan el péndulo simple, el movimiento armónico simple, el movimiento armónico amortiguado y el oscilador forzado. Los resultados obtenidos mediante simulaciones por computadora concuerdan con la teoría y muestran el comportamiento de las variables como la posición, velocidad y aceleración en función del tiempo para cada sistema.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    PENDULO

    Cargado por

    ribert
    84 vistas9 páginas
    Este documento presenta un estudio de sistemas oscilatorios lineales mediante el método de Euler. Se analizan el péndulo simple, el movimiento armónico simple, el movimiento armónico amortiguado y el oscilador forzado. Los resultados obtenidos mediante simulaciones por computadora concuerdan con la teoría y muestran el comportamiento de las variables como la posición, velocidad y aceleración en función del tiempo para cada sistema.

    Descripción original:

    FISICA

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta un estudio de sistemas oscilatorios lineales mediante el método de Euler. Se analizan el péndulo simple, el movimiento armónico simple, el movimiento armónico amortiguado y el oscilador forzado. Los resultados obtenidos mediante simulaciones por computadora concuerdan con la teoría y muestran el comportamiento de las variables como la posición, velocidad y aceleración en función del tiempo para cada sistema.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    84 vistas9 páginas

    PENDULO

    Cargado por

    ribert
    Este documento presenta un estudio de sistemas oscilatorios lineales mediante el método de Euler. Se analizan el péndulo simple, el movimiento armónico simple, el movimiento armónico amortiguado y el oscilador forzado. Los resultados obtenidos mediante simulaciones por computadora concuerdan con la teoría y muestran el comportamiento de las variables como la posición, velocidad y aceleración en función del tiempo para cada sistema.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    de 9

    Estudio de las soluciones para sistemas oscilatorios lineales

    mediante el método de Euler.


    Ribert Jhon Hancco Cuba1
    1
    Escuela Profesional de Física, Facultad de Ciencias Naturales y Formales, Universidad
    Nacional de San Agustín, Arequipa, Perú

    Abril 03, 2019

    En esta práctica se muestra la aplicación de las ecuaciones diferenciales ordinarias para


    predecir el movimiento de sistemas oscilatorios mediante métodos computacionales y
    explicar su comportamiento. Se estudió el dispositivo experimental del péndulo simple, el
    sistema masa resorte (MAS), el movimiento armónico amortiguado y las oscilaciones
    forzadas. Los resultados obtenidos concuerdan bastante bien, teniendo en cuenta los
    intervalos de error por redondeo o los errores globales y/o locales propias de los métodos
    computacionales.

    1. Introducción [1]
    El movimiento de Oscilación o vibración es una de las formas más comunes de movimiento
    de la naturaleza. Algunos ejemplos de esta situación se pueden apreciar cuando: una masa
    de péndulo es liberada desde una cierta altura, el extremo de un resorte se separa de su
    posición de equilibrio y luego se suelta, modelar sistemas cuánticos o predecir parámetros
    macroscópicos como el calor especifico de un sólido a parir de las vibraciones de la
    estructura cristalina, etc.

    Para la realización de esta experiencia se resolvió las ecuaciones diferenciales mediante el


    método de Euler que pertenece a las llamadas técnicas de Runge-Kutta. En esta primera
    aproximación se debe tomar en cuenta que los pasos para las iteraciones deben ser lo más
    cortos posibles para acercarnos a las predicciones teóricas, considerando que el tiempo de
    procesamiento también aumenta.

    2. Fundamento teórico [3].

    2.1 Movimiento armónico simple (MAS).


    La ecuación que gobierna el MAS está definido por:
    (1)

    Donde m es la masa del objeto k la constante del resorte, también podemos


    definir en función de la frecuencia angular del sistema (ω) y la aceleración
    lineal (α).
    (2)
    2.2 Péndulo simple.

    Para escribir la ecuación diferencial del péndulo se restringió el ángulo para


    dimensiones muy pequeñas, con la finalidad de obtener un movimiento armónico.

    (3)

    En términos de la aceleración angular:

    (4)

    Utilizando el método de Euler:


    (5)
    (6)

    2.3 Movimiento armónico amortiguado.

    La E.D. del sistema es:

    (7)

    Dónde: k es la constante del resorte y b constante de amortiguamiento. Se


    define de acuerdo a la ecuación diferencial tres tipos de soluciones según el
    valor de b2-4km.

     b2-4km>0 sobre amortiguado.


     b2-4km=0 críticamente amortiguado.
     b2-4km<0 subamortiguado.

    2.4 Oscilador forzado amortiguado.

    Para este caso, en la ecuación 7 hay que añadirle una fuerza periódica externa
    F=Fmcos(ωt).

    (8)

    Por lo que la ecuación diferencial se compone de la parte homogénea y la parte heterogenea


    que generara la solución particular del sistema.
    3. Implementación y Resultados:
    3.1 Péndulo Simple

    La Figura 1 muestra que la velocidad es máxima cuando el péndulo pasa por la posición de
    equilibrio además se observa que el sentido de la aceleración angular es opuesta a la
    velocidad angular y llega a su valor máximo cuando pasa por los extremos.
    Para el análisis se consideró los siguientes parámetros: g=9.8m/s2 l=1m y las condiciones
    iniciales: ω1=0, ϴ1=0.5. En el Anexo A1 se muestra el pseudocódigo para su
    implementación.

    Figura 1. Desplazamiento angular, velocidad y aceleración angular en función del tiempo.

    3.2 Movimiento armónico simple (mas).

    Análogamente al caso del péndulo, el comportamiento de la velocidad muestra un máximo


    cuando pasa por la posición de equilibrio, en tanto que la aceleración alcanza su máximo
    valor en los extremos del desplazamiento ver Figura 2. Las condiciones iniciales que se
    consideraron son: x1=1, v1=0 con m=1 y k=1 las cuales poseen unidades arbitrarias. En el
    Anexo A2 se muestra el pseudocódigo para su implementación en MatLab.
    Figura 2. Desplazamiento, velocidad y aceleración en función del tiempo.

    3.3 Movimiento Armónico Amortiguado

    3.3.1 Sobre amortiguado: (𝒃𝟐 > 𝟒𝒌𝒎)

    Para este caso consideramos que la constante de amortiguamiento b es relativamente


    grande, con una fuerte resistencia a la oscilación por lo que se deduce de la figura 3 que
    x(t)0 cuando t+∞ además se muestra la evolución del sistema para diferentes
    velocidades iniciales. En todos los casos las oscilaciones son eliminadas por el
    amortiguamiento.Se consideraron las siguientes condiciones: m=1kg, k=1N/m, b=3kg/s

    En el Anexo A3 se muestra el pseudocódigo para su implementación en MatLab.

    Figura 3. Grafica de la posición para diferentes velocidades en un sistema sobre amortiguado.


    3.3.2 Críticamente amortiguado(𝒃𝟐 = 𝟒𝒌𝒎)

    El comportamiento es similar al movimiento sobreamortiguado ver Figura 4, con la


    diferencia de que el amortiguamiento es lo suficientemente grande como para absorber las
    oscilaciones, pero incluso cualquier disminución en la resistencia nos llevaría al caso
    subamortiguado, además se muestra que el cuerpo pasa por el punto de equilibrio a lo más
    una vez.

    Se consideraron las siguientes condiciones:

    m=1kg, k=1N/m, b=2kg/s

    En el Anexo A4 se muestra el pseudocódigo para su implementación en MatLab.

    Figura 4. Grafica de la posición para diferentes velocidades en un sistema críticamente amortiguado

    3.3.3 Subamortiguado (𝒃𝟐 < 𝟒𝒌𝒎)

    Deacuerdo a la Figura 5 se observa que el cuerpo presenta oscilaciones “exponencialmente”


    amortiguadas alrededor de su posicion de equilibrio. Por lo tanto la accion de
    amortiguamiento tiene al menos tres efectos:

    1. Amortigua exponencialmente las oscilaciones, deacuerdo con la amplitud de


    tiempo variable.
    2. Retarda el movimiento; esto es el amortiguamiento retarda la frecuencia del
    movimiento.
    3. Retrasa el movimiento, deacuerdo al angulo de desfase.

    Se consideraron las siguientes condiciones: m=1kg, k=1N/m, b=0.2kg/s


    En el Anexo A5 se muestra el pseudocódigo para su implementación en MatLab.

    Figura 5: grafica de la posicion frente al tiempo en un sistema sub amortiguado.

    3.4 Oscilador Forzado

    La Figura 6 muestra el comportamiento del sistema para diferentes frecuencias de la fuerza


    externa pasando por la por la frecuencia de resonancia cuya amplitud de oscilación aumenta
    producto de un reforzamiento de las vibraciones naturales del sistema mediante un agente
    externo.Se consideraron las siguientes condiciones:

    m=1kg, k=1N/m, b=.1kg/s

    En el Anexo A6 se muestra el pseudocódigo para su implementación en MatLab.

    Figura 6: Grafica de la poscion para diferentes velocidades.


    4. Conclusiones.

    Se mostró que los métodos computacionales nos brinda una alternativa para
    resolver ecuaciones diferenciales y predecir la evolución del sistema.

    Se demostró que las condiciones iniciales del sistema son imprescindibles


    para las predicciones posteriores del sistema.

    5. Referencias.

    [1] ARDILA, William; CHAVES, José y CORTES, Jimy.Oscilaciones y Ondas en


    ciencias básicas para ingenierías.ECOE Ediciones, 2013. 219 páginas.
    [2] Chapra et al, Métodos numéricos para ingenieros, Mc Graw Hill, Quinta
    Edición, México,( 2006 )
    [3] Vizcardo Miguel, Curso: Fisica computacional 2

    ANEXO
    A1.PENDULO SIMPLE
    clear all
    l=1;
    th(1)=0.5;
    dt=0.0004;
    g=9.8;
    w(1)=0;
    t(1)=0;
    a(1)=-g/l*th(1);
    for i=1:15000
    th(i+1)=th(i)+w(i)*dt;
    w(i+1)=w(i)+a(i)*dt;
    a(i+1)=-g/l*th(i+1);
    t(i+1)=i;
    %fprintf('%f %f %f\n',th(i),w(i),t(i));
    end
    xx=plot(t,th,'r');
    hold on
    yy=plot(t,w,'g');
    zz=plot(t,a,'b');
    title(['\fontsize{14} Desplazamiento(\theta), Velocidad(\omega) y Aceleracion
    angular(\alpha)'])
    xlabel('\fontsize{16} Tiempo(s)')
    ylabel('\fontsize{16}(\theta, \omega, \alpha)')
    legend([xx,yy,zz],'\fontsize{16}(\theta)','\fontsize{16}(\omega)','(\fontsize{16}\alpha)')
    grid on

    A2. MOVIMIENTO ARMONICO SIMPLE

    clear all
    m=1;
    x(1)=1;
    dt=0.0004;
    k=1;
    v(1)=0;
    t(1)=0;
    a(1)=-k/m*x(1);
    for i=1:30000
    x(i+1)=x(i)+v(i)*dt;
    v(i+1)=v(i)+a(i)*dt;
    a(i+1)=-k/m*x(i+1);
    t(i+1)=i;
    %fprintf('%f %f %f\n',th(i),w(i),t(i));
    end
    %plot(t,th)
    %plot(t,w)
    xx=plot(t,x,'r');
    hold on
    yy=plot(t,v,'g');
    zz=plot(t,a,'b');
    title(['\fontsize{18} Desplazamiento(X), Velocidad(V) y Aceleracion(a)'])
    xlabel('\fontsize{16} Tiempo(s)')
    ylabel('\fontsize{16}(X, V, a)')
    legend([xx,yy,zz],'\fontsize{16}(X)','\fontsize{16}(V)','\fontsize{16}(a)')
    grid on

    A3. SOBRE AMORTIGUADO

    clear all
    m=1;
    v(1,1)=1;v(2,1)=0;v(3,1)=-1;v(4,1)=-2;v(5,1)=-3;
    dt=0.005;
    k=1;
    b=3;
    t(1)=0;
    color=[1 0 0;0 1 0;0 0 1;1 0 1;0 1 1];
    for i=1:5
    x(i,1)=1;
    a(i,1)=(-k*x(i,1)-b*v(i,1))/m;
    for j=1:2000
    x(i,j+1)=x(i,j)+v(i,j)*dt;
    v(i,j+1)=v(i,j)+a(i,j)*dt;
    a(i,j+1)=(-k*x(i,j+1)-b*v(i,j+1))/m;
    t(j+1)=j;
    end
    plot(t,x(i,:),'color',color(i,:));
    hold on
    end
    grid on

    title(['\fontsize{18} Posición(X) con diferentes velocidades iniciales(V)'])


    xlabel('\fontsize{16} Tiempo(s)')
    ylabel('\fontsize{16} Posición(m)')
    legend('\fontsize{16} V=1m/s','\fontsize{16} V=0m/s','\fontsize{16} V=-1m/s','\fontsize{16}
    V=-2m/s','\fontsize{16} V=-3m/s')
    grid on

    A4. CRÍTICAMENTE AMORTIGUADO

    clear all
    m=1;
    v(1,1)=1;v(2,1)=0;v(3,1)=-1;v(4,1)=-2;v(5,1)=-3;
    dt=0.005;
    k=1;
    b=2;
    t(1)=0;
    color=[1 0 0;0 1 0;0 0 1;1 0 1;0 1 1];
    for i=1:5
    x(i,1)=1;
    a(i,1)=(-k*x(i,1)-b*v(i,1))/m;
    for j=1:2000
    x(i,j+1)=x(i,j)+v(i,j)*dt;
    v(i,j+1)=v(i,j)+a(i,j)*dt;
    a(i,j+1)=(-k*x(i,j+1)-b*v(i,j+1))/m;
    t(j+1)=j;
    end
    plot(t,x(i,:),'color',color(i,:));
    hold on
    end
    grid on
    title(['\fontsize{18} Posición(X) con diferentes velocidades iniciales(V)'])
    xlabel('\fontsize{16} Tiempo(s)')
    ylabel('\fontsize{16} Posición(m)')
    legend('\fontsize{16} V=1m/s','\fontsize{16} V=0m/s','\fontsize{16} V=-1m/s','\fontsize{16}
    V=-2m/s','\fontsize{16} V=-3m/s')
    grid on

    A5. SUB AMORTIGUADO

    clear all
    m=1;
    v(1)=0;
    dt=0.005;
    k=1;
    b=.2;
    t(1)=0;
    x(1)=1;
    a(1)=(-k*x(1)-b*v(1))/m;
    for j=1:6000
    x(j+1)=x(j)+v(j)*dt;
    v(j+1)=v(j)+a(j)*dt;
    a(j+1)=(-k*x(j+1)-b*v(j+1))/m;
    t(j+1)=j;
    end
    plot(t,x);
    hold on
    grid on
    title(['\fontsize{18} Posición(X)'])
    xlabel('\fontsize{16} Tiempo(s)')
    ylabel('\fontsize{16} Posición(m)')
    legend('\fontsize{16} Velocidad inicial=0m/s')
    grid on

    A6. OSCILADOR FORZADO

    clear all
    m=1;
    w(1,1)=2;w(2,1)=1;w(3,1)=0.5;
    dt=0.005;
    k=1;
    b=.1;
    t(1)=0;
    color=[1 0 0;0 1 0;0 0 1;1 0 1;0 1 1];
    for i=1:3
    x(i,1)=1;
    f(i,1)=cos(w(i,1)*dt*t(1));
    v(i,1)=0;
    a(i,1)=(-k*x(i,1)-b*v(i,1)+f(i,1))/m;
    for j=1:6000
    f(i,j+1)=cos(w(i,1)*dt*(j+1));
    x(i,j+1)=x(i,j)+v(i,j)*dt;
    v(i,j+1)=v(i,j)+a(i,j)*dt;
    a(i,j+1)=(-k*x(i,j+1)-b*v(i,j+1)+f(i,j+1))/m;
    t(j+1)=j;
    end
    r(i)=plot(t,x(i,:),'color',color(i,:));
    hold on
    end
    grid on

    title(['\fontsize{18} Posición(X) con diferentes frecuencias'])


    xlabel('\fontsize{16} Tiempo(s)')
    ylabel('\fontsize{16} Posición(m)')
    legend([r(1),r(2),r(3)],'\fontsize{16}w^2>k/m','\fontsize{16}
    w^2=k/m(Resonancia)','\fontsize{16} w^2<k/m')
    grid on

    También podría gustarte