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
По окончании расчета иметь в виду:
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.