По методике приведенной в автореферате для расчета гамма-процентного ресурса была составлена программа в среде Visual Basic Excel.
Код программы
Sub project()
'Объявление переменных
Dim tolshina_n(11), tolshina_k(11), rasstoyanie_t(11), F_s_tochkoy, gamma, summa_alfa_1, alfa_1(11), alfa_0, t_rab, Q_d, s_s, Q, V_sr, Q_k(11), Q_sr, Qu_a, m, summa_Q_k, tochka_num, index_arr, Q_д, t_д, Q_sr_dop, F, G, gamma_rasch, t_ostatoch, N, Q_sr_priblizh, Q_д_priblizh, alfa_б, U_r, ostat_ressurs, summa_tolshina_n
'Ввод данных
Worksheets("Лист1").Select
S_r = Range("C3").Value
gamma = Range("C4").Value
m = Range("C5").Value
t_д= Range("C6").Value
alfa_0 = Range("C7").Value
U_q = Range("C8").Value
U_r = Range("C9").Value
F_s_tochkoy = Range("C10").Value
index_arr = 0
tochka_num = 3
summa_Q_k = 0
Do Until IsEmpty(Range("F" & tochka_num).Value)
index_arr = index_arr + 1
tolshina_n(index_arr) = Range("G" & tochka_num).Value
tolshina_k(index_arr) = Range("H" & tochka_num).Value
Q_k(index_arr) = 1 - (tolshina_k(index_arr) / tolshina_n(index_arr))
Range("I" & tochka_num).Value = Q_k(index_arr)
tochka_num = tochka_num + 1
Loop
summa_Q_k = 0
summa_tolshina_n = 0
For i = 0 To index_arr
summa_Q_k = summa_Q_k + Q_k(i)
summa_tolshina_n = summa_tolshina_n + tolshina_n(index_arr)
Next
Q_sr = summa_Q_k / index_arr
Range("I15").Value = Q_sr
tolshina_n_sr = summa_tolshina_n / index_arr
'Расчет параметров
V_sr = Q_sr / (t_д ^ m)
index_arr = 0
tochka_num = 3
Do Until IsEmpty(Range("F" & tochka_num).Value)
index_arr = index_arr + 1
Q_k(index_arr) = Range("I" & tochka_num).Value
alfa_1(index_arr) = ((Q_k(index_arr) ^ 2 - alfa_0 ^ 2) / t_a ^ (2 * m)) - V_sr ^ 2
Range("J" & tochka_num).Value = alfa_1(index_arr)
tochka_num = tochka_num + 1
Loop
summa_alfa_1 = 0
For i = 1 To index_arr
summa_alfa_1 = summa_alfa_1 + alfa_1(i)
Next
Range("J15").Value = summa_alfa_1
'Условие
alfa_б = Sqr(summa_alfa_1 / (index_arr - 1))
If alfa_б <= alfa_0 Then Q_d = 0 Else Q_d = Sqr(alfa_б^ 2 - alfa_0 ^ 2)
Q_sr_priblizh = Q_sr + (U_q * Q_d / Sqr(index_arr - 2))
Q_д_priblizh = Q_d + (U_q * Q_d / Sqr(2 * index_arr - 8))
Q_sr_dop = 1 - (S_r / tolshina_n_sr)
F = (Q_sr_dop - Q_sr_priblizh) / Sqr(alfa_0 ^ 2 + Q_д_priblizh ^ 2)
G = (gamma / 100) * F_s_tochkoy
gamma_rasch = (Q_sr_dop * Q_sr - U_r * Q_d * Q_sr_dop ^ 2 + alfa_0 ^ 2 * (Q_sr ^ 2 - U_r ^ 2 * Q_d)) / (Q_sr ^ 2 - U_r ^ 2 * Q_d)
ostat_ressurs = t_д* (gamma_rasch ^ (1 / m) - 1)
'Результат расчета
Range("C19").Value = alfa_б
Range("C20").Value = Q_sr_priblizh
Range("C21").Value = Q_д_priblizh
Range("C22").Value = Q_sr_dop
Range("C23").Value = F
Range("C24").Value = G
Range("C25").Value = gamma_rasch
Range("C26").Value = ostat_ressurs
End Sub
Использованная литература