
книги из ГПНТБ / Геллер Б. Импульсные процессы в электрических машинах
.pdfТогда остается |
выражение |
k |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
(тк)2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
пгк |
|
|
|
|
|
|
|
|
|
|
|
|
2п |
^ |
k |
|
|
|
|
(10-148) |
|||
|
|
|
|
|
~2п |
2п sin • 2п |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
ci2 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
По (10-148) были рассчитаны относительные откло |
||||||||||||||
нения |
величин |
амплитуд а= (Ат—Amw)/Am |
|
для значе- |
|||||||||||
% а |
|
ч |
|
|
|
1 |
|
% |
|
|
|
|
|
||
а |
|
|
|
// |
|
|
|
|
|
|
|||||
|
|
|
1 |
|
|
|
|
60 |
|
|
|
|
|
||
30 |
* і |
|
|
|
|
; |
|
|
|
|
|
|
|||
|
|
|
|
|
|
1 і |
|
|
|
|
|
|
|
||
|
|
ч! |
|
|
|
|
50 |
\ |
|
н |
1_ |
|
|||
|
|
|
|
|
11 |
|
|
|
|
|
|||||
20 - J zI |
|
|
|
11 |
|
W |
|
|
|
|
|
||||
|
|
|
11 |
|
|
\ |
|
cl2 |
100 |
|
|||||
|
|
|
hi |
|
|
11 |
|
|
|
|
|
||||
|
|
|
/ Л > |
|
|
|
|
|
|
|
В |
|
|||
70 Щ// hj |
|
|
|
|
|
|
|
|
|
||||||
|
|
J |
|
/% |
|
|
|
30 |
|
|
|
|
|
||
|
|
s |
|
|
|
m |
|
20 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
А |
£ |
|
|
0 |
|
2 |
|
4> |
|
6 |
8 |
10 |
|
10 |
в' \ |
|
|||
|
|
|
|
к |
|
||||||||||
Рис. |
Ю-17. |
Относительное от |
|
|
2 10 |
|
|||||||||
клонение |
амплитуд |
в зависи |
|
|
|
|
|
|
|
||||||
мости |
от |
номера |
простран |
|
О |
|
|
|
S |
10 |
|||||
ственной |
гармоники. |
|
|
|
Рис. 10-18. Максимально |
||||||||||
Обозначения — см. рис. 10-16. |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
возможное |
отклонение |
||||
ний |
/г/с/2 |
= 1/100 |
и |
/г/с/2 =1/10 |
амплитуд |
в |
зависимости |
||||||||
при |
числе |
элементов |
щ=2, |
3, |
от числа |
элементов в об |
|||||||||
мотках |
с |
заземленным |
|||||||||||||
5 и 10 для различных т. Было |
(A) |
и |
изолированным |
||||||||||||
также |
рассчитано |
максималь |
(B) |
концами. |
|
|
|||||||||
но |
возможное |
отклонение |
Д, |
|
|
|
|
|
|
||||||
полученное |
как арифметическая |
сумма |
амплитуд |
всех |
оставшихся высших гармоник схемы с распределенными
постоянными, не |
воспроизведенных |
в |
цепной |
схеме. |
|
На рис. 10-17 |
показана |
зависимость |
относительных |
||
отклонений от порядка гармоники. На |
рис. 10-18 |
пока |
|||
зана зависимость |
максимально возможного отклонения |
||||
от числа элементов |
(кривые |
А). |
|
|
|
2. Обмотка с изолированным |
концом |
|
|
|
Анализ, проведенный для обмотки с заземленным концом, применим с соответствующими изменениями и для обмотки с изолированным концом. На рис. 10-18
360
приведена зависимость максимально возможного откло нения для этого случая (кривые В) от числа элементов модели.
3.Выводы
Влияние числа элементов на точность моделирования видно из рис. 10-18.
Если допустить максимальное отклонение в резуль тирующей кривой напряжения 10%, то при &/с/2 = 1/100
потребуется |
примерно |
10 элементов, а при £/с/2 =1/10 — |
примерно 3 |
элемента |
(как при заземленном, так и при |
изолированном конце |
обмотки). Поскольку значение |
k/cl2, как правило, близко к 1/10, достаточно |
ограничить |
||
ся при моделировании пятью элементами. |
|
|
|
Следует |
отметить, что показанное на рис. 10-18 |
ма |
|
ксимальное |
отклонение — это теоретический |
предел, |
ко |
торый едва |
ли может быть достигнут, поскольку ампли |
тудные значения отдельных гармоник не появляются одновременно и высшие гармоники почти полностью за тухают, так что они незначительно влияют на резуль тирующее напряжение. Это отклонение следует учиты вать лишь при оценке напряжений обмотки относительно
земли, но |
не при оценке напряжений на участках |
обмотки. |
|
Глава |
одиннадцатая |
АНАЛИЗ ИМПУЛЬСНЫХ ПРОЦЕССОВ
В ОБМОТКАХ ТРАНСФОРМАТОРОВ
СПОМОЩЬЮ ЦИФРОВЫХ ВЫЧИСЛИТЕЛЬНЫХ МАШИН
11-1. Общие соображения
В последние годы большое внимание уделялось про блеме расчета импульсных перенапряжений в. обмотках трансформаторов с помощью цифровых вычислительных машин (ЦВМ), [Л. 11-1 —11-14], поскольку на стадии проектирования непосредственный расчет на основе за данных параметров обмотки является единственным (за исключением моделирования) практическим методом оценки импульсных характеристик обмотки. Кроме того, с помощью расчета могут быть проверены импульсные
361
характеристики существующих обмоток без воздействия на них относительно высоких испытательных импульсов.
Одной из основных причин, которая привела к исполь зованию ЦВМ, явилась большая трудоемкость обычных расчетов даже для простейших обмоток в тех случаях, когда требуется получить результаты, достаточно точные для практических целей. Вскоре, однако, стало очевид ным, что даже с помощью современных расчетных мето дов, использующих ЦВМ, не удается справиться со все ми стоящими проблемами, и некоторые из них до сих пор фактически не решены. Сложность физических про цессов, происходящих в обмотках при импульсном воз действии, большое число различных типов обмоток и их устройств, используемых в трансформаторах высокого и сверхвысокого напряжений, привели к таким трудностям расчета, которые не смогли быть преодолены даже при принятии существенных (однако, практически приемле мых) упрощений. Для того чтобы получить приемлемые для практики результаты, используя ЦВМ малой и сред ней мощности, требуется достаточно большое время.
Численный анализ распределения импульсного напря жения в обмотке трансформатора — физического процес са, происходящего в обмотке при распространении им пульсной волны, производится пока двумя, в принципе различными, математическими методами.
Первый метод исходит из уравнений Максвелла и описывает переходный процесс в трансформаторах при определенных упрощающих условиях в виде распростра
нения электромагнитных волн. |
|
|
||
Второй метод анализа |
импульсных процессов, |
имею |
||
щий больше |
преимуществ |
и более |
распространенный на |
|
практике, требует, чтобы |
обмотка |
трансформатора |
была |
|
представлена |
более или |
менее сложной цепной схемой |
с конечным или бесконечно большим числом элементов; соответственно этому параметры схемы рассматривают ся сосредоточенными или распределенными. Если число элементов конечно, импульсный процесс, происходящий в такой эквивалентной схеме, описывается системой со вместных обычных дифференциальных уравнений, в то время как при бесконечном числе элементов мы имеем систему уравнений в частных производных или систему интегральных уравнений.
Каждый из этих двух методов при использовании ЦВМ имеет свои преимущества и недостатки в зависимо-
362
сти от условий, при которых производится расчет, и от требований к расчету (тип вычислительной машины, главным образом ее память и быстродействие и т. п., время расчета и экономичность данной машины).
Описанные |
в |
статьях |
методы |
расчета |
распределе |
||
ния |
импульсного |
напряжения в |
обмотке |
трансформа |
|||
тора |
(во времени |
и пространстве) |
основаны |
на |
рассмот |
||
рении схемы захМещения с конечным числом |
элементов |
||||||
[Л. |
11-2]. Для |
этой цели |
обмотка |
трансформатора под |
разделяется на ряд частей, каждая из которых пред ставляется в схеме замещения одним элементом. Отдель ные элементы выбираются и соединяются таким обра зом, чтобы воспроизвести емкости и индуктивности реальных частей обмотки, а также емкостные и индуктив ные связи между ними. Их соединения точно соответст вуют фактическим связям между 'частями в реальной обмотке. Схема замещения может быть составлена с уче том затухания импульсного процесса вследствие поверх ностного эффекта и разного вида потерь. Переходный процесс в таких схемах описывается системой диффе ренциальных уравнений второго порядка. Затем эта си стема может быть решена или, после соответствующего преобразования, на аналоговой машине [Л. 11-2], или, при использовании обычных численных методов, на циф
ровой |
машине |
{Л. |
11-1], |
[Л. 11-3], или, наконец, с по |
мощью |
комбинации |
этих |
методов. |
|
Теоретически |
эти |
два |
метода могут быть использова |
ны для любых типов катушечных и слоевых обмоток трансформатора, даже неоднородных. Они без труда мо гут быть применены при любой форме воздействующей импульсной волны, а число элементов в эквивалентной схеме может быть увеличено в пределах памяти вычис
лительной машины для того, чтобы получить |
результаты |
|
с |
требуемой точностью. Следовательно, они |
пригодны, |
в |
частности, для слоевой обмотки сложной конструкции |
|
в |
тех случаях, когда требуется не детальное |
рассмотре |
ние распределения импульсного напряжения и градиен тов внутри слоев, а только форма импульсного напря жения в отдельных точках обмотки, где ожидаются наи большие воздействия.
Оба метода имеют недостатки. Во-первых, распреде ление импульсного напряжения вдоль обмотки в про странстве и времени должно быть рассчитано от точки к точке, т. е. дискретно, так что вычислительная машина
363
должна быть использована снова для каждого индиви дуального случая импульсных воздействий на анализи руемую обмотку, а также после каждого изменения ее конструктивных параметров. Во-вторых, более подроб ный расчет распределения импульсного напряжения и особенно градиентов внутри слоев обмотки требует схе мы замещения с большим числом элементов, а с уве личением числа элементов требуемое расчетное время и память вычислительной машины увеличиваются. Даль нейшие трудности возникают вследствие того, что, как правило, мы не знаем, каким образом и на сколько эле ментов следует разбить обмотку, чтобы получить резуль таты с высокой точностью.
В дальнейшем рассмотрим мегэд расчета |
импульс |
ных перенапряжений в катушечных и слоевых |
обмотках |
с помощью цифровой машины. Этот метод основан на анализе схемы замещения с бесконечно большим числом элементов, воспроизводящей реальную обмотку транс форматора. В данном случае возможно получить общее решение распределения импульсного напряжения вдоль обмотки в замкнутой форме в зависимости от парамет ров обмотки. Применение этого метода целесообразно для определения импульсных перенапряжений в одно родных катушечных обмотках вследствие существенного упрощения расчета в этом случае. Он также может быть использован при получении расчетных таблиц и всех параметров обмоток конкретных конструкций, необходи мых для простой и быстрой оценки импульсных напря жений с помощью вычислительной машины или для обычного расчета. Его использование также целесооб разно в случае слоевых обмоток (с числом слоев менее 5), особенно когда необходимо узнать распределение им пульсного напряжения относительно земли и градиенты внутри каждого слоя. Метод, описанный ниже, справед лив только для однородных обмоток при воздействии прямоугольного импульса; для неоднородных обмоток он должен быть видоизменен. Импульсные напряжения для других форм волн могут быть легко определены с помощью хорошо известной формулы исходя из реше ния, полученного при воздействии прямоугольного им пульса.
Заслуживают внимания также методы расчета на чального распределения напряжения, возникающего в об мотке сразу после воздействия импульса, перед тем как
364
начнет сказываться действие магнитных связей; оно дает нам сведения о максимальных напряжениях в некоторых частях обмотки, например вблизи линейного конца ка тушечной обмотки (с помощью вычислительных машин было получено начальное распределение для большин ства различных типов обмоток) [Л. 4-38 и 5-1].
11-2. Расчет собственных частот и амплитуд
импульсных напряжений в однослойной и катушечной обмотках
Рассмотрим применение вычислительной машины для расчета распределения импульсного напряжения в одно родной однослойной обмотке трансформатора с зазем ленной нейтралью при воздействии на ее линейный ко нец единичного импульса. Метод и результаты расчета могут (с упрощением) быть также применимы для ана лиза распределения импульсного напряжения в одно родной дисковой обмотке с заземленной нейтралью. Рас сматриваемая обмотка, находящаяся на магнитопрово де, была представлена схемой замещения с бесконечным числом элементов (см. гл. 2, рис. 2-1). Элементы соответ ствуют отдельным виткам в однослойной обмотке или отдельным катушкам в случае катушечной обмотки, раз меры которых пренебрежимо малы по сравнению с раз
мерами всей обмотки. Они |
характеризуются |
емкостями |
|
и индуктивностями. Первые |
представляются |
продольной |
|
емкостью и емкостью элементов реальной обмотки |
(т. е. |
||
витков и катушек) на землю, последние отражают |
вза |
имные магнитные связи в установившемся режиме меж ду различными элементами обмотки. Реальный электро магнитный процесс в обмотке при воздействии единич ного импульса аппроксимируется квазистационарными процессами, происходящими при одних и тех же началь ных и граничных условиях в схеме замещения. Затуха нием импульсного процесса вследствие поверхностного эффекта и потерь пренебрегаем. Реальная зависимость взаимной индуктивности от аксиального расстояния меж
ду двумя элементами обмотки для |
упрощения |
аппрокси |
|||
мируется следующими |
экспоненциальными |
функциями |
|||
(2-99) : |
|
|
|
|
|
Мх (X - |
S) = |
М0е-Цх~1) |
для 5 < |
х; |
|
М^-х) |
= |
Мйе~ш~х) |
для S > |
je, |
|
ЗѲ5
где М0—индуктивность |
элемента обмотки, зависящая от |
конструкции обмотки |
(см. гл. 2 и 3). Экспоненциальная |
функция воспроизводит действительную функцию взаим ной индуктивности только в случае, если вторичная об мотка трансформатора разомкнута или если магнитопровод полностью экранирован (например, медным листом). Если вторичная обмотка трансформатора замкнута на коротко или на небольшое сопротивление, экспоненциаль ная функция не отражает фактическую функцию взаим ной индуктивности [Л. 2-24], которая в этом случае луч ше воспроизводится экспонентой, смещенной на посто янную величину [Л. 2-31]. Однако экспоненциальная функция М(х—І) может быть использована для расчета низших гармоник свободных колебаний обмоток с зазем ленной нейтралью даже тогда, когда вторичная обмотка рассматриваемого трансформатора короткозамкнута при условии, что коэффициент трансформации большой.
Рассчитанное распределение импульсного напряжения относительно земли отражает соответственно зависимо сти напряжения между витками в однослойной катушке или между катушкой в дисковой обмотке и землей (если размеры этой катушки пренебрежимо малы по сравне нию с размерами обмотки) от времени. В результате с учетом всех сделанных предположений дифференци
альное |
и интегральное уравнения для напряжения |
и(х, t) |
и тока i(x, t) схемы замещения, показанной на |
рис. 2-2, могут быть получены из уравнения в частных производных четвертого порядка относительно напряже
ния и(х, |
t) аналогично |
(2-100). Решение этого |
уравнения |
|||
в данном случае имеет |
вид: |
|
|
|
||
|
u{sc,t) = |
|
'£un(x,t)-\-uPf(x), |
(11-1) |
||
|
|
|
п |
|
|
|
где |
|
|
|
|
|
|
ип(х, |
і) ={Ап |
sin ащх + Вп |
cos |
umx + Pn sh |
агпхЛ- |
|
|
|
+iRn |
ch aznx] cos |
wnt |
(11 -2) |
|
— составляющие |
напряжения, |
соответствующие собст |
венным колебаниям обмотки, содержащие для каждой собственной частоты обмотки четыре постоянные инте
грирования: Ап, |
Вп, |
Рп, Rn. |
аксиальное |
Переменной |
х обозначено относительное |
||
расстояние (х = х/1) |
элементов обмотки от |
линейного |
366
Конца |
( / — а к с и а л ь н а я |
длина |
обмотки). Между собствен |
ными |
частотами ш„ и |
а і п = |
аі„/_существует связь [см. |
(2-102)], зависящая от величины Х=Х1:
{„2 _ |
An ft1 + An) |
(11-3) |
|
|
|
и между а 1 П и аІП |
отношение |
|
|
|
(11-4) |
Upf{x) — так называемое квазистационарное распреде
ление |
напряжения, |
является частным решением |
задачи |
|||||||||||
и не зависит от времени. |
|
|
|
|
|
|
||||||||
|
Из |
исходной системы уравнений и из уравнения |
||||||||||||
(11-2) |
можно получить выражение для тока |
|
|
|||||||||||
i |
(X, t) = |
Ак |
(t) |
+ |
^ |
[ - |
Ап |
( 1 + |
- ^ * \ п |
j cos |
âinx |
4- |
||
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
+ |
BnSi. |
|
С1+ |
- f c M s i n â i n * + Я я |
-£i- |
X |
|
|||||
|
|
|
|
X |
( 1 - |
Tf^ri ) s h *»n* |
sin «„г1, |
|
(11-5) |
|||||
где |
Cg — общая |
емкость |
обмотки |
на |
землю |
(/=1), |
||||||||
a /(s — |
общая |
продольная |
емкость; |
An(t)—частное ре |
||||||||||
шение, зависящее только от времени. |
|
|
|
|
||||||||||
|
Частные решения для свободных колебаний напря |
|||||||||||||
жения |
ип(х, |
t) |
|
и тока |
in(x, |
t) обмотки |
с заземленной |
нейтралью должны удовлетворять для каждого n исход ной системе уравнений я граничным условиям:
М М = 0 ; }
ип(і,о=о. !
Их постоянные интегрирования должны удовлетво рять, следовательно, для каждого п однородной системе четырех линейных алгебраических уравнений: (2-117),
367
(2-118), (2-121) и (2-122), определитель которой имеет вид:
0 |
|
|
|
|
1 |
|
|
|
|
|
|
sina,„ |
|
|
|
|
cos<xm |
|
|
|
|
|
|
|
|
|
|
|
~2 |
|
|
|
|
|
|
аіп4п |
|
|
|
|
"2" |
|
|
|
|
|
|
-2 f Т |
- |
|
. - |
|
\ _=± |
К |
. - |
|
, |
|
|
a |
cosai n |
— sinai n |
- 2 |
=3—srnam -f-cosai n |
|||||||
...0sh а 2 П |
|
|
|
|
Ich а2 |
|
|
|
|||
"In «an |
|
|
|
|
|
|
An |
|
|
|
|
••• J |
— ( — |
chа 2 П + |
=4-sh а 2 П |
) -ü— ( - =ф - shа 2 П - f |
|||||||
|
|
|
|
|
|
|
|
+ -=4—cha2 |
|||
|
|
|
|
|
|
|
|
|
|
|
'(11-7) |
Его столбцы образованы из коэффициентов при по |
|||||||||||
стоянных Ап, |
Вп, Рп, Rn- Кроме |
того, |
частные |
решения |
|||||||
uPf(x), |
Ak(t) |
должны |
удовлетворять |
соотношению |
|||||||
и р / (х ) = |
л™. È£, |
|
2* - ± - e - r |
x + |
± e Î C x - l ) |
]+Dk, |
|||||
где Dk — постоянная |
интегрирования. |
|
|
(11-8) |
|||||||
|
|
|
|||||||||
Общее решение (11-1) должно удовлетворять гранич |
|||||||||||
ным условиям |
на линейном |
конце, где воздействует еди |
|||||||||
ничный |
импульс: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
« ( 0 , 0 = 1 - |
|
|
|
(П-9) |
||
и, следовательно, для любого t |
|
|
|
|
|
||||||
ц р / (0) = |
^ |
^ ( |
- |
1 + |
е ~ |
Г ) |
+ |
/Л = 1. |
(П-Ю) |
||
Кроме того, |
если |
|
нейтраль |
обмотки |
заземлена, мы |
||||||
для любого t имеем: |
|
|
|
|
|
|
|
|
|||
« р / |
І Ѵ Ш 0 |
OA*. |
, |
п - |
„-ХТ |
+ 1 ) |
+ |
/Л = 0. |
(Ц.11) |
||
( 1 ) = ^ - ^ ( - 2 Я - < Г |
368
Оба |
уравнения |
справедливы |
лишь в случае, |
если |
dAh/dt |
не зависит |
от времени, |
г. е. Ak{t)='Äht. |
Тогда |
постоянные Ah, Du могут быть определены из (11-10) и (11-11):
z _ J> |
1 |
. |
n — 2 І + е ~ Х + 1 |
"Т |
|
ОД/2 м |
— |
Т" ' |
ft |
— |
|
|
|
|
|
|
|
Таким образом, |
|
|
|
|
|
Ц р / а д = ^ |
+ |
|
21 ( ! - « ) - ! ^ |
( 1 М 2 ) |
|
|
|
2(Х + е - х — 1) |
|
|
Однородная система линейных уравнений с неизвест ными коэффициентами Ап, Вп, Рп, Rn имеет решение, отличное от 'Нуля лишь в случае, если определитель D (11-7) равен нулю
|
|
|
|
D = 0. |
|
|
|
|
|
(П-13) |
|||
Как известно, это условие удовлетворяется |
только |
||||||||||||
при определенных am, um |
{n=l, |
2, 3 ... ) и |
соответству |
||||||||||
ющих |
им сои, |
которые |
являются |
корнями |
уравнения |
||||||||
(11-13). Подставляя эти значения |
(так называемые |
соб |
|||||||||||
ственные частоты обмотки) |
в однородную систему |
урав |
|||||||||||
нений с коэффициентами Ап, |
Вп, |
|
Рп, |
Rn, |
можно |
полу |
|||||||
чить для каждого га отношения: |
|
|
|
|
|
|
|
||||||
|
Вп/Ап |
= аП2-, Рп/Ап |
= апз; |
RnlAn |
= ani. |
(11-14) |
|||||||
Подставляя |
в элементы |
определителя D значения |
|||||||||||
|
|
|
а т = і л + 8п , |
|
|
|
(11-15) |
||||||
где En — величина, зависящая |
от га, из |
(11-13) |
получаем, |
||||||||||
что en—>-0 при га—>-оо и, |
следовательно, |
a m — п р и |
|||||||||||
га—м». |
|
|
|
|
|
|
|
|
|
|
|
|
|
Подстановкой (11-15) в (11-13) можно получить при |
|||||||||||||
ближенную |
зависимость |
а і и ( я ) |
для больших га. Прене |
||||||||||
брегая для этих га величинами |
еп/ |
{пп + гп)^, |
8n/(ran-f-en)3, |
||||||||||
Eni(пп |
+ еп)2, |
1/(гал + 8п)3 |
по |
сравнению с |
1/(гая + е п ) и |
||||||||
En и |
полагая |
sine„ = en , |
c o s e n = l , |
получаем |
для |
доста |
|||||||
точно |
больших га уравнение относительно |
еп- |
|
|
|||||||||
|
|
|
s2 -f- Ппеп |
— 21 |
= |
0. |
|
|
|
|
|
24—8 |
369 |