Scribd
Upload
30 views3 pages

Material Properties

The document defines material properties and member properties for a structural analysis. It then calculates stiffness matrices for each element and assembles them into a global stiffness matrix for the overall structure. Forces are applied and the document solves for displacements.

Uploaded by

Mayank Agrawal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
30 views3 pages

Material Properties

The document defines material properties and member properties for a structural analysis. It then calculates stiffness matrices for each element and assembles them into a global stiffness matrix for the overall structure. Forces are applied and the document solves for displacements.

Uploaded by

Mayank Agrawal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

%Material Properties%

fck=30;
E=5000*fck^0.5*10^3;
% Member Properties%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L=[2,2,3,3,3,3];
W=[.23,.23,.3,.3,.3,.3];
d=[.45,.45,.5,.5,.5,.5];
inc=[180,180,90,90,-90,-90];
for j= 1:6
l1(1,j)=cosd(inc(1,j));
m1(1,j)=sind(inc(1,j));
end
I=W.*d.^3/12;
EI=E*I;
A=W.*d;
AE=E.*A;
ke=zeros(6,6,6);
T=zeros(6,6,6);
kge=zeros(6,6,6);
kgf=zeros(18,18);
F=zeros(18,1);
kgb=zeros(12,12);
Fb=zeros(12,1);
for i= 1:6
T(:,:,i)=[l1(1,i) m1(1,i) 0 0 0 0; -m1(1,i) l1(1,i) 0 0 0 0; 0 0 1 0 0
0; 0 0 0 l1(1,i) m1(1,i) 0; 0 0 0 -m1(1,i) l1(1,i) 0; 0 0 0 0 0 1];
ke(:,:,i)=[AE(1,i)/L(1,i) 0 0 -AE(1,i)/L(1,i) 0 0;0 12*EI(1,i)/L(1,i)^3
6*EI(1,i)/L(1,i)^2 0 -12*EI(1,i)/L(1,i)^3 6*EI(1,i)/L(1,i)^2;0
6*EI(1,i)/L(1,i)^2 4*EI(1,i)/L(1,i) 0 -6*EI(1,i)/L(1,i)^2 2*EI(1,i)/L(1,i);-
AE(1,i)/L(1,i) 0 0 AE(1,i)/L(1,i) 0 0;0 -12*EI(1,i)/L(1,i)^3 -
6*EI(1,i)/L(1,i)^2 0 12*EI(1,i)/L(1,i)^3 -6*EI(1,i)/L(1,i)^2;0
6*EI(1,i)/L(1,i)^2 2*EI(1,i)/L(1,i) 0 -6*EI(1,i)/L(1,i)^2 4*EI(1,i)/L(1,i)];
kge(:,:,i)=T(:,:,i)'*ke(:,:,i)*T(:,:,i);
end
kg=zeros(18,18,6);
% for element 1 %
for i=1:6
for j=1:6
kg(i,j,1)=kge(i,j,1);
end
end
%for element 2%
for i=6:11
for j=6:11
kg(i,j,2)=kge(i-5,j-5,2);
end
end
%for element 3%
for i=1:3
for j=1:3
kg(i,j,3)=kge(i,j,3);
end
end
for i=7:9
for j=7:9
kg(i,j,3)=kge(i-3,j-3,3);
end
end
for i=1:3
for j=7:9
kg(i,j,3)=kge(i,j-3,3);
end
end
for i=7:9
for j=1:3
kg(i,j,3)=kge(i-3,j,3);
end
end
%for element 4%
for i=7:9
for j=7:9
kg(i,j,4)=kge(i-6,j-6,4);
end
end
for i=13:15
for j=13:15
kg(i,j,4)=kge(i-9,j-9,4);
end
end
for i=7:9
for j=13:15
kg(i,j,4)=kge(i-6,j-9,4);
end
end
for i=13:15
for j=7:9
kg(i,j,4)=kge(i-9,j-6,4);
end
end
%for element 5%
for i=4:6
for j=4:6
kg(i,j,5)=kge(i-3,j-3,5);
end
end
for i=10:12
for j=10:12
kg(i,j,5)=kge(i-6,j-6,5);
end
end
for i=4:6
for j=10:12
kg(i,j,5)=kge(i-3,j-6,5);
end
end
for i=10:12
for j=4:6
kg(i,j,5)=kge(i-6,j-3,5);
end
end
%for element 6%
for i=10:12
for j=10:12
kg(i,j,6)=kge(i-9,j-9,6);
end
end
for i=16:18
for j=16:18
kg(i,j,6)=kge(i-12,j-12,6);
end
end
for i=10:12
for j=16:18
kg(i,j,6)=kge(i-9,j-12,6);
end
end
for i=16:18
for j=10:12
kg(i,j,6)=kge(i-12,j-9,6);
end
end
kgf=sum(kg,3);
F(1,1)=10;
F(7,1)=5;
for i=1:12
for j=1:12
kgb(i,j)=kgf(i,j);
end
end
for i=1:12
Fb(i,1)=F(i,1);
end
d=inv(kgb)*Fb;

You might also like