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 la estructura: '); 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 ESAS UNIONES 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 …ver más…

=MATUNIO(:,2:5); K(VRE(i,:),VRE(i,:)) = K(VRE(i,:),VRE(i,:))+MKe; end; % Matriz K de todo el sistema con base en la matriz de K % uniones MATUNIO y la Matriz K de cada elemento

% SE DEBE INGRESAR LOS APOYOS O RESTRICCIONES QUE TENGA LA CERCHA % este vector acomoda la matriz para los calculos Vector.restricciones=1:2*nodos; NR = input(' Indique numero de nodos reestringidos:'); P=0; AUX=0; for i=1:NR NNR = input('indique el numero del nodo restringido:'); RE=input('\n Que tipo de reestriccion tiene(1=en 2 direcciones,0=en una direccion):'); if RE==1 Vector.restricciones(2*nodos-2*P)=2*NNR; Vector.restricciones(2*nodos-1-2*P)=2*NNR-1; Vector.restricciones(2*NNR)=2*nodos-2*P; Vector.restricciones(2*NNR-1)=2*nodos-1-2*P; AUX=AUX+2; else REXY=input('\n En que direccion esta restringido(1=X, 0=Y):'); if REXY==1 Vector.restricciones(2*nodos-2*P)=2*NNR-1; Vector.restricciones(2*NNR-1)=2*nodos-2*P; else Vector.restricciones(2*nodos-2*P)=2*NNR; Vector.restricciones(2*NNR)=2*nodos-2*P; end AUX=AUX+1; end P=P+1; end

% se cre la matriz de rigidez del sistema organizada para efecto de los
% calculos KM0=zeros(2*nodos); for i=1:2*nodos