Научная статья УДК 621.314.225
DOI: https://doi.org/10.18721/JEST27405
Л.И. Сахно н, О.И. Сахно, М.С. Крылов
Санкт-Петербургский политехнический университет Петра Великого, Санкт-Петербург, Россия
н lsahno2010@yandex.ru
МЕТОДИКА РАСЧЕТА И ИЗМЕРЕНИЯ АКТИВНЫХ МОЩНОСТЕЙ И ПОТЕРЬ В ТРАНСФОРМАТОРАХ КОНТАКТНОЙ СТЫКОВОЙ СВАРКИ ОПЛАВЛЕНИЕМ ТРУБ
Аннотация. В работе рассматриваются трехобмоточные трансформаторы с двумя параллельно соединенными секциями первичной обмотки и одной вторичной обмоткой. Такая конструкция используется в трансформаторах для стыковой контактной сварки оплавлением труб. В статье дается методика расчета токов, активной мощности и потерь в обмотках, основанная на схеме замещения многообмоточного трансформатора. В основе схемы замещения трехобмоточного трансформатора лежит замена его двумя двухобмоточными трансформаторами. Все параметры схемы замещения могут быть измерены. Установлена связь между активной мощностью и потерями в обмотках, которая лежит в основе косвенного метода измерения потерь в обмотках. Анализ точности методики проведен на модельном трансформаторе МКЭ. Для пояснения методики построены векторные диаграммы, показывающие, в частности, что активная мощность в одной из первичных обмоток при коротком замыкании может быть отрицательной (угол между векторами тока и напряжения в обмотке превышает 90°). Получены соотношения между параметрами двухобмо-точных трансформаторов, при которых в одной из параллельно соединенных первичных обмоток возможно появление отрицательной активной мощности. Результаты расчета токов и активных мощностей в трехобмоточном трансформаторе сварочного комплекса КСС-04 для сварки труб хорошо согласуются с результатами измерений.
Ключевые слова: трехобмоточный трансформатор, схема замещения, активная мощность, магнитное поле рассеяния, взаимная индуктивность по потокам рассеяния.
Для цитирования:
Сахно Л.И., Сахно О.И., Крылов М.С. Методика расчета и измерения активных мощностей и потерь в трансформаторах контактной стыковой сварки оплавлением труб // Материаловедение. Энергетика. 2021. Т. 27, № 4. С. 68-84. DOI: https://doi.org/10.18721/ JEST.27405
Это статья открытого доступа, распространяемая по лицензии СС BY-NC 4.0 (Ы^://сгеа-tivecommons.Org/licenses/by-nc/4.0/)
© Сахно Л.И., Сахно О.И., Крылов М.С., 2021. Издатель: Санкт-Петербургский политехнический университет Петра Великого
Scientific article
DOI: https://doi.org/10.18721/JEST.27405
L.I. Sakhno H , O.I. Sakhno, M.S. Krylov
Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia
H lsahno2010@yandex.ru
METHOD FOR CALCULATING AND MEASURING ACTIVE POWERS AND LOSSES IN FLASH-BUTT WELDING TRANSFORMERS
Abstract. The paper considers three-winding transformers with two parallel-connected sections of the primary winding and one secondary winding. This design is used in flash-butt welding transformers. The article provides a method for calculating currents, active power and losses in the windings, based on the equivalent circuit of a multi-winding transformer. The equivalent circuit of a three-winding transformer is based on its replacement with two two-winding transformers. All parameters of the equivalent circuit can be measured. The relationship between active power and winding losses has been established, which underlies the indirect method of measuring winding losses. The analysis of the accuracy of the technique was carried out on a model FEM transformer. To clarify the technique, vector diagrams have been constructed, showing, in particular, that the active power in one of the primary windings during a short circuit can be negative (the angle between the current and voltage vectors in the winding exceeds 90°). Relationships between the parameters of two-winding transformers in this case are obtained. The results of calculating the currents and active powers in the three-winding transformer of the KSS-04 welding complex for pipe welding are in good agreement with the measurement results.
Keywords: three-winding transformer, equivalent circuit, active power, magnetic leakage flux, mutual inductance on leakage flux.
Citation:
L.I. Sakhno, O.I. Sakhno, M.S. Krylov, Method for calculating and measuring active powers and losses in flash-butt welding transformers, Materials Science. Power Engineering, 27 (04) (2021) 68-84, DOI: https://doi.org/10.18721/JEST.27405
This is an open access article under the CC BY-NC 4.0 license (https://creativecommons.org/ licenses/by-nc/4.0/)
Введение. Эффективным способом сварки труб при строительстве трубопроводов является применение контактной стыковой сварки оплавлением [1—2]. При сварке магистральных трубопроводов большого диаметра (800 — 1420 мм) с 70-х годов ХХ века успешно используются машины контактной сварки внутритрубного исполнения, при котором сварочная машина перемещается внутри трубопровода и поочередно приваривает к его концу следующую трубу.
В последние годы для увеличения пропускной способности газопроводов давление в них стали поднимать до 11,8 МПА. Для обеспечения прочности трубопроводов необходимо применять трубы из сталей повышенного качества, что требует увеличения мощности сварочных машин. Источником сварочного тока внутритрубных машин является сварочный трансформатор тороидальной конструкции, обеспечивающий равномерное распределение тока по периметру свариваемых труб. Увеличение мощности трансформатора требует увеличения сечения его первичной обмотки без увеличения наружного габарита трансформатора. Конструктивные и технологические сложности увеличения мощности трансформаторов были преодолены на ЗАО «Псковэлектросвар»1
1 ЗАО «Псковэлектросвар» https://pskovelektrosvar.ru (дата обращения 5.07.2021)
© L.I. Sakhno, O.I. Sakhno, M.S. Krylov, 2021. Published by Peter the Great St. Petersburg Polytechnic University
нетрадиционным путем — секционированием обмотки на параллельные секции, вмотанные в тороидальный магнитопровод одна поверх другой. При проектировании такого трансформатора необходимо рассчитать потери в секциях и температуру нагрева обмоток в номинальном режиме работы. Стандартные методики расчета трансформаторов контактной сварки [3] не позволяют найти потери в секциях параллельных обмоток для указанной новой конструкции.
Испытания трансформатора для сварки труб с двумя параллельными секциями первичной обмотки в режиме короткого замыкания показали, что активные мощности и токи в секциях не одинаковы, причем в одной из них ваттметр показывает отрицательное значение (угол между векторами тока и напряжения в секции превышает 90°). Поскольку тепловой режим работы трансформатора является весьма напряженным и в значительной степени определяет работоспособность и надежность машины, наиболее сложной и актуальной в настоящее время задачей, стоящей перед разработчиками источников питания внутритрубных машин, является задача расчета и измерения потерь в трансформаторе с несколькими параллельными секциями первичной обмотки.
Целью настоящей статьи является разработка методики расчета активных мощностей и потерь в обмотках тороидального сварочного трансформатора с двумя параллельно соединенными секциями первичной обмотки и установление связи между потерями в секциях и показаниями ваттметров, установленных в этих секциях. Кроме того, в статье должна быть исследована точность расчета с использованием разработанной методики на основе сравнения результатов расчета с результатами численного и натурного экспериментов.
Метод исследования
Для разработки методики расчета активных мощностей и потерь используется предложенная нами в [4] схема замещения многообмоточного трансформатора, поскольку эта схема в отличие от классической трехлучевой схемы замещения трехобмоточного трансформатора [5] и ее модификаций [6—10] позволяет установить связь между активной мощностью в обмотках и потерями в них. В оригинальных схемах замещения многообмоточных трансформаторов [11, 12] не рассмотрены тороидальные трансформаторы и не установлена связь между потерями и активной мощностью в обмотках. В данной работе установить эту связь удалось благодаря физичности параметров, входящих в схему замещения [4], и наличию в ней ветвей эквивалентных секциям первичной обмотки.
В [4] рассмотрен трансформатор с n первичными обмотками и одной вторичной обмоткой. Для удобства записи уравнений трансформатора в этой работе индексы первичных обмоток имеют две цифры 11, 12... 1n, вторичная одну цифру 2. Первая цифра 1 индекса обмотки с двумя цифрами означает, что обмотка первичная, вторая цифра — номер первичной обмотки. Первичные обмотки не соединены гальванически и подключены к источникам напряжения un, un... uin. Индексы источников показывают, к каким первичным обмоткам они присоединены. Вторичная обмотка 2 подключена к нагрузке zH (для синусоидального режима zH = rH + j^>LH, где rH — активное сопротивление нагрузки, LH индуктивность, ю = 2nf, f — частота). Этот трансформатор заменяется на n двухобмоточных трансформаторов с обмотками 11 и 2, 12 и 2, ., 1n и 2 (далее будем обозначать эти трансформаторы как 11-2, 12-2. 1n-2). На основе классических уравнений двухобмоточных трансформаторов и метода наложения в [4] получены n-уравнений для этого трансформатора:
UU = j (А;-2 + L'Hl ) 4 + (Г-2 + r'i ) 4 +
n n __(1)
+j®S( mHm ,2 + LH.) 4+Z( r i,1l ,2 + rHi ) Wu, j = i, n
i=i i=i
Ш Ш
где I — номер первичной обмотки, ] = ии - напряжение на зажимах обмотки 11, Ь1._2, г1._2 -
соответственно индуктивность и активное сопротивление короткого замыкания трансформатора 1.-2, приведенные к обмотке И, к'т, г'т — соответственно индуктивность и активное сопротивление нагрузки, приведенные к обмотке И, к. = w;i/w2, к1 = — коэффициенты трансформации двухобмоточных трансформаторов с обмотками 1.-2, 11-2, w1, w1, w2 — количество витков обмоток И, 11, 2, М., — взаимная индуктивность по потокам рассеяния трансформаторов 1.-2и 11-2. Взаимная индуктивностьМ. 1х 2и параметр г.112 определяются из опытов короткого замыкания трансформаторов 1.-2, 11-2, 1.-11:
Ми,„,2 =(к-2 + к-2 - к-11 )/2' (2)
ГИ,1Х,2 =(Гц-2 + Гц-2 - Г1г-1Х V2' (3)
где к. 2, г. 2; кп 2, г}[ 2; к. , г. п — индуктивности рассеяния и активные сопротивления трансформаторов 1.-2; 1Х-2; И-1[ в соответствующих опытах короткого замыкания [4].
Для трехобмоточного трансформатора (рис. 1а) с двумя параллельно соединенными первичными 11 и 12 обмотками и одной вторичной обмоткой 2, подключенной к нагрузке 1Н = ГН + у'юкд, уравнение (1) примет вид:
и1 = ,/Ц к11-2 + к') ¡11 +(Г11-2 + Г'н) ¡11 + 7Ю(М11,12,2 + кН) ¡12 +(Г
1,12,2 + гн ) кк ¡12 ,
(4)
и2 = М к12-2 + ) ¡12 + (Г12-2 + ГН ) ¡12 + -/ЦМ12,11,2 + ) ¡11 + (Г12,11,2 + ГН ) КК^Г
Для данной схемы и1 = и 2, и, так как нагрузка одна, второй индекс у нагрузки опущен.
Введем обозначение Мп 122 = М12112 = М, ги 122 = г12112 = г. С учетом введенных обозначений
для режима короткого замыкания (к'н = 0, г'т = 0) (4) перепишем в виде:
и1 = 7юк11-2¡11 + Г 1-2¡11 + >^12 + Гк1к2¡12 '
(5)
и1 = ¡Юк^к + Г 2-2 ¡12 + + ^2^11-
Схема замещения, соответствующая уравнениям (5), показана на рис. 1б. На этом рисунке зависимый источник тока определяет ток во вторичной обмотке. Параллельно соединенные ветви этой схемы замещают двухобмоточные трансформаторы 11-2 и 12-2 с учетом их взаимного влияния. Как будет показано и подтверждено далее численным и натурным экспериментами активная мощность в каждой параллельной ветви схемы замещения совпадает с активной мощностью, измеряемой ваттметрами, установленными в секциях 11 и 12 первичной обмотки. Здесь важно, что все параметры схемы замещения, представленной на рис. 1б, можно найти экспериментально. В зависимости от соотношения индуктивностей и взаимной индуктивности в схеме на рис. 1б активные мощности (показания ваттметров) в ветвях этой схемы могут принимать как положительные, так и отрицательные значения. Это явление рассмотрено для индуктивно связанных катушек в [13]. Однако непосредственно перенести результаты [13] на рассматриваемый трансформатор нельзя, поскольку в схеме замещения присутствуют не индуктивности катушек, а индуктивности рассеяния трансформаторов. Кроме того, в схеме замещения имеются зависимые источники напряжения, которые появились из-за связи трансформаторов 11-2 и 12-2 через общую вторичную обмотку 2.
Рис. 1. Электрическая схема (а) и схема замещения (б) трехобмоточного трансформатора Fig. 1. Electrical circuit (a) and equivalent circuit (b) of a three-winding transformer
Сопоставление результатов, полученных с использованием схемы замещения с результатами, полученными с использованием расчета трансформатора методом конечных элементов (МКЭ) с присоединённой цепью. В этом сопоставлении сравниваются результаты расчетов токов, активных мощностей и потерь в модельном трансформаторе, полученные на основе МКЭ для квазистационарного магнитного поля трансформатора совместно с уравнениями для присоединенной электрической цепи нагрузки (первый способ) с результатами, полученными по схеме замещения, рис. 1б (второй способ). Погрешность первого способа обусловлена только размерами ячеек сетки и положением границ, которые были определены предварительно, таким образом, чтобы их дальнейшее изменение (уменьшение размеров ячеек и увеличение расстояния от трансформатора до границ расчетной области) не влияли на точность расчета токов и активных мощностей в обмотках трансформатора. Этот способ соответствует натурному эксперименту, в котором амперметры и ваттметры, установленные в секциях первичной обмотки, показывают токи и активные мощности в секциях. Первый способ дает возможность найти и потери в каждой обмотке, что существенно расширяет возможности анализа точности расчета потерь в обмотках с помощью численного эксперимента по сравнению натурным опытом короткого замыкания.
В выполняемом сопоставлении параметры схемы замещения (для второго способа) находятся МКЭ из расчета режимов противовключения трансформаторов 11-2, 12-2 и 11-12. Дискретизация расчетной области и положение границы в первом и втором способе расчета одинаковы. Для расчета первым и вторым способом используется пакет программ ELCUT [14]. Выполненное сопоставление дает возможность показать справедливость предлагаемого в статье метода расчета активной мощности и потерь в параллельных секциях первичной обмотки в широком диапазоне изменения электромагнитных параметров трансформатора, в том числе при изменении этих параметров из-за влияния поверхностного эффекта. Для детального описания обоих подходов, участвующих в сопоставлении, выполним подробное описание объекта исследования.
Необходимый нам многовариантный расчет трехмерного поля даже простого по конструкции трансформатора требует значительных вычислительных ресурсов, поэтому для исследования выбрана осесимметричная конструкция трансформатора (рис. 2а), которая позволяет использовать двумерную формулировку задачи. Прототипом этого трансформатора является трансформатор, описанный в [15]. Этот трансформатор имеет повышенное рассеяние, что обеспечивает возможность тщательного исследование влияния взаимной индуктивности по потокам рассеяния на токи в первичных обмотках и активную мощность в них.
Также для упрощения расчетов принято, что все обмотки трансформатора имеют одинаковое количество витков 110 (коэффициенты трансформации трансформаторов 11-2 и 12-2 равны единице). Обмотки выполнены из медного провода сечением 5^5 мм2.
а)
б)
в)
Рис. 2. Модельный трансформатор (а), присоединенная электрическая цепь для режима короткого замыкания (б), электрическая схема режима противовключения (в)
Fig. 2. Model transformer (a), connected electrical circuit for short-circuit mode (b), electrical circuit for counter-switching mode (c)
При расчете магнитного поля трансформатора совместно с уравнениями токов и напряжений цепи нагрузки расчетная область представляет собой поперечное сечение трансформатора, проходящее через ось симметрии 00' (рис. 2а), ограниченное половиной окружности, радиус которой превышает размеры поперечного сечения трансформатора. Границы расчетной области отодвинуты достаточно далеко от трансформатора, поэтому на них задано условие отсутствия поля. Присоединенная электрическая цепь представлена в виде блоков на рис. 2б, блоки 11, 12, 2 соответствуют обмоткам 11, 12, 2 в расчетной модели поля на рис. 2а. Магнитная проницаемость магнитопровода принята постоянной. Результатом решения являются комплексные амплитуды токов первичных обмоток, активные мощности и потери в обмотках.
Для использования схемы замещения на рис. 1б необходимо рассчитать импедансы двухоб-моточных трансформаторов 11-2, 12-2, 11-12. Они находятся из режима противовключения, в котором вектора магнитодвижущих сил в первичной и вторичной обмотках равны и противоположно направлены [9]. При этом рассчитывается переменное магнитное поле трансформатора с электрической цепью, моделирующей режим противовключения. На рис. 2в показана электрическая схема для расчета импеданса двухобмоточного трансформатора 11-2. Согласно [9] источники тока Jn, J2 на рис. 2в имеют одинаковые произвольные модули, но противоположные фазы:
Ju = 1пМи , 32 = =-Ju , (6)
где 3 и 32 — токи обмоток 11 и 2 соответственно.
Результаты сопоставления способов 1 и 2. Сопоставление выполнялось для различных частот из диапазона 1—50 Гц. Широкий диапазон изменения частоты использован нами для изменения соотношений между параметрами схемы замещения трансформатора и повышения достоверности расчетов. На нижнем пределе 1 Гц все электромагнитные параметры схемы замещения исследуемого трансформатора соизмеримы, а при 50 Гц (сильный поверхностный эффект) реактивные сопротивления в схеме замещения на порядок выше активных. Амплитуда источника напряжения принята равной 8 В. Расчеты выполнялись при магнитной проницаемости ц = 1000 Дополнительное исследование показало, что изменение ц магнитопровода от 500 ц0 до 10000 ц0 меняет значения токов и активных мощностей обмоток не более, чем на 0,05%.
В табл. 1 приведены амплитуды и фазы токов в обмотках 11, 12 и амплитуда суммарного тока 11т, полученные двумя способами.
Таблица 1
Амплитуды и фазы токов обмоток
Table 1
Amplitudes and phases of winding currents
/, Гц Способ расчета I ,,, A mW ф„, град I 12, A ф12, град I ,, A m1 ф1, град
1 Первый 52,18 -35,6 73,78 -16,6 124,3 -24,4
Второй 52,24 -35,1 73,89 -16,3 124,5 -24,6
50 Первый 0,85 130,4 8,39 -80,8 7,68 -84,1
Второй 0,82 130,2 8,38 -80,2 7,69 -84,9
Результаты, представленные в табл. 1 показывают, что токи, рассчитанные двумя способами, отличаются менее, чем на 1%. Таким образом, схема замещения позволяет находить токи в обмотках с высокой точностью.
Выполним далее анализ активных мощностей и потерь, рассчитанных двумя способами. Для анализа активных мощностей первичных обмоток при коротком замыкании на рис. 3 представлены векторные диаграммы для амплитудных значений входного напряжения и токов в обмотках, полученных первым способом. Из этих диаграмм видно, что проекция OC вектора приложенного напряжения U1m на направление вектора тока In при частоте 1 Гц совпадает с направлением тока в этой обмотке, а при частоте 50 Гц противоположна этому направлению, поскольку в этом случае вектор тока в обмотке 11 сдвинут по отношению к питающему напряжению на угол, больше 90 градусов.
Эти проекции определяют активные мощности в обмотке 11 Pn = Ua11m-I11m /2 и в обмотке 12P12 = Ual2mIl2m/2 и суммарную активную мощность P в трансформаторе. Активные мощности Pn, P12, P равны активным мощностям, измеряемым ваттметрами, установленными в секциях первичной обмотки. Результаты расчета активных составляющих напряжения, активных мощностей, потерь в обмотках 11, 12, 2даны в табл. 2. В этой таблице Pn11, Pn12, Pn2 — потери в обмотках 11, 12 и 2.
Таблица 2
Амплитуды активных составляющих напряжения, активных мощностей и потерь (способ 1)
Table 2
Amplitudes of active voltage components, active powers and losses (method 1)
/, Гц Uu , В allrn1 U „ , В а 12m Pw Вт Pa, Вт P, Вт P Вт П11' P Вт П12 P Вт П2 Pw Вт
1 6,5 7,66 170 282 453 103 151 199 453
50 -5,19 1,27 -2,20 5,37 3,16 0,032 2,2 0,908 3,14
Суммарная активная мощность в модельном трансформаторе, найденная по формуле Р =
: Re (IV),
где * — символ комплексного сопряжения, для заданной частоты совпадает со значением суммарной активной мощности Р = Рп + Р12 и суммой потерь в обмотках Р = РП11 + РП12 + РП2. Из диаграмм на рис. 3 виден только результат — в обмотке 11 угол между током и напряжением больше 90° при частоте 50 Гц, что привело к отрицательной активной мощности в
этой обмотке. Однако не ясно, какие параметры трансформатора привели к получению такого результата.
Для анализа этого явления получим активные мощности в обмотках вторым способом, то есть используя схему замещения на рис. 1б. Параметры схемы замещения, найденные из режимов противовключения, а также активные сопротивления обмоток r11, r12, r2 приведены в табл. 3.
Таблица 3
Параметры схемы замещения токи и активные мощности (способ 2)
Table 3
Equivalent circuit parameters currents and active powers (method 2)
f, Гц ®L11-2, Ом ®L12_2, Ом rn_2, Ом rn, Ом Ar, Ом r12-2, Ом r12, Ом r2, Ом r, Ом юМ, Ом
1 0,0456 0,0208 0,101 0,0753 0.00001 0,0814 0,0555 0,0259 0,0259 0,0229
50 2,27 1,04 0,146 0,0892 0,0259 0,095 0,0654 0,0309 0,044 1,14
При 1 Гц (слабый поверхностный эффект) из (3) следует, что сопротивление вторичной обмотки совпадает с коэффициентом г. При 50 Гц (сильный поверхностный эффект) в сопротивление г11-2 входит добавочное сопротивление Аг, которое создается полем рассеяния в обмотке 12, расположенной между секцией 11 и обмоткой 2 (рис. 2а). Коэффициент г отличается от сопротивления вторичной обмотки на величину Аг/2, что следует из (3). Векторные диаграммы для амплитуд, построенные с использованием параметров схемы замещения в табл. 2 согласно (5), представлены на рис. 4.
Из диаграмм видно, что активные составляющие напряжения на каждой обмотке складываются не только из падения напряжения на активных сопротивлениях г11-2 или г12-2, но и из проек-
а) б)
Рис. 3. Векторные диаграммы для частоты 1 Гц (а) и 50 Гц (б) Fig. 3. Vector diagrams for a frequency of 1 Hz (a) and 50 Hz (b)
а)
б)
Рис. 4. Векторные диаграммы для схемы замещения при частотах 1 Гц (а) и 50 Гц (б) Fig. 4. Vector diagrams for the equivalent circuit at frequencies of 1 Hz (a) and 50 Hz (b)
ций двух векторов — вектора зависимого источника напряжения и вектора напряжения, зависящего от взаимной индуктивности. Как для частоты 1 Гц, так и для частоты 50 Гц проекция вектора &М112т на направление вектора 1Пт (отрезок АВ) противоположно направлению этого вектора, а проекция ®М1ит на направление вектора 112т (отрезок CD) совпадает с направлением этого вектора, поэтому эти составляющие имеют разные знаки. Активные составляющие напряжений на обмотках 11 и 12 на рис. 4 равны:
Ua11m = ru 2I11m + rI12m cos0_<MI„m sin0,
a11m 11-2 11m 12m 12m ?
Ua12m = Г12- 2I12m + rI11m COS 0 + Sin 0:>
(7)
(8)
где 0 — угол между векторами токов 111т, 112т. Активные мощности в обмотках:
P1 = UaUIu = Г11-21121 + rI11I12 COS 0 - ^MI„I12 Sin 0 = Pu_2 + Pr _ Pm ,
(9)
Pu = UaUIu = ru_2^ + rlulu cos 0 + (aMJn sm 0 = Pu_2 + Pr + PM, где I , I , U .., U ,, — соответственно действующие значения токов в обмотках 11, 12 и активных
4P "12' all' a12
напряжений на них.
В (9) слагаемые Р11-2 = 111г11-2, Р12-2 = 112г12-2 определяют мощности потерь энергии в трансформаторах 11-2,12-2, приведенные к своим первичным обмоткам. Выражение РМ = шМГп/п81п0
определяет мощность, передаваемую из одного трансформатора в другой вследствие индуктив-
ной связи между ними через магнитное поле рассеяния. Выражение Pr = /n/12rcos0 определяет мощность потерь энергии в одном трансформаторе за счет падения напряжения на общей для двух трансформаторов вторичной обмотке 2 (с учетом добавочных потерь). Полная активная мощность всей цепи из (9) равна сумме потерь в обмотках:
P = Pi +P2 =r11-2/121 +r12-2/122 + 2r/11/12 cos е = рП11 +рП12 +pn2. (10)
Полная мощность P не зависит от мощности PM, передаваемой из одного трансформатора в другой вследствие индуктивной связи между ними через магнитное поле рассеяния. Амплитуды активных составляющих напряжения, рассчитанные по (7), (8) представлены в табл. 4, а составляющие активных мощностей — в табл. 5. Активные составляющие напряжения для частоты 1 Гц в обмотке 11 определяются, в основном, падением напряжения на активных сопротивлениях трансформаторов 11-2 и 12-2. Отрицательная составляющая активного напряжения, определяемая взаимной индукцией, мала по сравнению с положительными, поэтому активная мощность в обмотке 11 положительна. Для 50 Гц проекция вектора ®MI12m на направление вектора IUm, направленная противоположно этому вектору, существенно превышает все положительные составляющие активного напряжения обмотки 11, что приводит к отрицательной активной мощности в этой обмотке.
Таблица 4
Амплитуды активных составляющих напряжения
Table 4
Amplitudes of active voltage components
f, Гц I r 11m 11-2 I12mrCOS0 юМ!,. sin0 12m U „ allm I r 12m 12-2 I11mrCOS0 «¡ML, sin0 11m U„ a 12m
1 5,31 1,81 0,545 6,57 6,0 1,29 0,386 7,67
50 0,124 -0,315 4,95 -5,14 0,797 -0,0319 0,50 1,26
Таблица 5
Cоставляющие активной мощности
Table 5
Active power components
f, Гц Трансформатор 11-2 Трансформатор 12-2
I2 r 11 11-2 IUI12 rCOs0 roMI11I12sin0 Pu I2 r 12 12-2 IUI12 rCOs0 roMI11I12sin0 P12
1 138,5 47,2 14,28 171 221 47,2 14,28 282
50 0,053 -0,134 2,10 -2,18 3,34 -0,13 2,09 5,29
Анализ результатов в табл. 4, 5, а также результатов , полученных для других частот указанного выше диапазона, показывает, что при положительном значении взаимной индуктивности по потокам рассеяния обмоток (М > 0) и Ьп > М > Ь12 энергия передается из двухобмоточного трансформатора с меньшей индуктивностью рассеяния в трансформатор с большей индуктивностью. Активная мощность в трансформаторе с большей индуктивностью будет отрицательной, если юМ/л 12 sin0 в этом трансформаторе превышает сумму двух активных мощностей 1121г11-2 и !п 112 соэб. Полученные соотношения отличаются от соотношений, полученных для двух индук-
тивно связанных катушек с индуктивностями L^ L2 и взаимной индуктивностью M, в которых при M > 0 отрицательная активная мощность появляется в катушке с большей индуктивностью при условии (üMIuI2 ■ sin 0 > r1/121 (I1, I2 — токи в катушках, 0 — угол между векторами токов).
Таким образом, нами сформулированы условия возможности получения отрицательной активной мощности в одной из параллельно соединенных первичных обмоток трансформатора. Кроме того, можно сделать вывод о том, что схема замещения позволяет с высокой точностью рассчитывать активную мощность в обмотках рассматриваемого трансформатора, поскольку результаты расчета активных мощностей в обмотках 11 и 12 в табл. 2 и 5 (первый и второй способ расчета) отличаются не более, чем на 1,5%.
Далее необходимо установить связь активной мощности в обмотках и потерями в них. Из (9) можно найти потери в трансформаторах 11-2 и 12-2:
г11-2 /121 = Р11 -г/п/12ео8 0 + шМ/п/128т 0, (11)
Г12-27П = Р2 - г/11/12 С°8 0 - ^М/11/12 ^ 0 (12)
Потери в этих трансформаторах состоят из потерь в первичных и вторичных обмотках:
11-241 = Гц!2! + ^2 /11, (I3) г /2 = г /2 + г /2 (14)
12-2 12 '1242 ^'242 • V1 V
При слабом поверхностном эффекте активное сопротивление Г2 совпадает с коэффициентом г, найденным из трех опытов короткого замыкания, что следует из (3). Тогда потери в обмотках 11 и 12 из (13), (14) равны:
Р = г /2 - г/2 (15)
-'пИ '11-2-41 11' V1-','
Р = г /2 - г/2 (16)
П12 12-2 12 12 *
Используя (10), находим потери во вторичной обмотке:
Р = Р - Р - Р (17)
1 П2 1 1 П11 1 П12 * V1 V
Для частоты 1 Гц, для которой г2 = г (табл. 2), потери в обмотках, рассчитанные по (15)-(17), отличаются от точных значений (первый способ расчета) не более, чем на 1,5%. Следовательно, схема замещения позволяет в этом случае с высокой точностью рассчитывать потери в обмотках.
При сильном поверхностном эффекте параметр г может отличаться от сопротивления вторичной обмотки г2 из-за добавочных сопротивлений, вызванных полем рассеяния. Например, в опыте короткого замыкания трансформатора 11-2 между обмотками 11 и 2 находится обмотка 12 (рис. 2а), в которой наводятся вихревые токи полем рассеяния. Добавочное сопротивление Дг, обусловленное вихревыми токами в обмотке 12, которое пренебрежимо мало при 1 Гц, при 50 Гц (табл. 2) составляет почти 17%. Согласно (3) коэффициент г будет отличаться от сопротивления вторичной обмотки в этом режиме на величину Дг/2, поэтому:
г2 = г -Дг/ 2. (18)
Результаты расчета потерь по (15)-(17) для 50 Гц даны в табл. 6.
Потери в обмотках Losses in windings
Таблица 6 Table 6
f Гц P Вт ± П11' P Вт ± П12' P Вт П2 P, Вт
1 103 150 200 453
50 0,032 2,24 0,90 3,17
Совпадение активных мощностей и потерь, рассчитанных двумя способами (табл. 2 и 6), подтверждает справедливость и высокую точность методики расчета потерь в обмотках с использованием схемы замещения.
Важно отметить, что активные сопротивления в (15)-(18) могут быть измерены в опытах короткого замыкания двухобмоточных трансформаторов 11-2, 12-2 и 11-12. Далее измерив токи в обмотках в опыте короткого замыкания трехобмоточного трансформатора, можно найти потери в каждой его обмотке. Таким образом, (11)-(18) лежат в основе косвенного измерения потерь в обмотках.
Результаты расчета активной мощности и потерь в трансформаторе комплекса КСС-04 для стыковой контактной сварки труб
Применим разработанную методику к расчету и анализу активной мощности и потерь в трансформаторе КСС-04 для контактной стыковой сварки оплавлением труб диаметром 1420 мм, который разработан и испытан в ЗАО «Псковэлектросвар». На рис. 5а показано поперечное сечение, проходящее через ось симметрии этого трансформатора».
На тороидальный магнитопровод 1 (рис. 5а) намотана первичная обмотка 2, состоящая из двух параллельно соединенных секций. Вторичная обмотка 3 имеет один виток, который представляет собой медный корпус, охватывающий первичную обмотку. Поперечное сечение вторичной обмотки заштриховано. В этом корпусе сделан разрез 4, к краям которого в определенных местах периметра разреза прикрепляются гибкие медные шины 5. В данном трансформаторе шины установлены на одинаковом расстоянии в 18 местах периметра разреза. Для пояснения конструкции на рис. 1б показан внешний вид трансформатора с гибкими шинами, на котором виден корпус (вторичная обмотка) и гибкие шины, расположенные по периметру разреза.
К гибким шинам крепятся массивные проводники 6. Эти проводники плотно примыкают к предварительно зачищенным внутренним поверхностям свариваемых труб 7 и 8. При сближении труб ток во вторичной обмотке замыкается по гибким шинам 5, массивным проводникам 6 и стыку свариваемых труб. Высокая плотность тока в стыке труб приводит к оплавлению кромок труб, после чего трубы с усилием прижимаются и происходит их соединение. Первичная обмотка трансформатора выполнена из медного провода ПМГ ЛБН с прямоугольным сечением 20*9 мм2, количество витков первичной обмотки — 52. Номинальное первичное напряжение = 380 В, частота / = 50 Гц. Номинальная эквивалентная длительная мощность трансформатора составляет 577 кВА, номинальный эквивалентный длительный первичный ток — 1442 А.
Схема опыта короткого замыкания показана на рис. 6, результаты измерений в табл. 7. Активные мощности Рп, Рп и токи /и, /12 в секциях первичной обмотки 1, 2 измерялись клещами-ваттметрами АКИП2.
Входная активная мощность Р и входной ток I измерялись с помощью комплекта измерительного К-5403.
Рис. 5. Конструкция трансформатора КСС-04 Fig. 5. KSS-04 transformer design
Рис. 6. Схема опыта короткого замыкания 11, 12-секции первичной обмотки Fig. 6. Short circuit experiment circuit 11, 12 are primary winding sections
Таблица 7
Результаты опыта короткого замыкания
Table 7
Results of the short circuit test
I, А U, B P, кВт Ilv А А P„, кВт P12, кВт
847,5 25,6 6,2 495 405 6,9 -0,7
Данные табл. 7 показывают, что измеренные в секциях первичной обмотки активные мощности отличаются приблизительно в десять раз. Отрицательное показание ваттметра в секции 12 означает, что угол ф2 между входным напряжением и и током / в секции 12 больше 90°.
Для определения параметров схемы замещения выполнен расчет 2Б магнитных полей рассеяния трансформаторов 11-2, 12-2 и 11-12. Результаты расчета и измерения параметров схемы замещения отличаются на 10—12%. Параметры схемы замещения, найденные из опытов, приведены в табл. 8.
Рис. 7. Векторная диаграмма для схемы замещения Fig. 7. Vector diagram for equivalent circuit
Параметры схемы замещения
Equivalent circuit parameters
Таблица 8 Table 8
<üL112, Ом rl1-2, 0M o>L12-2, Ом ri22 Ом шМ, Ом r, Ом
0,028 0,013 0,032 0,012 0,027 0,004
Для данного трансформатора взаимная индуктивность по потокам рассеяния положительна, следовательно, отрицательная активная мощность может появиться в трансформаторе 12-2, так как L12_2 > Ь112. На рис. 7 показана векторная диаграмма, из которой можно найти активные составляющие напряжений на секциях первичной обмотки и активные мощности в них (табл. 9). Из диаграммы видно, что угол между первичным напряжением и током в секции 12 больше 90°.
Таблица 9
Активные мощности
Table 9
Active powers
I2 r ll 11-2> Вт IuIj2 rcos0, Вт (oMI11I12 sin0, Вт рш Вт I2 r 12 12-2' Вт IuIj2 rCOS0, Вт (oMI11I12 sin0, Вт pxv Вт P, Вт
3,1 0,5 3,1 6,7 2,0 0,6 -3,24 -0,64 6,06
Поскольку в обмотке 12 составляющая напряжения юМ/и8т0, имеющая отрицательное значение, превышает сумму двух падений напряжения /12г122 и /пгео^0, результирующее активное напряжение в обмотке 12 отрицательно, следовательно, и активная мощность в обмотке, которую оно определяет, также имеет отрицательное значение. Разница между рассчитанными и измеренным значениями активных мощностей в обмотке 11 составляет 3%, а в обмотке 12—9%. В обмотке 12 погрешность больше, чем в 11, поскольку потери в ней представляют собой разницу близких чисел. Для увеличения точности необходимо увеличить точность расчета и измерения параметров схемы замещения.
Найдем потери в обмотках трансформатора, используя измеренные активные мощности в обмотках и параметры схемы замещения. Поскольку первичная обмотка выполнена из литцен-драта, пренебрежем добавочными потерями в обмотке 12 в опыте короткого замыкания (режиме противовключения) трансформатора 11-2. Тогда по (11)-(17) получаем: Р = 2,2 кВт, РП12 = = 1,4 кВт, РП11 = 2,5 кВт. Эти потери являются исходными данными для дальнейшего расчета теплового режима трансформатора.
Заключение
В статье разработана методика расчета активных мощностей и потерь в трехобмоточном трансформаторе с двумя параллельными первичными обмотками, основанная на схеме замещения трансформатора [4] и анализе векторных диаграмм токов и напряжений в трансформаторе. Новизна методики состоит в установлении связи между активными мощностями в обмотках и потерями в них. Эта связь необходима как для расчета потерь в обмотках, так и для косвенного измерения потерь в обмотках. Благодаря применению этой методики удалось объяснить появление отрицательной активной мощности в одной из секций первичной обмотки. Методика расчета позволила установить новые соотношения между параметрами двухобмоточных трансформаторов, входящих в трехобмоточный трансформатор, при которых в одной из параллельно соединенных первичных обмоток возможно появление отрицательной активной мощности. Эти соотношения подтверждены результатами численного и натурного эксперимента. Численный эксперимент на модельном трансформаторе показал высокую точность расчета с использованием разработанной методики.
СПИСОК ИСТОЧНИКОВ
[1] Хоменко В.И., Теренин А.Н., Кучук-Яценко С.И., Швец Ю.И. Контактная стыковая сварка оплавлением современных трубопроводов большой мощности // Журнал нефтегазового строительства. 2014, № 3. С. 52-58.
[2] Журавлёв С.И. Пути решения проблем внедрения контактной стыковой сварки оплавлением на магистральных трубопроводах больших диаметров (обзор) // Изв. Высших учебных заведений. Машиностроение. 2013, № 8, с. 64-70.
[3] Рыськова З.С., Федоров П.Д., Жимерова В.И. Трансформаторы для электрической контактной сварки. Л.: Энергоатомиздат, 1990. 124 с.
[4] Сахно Л.И., Сахно О.И., Крылов М.С. Методика расчета токов в обмотках трансформаторов для стыковой сварки газопроводов // Известия РАН. Энергетика. 2019. № 5. С. 88-97.
[5] Вольдек А.И. Электрические машины: Л.: Энергия, 1978. 832 с.
[6] Kant P., Singh B. Multiwinding transformer FED CHB inverter with on-line switching angle calculation based SHE technique for vector controlled induction motor drive, IEEE Transaction on Industry Application. 56 (3)(2020) 2807-2815.
[7] Alvarez-Mario C., Leon F., Lopez-Fernandez X.M. Equivalent Circuit for Leakage Inductance of Mul-tiwinding Transformers: Unification of terminal and duality models, IEEE transactions on power delivery. 27 (1) (2012) 353-361.
[8] Leon F., Martinez J.A. Dual Three-winding Transformer Equivalent Circuit matching Leakage measurements, IEEE transaction on power delivery. 24 (1) (2009) 160-168.
[9] Бутырин П.А., Алпатов М.Е. Аналитическая теория трансформаторов. М.: Национальный исследовательский университет МЭИ. 2019. C. 112.
[10] Миронюк Н.Е., Соболев А.С., Пудов В.И. Расчетная модель для оценки характеристик электромагнитных трансформаторов тока // Электричество. 2016. № 2. С. 19-28.
[11] Шакиров М.А. Вектор Пойнтинга и новая теория трансформаторов. Часть 3. «Анатомия» трансформатора // Электричество, 2017, № 3, С. 37-49.
[12] Шакиров М.А. Вектор Пойнтинга и новая теория трансформатора // Ч. 9. Размещение обмоток на разных стержнях магнитопровода. 2020, № 1, С. 39-50.
[13] Модеров А.А. Теоретические основы электротехники. Теория электрических цепей. Лабораторный практикум. 2003. С. 195.
[14] Дубицкий С.Д. Инженерное моделирование квазистатического электромагнитного поля в программе ELCUT для задач электроники // Силовая электроника. № 6. 2017. С. 64-68.
[15] Kharlamova E., Sakhno L., Sakhno O. Calculation and Measurement of the Magnetic Flux in the Magnetic Core of Welding Transformers International Journal of Applied Engineering Research (IJAER). 11 (22) (2016) 11055-11059.
СВЕДЕНИЯ ОБ АВТОРАХ
САХНО Людмила Ивановна — профессор, Санкт-Петербургский политехнический университет Петра Великого, д-р техн. наук. E-mail: lsahno2010@yandex.ru ORCID: https://orcid.org/0000-0002-7204-5449
САХНО Ольга Ивановна — доцент, Санкт-Петербургский политехнический университет
Петра Великого, канд. техн. наук.
E-mail: Olasahno@mail.ru
ORCID: https://orcid.org/0000-0002-8879-986X
КРЫЛОВ Михаил Сергеевич — аспирант, Санкт-Петербургский политехнический университет Петра Великого, без степени. E-mail: akpo6at11@yandex.ru ORCID: https://orcid.org/0000-0001-7547-7709
REFERENCES
[1] V.I. Khomenko, A.N. Terenin, S.I. Kuchuk-Yatsenko, Yu.I. Shvets, Kontaktnaya stykovaya svarka oplav-leniyem sovremennykh truboprovodov bolshoy moshchnosti // Zhurnal neftegazovogo stroitelstva. 2014, № 3. S. 52-58.
[2] S.I. Zhuravlev, Puti resheniya problem vnedreniya kontaktnoy stykovoy svarki oplavleniyem na magis-tralnykh truboprovodakh bolshikh diametrov (obzor) // Izv. Vysshikh uchebnykh zavedeniy. Mashinostroyeni-ye. 2013, № 8, s. 64-70.
[3] Z.S. Ryskova, P.D. Fedorov, V.I. Zhimerova, Transformatory dlya elektricheskoy kontaktnoy svarki: L: Energoatomizdat, 1990. 124 s.
[4] L.I. Sakhno, O.I. Sakhno, M.S. Krylov, Metodika rascheta tokov v obmotkakh transformatorov dlya stykovoy svarki gazoprovodov // Izvestiya RAN. Energetika. 2019. № 5. S. 88-97.
[5] A.I. Voldek, Elektricheskiye mashiny. L.: Energiya, 1978. 832 s.
[6] P. Kant, B. Singh, Multiwinding transformer FED CHB inverter with on-line switching angle calculation based SHE technique for vector controlled induction motor drive, IEEE Transaction on Industry Application. 56 (3)(2020) 2807-2815.
[7] C. Alvarez-Mario, F. Leon, X.M. Lopez-Fernandez, Equivalent Circuit for Leakage Inductance of Multiwinding Transformers: Unification of terminal and duality models, IEEE transactions on power delivery. 27 (1) (2012) 353-361.
[8] F. Leon, J.A. Martinez, Dual Three-winding Transformer Equivalent Circuit matching Leakage measurements, IEEE transaction on power delivery. 24 (1) (2009) 160-168.
[9] P.A. Butyrin, M.Ye. Alpatov, Analiticheskaya teoriya transformatorov. M.: Natsionalnyy issledovatelskiy universitet MEI. 2019. S. 112
[10] N.Ye. Mironyuk, A.S. Sobolev, V.I. Pudov, Raschetnaya model dlya otsenki kharakteristik elektromag-nitnykh transformatorov toka // Elektrichestvo. 2016. № 2. S. 19-28.
[11] M.A. Shakirov, Vektor Poyntinga i novaya teoriya transformatorov. Chast 3. «Anatomiya» transformato-ra // Elektrichestvo, 2017, № 3, S. 37-49.
[12] M.A. Shakirov, Vektor Poyntinga i novaya teoriya transformatora // Ch. 9. Razmeshcheniye obmotok na raznykh sterzhnyakh magnitoprovoda. 2020, № 1, S. 39-50.
[13] A.A. Moderov, Teoreticheskiye osnovy elektrotekhniki. Teoriya elektricheskikh tsepey. Laboratornyy praktikum. 2003. S. 195.
[14] S.D. Dubitskiy, Inzhenernoye modelirovaniye kvazistaticheskogo elektromagnitnogo polya v programme ELCUT dlya zadach elektroniki // Silovaya elektronika. № 6. 2017. S. 64-68.
[15] E. Kharlamova , L. Sakhno, O. Sakhno, Calculation and Measurement of the Magnetic Flux in the Magnetic Core of Welding Transformers International Journal of Applied Engineering Research (IJAER). 11 (22) (2016) 11055-11059.
THE AUTHORS
SAKHNO Liudmila I. — Peter the Great St. Petersburg Polytechnic University.
E-mail: lsahno2010@yandex.ru
ORCID: https://orcid.org/0000-0002-7204-5449
SAKHNO Olga I. — Peter the Great St. Petersburg Polytechnic University.
E-mail: Olasahno@mail.ru
ORCID: https://orcid.org/0000-0002-8879-986X
KRYLOV Mikhail S. — Peter the Great St. Petersburg Polytechnic University.
E-mail: akpo6at11@yandex.ru
ORCID: https://orcid.org/0000-0001-7547-7709
Статья поступила в редакцию 06.07.2021; одобрена после рецензирования 22.07.2021; принята к публикации 27.09.2021.
The article was submitted 06.07.2021; approved after reviewing 22.07.2021; accepted for publication 27.09.2021.