Programa Cerchas En Matlab

Gonzalo Díaz Trillos
Joyce María Rojas
Elementos finitos


%PROGRAMA CERCHA EN MATLAB

% eL USUARIO INGRESA EL NUMERO DE NODOS Y ELEMENTOS DE LA CERCHA
nodos=input('Cuantos nodos tiene laestructura: ');
numel=input('Cuantos elementos tiene la estructura: ');

% EL USUARION INGRESA LAS UNIONES QUE TIENE CADA ELEMENTO CON CADA NODO
% SE CREA UNA MATRIZ CON LA INFORMACION DE ESASUNIONES
for i=1:numel
MATUNIO(i,1)=i;
fprintf('\n Indique las uniones de los elementos %d:',i);
NodoI(i)=input('\n nodo inicial:');
MATUNIO(i,2)=2*NodoI(i)-1;MATUNIO(i,3)=2*NodoI(i);
NodoF(i)=input(' nodo final:');
MATUNIO(i,4)=2*NodoF(i)-1;
MATUNIO(i,5)=2*NodoF(i);
end

% EL USUARIO TIENE QUE INGRESAR LAS COORDENADAS DE CADA NODOfprintf('\n Indique coordenadas de los nodos:\n');
for i=1:nodos
fprintf('\n Nodo %d:',i);
x(i)=input('\n Coordenada x:');
y(i)=input(' Coordenada y:');
end

%SE DEBEN INGRESAR LAS PROPIEDADES DE CADA ELEMENTO
for i=1:numel
fprintf('\n Indique propiedades de los elementos %d:',i);
E(i)=input('\n Modulo de elasticidad (E):');AREA(i)=input('\n Area del elemento:');

end
PE(i,1)=E(i);
PE(i,2)=AREA(i); % se almacena la info de los elementos
Dx(i)=x(NodoF(i))-x(NodoI(i)); % se calculael angulo de estos
Dy(i)=y(NodoF(i))-y(NodoI(i)); %con base a las coordenadas
b(1,i)=Dx(i); % de los nodos y se calcula
b(2,i)=Dy(i); %la longitud L de cada elemento
L(i)=(Dx(i)^2+Dy(i)^2)^0.5;
SENO(i)=Dy(i)/L(i);
COSENO(i)=Dx(i)/L(i);


% SE CREA LA MATRIZ DE RIGIDEZ DE CADA ELEMENTO Y SE CREA TAMBIENLA
% MATRIZ DE RIGIDEZ DE TODO EL SISTEMA
K=zeros(2*nodos);
for i=1:numel

E=PE(i,1); A=PE(i,2);
MKle=E*A/L(i)*[ 1,-1;-1,1 ];
n=b(:,i)'/L(i);...