на щільність струму
* Визначаємо довжини кільця
dl = 0
do i = 1, n +1
dl = dl + l (i)
enddo
Write (6, *) 'Обьщяя довжина кільця А У З '
Write (6, *) 'dl =', dl, ' м '
Pp (1) = 0
Qp (1) = 0
* Потужність джерела харчування
do i = 1, n
do j = i +1, n +1
Pp (1) = Pp (1) + Pn (i) * l (j)/dl
Qp (1) = Qp (1) + Qn (i) * l (j)/dl
enddo
enddo
Write (6, *) 'Р джерела живлення Pp (1) =', Pp (1), 'кВт'
Write (6, *) 'Q джерела харчування Qp (1) = ', Qp (1),' кВАр '
* Визначимо перетоки потужності гілки
do i = 2, n +1
Pp (i) = Pp (i-1)-Pn (i-1)
Qp (i) = Qp (i-1)-Qn (i-1)
Write (6, *) 'Pp (', i, ') =', Pp (i)
Write (6, *) 'Qp (', i, ') =', Qp (i)
enddo
* Знайдемо точки Струмоподіл активних і реактивних потужностей
do i = 1, n
if (Pp (i +1). lt.0) then
Write (6, *) 'Точка Струмоподіл активних потужностей '
ia = i
Write (6, *) 'ia =', ia
Goto 1
endif
enddo
1do i = 1, n
if (Qp (i +1). lt.0) then
Write (6, *) 'Точка Струмоподіл реактивних потужностей '
ir = i
Write (6, *) 'ir =', ir p> Goto 2
endif
enddo
2do i = 1, n +1
Pp (i) = abs (Pp (i))
Qp (i) = abs (Qp (i))
enddo
* Визначення струмів на ділянках лінії
Write (6, *) 'Токі на ділянках лінії '
do i = 1, n +1
Car (i) = (sqrt (Pp (i) ** 2 + Qp (i) ** 2))/(Unom * sqrt (3))
Write (6, *) 'Car (', i, ') =', Car (i), ' A '
enddo
* Визначимо еквівалентний струм
Carl = 0
do i = 1, n +1
Carl = carl + (Car (i) ** 2 * l (i))/dl
enddo
Carl = sqrt (Carl)
Write (6, *) 'Еквівалентний струм дорівнює'
Write (6, *) 'Carl =', Carl, ' A '
* Вибираємо перетин кабелю
3 s = Carl/ej
Write (6, *) 'Попередній розрахунок перетину кабелю '
Write (6, *) 'S =', S, ' мм2 '
Write (*, *) 'Перетин кабелю'
Write (*, *) 'S =', S
Write (*, *) 'Введіть перетин кабелю по паспорту '
Write (*, *) 'S0 ='
Read (*, *) S0
Write (*, *) 'Введіть акт.сопротівленіе кабелю r0 Ом/км '
Write (*, *) 'r0 ='
Read (*, *) r0
Write (*, *) 'Введіть інд.сопротівленіе кабелю x0 Ом/км '
Read (*, *) x0
* Визначимо втрати напруги на ділянках лінії
Write (6, *) 'Втрати напруги на ділянках лінії '
do i = 1, n +1
Ul (i) = Ul (i) + Pp (i) * l (i) * r0/Unom + Qp (i) * l (i) * x0/Unom
Write (6, *) 'Ul (', i, ') =', Ul (i), ' B '
enddo
Визначимо втрати напруги до точок Струмоподіл
Ua = 0
if (ia.le.ir) then
do i = 1, ia
Ua = Ua + Ul (i)
enddo
Write (6, *) 'Втрата напруги до активної точки Струмоподіл '
Write (6, *) 'Ua =', Ua, 'B'
do i = ir +1, n +1
Ur = Ur + Ul (i)
enddo
Write (6, *) 'Втрата напруги до реактивної точки Струмоподіл '
Write (6, *) 'Ur =', Ur, 'B'
else
do i = ia +1, n +1
Ua = Ua + Ul (i)
enddo
Write (6, *) 'Втрата напруги до активної точки Струмоподіл '
Write (6, *) 'Ua =', Ua, 'B'
do i = 1, ir
Ur = Ur + Ul (i)
enddo
Write (6, *) 'Втрата напруги до реактивної точки Струмоподіл '
Write (6, *) 'Ur =', Ur, 'B'
endif
if (Ua.ge.Ur) then
UU = Ua
else
UU = Ur
endif
if (UU.le.Udop) then
Write (6, *) 'Вибрати провід перетином'
Write (6, *) 'S0 =', S0, ' мм2 '
else
Write (*, *) 'Вибрати провід більшого перетину '
Goto 3
endif
* Розрахунок для Аварійного режиму
* Аварія на початку лінії
UU = 0
do i = 1, n
do j = i +1, n +1
UU = UU + Pn (i) * l (j) * r0/Unom + Qn (i) * l (j) * x0/Unom
enddo
enddo
Write (6, *) 'При аварії в початку лінії: '
Write (6, *) 'UU =', UU, 'В'
if (UU.gt.Udopa) then p> Write (6, *) 'UU>', Udopa p> endif
* Аварія в кінці лінії
UU = 0
do i = +1, n
do j = 1, i
UU = UU + Pn (i) * l (j) * r0/Unom + Qn (i) * l (j) * x0/Unom
enddo
...
