Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,737

THE MATHEMATICAL MODEL OF TRANSFORMER WHEN ITS DIAGNISIS BY SHORT DURATION PULS TAKES PLACE

Isaev Yu.N. 1 Kolchanova V.A. 1 Elgina G.A. 1 Lavrinovich A.V. 1
1 National Research Tomsk Polytechnic University
Приводится оригинальный синтез модели обмоток трехфазного двухобмоточного трансформатора. Ис-ходный сигнал входного воздействия в виде импульса короткой длительности раскладывался на со-ставляющие ряда Фурье. Для каждой составляющей в соответствии с принципом суперпозиции рассчи-тывался переходный процесс. Для определения функций токов и напряжений использовался метод про-странства состояний с учётом распределённости параметров схемы замещения, а также влияния скин-эффекта. В работе приведен пример расчета индуктивности и сопротивления витков обмотки трансфор-матора круглого и прямоугольного сечения с учетом скин-эффекта на основе расчёта поля распределе-ния магнитного потенциала c использованием метода конечных элементов. Для реализации метода ко-нечных элементов был использован пакет COMSOL Multiphysics. Приведены функции переменных со-стояния при входном воздействии в виде зондирующего импульса короткой длительности. По данным сравнительного анализа, полученные нами характеристики на основе моделирования и характеристики натурного эксперимента совпадают.
The synthesis of the model of winding of a three-phase two-winding transformer has been considered. The original input signal in the form of a short duration pulse into components of the Fourier series has been decomposed. In accordance with the superposition principle for each component of decomposition the transient process has been calculated. To determine the transients of currents and voltages the state space method have been used, the distribution of parameters of the equivalent circuit as well as the influence of the skin effect have been taken in account. The examples of the calculation of inductance and resistance of the transformer windings with circular and rectangular cross section with the skin effect by calculating the distribution of the magnetic field potential using the finite element method have been given. To implement the finite element method the COMSOL Mul-tiphysics soft ware has been used. The functions of the state variables with input action as a probe pulse of short duration are given. According to a comparative analysis, the modeling characteristics obtained by us and the characteristics of a natural experiment have been coincided.
transients.
Maxwell equations
skin effect
state space method
wire cross-section
transformer diagnostics
Тransformer

Введение. Для уменьшения затрат на ревизию трансформаторов предлагается исследовать обмотки трансформаторов при подключении импульсного напряжения [1, 2, 5]. Длина импульса – теста составляет 400-500 нс.

Цель исследования создать математическую модель трансформатора. Модель можно использовать для моделирования дефектов обмоток, формировать эталонные сигналы при различных повреждениях исследуемого оборудования.

Вывод основных соотношений

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

Рис. 1. Обмотки трехфазного двухобмоточного трансформатора и поперечные сечения проводов обмоток, а – внутренняя обмотка трансформатора (НН); б – внешняя обмотка трансформатора (ВН); в – трехфазный трансформатор; г – поперечное сечение провода обмоток соответственно.

1. Электротехническая схема замещения одного витка обмотки

Виток обмотки заменяется электротехническим звеном, состоящим из элементов, приведенных на рисунке 2, а. Здесь – межвитковая емкость, – емкость между витком и заземленными частями трансформатора, – индуктивность рассеяния обмотки с учетом взаимной индуктивности между соседними элементами обмотки, – активное сопротивление витка соответственно.

а

б

Рис. 2. Электротехническая схема замещения: а – одного витка обмотки трансформатора; б – схема замещения обмотки трансформатора

Все приведенные величины являются погонными, т.е. значения величин отнесены к единице длины. Таким образом, обмотку трансформатора можно заменить последовательно соединенными звеньями (рисунок 2, б).

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

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

2. Частотное разложение зондирующего импульса.

В качестве модельного зондирующего импульса выберем функцию с длительностью мкс в виде:

а

б

в

Рис. 3. Зондирующий импульс (t в мкс): а – исходный импульс показан сплошной линией, пунктирной показан зондирующий импульс, полученный в среде MathCAD (t в мкс)); б - амплитудно-частотный спектр; в – фазо-частотный спектр

Определим частотный спектр импульса. Для этого разложим импульс в ряд Фурье, определив коэффициенты разложения . Для определения коэффициентов использовалась стандартная программа быстрого преобразования Фурье (FFT) программно-интегрированной среды MathCAD. Далее функция представлялась в виде ряда:

где

Спектр и результат восстановления функции по коэффициентам представлены ниже.

Число коэффициентов разложения бралось N=30. Из спектра импульса видно, что начиная с 10 гармоники, относительный вклад высших гармоник в импульсе не превышает 5%.

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

Проведем расчет плотности распределения тока в поперечном сечении проводника. Система уравнений Максвелла [3, 4] позволяют получить уравнение для векторного потенциала , позволяющего определить плотность распределения тока по сечению проводника:

, (1)

где – векторный магнитный потенциал, – напряжения проводника, диэлектрическая проницаемость, – частота, – магнитная проницаемость, – проводимость, мнимая единица.

Плотность тока и ток определяются выражениями [3]:

. (2)

Магнитная энергия определяется выражением [4]:

, (3)

из которого при известном токе определяется индуктивность с помощью соотношения [3]:

(4)

При известном напряжении и токе определяем сопротивление проводника, используя формулу:

(5)

Для расчета представленных уравнений использовался математический пакет COMSOL Multihpysics, основанный на методе конечных элементов (Finite Element Method (FEM)). Результаты расчетов приведены на рисунках 4, 5.

Рис. 4. Поле распределения магнитного потенциала кольцевых проводников: а – круглого и б – прямоугольного сечения при фиксированной частоте.
Распределение плотности тока в поперечном сечении проводника при частоте Гц, в – круглого; г – прямоугольного соответственно

а

б

Рис. 5. Частотные зависимости индуктивности L(f) и сопротивления R(f) провода: а – с круглым поперечным сечением; б – с прямоугольным поперечным сечением

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

Запишем уравнение дифференциальные уравнения для определения напряжения и токов в реактивных элементах схемы замещения (рисунок 6) для одного витка с учетом внутреннего сопротивления источника r. Для решения дифференциальных уравнений использовался метод пространства состояний, дифференциальное уравнение третьего порядка представлялось в виде системы 3-х уравнений первого порядка. Уравнения рассчитывались численно с использованием метода Рунге-Кутта четвертого порядка с фиксированным шагом. Число точек дискретизаций равнялось N=103 .

Рис. 6. Электротехническая схема замещения одного витка обмотки трансформатора

а

б

Рис.7. Напряжение на конденсаторе : а – в результате моделирования; б – экспериментальные осциллограммы напряжений на витке при зондировании пробным импульсом

Здесь – гармоническая составляющая входного импульса при фиксированной частоте, – ток индуктивности, – напряжение на первом и втором конденсаторах. Результирующие токи и напряжения представляются в виде ряда:

рад/сек.

Результаты расчетов для провода круглого сечения приведены на рисунке 7, а.

График напряжения на конденсаторе С1, рассчитанный теоретически, удовлетворительно совпадает с экспериментом. Что свидетельствует об адекватности модели.

Для провода с прямоугольным поперечным сечением графические зависимости выглядят аналогично.

Выводы

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

Рецензенты:

Усов Ю.П., д.т.н., профессор кафедры ЭСиЭ ЭНИН ФГБОУ ВПО «НИ ТПУ», Национальный исследовательский Томский политехнический университет, г. Томск.

Сивков А.А., д.т.н., профессор кафедры ЭПП ЭНИН ФГБОУ ВПО «НИ ТПУ», Национальный исследовательский Томский политехнический университет, г. Томск.