solve a base excited spring damper system with sinusoidal excitation.

27 views (last 30 days)
I'm modelling a spring damper system governed by the equation of motion m*x''+c*x'+kx=ky+cy''.
y is the input y=Y0sin(wt),
x is the output (motion of the mass)
k=spring constant
c=damping coefficient
m=mass
I need velocity and displacement profiles for the mass as outputs.
clear
t=0:0.1:15; %time peroid
Y0= input('wave amplitude ') ; %Wave amplitude
l= input('length of wave ') ; %length of the wave
u=input('speed of boat'); %Boat Velocity
w=u/l; %frequency of wave
y=Y0*sin(w*t); %wave height model
c=2*Dr*(k*m)^(0.5); %Damping coefficient
Dr=0.5; %Required Damping Ratio
k=17000; %Spring Constant of suspension
m=100;
wn=(k/m)^0.5;
wd=wn*(1-Dr^2)^0.5;
syms x(t)
ode=m*diff(x,t,2)+c*diff(x,t)+k*x==y+(2*Dr/wn)*diff(y,t);
xSol(t)=dsolve(ode);

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!