%zadanie 4 clc; clear; clf; %DANE E=200e9;%N/m2 ro=7850;%kg/m3 L=5;%m b=0.1;%m h=0.2;%m ksi=0.03; u0=0;%m v0=0;%m/s %WYKRES OBCI¥¯ENIA ts=0.04;%s tk=0.08;%s p0=2000;%N x1=0:0.001:ts; x2=ts:0.001:tk; x3=tk:0.001:1; y1=50000*x1; y2=4000-50000*x2; y3=0*x3; plot(x1,y1,'b',x2,y2,'b',x3,y3,'b') grid on; title('Wykres obci¹¿eñ'); ylabel('p(t) [N]'); xlabel('t [s]'); %OBLICZENIA Ix=b*h^3/12; %m4 k=48*E*Ix/L^3; %N/m m=0.5*L*b*h*ro; %kg wn=sqrt(k/m); %rad/s c_kr=2*m*wn; %kg/s c=c_kr*ksi; %kg/s dtkr=2/wn; dt=0.001; if dt>dtkr disp('Zmniejsz krok ca³kowania') end t=[0:dt:1]; p=[y1 y2 y3]; [u,v,a]=mrc(m,c,k,p,t,u0,v0); subplot(3,1,1); plot(t,u); xlabel('t [s]'); ylabel('u [m]'); grid on; subplot(3,1,2); plot(t,v); xlabel('t [s]'); ylabel('v [m]'); grid on; subplot(3,1,3); plot(t,a); xlabel('t [s]'); ylabel('a [m]'); grid on;
KaC2009