%--------------------------------------------------------------------------
%Archivo de coordenadas y comienzo de Dibema
%--------------------------------------------------------------------------
clear all
close all
%Coordenadas de los nodos
coordena=[];
nodos=65;%65 128
dist_nodos=1.27/(nodos-1); %distancia entre nodos
%---------------------%
%Para análisis donde se ha prescindido de alguna señal
%i=1;%número de nodos que se han quitado del análisis
%nodo_elim=33;%nodo eliminado
%nodos=nodos-i;
%dist_nodos=1/(nodos);
%for i=1:nodo_elim-1
%coordena=[coordena;i dist_nodos*(i-1) 0];
%end
%for i=nodo_elim+1:nodos+1
%coordena=[coordena;i dist_nodos*(i-1) 0];
%end
%---------------------%
m=5;%número de modos
n_elem=129;%129 255%número de elementos para la discretización de Stubbs

nbarra=1;
n_puntos_barra=[nodos]; %nºnodos en cada barra
inic_fin=[1 nodos];%nodos inicial y final en cada barra 
factor_escala=1;%factor de escala para escalar los modos de vibración
%------------------
%Para análisis donde todos las señales son válidas
for i=1:nodos
coordena=[coordena;i dist_nodos*(i-1) 0];
end
%-----------------------------------
coordenadas=coordena(:,2:3);
[n,nb]=size(coordenadas);

%Matriz de masa
masa=diag([0.408 3.315 6.732 6.477 6.018 6.732 7.956 6.834 6.018 6.834 3.366 0.306]);

ipe100sin;
%En esta línea, poner el nombre del fichero del estado que hará las veces
%de estado sin daño
phisin=matriz(:,1);
angsin=matriz(:,2);
clear matriz;
ipe100con;
%En esta línea, poner el nombre del fichero del estado que hará las veces
%de estado con daño
phicon=matriz(:,1);
angcon=matriz(:,2);


dibema;