This is the entire problem. I am looking for help with part 3.(adsbygoogle = window.adsbygoogle || []).push({});

For the mass-spring system with damping, we obtain the following DE

my′′ + gammay′ + ky = 0, or y′′ + cy′ + omega^2y = 0, where c = gamma/m and omega^2 = k/m. Given m = 1,k = 4, gamma= 1, and initial conditions y(0) = 0.1, y′(0) = 0.

(1) Plot solutions y and v = y′ = dy/dt as functions of time.

(2) Plot v vs y (phase plot). Comment on the behavior of the curve.

(3) Plot the quantity E = 1/2mv^2 + 1/2ky^2 as a function of time. What do you observe?

(4) Show that dE/dt < 0 for gamma > 0 while dE/dt > 0 for gamma < 0 (analytically or by graphic).

For part 1

function Q2part1

m = 1; % mass [kg]

k = 4; % spring constant [N/m]

gamma = 1; % friction coefficient [Ns/m]

omega = sqrt(k/m); c = gamma/m;

y0 = 0.1; v0 = 0; % initial conditions

[t,Y] = ode45(@f,[0,10],[y0,v0],[],omega,c); % solve for 0<t<10

y = Y(:,1); v = Y(:,2); % retrieve y, v from Y

figure(1); plot(t,y,'b+-',t,v,'ro-'); % time series for y and v

%-------------------------------------------

function dYdt = f(t,Y,omega,c)

y = Y(1); v = Y(2);

dYdt = [ v ;-c*(v)-omega^2*y];

for part 2

function Q2part2

m = 1;

k = 4;

gamma = 1;

omega = sqrt(k/m); c = gamma/m;

y0 = 0.1; v0 = 0;

[t,Y] = ode45(@f,[0,10],[y0,v0],[],omega,c);

y = Y(:,1); v = Y(:,2);

figure(1); plot(v,y);

function dYdt = f(t,Y,omega,c)

y = Y(1); v = Y(2);

dYdt = [ v ;-c*(v)-omega^2*y];

I know that i need to plug in formulas for y (and v as well?) so that I have E(t), but is it as simple as solving for and plugging in y? I am in need of direction here I suppose...a nice big fat hint to put me in the correct direction.

**Physics Forums - The Fusion of Science and Community**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# MatLab Spring Damping and Energy

Loading...

Similar Threads for MatLab Spring Damping | Date |
---|---|

MATLAB: forces acting on an object | Apr 6, 2018 |

MATLAB Matlab histogram | Apr 4, 2018 |

MATLAB Sovling Bernoulli's differential equation in matlab? | Apr 3, 2018 |

MATLAB Matlab-like software | Mar 27, 2018 |

MathCAD/MatLab program for non-linear spring-mass-damper system | Mar 19, 2011 |

**Physics Forums - The Fusion of Science and Community**