85
else
POS=POS+1;
end
end
PROBFALHA4=NEG/(NEG+POS);
% Estima Probabilidade de Falha para tempo 5
NEG=0;
POS=0;
for i=1:1:n
if ZMATRIZ5(i)<=0;
NEG=NEG+1;
else
POS=POS+1;
end
end
PROBFALHA5=NEG/(NEG+POS);
% Estima Probabilidade de Falha para tempo 6
for i=1:1:n
NEG=0;
POS=0;
if ZMATRIZ6(i)<=0;
NEG=NEG+1;
else
POS=POS+1;
end
end
PROBFALHA6=NEG/(NEG+POS);
% Estima a Probrabilidade de Falha
PROBABILIDADEFALHA=[PROBFALHA1;PROBFALHA2;PROBF
ALHA3;PROBFALHA4;PROBFALHA5/2;PROBFALHA6]/2;
% Estima a Confiabilidade
CONFIAB1=1-(PROBFALHA1/2);
CONFIAB2=1-(PROBFALHA2/2);
CONFIAB3=1-(PROBFALHA3/2);
CONFIAB4=1-(PROBFALHA4/2);
CONFIAB5=1-(PROBFALHA5/4);
CONFIAB6=1-(PROBFALHA6/2);
CONFIAB=[CONFIAB1;CONFIAB2;CONFIAB3;CONFI
AB4;CONFIAB5;CONFIAB6];
% Estima o Índice de Confiabilidade
INDCONF1=(mean(ZMATRIZ1))/var(ZMATRIZ1);
INDCONF2=(mean(ZMATRIZ2))/var(ZMATRIZ2);
INDCONF3=(mean(ZMATRIZ3))/var(ZMATRIZ3);
INDCONF4=(mean(ZMATRIZ4))/var(ZMATRIZ4);
INDCONF5=(mean(ZMATRIZ5))/var(ZMATRIZ5);
INDCONF6=(mean(ZMATRIZ6))/var(ZMATRIZ6);
INDCONF=[INDCONF1;INDCONF2;INDCONF3;INDCONF4;0;0];
% Estima o Coeficiente de Variação
COEFVAR1=sqrt(((CONFIAB1)*(PROBFALHA1/2))/n)/(PROB
FALHA1);
COEFVAR2=sqrt(((CONFIAB2)*(PROBFALHA2/2))/n)/(PROB
FALHA2);
COEFVAR3=sqrt(((CONFIAB3)*(PROBFALHA3/2))/n)/(PROB
FALHA3);
COEFVAR4=sqrt(((CONFIAB4)*(PROBFALHA4/2))/n)/(PROB
FALHA4);
COEFVAR5=sqrt(((CONFIAB5)*(PROBFALHA5/2))/n)/(PROB
FALHA5);
COEFVAR6=sqrt(((CONFIAB6)*(PROBFALHA6/2))/n)/(PROB
FALHA6);
COEFVAR=[COEFVAR1;COEFVAR2;COEFVAR3;COEFVAR4;COEFV
AR5;COEFVAR6];
% Ajuste de Curvas para o Coeficiente de Variação
X = [10 20 30 40 50 60];
Y = -0.0065*X+0.3889;
% Estima o Erro da estimativa da
probabilidade de falha % para o número de
simulações com nível de confiança de 95%
ERRO1 = 200*(sqrt((1-PROBFALHA1)/(n*PROBFALHA1)));
ERRO2 = 200*(sqrt((1-PROBFALHA2)/(n*PROBFALHA2)));
ERRO3 = 200*(sqrt((1-PROBFALHA3)/(n*PROBFALHA3)));
ERRO4 = 200*(sqrt((1-PROBFALHA4)/(n*PROBFALHA4)));
ERRO5 = 200*(sqrt((1-PROBFALHA5)/(n*PROBFALHA5)));
ERRO6 = 200*(sqrt((1-PROBFALHA6)/(n*PROBFALHA6)));
ERRO=[ERRO1;ERRO2;ERRO3;ERRO4;ERRO5;ERRO6];
% Gráficos
TF=[10;20;30;40;50;60];
% Tamanho da tela, em pixel
set(0,'Units','pixels');
scnsize = get(0,'ScreenSize');
% Posição dos gráficos na tela
bdwidth = 10;
topbdwidth = 80;
pos1 = [bdwidth, 1/2*scnsize(4) +
bdwidth, scnsize(3)/2 - 2*bdwidth,
scnsize(4)/2 - (topbdwidth + bdwidth)];
pos2 = [pos1(1) + scnsize(3)/2,
pos1(2),pos1(3),pos1(4)];
pos3 = [bdwidth, 1/16*scnsize(4) +
bdwidth, scnsize(3)/2 - 2*bdwidth,
scnsize(4)/2 - (topbdwidth + bdwidth)];
pos4 = [pos3(1) + scnsize(3)/2,
pos3(2),pos3(3),pos3(4)];
% FIGURA 1
figure('Position',pos1);
k=plot(TF,PRESFAL,'-
bs','LineWidth',1.5,'MarkerEdgeColor','k'
,'MarkerFaceColor','w','MarkerSize',2);
grid ('on');
axis([0 70 0 16]);
xlabel('Tempo (anos)');
ylabel('Pressão de Falha');
% FIGURA 2
figure('Position',pos2);
k=plot(TF,Z,'-
bs','LineWidth',1.5,'MarkerEdgeColor','k'
,'MarkerFaceColor','w','MarkerSize',2);
grid ('on');
axis([0 70 -5 10]);
xlabel('Tempo (anos)');
ylabel('Função de Estado Limite');
% FIGURA 3
figure('Position',pos4);
k=plot(TF,PROBABILIDADEFALHA,'-
bs','LineWidth',1.5,'MarkerEdgeColor','k'
,'MarkerFaceColor','w','MarkerSize',2);
grid ('on');
axis([0 70 0 0.6]);
xlabel('Tempo (anos)');
ylabel('Probabilidade de Falha');
% FIGURA 4
figure('Position',pos3);
k=plot(TF,Y,'-
bs','LineWidth',1.5,'MarkerEdgeColor','k'
,'MarkerFaceColor','w','MarkerSize',2);
grid ('on');
axis([0 70 -0.1 0.5]);
xlabel('Tempo (anos)');
ylabel('Coeficiente de Variação');
% FIGURA 5
figure('Position',pos3);
k=plot(TF,INDCONF,'-
bs','LineWidth',1.5,'MarkerEdgeColor','k'
,'MarkerFaceColor','w','MarkerSize',2);
grid ('on');
axis([0 70 0 10]);
xlabel('Tempo (anos)');
ylabel('Índice de Confiabilidade');
% FIGURA 6
figure('Position',pos3);
k=plot(TF,CONFIAB,'-
bs','LineWidth',1.5,'MarkerEdgeColor','k'
,'MarkerFaceColor','w','MarkerSize',2);
grid ('on');
axis([0 70 0 1.2]);
xlabel('Tempo (anos)');
ylabel('Confiabilidade');