%-------------------------------------------------------------------------- %Archivo de datos de las propiedades dinámicas de la estructura, %provenientes de Spice %-------------------------------------------------------------------------- %m=5;%número de modos mat=[]; opcion=menu('Cargar archivos con las propiedades dinámicas de la barra SIN daño','Buscar archivo con las frecuencia y amortiguamiento',... 'Salir'); %Cargar archivo con los vectores de frecuencia y amortiguación if opcion==1 uiload else %Salida del programa disp('Gracias por usar "DIbEMA".') break end opcion=menu('Cargar archivos con las propiedades dinámicas de la barra SIN daño','Buscar archivo con los modos reales',... 'Salir'); %Cargar archivo con la matriz de modos if opcion==1 uiload else %Salida del programa disp('Gracias por usar "DIbEMA".') break end %Ordenar datos en vector columna para introducirlo en Dibema for i=1:m mat(:,i)=[ i;freq_sel(i);damp_sel(i);mod_real(:,i)]; %nº modo; Frecuencia natural (Hz);Amortiguación,Desplazamiento de cada grado de libertad end [q,m]=size(mat); matriz(1:q,1)=mat(:,1); if m>1 for i=2:m matriz((i-1)*q+1:i*q,1)=mat(:,i); end end matriz(:,2)=zeros(length(matriz(:,1)),1); for i=1:m %m es el nº de modos matriz(((i-1)*n+3*i+1):i*(n+3),1)=matriz(((i-1)*n+3*i+1):i*(n+3),1)*... (1/sqrt((matriz(((i-1)*n+3*i+1):i*(n+3),1))'*(matriz(((i-1)*n+3*i+1):i*(n+3),1)))); end