Координация изоляции электрооборудования: Учебное пособие к практическим занятиям, страница 42

d12 = SQR((Hsr - Lprtr + Hsr1 - Lprpr) ^ 2 + b ^ 2)

d112 = SQR((Hsr - Lprtr - Hsr1 + Lprpr) ^ 2 + b ^ 2)

k21 = ABS(LOG(d12 / d112) / LOG(2 * Hsr / r1 / 1000))

Kpt = 1.34 * k21

PRINT "Коэфф. связи между тросом и ближ. проводом, k21="; k21

PRINT "С учетом влияния импульсной короны, Kpt ="; Kpt

Npum1 = 5 * Hop + Hsr ^ 2 / 30

PRINT "Удельное число ударов молнии в год на 100 км линии и на 100 грозовых часов, Nпум1 ="; Npum1

82

Npum = Npum1 * L * Ngr / 10000

PRINT "Среднее число ударов молнии по линии в год, Nпум ="; Npum

Npumt = .5 * Npum

PRINT "Среднегодовое число ПУМ в середину троса, Nпумт ="; Npumt

PRINT "Среднегод. число ПУМ в вершину опоры, Nпум оп ="; Npumt

Npr = Npum * 10 ^ (Aug * SQR(Hop) / 90 - 4)

PRINT "Среднегодовое число ударов молнии в провод при наличии тросовой защиты углом Aug, Nпр="; Npr

INPUT "Нажмите ENTER"; AAAA

z0 = 600 * Zp / (600 + Zp)

Imp08 = u05 * 15 / z0

PRINT "Уровень грозоупор. при ПУМ в провод, Iмп0="; Imp08; "кА"

INPUT "Введите вероятность попадания ПУМ в провод, P(Iмпо)="; P3

PRINT "Введите вероят. попадания молнии в трос в середине пролета,"

INPUT "ею можно пренебречь,P(Iмто)="; P2

nu = u1 / u2: t1 = 2: t2 = 20: Ufmax = u / SQR(3): tf(0) = 0

a = ABS((nu ^ 2 - 1) / (1 - nu ^ 2 * t1 / t2)): x = 1

r1 = Rzi * 600 / (Rzi + 600): Re = 550 * r1 / (550 + r1)

CLS

PRINT "Кривые опасных параметров линии"

PRINT "  Tф, мкс     Iмо, кА     Uт, кВ        I'мо, кА "

FOR i = 1 TO 20

tf(i) = tf(i - 1) + .5

c(i) = (Re * x + (.6 * Hop / tf(i)) * x ^ 2 + .2 * Hop / tf(i)) * (1 - Kpt)

Imo(i) = (ut - u * Hsr * (1 - Kpt) - Ufmax) / c(i)

Imoo(i) = Imo(i) / tf(i)

Utt(i) = u2 * SQR(1 + a * t1 / tf(i)) / (SQR(1 + a * t1 / t2))

PRINT TAB(3); tf(i); TAB(13); Imo(i); TAB(33); Utt(i); TAB(55); Imoo(i)

NEXT

INPUT "Нажмите ENTER "; AAAAA

CLS

X1 = -Imoo(11): x2 = Imoo(10): x3 = Imoo(9)

83

y1 = Imo(11): y2 = Imo(10): y3 = Imo(9)

dx = (x3 - x2) / (x2 - X1):  dy = (y2 - y3) / (y1 - y2)

Imomin = (X1 * dx - x3 * dy) / (dx - dy)

Imoomin = (y3 * dx - y1 * dy) / (dx - dy)

a1 = (x2 - Imomin) * (y2 - Imoomin)

IF Imoomin < 0 THEN Imomin = (X1 * y1 - x2 * y2) / (y1 - y2)

IF Imoomin < 0 THEN Imoomin = 0: a1 = ((x2 - X1) / (y1 - y2)) * y1 * y2

m = SQR((1 + 3.1415 * SQR(.038 * .08 * a1)))

P1 = EXP(-.038 * iommin - .08 * iommin - 2 * SQR(.038 * .08 * a1)) * m

PRINT "Двумеpный интеграл,P(im,i'm)="; P1

PRINT "Годовое число перекрытий изоляции ЛЭП:"

Npert = .5 * Npum * P2

PRINT "при ударе молнии в середину троса, Nперт="; Npert

Nperp = Npum * Npr * P3

PRINT "при прорыве молнии мимо тросов на провод, Nперп="; Nperp

Nperop = .5 * Npum * P1

PRINT "при ударе молнии в вершину опоры, Nпероп="; Nperop

Nper = Npert + Nperp + Nperop

PRINT "   суммарное число перекрытий, Nпер="; Nper

Pd = 1.6 * (u / Lper) - .06

IF Pd > 1 THEN Pd = 1

IF Pd < .1 THEN Pd = .1

Notkl = Nper * Pd

PRINT "Годовое число отключений ЛЭП в следст вие грозовых воздействий, Notkl="; Notkl

Naw = Notkl * (1 - kapv)

PRINT "Cреднегодовое число аварий на линии  вызванных грозовой деятельностью, Naw="; Naw

Nnad = 1 / Naw

PRINT "Среднегодовое число лет безаварийной работы, Nnad="; Nnad

Nonad = Nnad * 100 / L * (100 / Ngr)

PRINT "Удельное число лет безаварийной работы(на 100 км и 100 грозовых часов, N'над="; Nonad

PRINT " end"

END

84

По окончании расчета иметь в виду: