Единое окно доступа к образовательным ресурсам

Усталостная долговечность и повреждаемость авиационных конструкций

Голосов: 0

В учебном пособии приведены результаты исследований в области расчетов долговечности авиационных конструкций при циклическом нагружении на основе методов схематизации переменной нагруженности, моделирования накопления повреждений в материале, стабилизации рассеяния свойств материалов путем оптимизации статистических моделей, обоснования вероятностных распределений показателей надежности, обработки цензурированных выборок, возникающих при целевых осмотрах самолетов. Приведены методы моделирования вертикальных перегрузок, возникающих в опасных зонах планера самолета, расчета повреждаемости и эквивалентной наработки в этих зонах, статистического анализа разброса усталостных свойств авиационных материалов по данным испытаний конструктивно-подобных образцов, поддержания жизненного цикла изделий авиационной техники. Учебное пособие предназначено для студентов, обучающихся по направлениям подготовки бакалавров, магистров и специалистов «Авиастроение», «Машиностроение», «Прикладная механика», «Конструкторско-технологическое обеспечение машиностроительных производств», «Испытание летательных аппаратов».

Приведенный ниже текст получен путем автоматического извлечения из оригинального PDF-документа и предназначен для предварительного просмотра.
Изображения (картинки, формулы, графики) отсутствуют.
     122 
Приложение 4.1.1. 
Программа схематизации случайного процесса нагружения методом полных циклов 
 
Public Sub FullCycle() 
Dim i As Integer, NX As Integer 
Dim MAX As Double, MIN As Double, n As Integer, kxx As Integer 
Dim LL As Integer, YFF As Double, k As Integer, j As Integer 
Dim ii As Integer, XXX As Double, ks As Integer, Delta As Double 
Dim XMIN As Double, XMAX As Double, s As Double 
Dim ws As String, wss As String, SS As Double 
Dim A(100) As Double, x(100) As Double 
Dim H(100) As Double, y(100) As Double, B(100) As Double 
Dim kprint As Integer 
 
ws = ActiveSheet.Name 
wss = ws: 
NX = CountRows(ws, 2, 500, 2) 
 
For i = 1 To NX 
 A(i) = CDbl(ReadFromCell(ws, i + 1, 2))     ‘считывание пиков 
Next 
 
MAX = FMAX(A): MIN = FMIN(A): n = Int((NX - 1) / 2 + 1): ii = NX: x(1) = (MAX - MIN)/2 
For LL = 2 To n 
YFF = x(1) + 1 
For i = 2 To ii 
 If Abs(A(i) - A(i - 1)) / 2 <= YFF Then 
  k = i 
   YFF = Abs(A(i) - A(i - 1)) / 2 
 End If 
Next 
For j = k - 1 To ii - 2 
A(j) = A(j + 2) 
Next 
x(LL) = YFF 
 ii = ii - 2 
For i = 2 To ii 
  If A(i) = A(i - 1) Then 
   k = i 
    For j = k To ii - 1 
      A(j) = A(j + 1) 
    Next 
     ii = ii - 1 
  End If 
Next 
For i = 3 To ii 
 test = 0 
 If A(i) > A(i - 1) And A(i - 1) < A(i - 2) Then test = 1 
 If A(i) < A(i - 1) And A(i - 1) > A(i - 2) Then test = 1 
 Select Case test 
 Case 0 
 MsgBox ("Three peeks on the line")  

 123 
 k = i - 1 
  For j = k To ii - 1 
    A(j) = A(j + 1) 
  Next 
  ii = ii - 1 
  Case 1 
  End Select 
   
 Next 
If ii <= 3 Then Exit For 
Next 
 XXX = x(1) 
For i = 2 To n 
 x(i - 1) = x(i) 
Next 
x(n) = XXX 
 
kprint = 3 
For i = 1 To n 
 Call WriteToCell(wss, x(i), i + 1, kprint)         ‘вывод амплитуд полных циклов 
Next 
 
'k = Int(ReadFromCell(ws, 2, 2)) 
k = 300 
Delta = (x(n) - x(1)) / k 
i = 1 
 XMIN = x(1) 
  XMAX = XMIN 
H(1) = 0 
y(1) = XMIN 
Do 
 i = i + 1 
 If i > k Then Exit Do 
 XMAX = XMAX + Delta 
  For j = 1 To n 
    If x(j) >= XMIN And x(j) < XMAX Then H(i) = H(i) + 1 
  Next 
 y(i) = (XMIN + XMAX) / 2 
 If H(i) <> 0 Then XMIN = XMAX 
 If H(i) = 0 Then 
  i = i - 1 
   k = k - 1 
 End If 
