GEOMETRIA.m

function[V,C1,C2]=GEOMETRIA(Lt,Lc,Dnt,Dnc)
%beam
Yt=0:Dnt:Lt;
Xt=zeros(1,length(Yt));
Zt=ones(1,length(Yt))*640;
V=[Xt' Yt' Zt'];
%column1
Zc1=0:Dnc:Lc;
Xc1=zeros(1,length(Zc1));
Yc1=ones(1,length(Zc1))*2600;
C1=[Xc1' Yc1' Zc1'];
%column2
Zc2=0:Dnc:Lc;
Xc2=zeros(1,length(Zc2));
Yc2=ones(1,length(Zc1))*8300;
C2=[Xc2' Yc2' Zc2'];
return