Scribd
Upload
28 views

Stiffness Final

Uploaded by

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

Stiffness Final

Uploaded by

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

%MATLAB PROGRAM FOR 2D frame

clear all;clc;close all;%clear all existing variables

%Obtain input file(Excel sheet)


[num]=xlsread('C:\Users\Roopak\Desktop\SEM-1\ASA\project\stiffness_input_s_i.xlsx')
xy_node=num(:,1:2);%coordinates vector
xy_node(any(isnan(xy_node),2),:)=[]
mem=num(:,4:5);
mem(any(isnan(mem),2),:)=[]
E=num(:,6);
E(any(isnan(E),2),:)=[]
A=num(:,7);
A(any(isnan(A),2),:)=[]
I=num(:,8);
I(any(isnan(I),2),:)=[]
dof_r=num(:,10:12);
dof_r(any(isnan(dof_r),2),:)=[]
ext_load=num(:,14:16);
ext_load(any(isnan(ext_load),2),:)=[]
no_elem=size(mem,1);% Found number of elements
no_node=size(xy_node,1);% founf number of nodes

dof_all=1:3*no_node; %martix with all dofs


sms=zeros(3*no_node); %initialised global stiffness matrix with zeros
d=zeros(3*no_node,1); %initialised displacement matrix with zeros
a=zeros(3*no_node,1); %initialised force vector with zeros

dof_res=[]; %initialised an empty matrix


for j=1:size(dof_r,1) %to find the dof corr. to restrains
dof_flag=3*(dof_r(j,1)-1)+dof_r(j,2);
dof_res=[dof_res dof_flag]; %concatinating restrained dof
d(dof_flag)=dof_r(j,3); %setting the input displacement corr. to restrains
end
dof_free=dof_all;%duplicating all dof matrix
dof_free(dof_res)=[];%deletinf the restained dofs

for j=1:size(ext_load,1)% Extracting loads from input file


a(3*(ext_load(j,1)-1)+ext_load(j,2))=ext_load(j,3);
end

for i=1:no_elem
elem_nodes=mem(i,1:2);
elem_cord=xy_node(elem_nodes,:);

E1=[(elem_cord(2,1)-elem_cord(1,1)) (elem_cord(2,2)-elem_cord(1,2))];
len_elem=sqrt((E1(1)^2)+((E1(2))^2));
E1=E1/len_elem;
E2=[-E1(2) E1(1)];

sm_bend=[(((12*E(i)*I(i)))/(len_elem^3)) (((6*E(i)*I(i)))/(len_elem^2)) -
(((12*E(i)*I(i)))/(len_elem^3))
(((6*E(i)*I(i)))/(len_elem^2));(((6*E(i)*I(i)))/(len_elem^2))
(((4*E(i)*I(i))/len_elem)) -(((6*E(i)*I(i)))/(len_elem^2))
(((2*E(i)*I(i))/len_elem));-(((12*E(i)*I(i)))/(len_elem^3))
-(((6*E(i)*I(i)))/(len_elem^2)) (((12*E(i)*I(i)))/(len_elem^3))
-(((6*E(i)*I(i)))/(len_elem^2));(((6*E(i)*I(i)))/(len_elem^2))
(((2*E(i)*I(i))/len_elem)) -(((6*E(i)*I(i)))/(len_elem^2))
(((4*E(i)*I(i))/len_elem))];
sm=(((E(i)*A(i))/len_elem))*[1 -1;-1 1];

sm_loc([1,4],[1,4])=sm;
sm_loc([2,3,5,6],[2,3,5,6])=sm_bend;

trans=[E1; E2];
trans(3,3)=1;
R=[trans zeros(3);zeros(3) trans];

k_elem=R'*sm_loc*R;

