%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%MÉTODO 6%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%AF"%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%FUNCIÓN AF
F=matsin*diag(1./wsin.^2)*matsin';
%No importa como estén normalizados los modos, puesto que luego se
%deriva dos veces
Fc=matcon*diag(1./wcon.^2)*matcon';
[fil,col]=size(F);
for i=1:1:col
%for j=2:1:fil-1
% curvflexsin(j-1,i)=(1/(h(j-1))^2)*(F(j-1,i)-2*F(j,i)+F(j+1,i));
%Estamos derivando dos veces
% curvflexcon(j-1,i)=(1/(h(j-1))^2)*(Fc(j-1,i)-2*Fc(j,i)+Fc(j+1,i));
%end
for j=2:1:3
curvflexsin(j-1,i)=(1/(h(j-1))^2)*(F(j-1,i)-2*F(j,i)+F(j+1,i));
%Estamos derivando dos veces
curvflexcon(j-1,i)=(1/(h(j-1))^2)*(Fc(j-1,i)-2*Fc(j,i)+Fc(j+1,i));
end
for j=5:1:8
curvflexsin(j-2,i)=(1/(h(j-1))^2)*(F(j-1,i)-2*F(j,i)+F(j+1,i));
%Estamos derivando dos veces
curvflexcon(j-2,i)=(1/(h(j-1))^2)*(Fc(j-1,i)-2*Fc(j,i)+Fc(j+1,i));
end
for j=10:1:11
curvflexsin(j-3,i)=(1/(h(j-1))^2)*(F(j-1,i)-2*F(j,i)+F(j+1,i));
%Estamos derivando dos veces
curvflexcon(j-3,i)=(1/(h(j-1))^2)*(Fc(j-1,i)-2*Fc(j,i)+Fc(j+1,i));
end
end
disp('Variación de la curvatura de la flexibilidad.')
disp('Estado con daño respecto al no dañado')
disp(curvflexcon-curvflexsin)
max(curvflexcon-curvflexsin)
disp('El valor máximo de esta matriz denota el grado de libertad donde quizá haya un daño')
disp('puesto que un daño provoca una mayor deformacion, y por tanto, una mayor curvatura de los modos')
disp('Nótese que hemos representado la curvatura de F del estado dañado menos lo mismo pero del estado sin daño')
disp('Ojo, este metodo todavia no está bien desarrollado, por lo que no es muy fiable')
disp('Le recomiendo usar otro método')
disp('SUMAMOS PARA TODOS LOS MODOS')
disp('grados de libertad: del 2 al n-1')
disp(sum((curvflexcon-curvflexsin)'))
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
