э л е к т р о э н е р г е т и к а
УДК 621.315/316.351.064
РАСЧЕТ СИЛ СЖАТИЯ ДИСТАНЦИОННЫХ РАСПОРОК РАСЩЕПЛЕННЫХ ФАЗ ПРИ КОРОТКОМ ЗАМЫКАНИИ
Докт. техн. наук, проф. СЕРГЕЙ И. И., асп. ПОТАЧИЦ Я. В.
Белорусский национальный технический университет E-mail: [email protected]
CALCULATION OF AERIAL CABLE SPACER COMPRESSION FORCES IN BUNDLED PHASES DURING SHORT CIRCUIT FAULT SERGEY I. I., POTACHITS Ya. V.
Belarusian National Technical University
Излагается приближенный метод расчета параметров электродинамической стойкости гибкой ошиновки с расщепленными фазами. Приводится сопоставление результатов приближенного расчета с компьютерными расчетами по программе FAZA. Дается оценка допустимых сил сжатия распорок проектируемых современных типовых распределительных устройств 330 и 500 кВ. Предлагается алгоритм определения параметров динамики расщепленной фазы с учетом влияния упруго-вязких характеристик демпфирующих распорок.
Ключевые слова: динамика токоведущих конструкций, электродинамическая стойкость, максимальные тяжения, силы сжатия, математическое моделирование, расщепленные фазы, дистанционная распорка.
Ил. 3. Табл. 2. Библиогр.: 8 назв.
The paper presents an approximate method for calculation of electrodynamic resistance parameters of a flexible bus-bar with bundled phases. Comparison of the approximate calculation results with the computer calculations using a FAZA software tool has been given in the paper. An evaluation of permissible aerial cable spacer compression forces for designed modern standard distribution substations of 330 and 550 kV has been provided in the paper. The paper proposes an algorithm for determination of dynamics parameters in a bundled phase with due account of influence pertaining to viscoelastic characteristics of damping spacers.
Keywords: dynamics of current-carrying structure, electrodynamic resistance, maximum tension, compression forces, mathematical simulation, bundled phases, расщепленные фазы, aerial cable spacer.
Fig. 3. Tab. 2. Ref.: 8 titles.
В открытых распределительных устройствах (ОРУ) высокого и сверхвысокого напряжения широко применяется гибкая ошиновка с расщепленными фазами. Это позволяет ограничить коронирование проводов и увеличить их пропускную способность. Особенностью поведения расщепленных фаз при коротком замыкании (КЗ) является слипание составляющих фазу
проводов на большей части пролета, вызванное большими токами КЗ. Развитие современных энергосистем сопровождается ростом токов КЗ. На ряде действующих электростанций и подстанций в России и Беларуси они перешагнули рубеж в 40 кА (действующее значение) на стороне 110-330 кВ. При разработке типовых конструкций ОРУ 330 кВ в проектных организациях Российской Федерации в качестве расчетного уже принимается ток КЗ величиной 63 кА. Ожидается рост токов КЗ и в Республике Беларусь в связи со строительством АЭС мощностью 2300 МВт, выдача мощности которой будет производиться на напряжении 330 кВ [1].
При токах КЗ более 40 кА механические усилия и смещения проводов могут оказывать решающее влияние на конструктивное выполнение гибкой ошиновки с расщепленными фазами. В связи с этим возникла необходимость в исследовании динамики расщепленной фазы при КЗ. Указанные исследования в международном формате проводились в рамках СИГРЭ и МЭК по двум направлениям: экспериментальные в опытных пролетах и теоретические. Их результаты обобщены в публикациях рабочей группы В2 «Воздушные линии» СИГРЭ и проекте международного стандарта [2].
Ключевым допущением при оценке динамики расщепленных фаз является ее представление одиночным эквивалентным проводом [3]. Указанная модель позволяет определить размах колебаний расщепленных фаз при междуфазных КЗ, а также второй и третий максимумы тяжения. Оценка первого максимума тяжения производится для «замороженной» фазы в момент максимального стягивания проводов в пучке. Достоверность этого статического метода подтверждена экспериментально.
При больших токах КЗ следует учитывать движение проводов расщепленных фаз, которые оказывают динамическое воздействие на распорки. Приложенные к распоркам силы раскладываются на две составляющие. Одна из них вызывает перемещение распорки как твердого тела, другая направлена на ее деформацию. В настоящее время указанная проблема стала актуальной для проектных организаций, разрабатывающих типовые проекты гибкой ошиновки с расщепленными проводами ОРУ 330 и 500 кВ. Важным конструктивным параметром расщепленной фазы является расстояние между распорками (длина подпролета). Его величина в значительной степени влияет на силы сжатия распорки при КЗ.
В БНТУ для математического описания динамики расщепленных проводов используются уравнения движения гибкой упругой нити [3]. Они справедливы для участков проводов между распорками, на которые воздействуют при КЗ только распределенные электродинамические усилия (ЭДУ). Из-за фиксации проводов в зажимах распорок при КЗ на провода в этих точках воздействуют также сосредоточенные силы, величины которых могут быть определены лишь с учетом упругих и инерционных характеристик распорок. При многофазных КЗ из-за разницы в расстояниях между составляющими различных фаз траектории движения проводов становятся несимметричными. Это ведет к повороту распорок и закручиванию фазы. Математическое описание динамики расщепленных фаз на основе принятой модели включает в себя уравнения поступательных перемещений и закручивания проводов и распорок, которые решаются совместно [3].
Жесткая распорка как неизменяемая система, с точки зрения теоретической механики, является абсолютно твердым телом. По определению, абсолютно твердое тело не может подвергаться никаким деформациям и представляет собой идеальный образ, который тем ближе подходит к реальной распорке, чем меньше последняя способна деформироваться под действием сил, приложенных от проводов. Подавляющее большинство распорок, установленных еще в энергосистемах бывшего СССР, можно отнести к распоркам жесткого типа [3].
С учетом сказанного в статье проведены исследования с использованием расчетной модели фазы в виде пучка гибких упругих нитей, которые в заданных точках пролета соединяются между собой жесткими дистанционными распорками. Такая модель учитывает геометрические и физико-механические характеристики расщепленных фаз в целом пролете. В момент максимального стягивания проводов в пучке под действием внутри-фазных ЭДУ происходит резкое увеличение тяжения фаз. Указанный пик тяжения был назван первым максимумом [3].
Рассматриваемый численный метод расчета параметров электродинамической стойкости расщепленной фазы пригоден для оценки электродинамической стойкости расщепленных фаз при любом взаимном расположении составных проводников для различных, в том числе и неравномерно установленных, распорок [4]. Однако расчет по вышеуказанному методу трудоемок и требует высокой квалификации пользователя. Поэтому для инженерных расчетов был модифицирован упрощенный метод расчета, основанный на рассмотрении расщепленной фазы как статической системы, к которой прикладываются максимальные электродинамические усилия при КЗ [5].
В статье рассматривается расчет сил сжатия, действующих на распорку, для фазы, состоящей из двух проводов. Расчетные формулы выводятся для положения равновесия провода. На основании условия упругой деформации растяжения прилегающих к распоркам участков (рис. 1) провода по закону Гука получена следующая формула для определения Нтах [5]:
Н - Н + ЕА
Н тах - Н 0 + —
5 (1 - ШЭ Ртах ) Ч3
8Ш Ртах 24Н2
(1)
где Нтах - максимальное тяжение проводов при их схлестывании; Н0 - начальное тяжение провода; 5 - шаг расщепления; I - длина подпролета, м; Е - модуль упругости, Н/мм2; А - поперечное сечение провода, мм2.
Расчет силы сжатия распорки производили по очевидному выражению
Ртах - 2Нтах18Ртах, (2)
где Ртах - максимальная сила сжатия, действующая на распорку, Н.
При допущении, что в момент максимального стягивания проводов в фазе имеет место равновесие приложенных к ним сил и моментов этих сил, получено следующее уравнение:
я 0 + f
S (1 - cos Pmax ) q_l3
sin P„
24Я2
tgPmax -
- i__ (cos Pmax )_ ^ [ K + S - d cos Pmax ][ K + S COS_Pmax - d ] = 0 10 Sin Pmax П Sd (_Sin_Pmax + 1 - COS_P,ax )_ '
где K = y]d2 + s2 -2dscos2pmax; d - удвоенный диаметр провода, м.
O
я
E
Рис. 1. К определению суммарного момента при предельном стягивании проводов расщепленной фазы
Трансцендентное алгебраическое уравнение с одним неизвестным ртах представляет собой (3). В [5] определение угла производится без решения (3) путем построения расчетных кривых, обобщающих расчеты по компьютерной программе, в зависимости от сочетания нескольких параметров: ток КЗ, шаг расщепления и удельный вес провода на единицу дли-
( т \
ны Pmax = f
4sq
что является частным случаем поставленной задачи.
Угол Ртах определяли из уравнения (3), которое решали методом половинного деления [6], сводящим процесс нахождения угла Ртах к выполнению конечного числа итераций и дающим результат в виде числового значения с погрешностью, приемлемой для решаемой задачи. Данный подход основан на последовательном сужении интервала, содержащего единственный корень (3) до того времени, пока не будет достигнута заданная точность в. Предварительно необходимо определить область локализации корня данного уравнения. Из анализа опытных данных установлено, что корень (3) находится на отрезке [0; п/4]. Разобьем отрезок [0; п/4] пополам. Определим новое приближение корня ртах в середине отрезка [0; п/4]: [0; П/4]
Ртах = —^—. Найдем значение, которое примет (3) в точках 0 и ртах. Если
произведение полученных чисел меньше нуля, то корень (3) находится на отрезке [0; п/8], если больше нуля, то на отрезке [п/8; п/4]. Область с корнем вновь делим пополам. Алгоритм выполняется до тех пор, пока значение ртах не будет найдено с необходимой точностью (рис. 2). При изменении конструктивных параметров расщепленных фаз потребуется корректировка угла ртах, что не представляет трудностей для проведения расчетов.
Начало
Конец
Рис. 2. Алгоритм определения ртах
Разработанный алгоритм был реализован в компьютерной программе. После нахождения ртах выполняли расчет Нтах, Ттах и Ртах. Оценка достоверности расчета Нтах и Ртах была выполнена при анализе параметров электродинамической стойкости гибкой ошиновки с расщепленными фазами при токах КЗ от 20 до 63 кА. Сопоставление результатов расчета, полученных по разработанному приближенному алгоритму и в компьютерной программе БЛ2Л, с опытными данными приводится в табл. 1.
Таблица 1
Сравнение результатов расчета Т1тах по предложенному алгоритму и программе ^АХА с опытными данными
/к, кА Марка провода Т1тгве кН опыт Т1тгве кН ЕЛ7Л ^тое кН метод Погрешность %
Дметод ДГЛ7Л
22,5 2хАБТЕЯ570 63,0 65,0 62,1 1,42 3,14
41,0 2хАБТЕЯ570 66,4 62,7 61,3 7,68 5,52
59,8 2хАБТЕЯ570 69,1 71,6 67,7 2,01 3,64
60,2 2хАБТЕЯ570 75,4 81,8 73,8 2,12 8,53
60,2 2хРЕТШ1Л612 81,8 105,6 79,4 2,93 29,1
60,2 2хРЕТШ1Л612 81,8 85,0 80,3 1,83 3,87
56,9 2хРЕТШ1Л612 78,2 81,9 77,7 1,28 4,71
Разработанный метод использован для расчета параметров электродинамической стойкости расщепленных фаз гибкой ошиновки проектируемых типовых открытых распределительных устройств напряжением 330-500 кВ в Российской Федерации. В результате этого расчета определялись силы сжатия дистанционных распорок. При сравнении результатов находили силы сжатия распорок по формуле Манузо, широко используемой в проектной практике. Результаты этого расчета представлены в табл. 2.
Таблица 2
Сравнение результатов расчета сил сжатия распорок Ртлх расч
I (3) кА Марка провода Р Н тах расч; ^^ ЕЛ7Л Тип распорки Р Н тах расч; ^^ метод 1доп, м Р Н тах; и (по Манузо)
20 2хАС 300/39 1920 РГ-2-400 1400 9,830 2574
20 2хАС 400/22 1800 РГ-2-400 1760 9,830 2635
20 2хАС 500/27 1600 РГ-2-400 1680 9,830 2659
20 2хПА 500 1720 РГ-2-400 1650 9,830 2929
40 2хАС 300/39 1920 РГ-2-400 1850 1,300 3022
40 2хАС 400/22 1760 РГ-2-400 1890 1,300 3076
40 2хАС 500/27 2100 РГ-2-400 1950 1,300 3100
40 2хПА 500 2120 РГ-2-400 2250 1,600 3466
63 2хАС 300/39 2200 РГ-2-400 2340 0,700 3490
63 2хАС 400/22 2350 РГ-2-400 2420 0,700 3657
63 2хАС 500/27 2150 РГ-2-400 2500 0,700 3759
63 2хПА 500 2520 РГ-2-400 2670 0,756 4654
Как видно из табл. 1, 2, приближенный метод расчета сил сжатия распорок и максимальных тяжений, действующих на дистанционную распорку, для фазы, состоящей из двух проводов, дает хорошие результаты, которые подтверждаются численными расчетами по компьютерной программе БЛ2Л. При токах КЗ 40-60 кА электродинамическая стойкость расщепленных фаз обеспечивается только при использовании демпфирующих распорок с повышенной механической прочностью. Следует отметить, что в компьютерной программе БЛ2Л и в приближенном методе не учитывали подвижность лучей демпфирующей распорки относительно рамы, что давало погрешность в результатах расчета сил сжатия распорки и макси-
мальных тяжений. Для оценки величины этой погрешности необходимо усовершенствовать компьютерную программу FAZA и учесть в расчетах особенности конструкции демпфирующей распорки. Демпфирующая распорка (рис. 3) выполняется в виде корпуса 1, подвижных лучей 2, которые с одного конца заканчиваются плашками-зажимами 3, а с другого - закреплены в демпферных узлах 4. Они состоят из эластомерных дисков, стянутых с определенным усилием болтом с гайкой. Эластомерные диски изготовлены из токопроводящей резины с высоким содержанием углерода [7].
Рис. 3. Демпфирующая распорка типа РД
Наличие в конструкции распорки-гасителя демпферного узла с эласто-мерными дисками, способными к энергопоглощению, не позволяет представлять ее одним абсолютно твердым телом. С учетом этого лучи распорки при математическом моделировании могут быть представлены абсолютно твердым телом, совершающим только вращательное движение вокруг точки упруго-вязкого шарнирного опирания на раме. Твердое тело рамы производит сложные вращательные и поступательные перемещения в пространстве[8].
Вращение /-го луча распорки вокруг точки его упруго-вязкого шарнирного опирания в плоскости, перпендикулярной пролету, приближенно можно описать следующим уравнением [8]:
(2а1 ( (а/ - аъ )
^+/-4^+с (а'- аъ У- М • (4)
где 3 - момент инерции луча относительно оси х; а■ - угол поворота/-го луча в упруго-вязком шарнире относительно рамы; с, / - угловая упругая и вязкостная жесткости узла демпфирования; аъ - динамический угол поворота рамы.
На основе предложенного алгоритма разработан модифицированный вариант компьютерной программы FAZA, учитывающий влияние подвижности лучей распорок.
В Ы В О Д
Разработан приближенный метод расчета максимальных тяжений и сил сжатия распорок в расщепленной фазе, который базируется на решении уравнения относительно угла ртах методом половинного деления. Выполнена оценка достоверности расчета максимального тяжения и сил сжатия
распорок при коротком замыкании с использованием компьютерной программы FAZA. Определены параметры электродинамической стойкости расщепленных фаз проектируемых типовых открытых распределительных устройств напряжениями 330 и 500 кВ. Предложен алгоритм определения параметров динамики с учетом влияния упруго-вязких характеристик демпфирующих распорок.
Л И Т Е Р А Т У Р А
1. Н а ц и о н а л ь н ы й интернет-портал Республики Беларусь [Электронный ресурс] / Официальный сайт Министерства энергетики Республики Беларусь. - Режим доступа: http: //www.minenergo.gov.by/ru/at_en/r_aes. - Дата доступа: 25.01.2014.
2. C I G R E Study Committee 23. Working group 23 - 11 ESCC Task Force. - Paris: CIGRE Central Office, 1996.
3. С е р г е й, И. И. Динамика проводов электроустановок энергосистем при коротких замыканиях: теория и вычислительный эксперимент / И. И. Сергей, М. И. Стрелюк. -М.: ВУЗ-ЮНИТИ, 2001. - 151 с.
4. С т р е л ю к, М. И. Расчет электродинамической стойкости гибкой ошиновки открытых распределительных устройств высокого напряжения / М. И. Стрелюк, И. И. Сергей // Электричество. - 1984. - № 1. - С. 10-14.
5. С е р г е й, И. И. Упрощенный расчет максимальных тяжений расщепленных проводов при коротком замыкании / И. И. Сергей, А. П. Андрукевич // Энергетика... (Изв. высш. учеб. заведений и энерг. объединений СНГ). - 2008. - № 6. - С. 18-23.
6. Д е м и д о в и ч, Б. П. Основы вычислительной математики / Б. П. Демидович, И. А. Марон. - М.: Наука, 1968. - 664 с.
7. Д и с т а н ц и о н н а я распорка-гаситель колебаний: а. с. 1651335 СССР, МКИ 5Н02Н3/16 / А. А. Виноградов, А. Ф. Голубев, Я. А. Кловский, С. Э. Шиков, Е. И. Виноградова, И. И. Сергей; МО СКТБ по изоляторам и арматуре ВПО «Союзэлектросетьизо-ляция». - № 435 3987/07; заявл. 31.12.1987; опубл. 23.05.1991 // Открытия. Изобретения. -1991. - № 19. - С. 207.
8. D y n a m i c Loads on Spacer Damfer at Short - Circuit / M. I. Streliuk [et al.] // Short -Circuit Currents in Power Systems: Proceedings of the 7th International Symp. 10-12 Sept. 1996, Technical University of Lodz. - Warsaw, 1996. - № 2.7. - P. 2.7.1-2.7.5.
R E F E R E N C E S
1. N a t i o n a l Legal Internet Portal of the Republic of Belarus. Official site ofMinistry of Energy of the Republic of Belarus. Available at: http://www.minenergo.gov.by/ru/at_en/r_aes. (accessed 25 January 2014).
2. CIG R E Study Committee 23. Working group 23 - 11 ESCC Task Force. Paris: CIGRE Central Office, 1996.
3. S e r g e y, I. I., & Streliuk, M. I. (2001) Dynamics of Electrical Installation Wires in Electric Power Systems During Short Circuit Faults: Theory and Calculative Experiment. Moscow: VUZ-UNITI.
4. S t r e l i u k, M. I, & Sergey, I. I. (1984) Calculation of Electrodynamic Resistance of Flexible Bus-Bar in Outdoor Switch-Gears of High Voltage. Elektrichestvo [Electricity], 1, 10-14.
5. S e r g e y, I. I., & Andrukevich, A. P. (2008) Simplified Calculation of Maximum Tensions of Split Conductors at Short Circuit. Izvestiia Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii - Energetika. [Proceedings of the Higher Education Institutions and Power Engineering Associations - Power Engineering], 6, 18-23.
6. D e m i d o v i c h, B. P., & Maron, I. A. (1968) Fundamental of Numerical Mathematics. Moscow: Nauka.
7. V i n o g r a d o v, A. A., Golubev, A. F., Klovskii, Ia. A., Shikov, S. E., Vinogradova, E. I., & Sergei, I. I. (1991) Distantsionnaia Rasporka - Gasitel' Kolebanii [Spacer-Damper]. Certificate of authorship no. 1651335 USSR. Otkrytiia. Izobreteniia [Discoveries. Inventions], (19), 207.
8. S t r e l i u k, M. I., Sergei, I. I., Bladyko, Iu. V., & Vinogradov, A. A. (1996) Dynamic Loads on Spacer Damper at Short - Circuit. Short - Circuit Currents in Power Systems. Proceedings of the 7th International Symp. 10-12 Sept. 1996. Technical University of Lodz. Warsaw, (2.7), 2.7.1-2.7.5.
Представлена кафедрой
электрических станций Поступила 14.02.2014