Программа расчета

CLS

DIM c(100), U_tt(100), I_mo(100), I_moo(100), tf(100)

PRINT " Введите исходные данные"

INPUT " Номинальное напряжение линии, кВ, U="; u

INPUT " Длина линии, км, L="; L

INPUT " Число цепей на опоре, N_cep="; N_cep

INPUT " Число грозозащитных тросов, N_tr="; N_tr

INPUT " Импульсное сопротивление заземление опоры, Ом, R_zi="; R_zi

INPUT " Среднее число гроз. часов в году на трассе ЛЭП, N_gr="; N_gr

INPUT "Коэффициент успешного действия АПВ, K_apv="; K_apv

INPUT "Длина пролета, м, L_prol="; L_prol

INPUT "Стрела провеса троса, м, L_prtr="; L_prtr

INPUT "Стрела провеса провода, м, L_prpr="; L_prpr

INPUT "Импульсное напряжение перекрытия, кВ, U_t="; u_t

INPUT "Импульсное напряжение перекрытия, кВ, U₀₅="; u₀₅

INPUT "Длина пути перекрытия, см, L_per="; L_per

INPUT "Введите расст. между трос. и провод. по горизонтали, м, b=";b

INPUT "Введите расст. между трос. и провод. по вертикали, м, h="; H

INPUT "Введите высоту подвеса троса, м, H_op="; H_op

INPUT "Введите волновое сопротивление провода, Ом, Z_р="; Z_p

INPUT "Введите импульсное разрядное напряжение при предразрядном времени t₁, кВ, u₁="; u₁

INPUT "Введите импульсное разрядное напряжение при предразрядном времени t₂, кВ, u₂="; u₂

INPUT "Введите радиус троса, мм, r₁="; r₁

CLS

A_ug = ATN(b / H) * 180 / 3.14

PRINT "Величина защитного угла, A_ug="; A_ug; "град"

H_sr = ABS(H_op - 2 * L_prtr / 3)

PRINT "Средняя высота троса над поверхн. земли, H_sr="; H_sr; "м"

H_op1 = H_op - H

H_sr1 = ABS(H_op1 - 2 * L_prpr / 3)

PRINT "Средняя высота провода над поверхн. земли,H_sr1="; H_sr1; "м"

d₁₂ = SQR((H_sr - L_prtr + H_sr1 - L_prpr) ^ 2 + b ^ 2)

d¹₁₂ = SQR((H_sr - L_prtr - H_sr1 + L_prpr) ^ 2 + b ^ 2)

k₂₁ = ABS(LOG(d₁₂ / d¹₁₂) / LOG(2 * H_sr / r1 / 1000))

K_pt = 1.34 * k₂₁

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

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

N_pum1 = 5 * H_op + H_sr ^ 2 / 30

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

N_pum = N_pum1 * L * N_gr / 10000

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

N_pumt =.5 * N_pum

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

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

N_pr = N_pum * 10 ^ (A_ug * SQR(H_op) / 90 - 4)

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

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

z₀ = 600 * Z_p / (600 + Zp) I_mp08 = u₀₅ * 15 / z₀

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

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

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

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

nu = u₁ / u₂: t₁ = 2: t₂ = 20: Uf_max = u / SQR(3): tf(0) = 0

a = ABS((nu ^ 2 - 1) / (1 - nu ^ 2 * t₁ / t₂)): x = 1

r₁ = Rz_i * 600 / (Rz_i + 600): Re = 550 * r₁ / (550 + r₁)

CLS

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

PRINT " T_ф, мкс I_мо, кА U_т, кВ I'_мо, кА "

FOR i = 1 T_O 20

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

c(i) = (R_e * x + (.6 * H_op / tf(i)) * x ^ 2 +.2 * H_op / tf(i)) * (1 - K_pt)

I_mo(i) = (u_t - u * H_sr * (1 - K_pt) - Uf_max) / c(i)

I_moo(i) = I_mo(i) / tf(i)

U_tt(i) = u₂ * SQR(1 + a * t₁ / tf(i)) / (SQR(1 + a * t₁ / t₂))

PRINT TAB(3); tf(i); TAB(13); I_mo(i); TAB(33); U_tt(i); TAB(55); I_moo(i)

NEXT

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

CLS

X₁ = -I_moo(11): x₂ = I_moo(10): x₃ = I_moo(9)

y₁ = I_mo(11): y₂ = I_mo(10): y₃ = I_mo(9)

dx = (x₃ - x₂) / (x₂ - X₁): dy = (y₂ - y₃) / (y₁ - y₂)

I_momin = (X₁ * dx - x₃ * dy) / (dx - dy)

I_moomin = (y₃ * dx - y₁ * dy) / (dx - dy) a1 = (x₂ - I_momin) * (y₂ - I_moomin)

IF I_moomin < 0 THEN I_momin = (X₁ * y₁ - x₂ * y₂) / (y₁ - y₂)

IF I_moomin < 0 THEN I_moomin = 0: a1 = ((x₂ - X₁) / (y₁ - y₂)) * y₁ * _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)="; P₁

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

N_pert =.5 * N_pum * P₂

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

N_perp = N_pum * N_pr * P₃

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

N_perop =.5 * N_pum * P₁

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

N_per = N_pert + N_perp + N_perop

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

P_d = 1.6 * (u / L_per) -.06

IF P_d > 1 THEN P_d = 1

IF P_d <.1 THEN P_d =.1

N_otkl = N_per * P_d

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

N_aw = N_otkl * (1 - k_apv)

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

N_nad = 1 / N_aw

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

N_onad = N_nad * 100 / L * (100 / N_gr)

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

PRINT " end"

END

ЛИТЕРАТУРА.

1. Лифанов В.Н. Электроизоляция и перенапряжения / Учеб. пособие. - Владивосток: Изд-во ДВГТУ, 2003. - С. 51-73.

4. Лифанов В.Н. Координация изоляции электрооборудования /Учеб. пособие. – Владивосток: Изд-во ДВГТУ, 2003. – 120 с.

3. Инструкция по устройству молниезащиты зданий, сооружений и промышленных коммуникаций.- СПб,: Издательство ДЕАН, 2005,- 64 с.

4. Карякин Р.Н. Заземляющие устройства электроустановок. Справочник. М.: ЗАО ЭНЕРГОСЕРВИС, 2000.- 376 с.

5. Лифанов В.Н. Программное обеспечение для расчета заземлителя молниеотводов ОРУ электростанций. - Владивосток.: Изд-во ДВГТУ, Вологдинские чтения, 2000.

6. Лифанов В.Н., Шайдуров И.Г. Электрические испытания изоляции крупных электрических машин. – Владивосток: Изд-во ДВГТУ, Вологдинские чтения, 2003.

Дата добавления: 2014-12-23