Scribd Logo
Upload

Welcome to Scribd!

  • 14 views

    Clear M 1 K 100 F 100 C 2 W 1:0.1:20 Sim V Simout (40:end) A Max (V) Plot (W, A,) Drawnow Hold

    Uploaded by

    irqovi
    This document contains MATLAB code that simulates and plots the response of a harmonic oscillator with and without an absorber. Without an absorber, the code simulates and plots the maximum response of the oscillator versus driving frequency. With an absorber, the code defines parameters for a two degree-of-freedom system with a primary oscillator and absorber, integrates the equations of motion, extracts and plots the maximum response versus driving frequency. A function defines the equations of motion for the two degree-of-freedom system with an absorber.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Clear M 1 K 100 F 100 C 2 W 1:0.1:20 Sim V Simout (40:end) A Max (V) Plot (W, A,) Drawnow Hold

    Uploaded by

    irqovi
    14 views3 pages
    This document contains MATLAB code that simulates and plots the response of a harmonic oscillator with and without an absorber. Without an absorber, the code simulates and plots the maximum response of the oscillator versus driving frequency. With an absorber, the code defines parameters for a two degree-of-freedom system with a primary oscillator and absorber, integrates the equations of motion, extracts and plots the maximum response versus driving frequency. A function defines the equations of motion for the two degree-of-freedom system with an absorber.

    Original Description:

    ggtg

    Original Title

    Absorber

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document contains MATLAB code that simulates and plots the response of a harmonic oscillator with and without an absorber. Without an absorber, the code simulates and plots the maximum response of the oscillator versus driving frequency. With an absorber, the code defines parameters for a two degree-of-freedom system with a primary oscillator and absorber, integrates the equations of motion, extracts and plots the maximum response versus driving frequency. A function defines the equations of motion for the two degree-of-freedom system with an absorber.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    14 views3 pages

    Clear M 1 K 100 F 100 C 2 W 1:0.1:20 Sim V Simout (40:end) A Max (V) Plot (W, A,) Drawnow Hold

    Uploaded by

    irqovi
    This document contains MATLAB code that simulates and plots the response of a harmonic oscillator with and without an absorber. Without an absorber, the code simulates and plots the maximum response of the oscillator versus driving frequency. With an absorber, the code defines parameters for a two degree-of-freedom system with a primary oscillator and absorber, integrates the equations of motion, extracts and plots the maximum response versus driving frequency. A function defines the equations of motion for the two degree-of-freedom system with an absorber.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    You are on page 1of 3

    clear

    m=1;k=100;F=100;c=2;
    for w=1:0.1:20
    sim('second_order_harmonic_response');
    v=simout(40:end);
    A=max(v);
    plot(w,A,'b--o')
    drawnow
    hold on
    end
    clear
    global f m1 w k1 c2 k2 m2
    m1=200;m2=10;f=100;k1=20000;k2=2000;c2=0;
    tspan=[0 20];
    y0=[0;0;0;0];
    for w=1:0.01:20
    [t,y]=ode45('absorber',tspan,y0);
    x=y(:,1);
    v=x(30:end);
    A=max(v);
    plot(w,A,'b--o')
    drawnow
    hold on
    end

    function yp = absorber(t,y)
    global f m1 w k1 k2 m2 c2
    yp=zeros(4,1);
    yp(1)=y(2);
    yp(2)=(((f/m1)*sin(w*t))-((k1/m1)*y(1))-((k2/m1)*(y(1)-
    y(3)))-((c2/m1)*(y(2)-y(4))));
    yp(3)=y(4);
    yp(4)=(-((k2/m2)*(y(3)-y(1)))-((c2/m2)*(y(4)-y(2))));

    You might also like