Метод расчета электромагнитных характеристик вертикальных стержневых заземлителей при стекании импульсных токов молнии Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.316.99

МЕТОД РАСЧЕТА ЭЛЕКТРОМАГНИТНЫХ ХАРАКТЕРИСТИК ВЕРТИКАЛЬНЫХ СТЕРЖНЕВЫХ ЗАЗЕМЛИТЕЛЕЙ ПРИ СТЕКАНИИ ИМПУЛЬСНЫХ ТОКОВ МОЛНИИ

Канд. техн. наук ГЕРАСИМОВИЧ Д. А., магистр техн. наук ДЕРЮГИНА Е. А.

ООО «Электросити», Белорусский национальный технический университет

Исследованию импульсных характеристик заземляющих устройств в последнее время уделяется большое внимание. Связано это, в первую очередь, с проблемой электромагнитной совместимости современных электроэнергетических объектов. Здесь, наряду с традиционными требованиями к защитному заземлению и заземлению грозозащиты, возникают вопросы выноса потенциала на цепи вторичной коммуникации.

При исследовании импульсных характеристик заземляющих устройств необходимо принимать во внимание ряд дополнительных факторов, связанных с большими величинами токов молнии и их быстрым изменением во времени. Большие величины токов требуют учета нелинейности вольт-амперных характеристик грунтов и пробоя грунта при достижении критической напряженности электрического поля в земле, а быстрое изменение токов молнии - электромагнитной составляющей и токов смещения в уравнениях поля (индуктивных и емкостных параметров заземлите лей).

Одновременный учет перечисленных выше факторов является сложной задачей. Следует выделить исследования [1-3] импульсных характеристик простых заземлителей с учетом пробоя и нелинейности характеристик грунтов. Для анализа импульсных характеристик в них предлагается использовать фиктивный радиус заземлителя с учетом зоны искрообразова-ния в грунте и импульсные коэффициенты для расчета сопротивлений. Исследования, учитывающие электромагнитные процессы в земле, базируются в основном на полевой модели [4-6] и модели линии с распределенными параметрами [1-3, 7-9]. Методы расчета импульсных характеристик на основе данных моделей разработаны для частных случаев протяженных заземлителей [1-3, 7], горизонтальных [8, 9] и вертикальных [4, 6] заземли-телей конечной длины, а также для сетчатых заземлителей [5, 8]. Удельные параметры схемы замещения (индуктивность и емкость) в моделях линий с распределенными параметрами при этом взяты как для поля постоянного тока, что некорректно для быстро изменяющегося во времени электромагнитного поля в проводящей среде. Для учета пробоя и нелинейности вольт-амперных характеристик земли в этих моделях предлагается использовать импульсные параметры земли с учетом нелинейности и радиуса искрооб-разования [1-3, 8, 9]. Соответствующие уравнения при этом становятся нелинейными. Их решение трудноосуществимо и выполняется с рядом допущений.

Наряду с указанными расчетными методами в последнее время получили распространение экспериментальные исследования [10, 11]. Для генерации импульсных токов используют генераторы апериодических импульсов

с заданными параметрами, а измерение напряжения на заземляющем устройстве выполняют импульсным вольтметром или цифровым осциллографом. Однако мощности генераторов импульсов не всегда достаточно для воспроизведения пробоя грунта, что может иметь место при реальном ударе молнии.

Анализ существующих методов расчета импульсных характеристик заземляющих устройств показывает, что универсальные и в то же время практически пригодные методы, учитывающие всю сложность протекающих в земле электромагнитных процессов, отсутствуют. В практике эксплуатации заземляющих устройств при этом необходимо использовать экспериментальные методы оценки импульсных характеристик. Дальнейшее развитие расчетных методов необходимо для качественной оценки импульсных характеристик заземляющих устройств и разработки практических рекомендаций для выбора конструкций заземляющих устройств на стадии проектирования.

Ниже предлагается метод расчета импульсных характеристик вертикального стержневого заземлителя конечной длины (рис. 1).

Рис. 1. Геометрическая модель вертикального стержневого заземлителя конечной длины

