Scribd Logo
Upload

Welcome to Scribd!

  • 87 views

    Exercise 8 SOLUTION

    Uploaded by

    Leng Chai
    This document is an exercise problem from a mechanical vibrations course. It provides instructions to numerically solve for the displacement of an oscillator with damping over time based on given parameters, and to plot the results. The student included their MATLAB code and output showing the time-varying displacement, velocity, acceleration, and applied force values from solving the problem.

    Copyright:

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

    Exercise 8 SOLUTION

    Uploaded by

    Leng Chai
    87 views6 pages
    This document is an exercise problem from a mechanical vibrations course. It provides instructions to numerically solve for the displacement of an oscillator with damping over time based on given parameters, and to plot the results. The student included their MATLAB code and output showing the time-varying displacement, velocity, acceleration, and applied force values from solving the problem.

    Original Description:

    MECHANICAL VIBRATION EXTERNAL FORCE HARMONIC VIBRATION

    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 is an exercise problem from a mechanical vibrations course. It provides instructions to numerically solve for the displacement of an oscillator with damping over time based on given parameters, and to plot the results. The student included their MATLAB code and output showing the time-varying displacement, velocity, acceleration, and applied force values from solving the problem.

    Copyright:

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

    Exercise 8 SOLUTION

    Uploaded by

    Leng Chai
    This document is an exercise problem from a mechanical vibrations course. It provides instructions to numerically solve for the displacement of an oscillator with damping over time based on given parameters, and to plot the results. The student included their MATLAB code and output showing the time-varying displacement, velocity, acceleration, and applied force values from solving the problem.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 6

    BYG550 MECHANICAL VIBRATIONS

    COORDINATOR ASSC. PROF. TERJE LEENHART ANDERSEN

    EXERCISE 8

    00 IZWAN BIN AHMAD 261972 24 October, 2020

    REV. NO NAME STUDENT ID DATE


    PROBLEM 1

    2
    3
    PROBLEM 2

    For the problem presented above, determine the displacement x(t) numerically, by the constant
    average acceleration method. Use the following data: F0=10kN; t0=0.4s; Tn=1.6s ; k=500kN/m. Plot
    x(t ) as function of time for t values from t=0 to t=10s. (See the recipe, Eq. (14) in lecture note #15.
    Use time-step h=0.01s. Use MATLAB, Octave, Scilab, Matchad, Excel or similar).

    Answer:

    The solution was done using MATLAB.


    NumericalSolution.m

    clear all;
    close all;
    clc;

    Tn = 1.6; % second
    wn = 2*pi/Tn; % Natural Eigen Frequency
    k = 500; % kN/m
    m = k/(wn^2); % 1e6 kg
    zeta = 0.0; % relative damping
    c = 2*m*wn*zeta;

    t(1) = 0; % time is zero at time step index 1


    x(1) = 0; % initial displacement is zero
    v(1) = 0; % initial velocity is zero
    h = 0.01;% time-step in seconds
    F(1) = InverselyLinearForce(t(1));
    a(1) = F(1)/m; % initial acceleration

    for I = 2:1001
    t(i) = t(i-1)+h;
    F(i) = InverselyLinearForce(t(i));
    x(i) = (F(i) + m*a(i-1) + (c+4*m/h)*v(i-1) + (4*m/h^2+2*c/h)*x(i-1))/(4*m/h^2 + 2*c/h + k);
    v(i) = 2/h*x(i) - 2/h*x(i-1) - v(i-1);
    a(i) = 4/h^2*x(i) - 4/h^2*x(i-1) - 4/h*v(i-1) - a(i-1);
    end

    disp('time [s] | Force [kN] | x [m] | v [m/s] | a [m/s/s]')


    [t' F' x' v' a']

    figure(1); clf;
    plot(t,x,'ok-'); ylabel('x [m]')
    title(['Solution with damping: \zeta = ',num2str(zeta)])
    grid on; grid minor;

    figure(2); clf;
    subplot(311); plot(t,F,'xr-'); ylabel('F [kN]')
    title(['Solution with damping: \zeta = ',num2str(zeta)])
    grid on; grid minor;

    subplot(312); plot(t,v,'*m-'); ylabel('v [m/s]')


    grid on; grid minor;

    subplot(313); plot(t,a,'db-'); ylabel('a [m/s/s]'); xlabel('time, t [s]');


    grid on; grid minor;

    InverselyLinearForce.m

    function F = InverselyLinearForce (t)


    ta = 0;
    tb = 0.4;
    to = 0.4;
    F = 10*(1-(t/to))*(ta<=t).*(t<=tb);
    end

    4
    5
    𝑇𝑖𝑚𝑒 𝐹𝑜𝑟𝑐𝑒 𝐷𝑖𝑠𝑝𝑙𝑎𝑐𝑒𝑚𝑒𝑛𝑡 𝑉𝑒𝑙𝑜𝑐𝑖𝑡𝑦 𝐴𝑐𝑐𝑒𝑙𝑒𝑟𝑎𝑡𝑖𝑜𝑛
    (s) (kN) (m) (m/s) (m/s²)
    0.0000 10.0000 0.0000 0.0000 0.3084
    0.0100 9.7500 0.0000 0.0030 0.3005
    0.0200 9.5000 0.0001 0.0060 0.2921
    0.0300 9.2500 0.0001 0.0089 0.2832
    0.0400 9.0000 0.0002 0.0117 0.2739
    0.0500 8.7500 0.0004 0.0144 0.2642
    0.0600 8.5000 0.0005 0.0170 0.2541
    0.0700 8.2500 0.0007 0.0194 0.2436
    0.0800 8.0000 0.0009 0.0218 0.2327
    0.0900 7.7500 0.0011 0.0241 0.2214
    0.1000 7.5000 0.0014 0.0262 0.2098
    0.1100 7.2500 0.0017 0.0283 0.1979
    0.1200 7.0000 0.0020 0.0302 0.1857
    0.1300 6.7500 0.0023 0.0320 0.1732
    0.1400 6.5000 0.0026 0.0337 0.1604
    0.1500 6.2500 0.0029 0.0352 0.1474
    0.1600 6.0000 0.0033 0.0366 0.1341
    0.1700 5.7500 0.0037 0.0379 0.1207
    0.1800 5.5000 0.0041 0.0390 0.1070
    0.1900 5.2500 0.0045 0.0400 0.0932
    0.2000 5.0000 0.0049 0.0409 0.0793
    0.2100 4.7500 0.0053 0.0416 0.0652
    0.2200 4.5000 0.0057 0.0422 0.0510
    0.2300 4.2500 0.0061 0.0426 0.0368
    0.2400 4.0000 0.0065 0.0429 0.0225
    0.2500 3.7500 0.0070 0.0431 0.0081
    0.2600 3.5000 0.0074 0.0431 -0.0062
    0.2700 3.2500 0.0078 0.0430 -0.0206
    0.2800 3.0000 0.0083 0.0427 -0.0349
    0.2900 2.7500 0.0087 0.0423 -0.0491
    0.3000 2.5000 0.0091 0.0417 -0.0633
    0.3100 2.2500 0.0095 0.0410 -0.0774
    0.3200 2.0000 0.0099 0.0402 -0.0914
    0.3300 1.7500 0.0103 0.0392 -0.1052
    0.3400 1.5000 0.0107 0.0380 -0.1189
    0.3500 1.2500 0.0111 0.0368 -0.1323
    0.3600 1.0000 0.0114 0.0354 -0.1456
    0.3700 0.7500 0.0118 0.0339 -0.1587
    0.3800 0.5000 0.0121 0.0322 -0.1715
    0.3900 0.2500 0.0124 0.0305 -0.1840
    0.4000 0.0000 0.0127 0.0285 -0.1963

    You might also like