%FUNCIÓN CAMBIOS_CURVATURA
format short
[fil,col]=size(matsin);
tipo_interp=menu('¿Qué tipo de interpolación desea?','Cúbica','Splines(recomendada)');
[curvsin_,curvcon_,h,coordcurv_,coord_int]=curv_todo_punto(tipo_interp,matsin,matcon,coord,n_elem,nbarra,n_puntos_barra);
disp('Decremento de la curvatura de los modos')
disp('curvcon_-curvsin_')
disp(curvcon_-curvsin_)
max(curvcon_-curvsin_)
disp('Cuanto mayor sea ese valor, mayor es la probabilidad de que exista un daño')
disp('SUMAMOS PARA TODOS LOS MODOS')
disp('grados de libertad: del 2 al n-1')
disp(sum((curvcon_-curvsin_)'))
xc=sum((curvcon_-curvsin_)');
for i=1:1:m
%%%%figure(i)
xcw=(abs(curvcon_(:,i))-abs(curvsin_(:,i)));
avexcw = mean(xcw(:)); %Media
stdxcw = std(xcw(:)); %Standard desviation
DIxcw = (xcw-avexcw)/stdxcw;%./abs(curvsin_(:,i));
%%%%xcw=(curvcon_(:,i)-abs(0.5*(curvcon_(:,i)+curvsin_(:,i))));
%%%%plot([2:1:n_elem n_elem+3:1:2*n_elem+1 2*n_elem+4:1:3*n_elem+2],xcw,'ro')
[kl,km]=max(DIxcw);
resultado_curvatura(i,:)=[i max(xcw) km coord_int(km)];
end
%figure(i+1)
%plot([2:1:n_elem n_elem+3:1:2*n_elem+1 2*n_elem+4:1:3*n_elem+2],xc,'r.-')
resultado_curvatura
%Tabla con curvatura en cada posición y diferencia de curvaturas
dist_elementos=1.27/(n_elem-2);
coordena_discret=[];
for i=1:n_elem-1
coordena_discret=[coordena_discret;i dist_elementos*(i-1)];
end
for i=1:m
xcw=(abs(curvcon_(:,i))-abs(curvsin_(:,i)));
avexcw = mean(xcw(:)); %Media
stdxcw = std(xcw(:)); %Standard desviation
DIxcw = (xcw-avexcw)/stdxcw;
tabla_curv_cadapunto(:,:,i)=[ coordena_discret(:,2) curvcon_(:,i) curvsin_(:,i) DIxcw];
figure(i)
plot(coordena_discret(:,2),curvcon_(:,i),coordena_discret(:,2),curvsin_(:,i),coordena_discret(:,2),DIxcw)
legend('Curvatura con daño','Curvatura sin daño','Variación de curvaturas');title('Curvaturas de los modos de vibración')
xlabel('Coordenada (m)');ylabel('Curvatura')
end
format short
curv_modo1=tabla_curv_cadapunto(:,:,1);
curv_modo2=tabla_curv_cadapunto(:,:,2);
curv_modo3=tabla_curv_cadapunto(:,:,3);
%curv_modo4=tabla_curv_cadapunto(:,:,4);
% curv_modo5=tabla_curv_cadapunto(:,:,5);
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