Loop While i <= k 
k = k + 1 
H(k) = 1 
y(k) = x(n) 
s = 0 
For i = 2 To k - 1 
 s = s + H(i) 
  H(i) = s / n 
Next  

 124 
‘вывод эмпирической функции распределения амплитуд полных циклов 
For i = 1 To k 
  Call WriteToCell(wss, y(i), i + 1, kprint + 1) 
  Call WriteToCell(wss, H(i), i + 1, kprint + 2) 
Next 
End Sub 
‘***********Вспомогательные функции********************* 
Public Sub WriteToCell(SheetName As String, txt As Variant, irow As Integer, icol As Integer) 
    ThisWorkbook.Worksheets(SheetName).Cells(irow, icol) = txt 
   End Sub 
Public Function ReadFromCell(SheetName As String, irow As Integer, icol As Integer) 
  ReadFromCell = ThisWorkbook.Worksheets(SheetName).Cells(irow, icol) 
End Function 
Public Function CountRows(SheetName As String, startrow As Integer, endrow As Integer, icol 
As Integer) As Integer 
 Dim i As Integer, txt As String 
  For i = startrow To endrow 
  txt = ReadFromCell(SheetName, i, icol) 
  If Trim(txt) = "" Then Exit For 
 Next 
 CountRows = i - startrow 
End Function 
Public Function FMIN(m As Variant) 
 Dim i As Integer, tm As Variant 
  tm = MaxRealNumber 
 For i = 1 To UBound(m) 
  If m(i) < tm Then tm = m(i) 
 Next 
  FMIN = tm 
End Function 
Public Function FMAX(m As Variant) 
Dim i As Integer, tm As Variant 
  tm = MinRealNumber 
 For i = 1 To UBound(m) 
  If m(i) > tm Then tm = m(i) 
 Next 
 FMAX = tm 
End Function 
‘*********************************************** 
   

 125 
Приложение 4.2.1. 
Программа расчета эмпирической функции распределения долговечности при 
эксплуатационном нагружении в среде VBA Excel 
 
 
 
Option Explicit 
Dim m As Double, cn As Double 
 Public Sub Resurs() 
    Dim i As Integer, ap As Double, A As Double 
    Dim P(9) As Double, L As Double, sigma As Double, zpx As Double 
    Dim ax As Double, bx As Double, f0 As Double, nmax As Integer, slgN As Double 
    Dim lgN As Double, teta As Double, xu As Double 
    Dim k As Integer, y() As Double, x() As Double 
    Dim s1 As Double, s2 As Double, s3 As Double, s4 As Double, kln As Double 
     
      slgN = 0.15 
      k = 7 
      kln = Log(10) 
      ReDim y(k), x(k) 
      s1 = 0: s2 = 0: s3 = 0: s4 = 0 
      For i = 1 To k 
       x(i) = Cells(i + 1, 2).Value 
       y(i) = kln * Cells(i + 1, 3).Value 
       Cells(i + 1, 4) = y(i) 
       s1 = s1 + x(i) 
       s2 = s2 + y(i) 
      Next 
       s1 = s1 / k: s2 = s2 / k 
       For i = 1 To k 
        s3 = s3 + (x(i) - s1) ^ 2 
        s4 = s4 + y(i) * (x(i) - s1) 
      Next 
       ax = s2: bx = s4 / s3: cn = -1 / bx: f0 = Exp(s2 - bx * s1)  

 126 
      
      Cells(2, 6) = cn 
      Cells(2, 7) = f0 
       For i = 1 To k 
        Cells(i + 1, 5) = ax + bx * (x(i) - s1) 'Log(f0 * Exp(-x(i) / cn)) 
      Next 
       xu = x(k) 
       nmax = 100 
       L = Cells(2, 8).Value 
       A = Cells(2, 9).Value 
       m = Cells(2, 10).Value 
       sigma = Cells(2, 11).Value 
       ap = 1 
       Cells(2, 12) = ap 
       P(1) = 0.01: P(2) = 0.05: P(3) = 0.1: P(4) = 0.3: P(5) = 0.5: P(6) = 0.7: P(7) = 0.9 
       P(8) = 0.95: P(9) = 0.99 
      For i = 1 To 9 
       zpx = WorksheetFunction.NormInv(P(i), 0, 1) 
       s3 = simpsonx(1, 0, xu, nmax) 
       teta = 10 ^ (zpx * slgN + A) * ap * cn / (L * f0 * sigma ^ m) 
       lgN = (Log(teta / s3)) / kln 
       Cells(i + 1, 13) = lgN 
       Cells(i + 1, 14) = zpx + 5 
       Cells(i + 1, 15) = P(i) 
     Next 
     
    End Sub 
     Public Function f(test As Integer, x As Double) 
      Select Case test 
      Case 1 
       f = (2 * x * (1 + x)) ^ (m / 2) * Exp(-x / cn) 
      Case 2 
       f = x * Exp(-x / cn) 
      Case 3 
       f = Exp(-x / cn) 
      End Select 
     End Function 