dof_elem=3*(elem_nodes(1)-1)+1:3*(elem_nodes(1));
dof_elem=[dof_elem 3*(elem_nodes(2)-1)+1:3*(elem_nodes(2))];
sms(dof_elem,dof_elem)=sms(dof_elem,dof_elem)+k_elem;
end
d(dof_free)=((sms(dof_free,dof_free))^-1)*(a(dof_free)-
(sms(dof_free,dof_res)*d(dof_res)))
a(dof_res)=sms(dof_res,:)*d

for i=1:no_elem
elem_nodes=mem(i,1:2);
elem_cord=xy_node(elem_nodes,:);

E1=[(elem_cord(2,1)-elem_cord(1,1)) (elem_cord(2,2)-elem_cord(1,2))]
len_elem=sqrt((E1(1)^2)+((E1(2))^2))
E1=E1/len_elem
E2=[-E1(2) E1(1)]

sm_bend1=[(((12*E(i)*I(i)))/(len_elem^3)) (((6*E(i)*I(i)))/(len_elem^2)) -
(((12*E(i)*I(i)))/(len_elem^3))
(((6*E(i)*I(i)))/(len_elem^2));(((6*E(i)*I(i)))/(len_elem^2))
(((4*E(i)*I(i))/len_elem)) -(((6*E(i)*I(i)))/(len_elem^2))
(((2*E(i)*I(i))/len_elem));-(((12*E(i)*I(i)))/(len_elem^3))
-(((6*E(i)*I(i)))/(len_elem^2)) (((12*E(i)*I(i)))/(len_elem^3))
-(((6*E(i)*I(i)))/(len_elem^2));(((6*E(i)*I(i)))/(len_elem^2))
(((2*E(i)*I(i))/len_elem)) -(((6*E(i)*I(i)))/(len_elem^2))
(((4*E(i)*I(i))/len_elem))]

sm1=(((E(i)*A(i))/len_elem))*[1 -1;-1 1]

sm_loc1([1,4],[1,4])=sm1
sm_loc1([2,3,5,6],[2,3,5,6])=sm_bend1
trans=[E1; E2]
trans(3,3)=1
n1=3*(elem_nodes(1)-1)+1:3*(elem_nodes(1))
n2=3*(elem_nodes(2)-1)+1:3*(elem_nodes(2))
mem_f([1,2,3],i)=((sm_loc1([1,2,3],[1,2,3])*trans*d(n1))+(sm_loc1([1,2,3],
[4,5,6])*trans*d(n2)))
mem_f([4,5,6],i)=((sm_loc1([4,5,6],[1,2,3])*trans*d(n1))+(sm_loc1([4,5,6],
[4,5,6])*trans*d(n2)))
end

f_temp([1,2,3],1)=mem_f([1,2,3],1);
f_temp([1,2,3],2)=mem_f([4,5,6],20);
f_temp([1,2,3],3)=mem_f([4,5,6],21);
f_temp([1,2,3],4)=mem_f([4,5,6],35);
p=[0,0,0,0;-14.6261,-14.6261,-14.6261,-14.6261;0,0,0,0];
tij=[0,1,0;-1,0,0;0,0,1];
reac(3,4)=zeros;
for i=1:4
reac(:,i)=(tij*f_temp([1,2,3],i))-p(:,i);
end
reac
mulim=670769
vulim=694.29
pulim=4993.832
dup=1:no_elem
dup([3,5,7,9,10,19,11,13,15,17,30,29,27,25,24])=[]
for i=1:size(dup',1)
r(dup(i))=(abs(mem_f(1,dup(i)))/vulim)+(abs(mem_f(2,dup(i)))/pulim)+
(abs(mem_f(3,dup(i)))/mulim);
end
qwe=[3,5,7,9,10,19,11,13,15,17,30,29,27,25,24];
for i=1:size(qwe',1)
r(qwe(i))=(abs(mem_f(1,dup(i)))/pulim)+(abs(mem_f(2,dup(i)))/vulim)+
(abs(mem_f(3,dup(i)))/mulim);
end
true=mem_f'
true1=r'

You might also like