В [12] получены выражения для распределения параметров электромагнитного поля в земле и заземлителе в первом приближении при гармоническом изменении тока с угловой частотой ю. Распределение комплексной амплитуды радиальной составляющей плотности тока в земле описывается выражением [12]

62тг(Г,2) =

а21т 2%Г

е 2 +

I ( еЬ(а2 (2 +/ )^8Ь(а2/) - еа22)

л/7 2 + (г - *2)2

, г е[ ^2,да); 2 е[0, -/ ];

(1)

1 -

ф2 + (г - Я2)2

г е [Я2,да), 2 е (-/,-да),

где а2 = ^7'юу2ц2; у2, ц2 - электрическая проводимость и магнитная проницаемость земли соответственно; I т - комплексная амплитуда тока, стекающего в заземлитель; / - длина заземлителя; Я2 - внешний радиус за-землителя.

Аналогичное выражение справедливо и для комплексной амплитуды угловой составляющей напряженности магнитного поля в земле

Н2тф^ 2) =

2%г

е"2' +

/Цех2(2 + /)у5Ь(а2/) -е"22), г е [^да), ^ е [0,-7]; ф2 + (г - Д2)2

(2)

1 -

V/2 + (г - Д2)2

г е [Л2, да), 2 е (- /,-да).

На основе (1), (2), а также аналогичных выражений для распределения параметров поля заземлителя в [12] получены приближенные формулы для расчета сопротивления растеканию тока в земле и сопротивления стержня заземлителя. Исследования показали, что первое приближение на основе (1), (2) дает незначительные отличия в сравнении с более точным численным методом расчета параметров поля для промышленной частоты. Зависимости сопротивлений при более высоких частотах, рассчитанные по приближенным формулам и более точно, иллюстрирует рис. 2. Их анализ показывает, что даже для частот порядка 107 Гц отличия не превышают 2 % (кривые 1-2 практически совпадают).

Таким образом, для дальнейшего расчета параметров импульсного электромагнитного поля можно основываться на первом приближении распределения параметров поля. Расчет параметров поля в импульсном режиме при этом целесообразно выполнить операторным методом. Для этого запишем (1), (2) в операторной форме:

СТ2Г (г, 2, р) =

а2 ( Р)1 (Р)

2%г (1 + 82р/У2 )

еа2( р)г +_

/1 сЬ(а2(р)(2 + 7)) _

_ еа2(Р)2

8Ь(а2( р)/)

ф2 + (г - Л2)2

г е [ ^2, да);

2 е

[0,-/]; (3)

,а2( р)г

1-

ф2 + (г - Л2)2

г е [^2,да); 2 е (-/,-да);

/

е

т

/

а^ 2

е

/

Н2Ф (г, г, р) =

I (Р)

2%г

/ \ 5Ка2(р)(г + 1)) _ _ „а2(р)г ,02(Р)г , I «Ь(а2(Р)1)_

