11 0 352 KB
Tugas Dinamika Struktur Moh Ra’afi Iskandar Harjanegara 41155020150046 Teknik Sipil A No.1 clear all; clc; v0 = 1; vdot0 =1; T =1; omega =2*pi/T; t=linspace(0,5,500); v =(vdot0/omega)*sin(omega*t)+v0*cos(omega*t); plot(t,v); xlabel('waktu'); ylabel('besar respon'); title('tugas no 1') grid on
Grafik no.1 No.2 %free-vibration response of single - degree - of - freedom system % with damping(redaman subkritis underdamped system) clear all clc zeta = input('fraction of critcal damping = ') v0=1 T=1 vdot0=0
omega =2*pi/T omegad=omega*sqrt(1-zeta.^2) t=linspace(0,5,100) % argl=exp(-zeta*omega*t) R=sqrt((vdot0+v0*zeta*omega)/omegad)^2 phi=atan((vdot0+(v0*zeta*omega))/(omegad*v0)) Env=R.*argl v=Env.*cos(omegad*t-phi) plot1=plot(t,v,t,Env,t,-Env,'linewidth',1,'color',[1 0 0]','linestyle','--'); set(plot1(1),'linewidth',2,'color',[ 0 0 0 ],'linestyle','-'); xlabel=('time') ylabel=('displacment') grid on
Grafik no.2 No.3 close all; clear; clc; zeta = [0.0, 0.1, 0.4, 1.0, 2.5]; tspan=linspace( 0 , 1 , 100 ); for n = 1:5 [t, v]=ode45(@DLSDOF , tspan, [1 0]', [], zeta(n)) plot(t, v(:,1)); hold on end xlabel('time (secon)') ylabel('displacement(in)') axis([0 , 1.0 ,-1.5 ,1.5]); plot([0,1.0],[0,0],'k-'); grid on
function v =DLSDOF (t ,v, zeta); v=[v(2); -2*27.8*zeta*v(2)-27.8*27.8*v(1)];
Grafik no.3