Time  Dependent  Cantilever Beams 

clc
clear
E=69e12;
D=0.005;
A=0.25*(pi)*(D)^2;
I=(pi/4)*(D)^4;
L=1;
M=100;
N=10;
DX=L/M;
Time=1;
Time0=0;
dt=(Time-Time0)/N;
for t=1:N;
for i=1:M;
    P(i,t)=-10000*t;
    PP(i,t)=(P(i,t)/(6*E*I));
    X(i,t)=i*DX;
    Y(i,t)=PP(i,t)*((X(i,t))^2)*(3*L-X(i,t));
end
end
for t=1:N;
plot(X(:,t),Y(:,t))
hold on
end
grid on
xlabel('X Distance (meter)')
ylabel('Defliction')
axis equal

---

Unless otherwise noted, all content on this site is @Copyright by Ahmed Al Makky 2012-2014 - http://cfd2012.com