V/2 + (г - П2у

г е[ Я2,да); г е[0, -/];

(4)

,а2( Р) г

1 _

V/2 + (Г -Я:)2

г е [ К2, да); г е (-/, -да),

где а2г (г, г, р), Н2ф(г, г, р) - операторные изображения радиальной составляющей плотности тока и угловой составляющей напряженности магнитного поля в земле соответственно; а2(р) = -^у2ц2р + в2ц2р2; I(р) -операторное изображение тока, стекающего в заземлитель; в2 - диэлектрическая проницаемость земли.

40 п

К,

Ом 30 -

25 -

/ = 30м / = 10м

10000 -Ки, Ом " 8000 -

7000 -

6000 -

5000 -

4000 -

3000 -

2000 -

1000 -

0 -

а

у 2 = 0,1 См/м

/ .Гц

б

— п ппги гчт/л..

/, Гц

Рис. 2. Зависимость сопротивления Ки медного заземлителя радиуса К2 = 6 мм от частоты / стекающего тока для: а - большой; б - малой проводимости земли; 1 - численный метод [12]; 2 - приближенные формулы; 3 - импульсный метод без учета токов смещения; 4 - импульсный метод с учетом токов смещения

Выражение (3) записано с учетом закона полного тока при наличии токов смещения в земле (сНф/& = (1 + в2р/у2 )сг).

/

0

Нахождение мгновенных значений плотности тока и напряженности магнитного поля в земле сводится к построению соответствующих переходных функций и последующему использованию интеграла Дюамеля. Изображения переходных функций получаются при подстановке в (3), (4) изображения единичной ступени тока I(р) = 1/р [13]. Тогда оригиналы можно записать следующим образом:

&, (г, гГ) = ^ 1

( Г) , 1 (^ (г, г) - (г, г)) г _ ) [п Я( > .

Я1аг (г,г) + ^—, 2 2—г е [Я2,да); г е [0,-I], Г > 0;

(Г)

1 -

ф2 + (г - И2у

ф2 + (г - Я2)2

, г е [ ^2,«); г е (-/, -да);

(5)

ЯЯ2Ф (г,0=

. , I (&яф (г,Г) - Я1Щ(г,Г)) г ) Гп , (г,Г) + ^---г е [Я2,да); г е [0,-I], Г > 0;

ф2 + (г - Д2)2

(6)

Я1Нф(гГ)

1 -

ф2 + (г - Д2)2

г е[Я2,да); 7 е(-1,-да),

где оригиналу (г, Г) соответствует изображение (г, р) =

' ' р{\ + е2р/у2)8Ца2(р)1)

«2( р)е'

«2( р) г

р(1 + В2 р/У2 )'

«2(р)г

Оригиналы я1аг (г, г) и я1Н (г, Г), учитывая, что а2(р) = ^/у2ц2р + в2ц2р2 определяются по формулам [13]:

Я1аг (г, г) =<

(

Уи1 ^

у2фГ -В2Ц2г' 2в2

2 Л У2

г > -гу]в2ц2;

(7)

Я1Нф(г Г) = <

0, г < - г у/в2Ц2 ;

^ Л2 2\ Л

У2 (г -В2Ц2г )

11

У2 г М2 " В2

Г /

2в2

2в2

- ^ В2Ц2

72 2 г -В2Ц2г

х + е

У2 г Шя 2\ В2

Г > - ^ В2ц2;

(8)

0, Г < в2Ц2 ,

I

I

2

е

е

где 10, II - модифицированные функции Бесселя нулевого и первого порядков соответственно.

Для нахождения оригиналов g2a (г, t) и g2н (г, t) используем теорему

разложения [13]. Итоговые выражения для их определения имеют вид:

g2a, (г, t) =<

(

У2,

1 - *52 +1

2с08 ( к №¡1 )(

е — е т<:

\

к=1

— 4к2 я^/ (^2 у2/2)

t > — г-у/82^2 ;

(9)

g2Hф (г,0 =

0, t < — г78гц7;

(

г +1

I

к=1

(—1)к

(—1)к

+ ^ кп + Щ^ Р2к 2к% 'Р1к 2к% р2

81П

( к%( г +1 ^ .-

4 / t > — г^82Ц2;

I

(10)

0, t < — гу!82Ц2,

где Р1к =

_У2_ 28 2

1 —

II

, 4к2Л287 1--2

^2у21 2

; Р2к =

_У2_ 282

1 +

, 4к2Л287

1--2

2

Зная переходные функции, находим мгновенные значения проекций векторов плотности тока и напряженности магнитного поля в земле с использованием интеграла Дюамеля [13]:

I

, (г, г, t) = | /'(т)ga2г (г, г, t — т^т;

(11)

Н 2ф (г, г, t) = |/'(т) gн 2ф (г, г, t — т)^т

(12)

где /'(т) - первая производная по времени тока, стекающего в заземлитель.

Интегралы (11), (12) записаны с учетом равенства нулю тока в начальный момент времени, /(0) = 0.

Таким образом, на основании выражений (5)-(12) рассчитываются параметры импульсного электромагнитного поля в земле.

Для расчета поля самого заземлителя учтем, что поверхностный эффект в нем при стекании импульсного тока ярко выражен, а токами смещения в металле можно пренебречь. Тогда электромагнитное поле заземлителя можно рассматривать как для полубесконечного проводящего пространства и операторные изображения для вертикальной проекции вектора плотности тока и угловой проекции вектора напряженности магнитного поля примут вид:

в1г(г,г,р) = аДр)Н2фЕг,р)еа'(р)(г—Е), ге[0, —I]; ге[Д,Е2]; (13)

Н\ф(г,г,р) = Н2ф(Е2,г,р)еа'(р)(г—Е), ге[0, —I]; ге[Д,Е], (14)

)

I

где а

1(р) = у1умр; Ух,

- электрическая проводимость и магнитная

проницаемость стержня заземлителя соответственно.

Переходные функции для этих проекций определяются по формулам [13]:

Лг,2,*)= <

gHгv (Кг,2,*Ы^е

У1№ (Кг-г)

4*

Л*

* > -2^;

(15)

0, * < -2у18гЦг;

gнlф (г,2, *) = <

gHгФ (Кг, 2,*)

( (

1 -ф

V V

Ут1м1 ( Кг - г ) г41

\\

, * > -2у18гЦг;

(16)

0, * < -ВгЦг ■

где Ф - функция Лапласа.

Мгновенные значения проекций векторов, как и для земли, определяем на основе интеграла Дюамеля:

I

Т1г (г, 2, *) = | /'(т)galz (г, 2, * - т)¿т;

(17)

Мф

н ф (г, 2, *) = | /'(т)gнlф (г, 2, * - т)Л.

(18)

По известным мгновенным значениям проекций векторов электромагнитного поля в земле и заземлителе находим импульсное сопротивление и напряжение на заземлителе. Как и для стационарного режима, расчет импульсных сопротивлений основан на применении вектора Пойнтинга. Импульсное сопротивление растеканию тока в земле рассчитываем по формуле

Кз.и (*) = Г гагг (г, 0, *)Нгф (г, 0, *)йг,

Уг/ (*) Кг

где К - расстояние до заземлителя с обратным током. Импульсное сопротивление стержня заземлителя

Кс.и (*) = ] (Кг, 2, * )Н1ф (Яг, 2, * )а2.

Общее сопротивление

Ки(*) = Кз.и(*) + Кс.и(*). Напряжение на заземлителе в импульсном режиме

К

и) = — |стгг (г,0, *)^г.

Уг Кг

(19)

(г0)

(г1)

(гг)

На основании приведенных выражений в среде программирования MatLab разработана программа расчета импульсных характеристик вертикальных стержневых заземлителей. Программа позволяет выполнять расчеты при любой зависимости тока, стекающего в заземлитель. Для ее апробации были проведены расчеты установившихся значений сопротивлений при стекании синусоидального тока различной частоты f = ю/2п (кривые 3, 4, рис. 2). Свойства земли и заземлителя приняты, как и в [12]. Диэлектрическая проницаемость земли s2 = 10s0 = 8,86 • 1011 Ф/м. Из приведенных зависимостей следует, что значения сопротивлений, рассчитанные методом [12] для установившегося режима, по приближенным формулам и разработанной программе для импульсного режима без учета токов смещения практически совпадают (отличия не превышают 5 %). Аналогичные результаты были получены для заземлителей других радиусов, а также для стальных заземлителей.

Следует отметить влияние токов смещения на частотную характеристику сопротивлений заземлителей. При большой проводимости земли влияние токов смещения несущественно (кривые 3, 4 на рис. 2а совпадают). С увеличением частоты сопротивление значительно возрастает за счет электромагнитной (индуктивной) составляющей поля в земле. При малой проводимости земли сопротивление, наоборот, резко падает с увеличением частоты за счет токов смещения (кривые 4, рис. 2б). Кроме того, для больших частот характерно колебательное изменение сопротивления. Аналогичные результаты были получены в [4] для вертикального заземлителя на основе модели антенны, излучающей в диссипативной среде (земле).

В Ы В О Д Ы

1. Предложен метод расчета импульсных характеристик вертикального стержневого заземлителя, основанный на решении уравнений переходного режима операторным методом при приближенном распределении параметров электромагнитного поля.

2. Разработана программа расчета импульсных характеристик и произведена ее апробация применительно к синусоидальным токам различной частоты.

Л И Т Е Р А Т У Р А

1. Р я б к о в а, Е. Я. Заземление в установках высокого напряжения / Е. Я. Рябкова. -М.: Энергия, 1978.-224 с.

2. Б у р г с д о р ф, В. В. Заземляющие устройства электроустановок / В. В. Бургсдорф, А. И. Якобс. - М.: Энергоатомиздат, 1987. - 400 с.

3. А н н е н к о в, В. З. Протяженные заземлители молниезащиты в грунтах с нелинейными вольт-амперными характеристиками / В. З. Анненков // Электричество. - 2001. -№ 7. - С. 22-29.

4. P o l j a k, D. Wire antenna model for transient analysis of simple grounding systems, part I: the vertical grounding electrode / D. Poljak, V. Doric // Progress in Electromagnetic Research. - 2006. - Vol. PIER 64. - P. 149-166.

5. G r c e v, L. D. Frequency dependent and transient characteristics of substation grounding systems / L. D. Grcev, M. Heimbach // IEEE Trans. Power Del. - 1997. - Vol. 12, № 1. - P. 172-178.

6. O l s e n, R. G. A comparison of exact and quasi-static methods for evaluating grounding systems at high frequencies / R. G. Olsen, M. C. Willis // IEEE Trans. Power Del. - 1996. -Vol. 11, № 2.-P. 1071-1081.

7. К а р я к и н, Р. Н. Удар молнии в протяженный неоднородный заземлитель / Р. Н. Карякин, В. Г. Лосев // Электричество. - 2009. - № 4. - С. 12-18.

8. L i u, Y. An engineering model for transient analysis of grounding system under lighting strikes : nonuniform transmission-line approach / Y. Liu, N. Theethayi, R. Thottappillil // IEEE Trans. Power Del. - 2005. - Vol. 20, № 2. - P. 722-730.

9. V e l a z q u e s, R. Analytical modeling of grounding electrodes transient behavior / R. Ve-lazques, D. Mukhedkar // IEEE Trans. Power Apparatus and System. - 1984. - Vol. PAS-103, № 6.-P. 1314-1322.

10. Э к с п е р и м е н т а ль н ы е исследования заземлителей молнезащиты при воздействии импульсных токов / Р. К. Борисов [и др.] // Электро. - 2004. - № 1. - С. 13-16.

11. Э л е к т р о м а г н и т н ы е характеристики заземляющих устройств опор высоковольтных линий при измерениях импульсным методом / В. Н. Боронин [и др.] // Известия РАН. Энергетика. - 2009. - № 2. - С. 11-20.

12. Г е р а с и м о в и ч, Д. А. Расчет электромагнитных характеристик уединенных вертикальных стержневых заземлителей конечной длины при стекании переменного тока промышленной частоты / Д. А. Герасимович, Е. А. Дерюгина // Энергетика... (Изв. высш. учеб. заведений и энерг. объединений СНГ). - 2011. - № 2. - С. 22-29.

13. А н г о, А. Математика для электро- и радиоинженеров / А. Анго; пер. с фр. -М.: Наука, 1965.-780 с.

Представлена кафедрой

электрических станций Поступила 30.05.2011

УДК 621.313.333

АСИНХРОННЫЙ ЭЛЕКТРОДВИГАТЕЛЬ СО ВСТРОЕННЫМ КОМБИНИРОВАННЫМ ТОРМОЗНЫМ УСТРОЙСТВОМ НА БАЗЕ ЭЛЕКТРОМЕХАНИЧЕСКОГО ТОРМОЗА И ЭЛЕКТРОМАГНИТНОЙ МУФТЫ

Канд. техн. наук, доц. СОЛЕНКОВ В. В., канд. техн. наук БРЕЛЬ В. В.

Гомельский государственный технический университет имени П. О. Сухого

Электродвигатели с тормозными устройствами широко используются в крановом электроприводе, конвейерах, лифтах, эскалаторах, электроприводе подач и других электроприводах, требующих быстрой и точной остановки. При этом электромеханические тормоза является относительно простым и дешевым средством для торможения электроприводов на базе асинхронных двигателей (АД) [1-3].

Они обеспечивают:

• наиболее быстрое торможение по сравнению с другими известными способами;