Упрощенная SPICE-модель контура индукционного нагревателя
В статье описан способ создания упрощенной модели существующего контура индукционного нагревателя.
Контур заменен эквивалентной схемой, состоящей из согласующего трансформатора, конденсатора, эквивалентных R и L элементов,
величины которых определены по затухающим колебаниям, вызванным генератором ударного возбуждения.
Измерения производятся при помещении разных объектов в индуктор после чего формируются наборы измеренных эквивалентных параметров для каждого объекта.
Полученные параметры заносятся в LTSpice-модель. В итоговой модели не предполагается учет динамики изменения параметров в процессе нагрева объектов.
Исходные данные
Рисунок 1 - Вид контура индукционного нагревателя
Рисунок 2 - Геометрия индуктора
Dind | 40 мм | - внутренний диаметр индуктора |
hind | 65 мм | - высота индуктора |
Nind | 10 вит | - число витков индуктора |
Dwire | 3/16'' | - внешний диаметр проводника |
dwire | 2 мм | - внутренний диаметр проводника |
C | 2,0 мкФ | - емкость конденсаторной батареи |
Nprim | 26 | - число витков первичной обмотки трансформатора |
Nsec | 1 | - число витков вторичной обмотки трансформатора |
Lprim | 5,07 мГн | - индуктивность первичной обмотки трансформатора (измеренная вне контура) |
Rprim | 81 мОм | - активное сопротивление первичной обмотки трансформатора (измеренное при постоянном токе) |
При построении модели описанным способом будут использованы не все данные, здесь они приведены для общего представления.
Эквивалентная схема
Рисунок 3 - Эквивалентная схема контура
Эквивалентная схема содержит следующие элементы:
Rprim – сопротивление постоянному току первичной обмотки (исходные данные);
Lprim – индуктивность первичной обмотки (исходные данные);
Lsec – индуктивность вторичной обмотки (расчет);
C – контурный конденсатор (исходные данные);
Requ – эквивалентное сопротивление контура (измерения);
Lequ – эквивалентная индуктивность контура (измерения).
Для определения индуктивности вторичной обмотки (витка) воспользуемся измеренной величиной индуктивности первичной обмотки
и соотношением витков:
Эквивалентное сопротивление Requ и
эквивалентная индуктивность контура Lequ будут рассчитаны по
параметрам свободных колебаний, вызванных в контуре генератором ударного возбуждения. Индуктивность определяется
по резонансной частоте и известной емкости в контуре, а сопротивление - по добротности контура (подробнее об измерении добротности).
Для измерения параметров контура нужно исключить влияние согласующего трансформатора. Поскольку отделить трансформатор от контура, не нарушив целостность контура, не представляется возможным,
измерения будут проводится при коротком замыкании первичной обмотки трансформатора. При этом полное сопротивление витка вторичной обмотки в контуре станет настолько малым,
что в данном случае её влиянием можно пренебречь.
Рисунок 4 - Схема экспериментальной установки
Рисунок 5 - Осциллограмма свободных колебаний с пустым индуктором (пример)
С применением генератора ударного возбуждения снимаются осциллограммы затухающих колебаний при помещении в индуктор различных объектов.
Параметры осциллограмм заносятся в таблицу 2, к ним относятся: F0 - частота свободных колебаний;
VA - амплитуда первого или второго колебания;
VB - амплитуда колебания, следующего за предыдущим измеренным (VA).
По измеренным параметрам рассчитываются промежуточные (Q - добротность) и недостающие для модели параметры (Lequ, Requ):
Таблица 2 - результаты измерений параметров контура
Объект | VA, мВ | VB, мВ | F0, кГц | Q | Lequ, мкГн | Requ, мОм |
Воздух | 69,2 | 65 | 77,64 | 50,2 | 2,10106 | 20,417 |
Латунь | 65 | 60 | 79,61 | 39,2 | 1,99837 | 25,5 |
Болт М10 | 43,8 | 33,6 | 78,12 | 11,9 | 2,07532 | 85,601 |
Болт М14 | 44,8 | 34,8 | 80,12 | 12,4 | 1,97301 | 80,099 |
Тигель | 19,2 | 10 | 83,33 | 4,8 | 1,82393 | 198,952 |
Тигель + болт М10 | 18,4 | 9,8 | 84,45 | 5 | 1,77587 | 188,461 |
Рисунок 6 - Модель для сверки с экспериментом (inductor_1_model_8.1.zip)
Таблица 3 - сравнение осциллограмм свободных колебаний в модели и в эксперименте
Объект | Эксперимент | Модель |
Воздух | | |
Латунь | | |
Болт М10 | | |
Болт М14 | | |
Тигель | | |
Тигель + болт М10 | | |
В результате сравнения модели и эксперимента видно сходство осциллограмм, поэтому принимаем полученные эквивалентные параметы контура для применения в последующих моделях.
Далее производится эксперимент с генератором ударного возбуждения, подключенным уже к первичной обмотке трансформатора,
как изображено в модели на рисунке 7.
Рисунок 7 - Модель с подачей импульса на первичную обмотку трансформатора (inductor_1_model_8.2.zip)
На рисунке 8 изображены осциллограммы напряжения на контурном конденсаторе в эксперименте и в модели.
Видно различие скорости затухания. К тому же кривая затухания колебаний в эксперименте не описывается экспонентой (переменный коэффициент затухания),
что может быть следствием того, что эквивалентное активное сопротивление в системе, влияющее на скорость затухания, зависит от амплитуды колебаний.
Рисунок 8 - Осциллограмма напряжения на C1 при подаче импульса на первичную обмотку трансформатора
Различия осциллограмм колебаний, проявившиеся в присутствии согласующего трансформатора, могут говорить
о том, что в модели трансформатора не всё учтено. Подробно методика моделирования трансформатора описана в [1].
В следующей модели добавлена нелинейная индуктивность L4, которая учитывает гистерезисные свойства феррита
(Hc - коэрцитивная сила, А*вит/м; Bs - индукция насыщения, Тл; Br - остаточная индукция, Тл),
а также площадь поперечного сечения A,
длину средней линии Lm и число витков первичной обмотки N.
Параметры ферритов можно посмотреть в [2].
Трансформатор намотан на ферритовом сердечнике из 3-х колец
М2000НМ1 45х28х12, тогда:
В модель также могли быть добавлены индуктивности рассеяния первичной и вторичной обмоток,
но по результатам измерений они составили 50мкГн и 74нГн соответственно,
поэтому их решено не добавлять в модель из-за малой величины, в данном случае соизмеримой с погрешностью измерений.
Рисунок 9 - Модель с добавлением нелинейной индуктивности трансформатора (inductor_1_model_9.zip)
Рисунок 10 - Осциллограмма напряжения на C1 при подаче импульса на первичную обмотку трансформатора
В последнем случае осциллограммы имеют лучшее сходство, поэтому в качестве модели контура решено оставить вариант, изображенный на рисунке 9.
Далее производятся эксперименты с подключением испытуемого контура к макету мостового инвертора. Инвертор будет управляться генератором,
с регулируемой частотой. Будут производиться сравнения осциллограмм напряжения и тока на выходе инвертора
при помещении в индуктор объектов, под которые была ранее рассчитана модель.
Рисунок 11 - Модель с подачей импульса на первичную обмотку трансформатора
Рисунок 12 - Модель с подачей импульса на первичную обмотку трансформатора (inductor_1_model_9_+_small_bridge.zip)
Таблица 4 - осциллограммы тока и напряжения на выходе инвертора при пустом индукторе
Частота | Эксперимент | Модель |
Fgen=F0=77,7кГц | | |
Fgen=F0+1 кГц | | |
Fgen=F0-1 кГц | | |
Таблица 5 - осциллограммы тока и напряжения на выходе инвертора при наличии болта М10 в индукторе
Частота | Эксперимент | Модель |
Fgen=F0=78,12кГц | | |
Fgen=F0+1 кГц | | |
Fgen=F0-0,6 кГц | | |
Из сравнения осциллограмм модели и эксперимента в целом видно сходство как по амплитудам, так и по форме.
Поэтому полученная модель контура может быть применена для моделирования некоторых режимов работы установки.
Источники информации:
- Володин В.Я. Создание модели трансформатора в симуляторе LTspice
http://valvolodin.narod.ru/articles/nransformer.pdf
- Магнитомягкие ферриты для радиоэлектронной аппаратуры.
Справочник. М.М. Михайлова, В.В. Филиппов, В.П. Муслаков, Москва “Радио и связь”, 1983.
Комментарии
|