‘************************************************************************* 
&#39;**********Вспомогательная функция численного интегрирования методом Симпсона 
   Public Function simpsonx(test As Integer, xl As Double, xu As Double, nk As Integer) 
    Dim dm As Double, ii As Integer, jj As Integer, sum As Double, sm(3) As Double 
  
    sm(1) = 1#: sm(2) = 4#: sm(3) = 1# 
    sum = 0#: dm = (xu - xl) / (nk) 
    For ii = 1 To nk 
     For jj = 1 To 3 
       sum = sum + sm(jj) * f(test, xl + (ii - 1) * dm + (jj - 1) * dm / 2#) 
     Next 
    Next 
      simpsonx = sum * dm / 6# 
End Function 
   

 127 
Приложение 4.4.1. 
Программа расчета предполетной условной часовой повреждаемости по формуле (4.4.16) 
в среде VBA 
 
 
 
Public Sub predpolet() 
   Dim T As Double, m As Double, teta As Double, s As Double,povr As Double, k As Integer, i 
As Integer 
   Dim , nycp As Double, r As Integer, nymax As Double, nymin As Double, HH As Double 
   Dim nyk As Double, deltanyk As Double, Group As String, Sny As Double, zp As Double 
   Dim C(8) As Double, H() As Double, xs() As Double, ny() As Double 
   Group = Cells(2, 4).Value: Group = Trim(Group) 
   T = Cells(2, 5).Value: m = Cells(2, 6).Value: r = Cells(2, 7).Value: Sny = Cells(2, 8).Value 
   T = T / 60 
      ReDim H(r), xs(r), ny(r) 
   For i = 1 To 8 
    If Group = "A" Then C(i) = Cells(i + 2, 2).Value 
    If Group = "B" Then C(i) = Cells(i + 2, 3).Value 
   Next 
   If Group = "A" And nymax >= 6 Then 
     C(7) = 80000000: C(8) = -3.19 
   End If 
   If Group = "B" And nymax >= 8 Then 
     C(7) = 100000000000000#: C(8) = -4.25 
   End If 
   deltanyk = (Log(T) * (C(3) - C(5)) + C(4) - C(6)) / r 
   nymin = C(5) * Log(T) + C(6): nymax = C(3) * Log(T) + C(4):  nycp = C(1) * Log(T) + C(2) 
   HH = C(7) * Exp(C(8) * nycp):  teta = C(7) * Exp(C(8) * nymin) * (1 - Exp(C(8) * deltanyk)) 
   For k = 0 To r 
       ny(k) = nymin + deltanyk * k:   H(k) = C(7) * Exp(C(8) * ny(k)) 
       Cells(2 + k, 15) = k 
       Cells(2 + k, 16) = ny(k) 
       Cells(2 + k, 17) = H(k) 
    Next 
    s = 0: povr = 0 
   For k = 0 To r - 1 
    xs(k) = teta * Exp(C(8) * deltanyk * k) * (nymin + deltanyk * (k + 0.5)) ^ m 
    povr = povr + xs(k) 
    s = s + (H(k) - H(k + 1)) * ((ny(k) + ny(k + 1)) / 2) ^ m 
    Cells(2 + k, 18) = xs(k)  

 128 
   Next 
   Cells(2, 9) = nymin 
   Cells(2, 10) = nycp 
   Cells(2, 11) = nymax 
   Cells(2, 12) = HH 
   Cells(2, 13) = povr 
   Cells(2, 14) = deltanyk 
   Cells(2, 19) = s 
End Sub 
Приложение 5.4.1. 
 
Зависимость повреждаемостей  от изгибающего момента Mи1 и вертикальных перегрузок 
ny для m=4,5,6 
Вертикальн
ая 
перегрузка=
ny 
Gny Mи1 
Повреждаемость=
      
1,32 32985 18,07535 1,067E+05 1,929E+06 3,488E+07 0,010486 0,003355 0,001074 
1,45 36062 19,09076 1,328E+05 2,536E+06 4,841E+07 0,041006 0,018453 0,008304 
1,65 40746 20,63648 1,814E+05 3,743E+06 7,723E+07 0,178506 0,116029 0,075419 
1,78 43747 21,62681 2,188E+05 4,731E+06 1,023E+08 0,370151 0,288717 0,2252 
1,58 38273 19,82039 1,543E+05 3,059E+06 6,063E+07 0,113165 0,065636 0,038069 
1,65 39872 20,34806 1,714E+05 3,488E+06 7,098E+07 0,178506 0,116029 0,075419 
2,69 63898 28,27664 6,393E+05 1,808E+07 5,112E+08 8,157307 13,78585 23,29809 
1,58 37251 19,48313 1,441E+05 2,807E+06 5,470E+07 0,113165 0,065636 0,038069 
1,84 42949 21,36347 2,083E+05 4,450E+06 9,507E+07 0,497871 0,418212 0,351298 
1,52 35122 18,78056 1,244E+05 2,336E+06 4,388E+07 0,073116 0,03802 0,019771 
2,50 57620 26,2049 4,716E+05 1,236E+07 3,238E+08 5,0625 7,59375 11,39063 
1,39 31954 17,73512 9,893E+04 1,755E+06 3,112E+07 0,023134 0,009022 0,003519 
2,69 61522 27,49256 5,713E+05 1,571E+07 4,318E+08 8,157307 13,78585 23,29809 
2,36 53836 24,95618 3,879E+05 9,680E+06 2,416E+08 3,42102 4,652587 6,327519 
2,10 47905 22,99895 2,798E+05 6,435E+06 1,480E+08 1,4641 1,61051 1,771561 
4,32 98547 39,71081 2,487E+06 9,875E+07 3,922E+09 121,4933 403,3578 1339,148 
3,41 77788 32,86034 1,166E+06 3,831E+07 1,259E+09 33,73403 81,299 195,9306 
3,54 80754 33,83912 1,311E+06 4,437E+07 1,501E+09 41,62314 105,7228 268,5359 
3,60 82123 34,29089 1,383E+06 4,741E+07 1,626E+09 45,6976 118,8138 308,9158 
4,32 98547 39,71081 2,487E+06 9,875E+07 3,922E+09 121,4933 403,3578 1339,148 
3,02 68892 29,92466 8,019E+05 2,400E+07 7,181E+08 16,64966 33,63232 67,93729 
2,36 53699 24,91097 3,851E+05 9,593E+06 2,390E+08 3,42102 4,652587 6,327519 
2,63 59843 26,93849 5,266E+05 1,419E+07 3,822E+08 7,059118 11,50636 18,75537 
2,95 66950 29,2838 7,354E+05 2,153E+07 6,306E+08 14,45901 28,19506 54,98037 
2,76 62638 27,86084 6,025E+05 1,679E+07 4,677E+08 9,595126 16,88742 29,72186 
2,63 59532 26,83586 5,186E+05 1,392E+07 3,735E+08 7,059118 11,50636 18,75537 
3,08 69718 30,19724 8,315E+05 2,511E+07 7,582E+08 18,71774 38,93289 80,98042 
2,82 63667 28,20041 6,324E+05 1,784E+07 5,030E+08 10,97199 19,96903 36,34363 
3,15 71117 30,65891 8,835E+05 2,709E+07 8,305E+08 21,36751 45,94014 98,7713 
2,63 59222 26,73356 5,108E+05 1,365E+07 3,650E+08 7,059118 11,50636 18,75537 
3,54 79508 33,42794 1,249E+06 4,174E+07 1,395E+09 41,62314 105,7228 268,5359 
2,17 48610 23,2316 2,913E+05 6,767E+06 1,572E+08 1,873887 2,192448 2,565164 
2,63 58914 26,63192 5,030E+05 1,340E+07 3,568E+08 7,059118 11,50636 18,75537 
3,54 79090 33,29 1,228E+06 4,089E+07 1,361E+09 41,62314 105,7228 268,5359  41uM 51uM 61uM 41yn 51yn 61yn  

 129 
Вертикальн
ая 
перегрузка=
ny 
Gny Mи1 
Повреждаемость=
      
2,82 63004 27,98162 6,130E+05 1,715E+07 4,800E+08 10,97199 19,96903 36,34363 
5,04 
11171
1 44,05493 3,767E+06 1,659E+08 7,311E+09 266,3946 1076,234 4347,987 
2,23 49427 23,50121 3,050E+05 7,169E+06 1,685E+08 2,288866 2,815306 3,462826 
1,78 39350 20,1758 1,657E+05 3,343E+06 6,745E+07 0,370151 0,288717 0,2252 
2,17 47972 23,02106 2,809E+05 6,466E+06 1,489E+08 1,873887 2,192448 2,565164 
2,56 56593 25,86599 4,476E+05 1,158E+07 2,995E+08 5,922409 9,238958 14,41277 
2,17 47844 22,97882 2,788E+05 6,407E+06 1,472E+08 1,873887 2,192448 2,565164 
2,95 65041 28,65383 6,741E+05 1,932E+07 5,535E+08 14,45901 28,19506 54,98037 
4,32 94737 38,45351 2,186E+06 8,408E+07 3,233E+09 121,4933 403,3578 1339,148 
3,08 67362 29,41976 7,491E+05 2,204E+07 6,484E+08 18,71774 38,93289 80,98042 
4,00 87016 35,90558 1,662E+06 5,968E+07 2,143E+09 81,00000 243,0000 729,0000 
3,86 83742 34,82516 1,471E+06 5,122E+07 1,784E+09 66,90586 191,3507 547,2631 
3,86 83742 34,82516 1,471E+06 5,122E+07 1,784E+09 66,90586 191,3507 547,2631 
3,80 82441 34,39583 1,400E+06 4,814E+07 1,656E+09 61,4656 172,1037 481,8903 
3,86 83742 34,82516 1470863,9 51223070 1,784E+09 66,90586 191,3507 547,2631 
4,32 93722 38,11856 2111280,5 80478974 3,068E+09 121,4933 403,3578 1339,148 
4,13 89600 36,7583 1825667,4 67108432 2,467E+09 95,97925 300,4151 940,2991 
1,97 42506 21,21728 202655,7 4299802,8 91230120 0,885293 0,858734 0,832972 
1,91 41212 20,79026 186827,02 3884182,3 80753160 0,68575 0,624032 0,567869 
1,58 33998 18,40964 114863,27 2114591,5 38928868 0,113165 0,065636 0,038069 
3,08 66096 29,00198 707474,18 20518152 595067035 18,71774 38,93289 80,98042 
2,95 63307 28,08161 621853,4 17462645 490379178 14,45901 28,19506 54,98037 
2,76 59055 26,67845 506573,39 13514593 360548392 9,595126 16,88742 29,72186 
2,30 49213 23,43059 301392,82 7061811,5 165462411 2,8561 3,71293 4,826809 
3,34 71465 30,77375 896853,87 27599557 849341858 29,9822 70,15834 164,1705 
3,21 68684 29,85602 794561,75 23722451 708257982 23,85443 52,7183 116,5074 
1,71 36482 19,22936 136728,59 2629203,2 50557895 0,254117 0,180423 0,1281 
1,91 40749 20,63747 181395,21 3743538,2 77257158 0,68575 0,624032 0,567869 
1,52 32333 17,86019 101752,31 1817315,7 32457603 0,073116 0,03802 0,019771 
1,58 33511 18,24893 110904,62 2023890,6 36933837 0,113165 0,065636 0,038069 
1,91 40511 20,55893 178649,59 3672844,5 75509753 0,68575 0,624032 0,567869 
3,15 66811 29,23793 730779,72 21366486 624711833 21,36751 45,94014 98,7713 
1,52 32239 17,82917 101047,25 1801588,6 32120830 0,073116 0,03802 0,019771 
2,04 43268 21,46874 212435,08 4560713,4 97912770 1,169859 1,216653 1,265319 
1,58 33017 18,08591 106994,5 1935092,9 34997917 0,113165 0,065636 0,038069 
Суммарные значения за 
полет 4,950E+07 1,623E+09 5,523E+10 1731,147 5140,283 16077,16 
 
Приложение 5.6.1. 
 
Программа расчета эквивалентного налета и текущей долговечности 
маневренного самолета в условиях эксплуатации 
 
Dim A(30000, 2)               &#39; Объявление массива для хранения Ny и GNy 
Const P4Norm = 48850000#      &#39; «Нормативная» повреждаемость для степени 4 
Const P5Norm = 2419000000#    &#39; «Нормативная» повреждаемость для степени 5 41uM 51uM 61uM 41yn 51yn 61yn  

 130 
Const P6Norm = 132400000000#  &#39; «Нормативная» повреждаемость для степени 6 
 
Sub Расчет_эквивалентного_налета() 
 
&#39; Расчет времени полета 
    n = 0 
    x = ActiveCell.Column 
    y = ActiveCell.Row 
    TT = 0 
    While Cells(y + TT, x).Value <> "" 
        TT = TT + 1 
    Wend 
    Tflight = (TT - 8) / 8 / 3600   &#39;Значение времени полета в часах 
     
&#39; Считывание последовательности значений Ny, GNy и запись экстремумов в массив 
     
    arow = y 
    acolumn = x + 4 
    k = 0 
    Max = -100 
    Min = 100 
    While Cells(y + k, x).Value <> "" 
            If (n = 0) Then 
                A(n, 1) = Cells(y + k, x).Value 
                A(n, 2) = Cells(y + k, x + 1).Value 
                n = n + 1 
                Delta = 0 
            Else 
                If (Delta = 0) Then 
                    If ((Cells(y + k, x).Value) > A((n - 1), 1)) Then 
                        Delta = 1 
                        Max = Cells(y + k, x).Value 
                        MaxT = Cells(y + k, x + 1).Value 
                        Min = 100 
                    End If 
                    If ((Cells(y + k, x).Value) < A((n - 1), 1)) Then 
                        Delta = -1 
                        Min = Cells(y + k, x).Value 
                        MinT = Cells(y + k, x + 1).Value 
                        Max = -100 
                    End If 
                Else 
                    If (Delta = 1) Then 
                        If (((Cells(y + k, x).Value) > Max) And ((Cells(y + k, x).Value) - Max) <= 1.5) 
Then 
                            Max = Cells(y + k, x).Value 
                            MaxT = Cells(y + k, x + 1).Value 
                        Else 
                            If (((Cells(y + k, x).Value) < (Max)) And (Max - (Cells(y + k, x).Value)) <= 
1.5) Then 
                                Delta = -1 
                                A(n, 1) = Max  

 131 
                                A(n, 2) = MaxT 
                                Max = -100 
                                Min = Cells(y + k, x).Value 
                                MinT = Cells(y + k, x + 1).Value 
                                n = n + 1 
                            End If 
                        End If 
                    Else 'Delta=-1 
                        If (((Cells(y + k, x).Value) < Min) And (Min - (Cells(y + k, x).Value)) <= 1.5) 
Then 
                            Min = Cells(y + k, x).Value 
                            MinT = Cells(y + k, x + 1).Value 
                        Else 
                            If (((Cells(y + k, x).Value) > Min) And ((Cells(y + k, x).Value) - Min) <= 1.5) 
Then 
                                Delta = 1 
                                A(n, 1) = Min 
                                A(n, 2) = MinT 
                                Max = Cells(y + k, x).Value 
                                MaxT = Cells(y + k, x + 1).Value 
                                Min = 100 
                                n = n + 1 
                            End If 
                        End If 
                    End If 
                End If 
            End If 
             
        k = k + 1 
    Wend 
    A(n, 1) = Cells(y + k - 1, x).Value 
    A(n, 2) = Cells(y + k - 1, x + 1).Value 
     
&#39; Фильтрация полученного массива экстремумов Ny и Gny (амплитуда фильтрации – 0,5) 
     
    For Delta = 5 To 50 Step 5 
    i1 = 1 
    While i1 < n 
        If (Abs(A(i1, 1) - A((i1 + 1), 1)) < (Delta / 100)) Then 
            For i = i1 To n - 1 
                A(i, 1) = A((i + 2), 1) 
                A(i, 2) = A((i + 2), 2) 
            Next i 
            A((n - 1), 1) = 0 
            A((n - 1), 2) = 0 
            A(n, 1) = 0 
            A(n, 2) = 0 
            n = n - 2 
            If i1 > 2 Then 
                i1 = i1 - 1 
            Else 
                i1 = 1  


    
Яндекс цитирования Яндекс.Метрика