Оценка параметрической надежности РЭС с использованием моделирования на ЭВМ постепенных отказов

write ('Верхняя граница температурного диапазона : '); read(Tv);Writeln(' -----------------------------------------------------------');Writeln(' Моделирование и рассчет займут некоторое время. ');Readkey;TextColor(13+Blink);Writeln(' ПРОИЗВОДИТСЯ МОДЕЛИРОВАHИЕ И РАССЧЕТ: ');TextColor(15);

sum:=0;sum1:=0;sum2:=0;sum3:=0;sum4:=0;sum5:=0;sum6:=0;sum7:=0;

sum8:=0;sum9:=0; sum10:=0;sum11:=0

;sum12:=0;sum13:=0;FOR n:=1 TO num DOBEGINif dR1<=5 thenR1:=Generator2(SR1-(SR1*dR1/100),SR1+(SR1*dR1/100))elseR1:=Generator(SR1,(SR1*dR1/300));if dR2<=5 thenR2:=Generator2(SR2-(SR2*dR2/100),SR2+(SR2*dR2/100))elseR2:=Generator(SR1,(SR1*dR1/300));if dR3<=5 thenR3:=Generator2(SR3-(SR3*dR3/100),SR3+(SR3*dR3/100))elseR3:=Generator(SR3,(SR3*dR3/300));RW:=Generator(SRW,(SRW*dRW/300));Corr(RW,SRW,SKOU,(SRW*dRW/300),(SKOU*dKOU/300),mzx,szx);KOU:=Generator(mzx,szx);Kexit:=(-R2/R1)*(1/(1+(1+R3/R1+2*R3/RW)/KOU));sum:=sum+(Kexit);sum1:=sum1+sqr(Kexit);temp:=Generator2(Tn,Tv);if (temp>=20) thenbegina:=a+1;dx1:=Generator(0,(Rtpol/300));R1:=R1+R1*Abs(20-Tv)*dx1;Rtemp1:=SR1+SR1*Abs(20-Tv)*dx1;dx2:=Generator(0,(Rtpol/300));R2:=R2+R2*Abs(20-Tv)*dx2;Rtemp2:=SR2+SR2*Abs(20-Tv)*dx2;dx3:=Generator(0,(Rtpol/300));R3:=R3+R3*Abs(20-Tv)*dx3;Rtemp3:=SR3+SR3*Abs(20-Tv)*dx3;dx4:=Generator(0,RWt/300);RW:=RW+RW*Abs(20-Tv)*dx4;RWtemp:=SRW+SRW*Abs(20-Tv)*dx4;Corr(dx4,0,0,RWt/300,KOUt/300,mzx,szx);dx5:=Generator(mzx,szx);KOU:=KOU+KOU*Abs(20-Tv)*dx5;KOUtemp:=SKOU+SKOU*Abs(20-Tv)*dx5;Kexit:=(-Rtemp2/Rtemp1)*(1/(1+(1+Rtemp3/Rtemp1+2*Rtemp3/RWtemp)/KOUtemp));sum2:=sum2+(Kexit);sum3:=sum3+sqr(Kexit);dx1:=Generator(0,(Rct/300));R1:=R1+R1*time*dx1;Rtime1:=SR1+SR1*time*dx1;dx2:=Generator(0,(Rct/300));R2:=R2+R2*time*dx2;Rtime2:=SR2+SR2*time*dx2;dx3:=Generator(0,(Rct/300));R3:=R3+R3*time*dx3;Rtime3:=SR3+SR3*time*dx3;dx4:=Generator(0,(Rct/300));RW:=RW+RW*time*dx4;RWtime:=SRW+SRW*time*dx4;Corr(dx4,0,0,RWct/300,KOUct/300,mzx,szx);dx5:=Generator(mzx,szx);KOU:=KOU+KOU*time*dx5;KOUtime:=SKOU+SKOU*time*dx5;Kexit:=(-Rtime2/Rtime1)*(1/(1+(1+Rtime3/Rtime1+2*Rtime3/RWtime)/KOUtime));sum4:=sum4+(Kexit); sum5:=sum5+sqr(Kexit);Kexit:=(-R2/R1)*(1/(1+(1+R3/R1+2*R3/RW)/KOU));sum6:=sum6+(Kexit);sum7:=sum7+sqr(Kexit);IF Kexit<(Kideal-Kideal*dKideal/100) THENIF Kexit>(Kideal+Kideal*dKideal/100) THEN P1:=P1+1;end;if (temp<20) thenbeginb:=b+1;dx1:=Generator(0,(Rtotr/300));R1:=R1+R1*Abs(20-Tn)*dx1;Rtemp1:=SR1+SR1*Abs(20-Tn)*dx1;dx2:=Generator(0,(Rtotr/300));R2:=R2+R2*Abs(20-Tn)*dx2;Rtemp2:=SR2+SR2*Abs(20-Tn)*dx2;dx3:=Generator(0,(Rtotr/300));R3:=R3+R3*Abs(20-Tn)*dx3;Rtemp3:=SR3+SR3*Abs(20-Tn)*dx3;dx4:=Generator(0,RWt/300);RW:=RW+RW*Abs(20-Tn)*dx4;RWtemp:=SRW+SRW*Abs(20-Tn)*dx4;Corr(dx4,0,0,RWt/300,KOUt/300,mzx,szx);dx5:=Generator(mzx,szx);KOU:=KOU+KOU*Abs(20-Tn)*dx5;KOUtemp:=SKOU+SKOU*Abs(20-Tn)*dx5;Kexit:=(-temp2/Rtemp1)*(1/(1+(1+Rtemp3/Rtemp1+2*Rtemp3/RWtemp)/KOUtemp));sum8:=sum8+(Kexit); sum9:=sum9+sqr(Kexit);dx1:=Generator(0,(Rct/300));R1:=R1+R1*time*dx1;Rtime1:=SR1+SR1*time*dx1;dx2:=Generator(0,(Rct/300));R2:=R2+R2*time*dx2;Rtime2:=SR2+SR2*time*dx2;dx3:=Generator(0,(Rct/300));R3:=R3+R3*time*dx3;Rtime3:=SR3+SR3*time*dx3;dx4:=Generator(0,RWct/300);RW:=RW+RW*time*dx4;RWtime:=SRW+SRW*time*dx4;Corr(dx4,0,0,RWct/300,KOUct/300,mzx,szx);dx5:=Generator(mzx,szx);KOU:=KOU+KOU*time*dx5;KOUtime:=SKOU+SKOU*time*dx5;Kexit:=(-Rtime2/Rtime1)*(1/(1+(1+Rtime3/Rtime1+2*Rtime3/RWtime)/KOUtime));sum10:=sum10+(Kexit); sum11:=sum11+sqr(Kexit);Kexit:=(-R2/R1)*(1/(1+(1+R3/R1+2*R3/RW)/KOU));sum12:=sum12+(Kexit); sum13:=sum13+sqr(Kexit);IF Kexit<(Kideal-Kideal*dKideal/100) THENIF Kexit>(Kideal+Kideal*dKideal/100) THEN P2:=P2+1;end;END;P1:=P1/a;P2:=P2/b;IF P2>P1 thenbegin P:=P1;mo1:=sum/num;mo2:=sum2/a;mo3:=sum4/a;mo4:=sum6/a;s1:=sqrt((sum1-sqr(sum)/num)/(num-1));s2:=sqrt((sum3-sqr(sum2)/a)/(a-1));s3:=sqrt((sum5-sqr(sum4)/a)/(a-1));s4:=sqrt((sum7-sqr(sum6)/a)/(a-1));end;if P2<P1 thenbeginP:=P2;mo1:=sum/num;mo2:=sum8/b;mo3:=sum10/b;mo4:=sum12/b;s1:=sqrt((sum1-sqr(sum)/num)/(num-1));s2:=sqrt((sum9-sqr(sum8)/b)/(b-1));s3:=sqrt((sum11-sqr(sum10)/b)/(b-1));s4:=sqrt((sum13-sqr(sum12)/b)/(b-1));end;

ClrScr;WriteLn('Коэффициент передачи: ',Kideal:6:3);WriteLn('Математическое ожидание, учитывая производственный допуск:',mo1:6:3);WriteLn('Среднеквадратичное отклоненение: ',s1:6:4);WriteLn('Математическое ожидание, учитывая температурный допуск: ' ,mo2:6:3);WriteLn('Среднеквадратичное отклоненение: ',s2:6:4);WriteLn('Математическое ожидание, учитывая старение: ',mo3:6:3);WriteLn('Среднеквадратичное отклоненение: ',s3:6:4);WriteLn('Математическое ожидание, учитывая все факторы: ',mo4:6:3);WriteLn('Среднеквадратичное отклоненение: ',s4:6:4);Writeln('-------------------------------------------------------------------------');WriteLn('Вероятность отсутствия параметрического отказа: ');WriteLn('P=',P:6:4);if num<4*Sqr(s4)/Sqr(0.01) thenBeginwriteln('Не достигнута заданная точность !');writeln('Следует сделать число реализаций процесса сделать>',num,'!');end;REPEAT UNTIL KeyPressed;END.

 Скачать реферат