Определение показателей безотказности РЭУ с учётом внезапных отказов и разных законов распределения

y">  w[3]:=StrToInt(Edit4.Text);

  w[4]:=StrToInt(Edit4.Text);

  w[5]:=StrToInt(Edit8.Text); ro:=StrToFloat(Edit9.Text); beta:=StrToFloat(Edit10.Text);

  w[6]:=StrToInt(Edit11.Text); MO:=StrToFloat(Edit12.Text); SKO:=StrToFloat(Edit13.Text);

  {организация цикла по индексу j}

  j:=1; T0:=0; Pt:=0; Tgamma:=0;

  while j<=N1 do

     begin

       {Получение случайных времен отказов  для каждого элемента РЭУ исходя  из его закона}

       for i:=1 to n do

         case w[i] of

         1: toe[i]:=exp(lambda);

         2: toe[i]:=veibull(ro,beta);

         3: toe[i]:=normal(MO,SKO);

         end;

       {Определение времени отказа всего  РЭУ}

       Toreu[j]:=toe[1];

       for i:=2 to n do

         if toe[i]<Toreu[j] then Toreu[j]:=toe[i];

       {Вывод нескольких реализаций на экран}

       if (j=1) or (j=111) or (j=222) or (j=333) or (j=444)

       or (j=555) or (j=666) or (j=777) or (j=888) or (j=999) or (j=1000) then

       Memo1.Lines.Add(IntToStr(j)+'-ая реализация t[R1]='+FloatToStrF(toe[1],fffixed,5,2)+

       ' t[R2]='+FloatToStrF(toe[2],fffixed,5,2)+' t[R3]='+FloatToStrF(toe[3],fffixed,5,2)+

       ' t[R4]='+FloatToStrF(toe[4],fffixed,5,2)+' t[DA1]='+FloatToStrF(toe[5],fffixed,5,2)+

       ' t[XP1]='+FloatToStrF(toe[6],fffixed,5,2)+' Toreu='+FloatToStrF(Toreu[j],fffixed,5,2));

       {Статистическая обработка результатов}

       if j=N1 then

         begin

           {Нахождение вероятности безотказной  работы РЭУ}

           Nt:=0;

           for i:=1 to N1 do

             if Toreu[i]<t then Nt:=Nt+1;

           Pt:=(N1-Nt)/N1;

           {Нахождение наработки на отказ}

           T0:=0;

           for i:=1 to N1 do T0:=T0+Toreu[i];

           T0:=T0/N1;

           {Проверка условия на необходимое  число реализаций}

           if N1<4*Pt*(1-Pt)/sqr(0.005) then N1:=N1+1000

              else

                begin

                  {Нахождение гамма-процентной наработки  до отказа}

                   QuickSort(Toreu,N1);

                   Tgamma:=Toreu[StrToInt(FloatToStrF(0.95*N1,fffixed,10,0))];

                end;

        end;

       j:=j+1;

     end;

   {Вывод результатов на экран}

    Memo1.Lines.Add('Вероятность безотказной  работы за заданное время tз  Pt='+FloatToStrF(Pt,fffixed,5,5));

    Memo1.Lines.Add('Наработка на отказ T0='+FloatToStrF(T0,fffixed,5,5));

    Memo1.Lines.Add('Гамма-процентная наработка  до отказа Tgamma='+FloatToStrF(Tgamma,fffixed,5,5));

    Memo1.Lines.Add('Число реализаций РЭУ  N1='+IntToStr(N1));

end;

end.

Приложение Б 

Результаты  работы программы 

1-ая реализация t[R1]=7524,10 t[R2]=6248,70 t[R3]=2553,80 t[R4]=6195,40 t[DA1]=1,0863E07 t[XP1]=6009,00 Toreu=2553,80

111-ая реализация t[R1]=6168,80 t[R2]=5109,00 t[R3]=5997,70 t[R4]=5974,70 t[DA1]=3,282E08 t[XP1]=1758,00 Toreu=1758,00

222-ая реализация t[R1]=15470,00 t[R2]=655,97 t[R3]=18714,00 t[R4]=17502,00 t[DA1]=2,7793E07 t[XP1]=4643,60 Toreu=655,97

333-ая реализация t[R1]=11661,00 t[R2]=6901,80 t[R3]=20615,00 t[R4]=731,05 t[DA1]=6,8772E07 t[XP1]=3928,10 Toreu=731,05

444-ая реализация t[R1]=21992,00 t[R2]=5414,00 t[R3]=6677,40 t[R4]=11061,00 t[DA1]=4,1106E05 t[XP1]=3375,20 Toreu=3375,20

555-ая реализация t[R1]=143,26 t[R2]=11626,00 t[R3]=2983,80 t[R4]=5679,70 t[DA1]=6,8487E06 t[XP1]=5111,90 Toreu=143,26

666-ая реализация t[R1]=1881,40 t[R2]=4151,00 t[R3]=3931,00 t[R4]=10649,00 t[DA1]=1,9709E08 t[XP1]=5339,50 Toreu=1881,40

777-ая реализация t[R1]=2138,70 t[R2]=8771,90 t[R3]=4915,80 t[R4]=994,71 t[DA1]=3,0661E07 t[XP1]=3518,20 Toreu=994,71

888-ая реализация t[R1]=170,83 t[R2]=3052,70 t[R3]=2347,00 t[R4]=8009,80 t[DA1]=1,0414E05 t[XP1]=3393,80 Toreu=170,83

999-ая реализация t[R1]=4244,30 t[R2]=6199,90 t[R3]=3059,30 t[R4]=3077,80 t[DA1]=4,2318E06 t[XP1]=4141,90 Toreu=3059,30

1000-ая реализация t[R1]=804,16 t[R2]=10347,00 t[R3]=5098,00 t[R4]=9210,50 t[DA1]=2,9218E07 t[XP1]=2263,60 Toreu=804,16

Вероятность  безотказной работы за заданное время tз Pt=0,54450

Наработка на отказ T0=1484,90000

Гамма-процентная наработка до отказа Tgamma=81,60800

Число реализаций РЭУ N1=40000