CC BY
210Моделирование сложных электромагнитных компонентов
при помощи Spice-симулятора LTspice/SwCAD III
Валентин ВОЛОДИН
В статье описывается моделирование трансформаторов и индукторов, имеющих сложные сердечники, средствами бесплатного симулятора LTspice/SwCAD III. Описанная в статье методика не ориентирована на специфические особенности LTspice/SwCAD III и поэтому может быть использована и в других подобных симуляторах.
Основные типы магнитопроводов
Практически все известные виды вторичных источников питания содержат электромагнитные компоненты, такие как трансформаторы и индуктивности (дроссели, индукторы, реакторы). Обычно эти компоненты имеют одну или несколько обмоток, намотанных поверх магнитопровода из ферромагнитного материала. Использование ферромагнитных материалов позволяет улучшить электрические параметры электромагнитных компонентов, а также уменьшить их размеры и массу.
По аналогии с электрическими цепями, маг-нитопроводы электромагнитных компонентов иногда называют магнитными цепями. В свою очередь, магнитные цепи, как и электрические, могут быть неразветвленными и разветвленными. На рис. 1 показаны два типа трансформаторов, один из которых имеет неразветвленный (рис. 1а), а другой разветвленный (рис. 1б) магнитопроводы. Здесь также указаны основные уравнения, соответствующие этим магнитопроводам.
На практике трансформаторы, используемые в источниках питания, имеют в основном неразветвленный магнитопровод или с небольшими допущениями могут таковыми считаться.
Трансформаторы с разветвленным магни-топроводом обычно используются в применениях с количеством фаз 3 и более.
Существуют также электромагнитные компоненты, имеющие сложный магнитопро-вод, являющийся комбинацией неразветв-ленных и разветвленных участков. Подобные магнитопроводы имеют различные интегрированные электромагнитные компоненты, совмещающие в себе несколько электромагнитных узлов, например, трансформатор и дроссель.
Встроенные средства SPICE-симуляторов, предназначенные для моделирования электромагнитных компонентов
Для моделирования электромагнитных компонентов во всех SPICE-симуляторах есть
стандартные библиотечные модели. Например, в симуляторе LTspice/SwCAD III такими моделями являются:
1. Модель линейной индуктивности.
2. Модель нелинейной индуктивности.
3. Модель магнитной связи.
Модель линейной индуктивности
Модель линейной индуктивности в LTspice/ SwCAD III имеет вид, представленный на рис. 2. Здесь, кроме индуктивности (Inductance), учитывается масса паразитных параметров, таких как последовательное сопротивление обмотки Rser (Series Resistance), сопротивление потерь между выводами обмотки Rpar (Parallel Resistance) и межвитковая емкость Cpar (Parallel Capacitance). Настроить основные и паразитные параметры модели можно в окне настройки (рис. 3), которое вызывается щелчком правой кнопкой мышки по индуктивности L1 (или L2). Здесь индуктивность измеряется в генри, емкость в фарадах, а сопротивление в омах. Опция Show Phase Dot позволяет делать видимой и скрывать индикацию начала обмотки.
о
[ 1 Rpar ] Rser Cpar <lnductance>
о
Рис. 2. Модель линейной индуктивности, используемая в программе LTspice/SwCAD III
Таблица. 2. Линейные размеры сердечника
Параметр Описание Единицы измерения
Hc Коэрцитивная сила А/м
Br Остаточная индукция Тл
Bs Индукция насыщения Тл
Параметр Описание Единицы измерения
Lm Средняя длина магнитной линии (исключая немагнитный зазор) м
Lg Длина немагнитного зазора м
A Сечение сердечника м2
N Количество витков -
Рис. 3. Меню настройки индуктивности
Указание паразитных параметров элементов позволяет получить более реалистичную картину моделирования. Некоторые паразитные параметры ЬТбрке/З'^АО III определяет сам, если пользователь не указывает их конкретную величину. Например, по умолчанию, LTspice/SwCAD III считает, что К^ег = 0,001 Ом.
Модель нелинейной индуктивност
В Spice-симуляторах для моделирования нелинейной индуктивности с ферромагнитным сердечником (далее по тексту индуктор) чаще всего используется модель Джилса-Атертона.
В LTspice/SwCAD III используется упрощенная гистерезисная модель магнитного сердечника, которую в 1991 году предложил Джон Чан [1]. Модель использует только основные параметры петли гистерезиса, перечисленные в таблице 1 и показанные на рис. 4.
Наряду с параметрами петли гистерезиса, при описании индуктора учитываются ли-
Таблица 1. Основные параметры петли гистерезиса
Модель магнитной связи
Модель магнитной связи позволяет создавать идеальный трансформатор, используя модели нескольких линейных индуктивностей. На рис. 5 изображен стандартный способ задания идеального трансформатора.
нейные размеры сердечника, а также количество витков катушки (табл. 2).
Допустим, нам нужно моделировать индуктивность с обмоткой из N = 21 витков, которая намотана на магнитопровод из феррита М3000НМС1, состоящий из четырех комплектов ПК40х18. Магнитопровод ПК40х18 имеет площадь сечения А = 2,2 см2 и среднюю длину магнитной силовой линии Lm я 200 мм. В отличие от обычных ферритов, М3000НМС1 может использоваться в сравнительно сильных магнитных полях и имеет следующие параметры петли гистерезиса: Б8 = 0,45 Тл (при Н = 800 А/м), Вт я 0,33 Тл (при Н = 100 А/м и Т = 60 °С), Бг = 0,1 Тл и Нс = 12 А/м. Если в магнитопроводе предусмотрен немагнитный зазор длиной Lg = 0,07 мм, то строка настройки модели будет выглядеть следующим образом:
Нс=12 Б8=0,45 Вг=0,10 А=0.00088 Lm=0,2 Lg=0,00007 N=21.
Эту строчку нужно ввести в качестве параметра Inductance (H) в меню настройки индуктивности (рис. 3).
Способ задания индуктора, используемый в LTspice/SwCAD III, достаточно привлекателен, так как при описании сердечника использует доступные справочные данные, а также понятные и легко контролируемые параметры самого сердечника и катушки.
Сравнение результатов моделирования с практическими результатами показало высокую достоверность используемой модели, которая в некоторых случаях даже превосходит достоверность модели Джилса-Атертона.
В данном случае трансформатор имеет две обмотки— Ь1 и Ь2. Индуктивности этих обмоток имитируют соответствующие индуктивности намагничивания трансформатора. Для создания магнитной связи между обмотками трансформатора служит фиктивный компонент К, в качестве параметров которого указаны индуктивности, между которыми нужно осуществить связь, и соответствующий коэффициент связи, который может принимать значения от 0 до 1. Для реальных трансформаторов этот коэффициент связи обычно находится в диапазоне 0,9-0,999.
Проблемы корректного моделирования трансформатора с разветвленным магнитопроводом
Идеальный трансформатор позволяет моделировать трансформатор с неразветвлен-ным магнитопроводом, работающим на линейном участке. Попытка таким же образом моделировать трансформатор с разветвленным магнитопроводом требует большой осторожности, так как может привести к серьезным ошибкам, одна из которых показана на рис. 6.
Согласно рис. 6, взаимная фазировка обмоток 3-обмоточного и трехстержневого трансформатора зависит от того, какая об-
Рис. 4. Параметры петли гистерезиса магнитного материала
мотка в данный момент подключена к источнику напряжения. Например, если источник напряжения подключен к обмотке N2 на среднем стержне (рис. 3а), то в рассматриваемый момент времени на выводе 1-Ш будет «-», а на выводе 2-Ш будет «+». Если источник напряжения подключить к обмотке N3 на крайнем стержне (рис. 3б), то фаза напряжения на выводах N1 будет сдвинута на 180°, то есть теперь на 1-Ш будет «+», а на 2-Ш будет «-». Если же этот трансформатор моделировать как неразветвленный, то ничего подобного происходить не будет. В этом случае взаимная фазировка обмоток будет сохраняться независимо от того, к какой из обмоток подключен источник напряжения.
Следовательно, для корректного моделирования электромагнитных компонентов с разветвленным или сложным магнитопро-водом необходимо создавать специальные модели, способные имитировать их физическую структуру. Наиболее распространенным является метод создания модели, использующий модель магнитного сопротивления. Этот метод имеет преимущество благодаря простоте генерации модели. При этом между имитационной моделью и физическим компонентом сохраняется вполне очевидная взаимосвязь, что упрощает процесс создания модели, а также позволяет избежать многих ошибок.
Принцип подобия электрических и магнитных цепей
3. Электрический ток (I) и магнитный поток (Ф).
4. Плотность электрического тока (I) и плотность магнитного потока (Б).
5. Электрическое сопротивление (Я) и магнитное сопротивление (Я).
6. Проводимость (8) и проницаемость (ц). Сравнение электрических и магнитных величин произведено в таблице 3.
Таблица 3. Сравнение между электрическими и магнитными параметрами цепей
Электричество Магнетизм
V = электрическое напряжение (электродвижущая сила); Е = напряженность электрического поля V = -|Ехс11с= Ех1с; Е = V/!, р = N1 = магнитное напряжение (магнитодвижущая сила); Н = напряженность магнитного поля Р=|Нхс1и=Нх|т; Н = Р/!т = М/!т
J = плотность тока; J = аЕ; а = удельная электропроводность В = плотность магнитного потока; В = |Н, где | = |0|г — абсолютная проницаемость; | — относительная проницаемость; |0 — проницаемость вакуума; |0= 4пх10-7 Н/т
I = электрический ток; 1=рхс15=^Ас, где А, — сечение проводника ф = магнитный поток; Ф = | Вхс15= ВхАт, где Ас — сечение магнитопровода
R = сопротивление; = V/! = !с/(аАс); Э = 1^ = проводимость R = магнитное сопротивление; R = Р/ф = !т/(|хАт) = N'/1; Р = 1^ = проницаемость
Индуктивность магнитной цепи непосредственно связана с магнитным сопротивлением Я и количеством витков в обмотке N
L = №/Я = №Р. (1)
Процесс создания модели связан с трансляцией физической структуры магнитного компонента в электрическую схему. Подобное преобразование возможно благодаря принципу подобия, согласно которому к расчету магнитных цепей применимы все приемы расчета электрических цепей с нелинейными сопротивлениями [2]. В этом случае используются следующие аналогии между электрическими и магнитными величинами:
1. Электродвижущая сила — ЭДС (V) и магнитодвижущая сила — МДС (Б).
2. Напряженность электрического поля (Е) и напряженность магнитного поля (Н).
Аналогично для взаимной индуктивности:
Ми= да2)/Д12 = N^2.
Воспользовавшись принципом подобия, создадим модель индуктора с неразветвленным магнитопроводом различного сечения и немагнитным зазором, общий вид которого изображен на рис. 7а. Обычно модель разрабатывается в следующей последовательности:
1. Магнитопровод делится на части, включая немагнитный зазор, и каждой части присваивается соответствующее магнитное сопротивление (рис. 7б).
Рис. 7. Моделирование индуктора с неразветвленным магнитопроводом и немагнитным зазором
2. Вычисляется магнитное сопротивление каждой части. Для участков магнитопровода:
Я1 = 11/(цеИгАт1); К2 = 12/(ц0 Цг Ат2)’ Я3 = 13/(ЦеИгАт3); К5 = 15/(ЦеИгАт5); Я6 = 16/(Ц0ИгАтб)-
Для немагнитного зазора:
К4 = 14/(Ц0Ат6)-
3. Присваивается магнитное напряжение источнику с ампервитками Б = №.
4. Рисуем эквивалентную цепь, показанную на рис. 8.
Рис. В. Результирующая модель магнитной цепи индуктора
Некоторые симуляторы, такие как SIMPLORER, используют именно этот принцип для моделирования сложных магнитных цепей. Для этого в библиотеке симулятора есть соответствующие элементы, такие как источник магнито-движущей силы (MMF), источник магнитного потока (FLUX), магниторезистор (RMAG), обмотка (ELTOMAG) и т. д. И все же, если вам понадобилось создать модель сложного трансформатора или дросселя, то можете пока не торопиться с покупкой и изучением достаточно дорогого и сложного симулятора. С таким же успехом подобное моделирование можно произвести и в понравившемся нам LTspice/SwCAD III. Для создания реалистичной модели разветвленного магни-топровода нам будет достаточно всего лишь его модели нелинейного индуктора.
Дуальность физических цепей
Под дуальностью или двойственностью обычно понимают аналогичность математи-
Таблица 4. Дуальные пары
Величина ^ Дуальная величина
V
V
Р Э = 1Д
Э Р = 1/Э
С 1_
і_ С
Разомкнутая цепь Замкнутая цепь
Замкнутая цепь Разомкнутая цепь
Генератор напряжения Генератор тока
Генератор тока Генератор напряжения
Ветвь Узел
Узел Ветвь
ческого описания различных процессов, происходящих в различных физических цепях. В таблице 4 указаны типичные дуальные пары для электрических цепей.
Модель неразветвленной магнитной цепи
Чтобы преобразовать модель магнитной цепи в соответствующую ей электрическую схему, необходимо проделать следующие шаги [3]:
1. Исходя из структуры моделируемого электромагнитного устройства, начертить соответствующую схему магнитной цепи (рис. 7, 8).
2. Используя двойственность, конвертировать модель сопротивлений (Я) в модель проводимостей (Р). В этом случае последовательные цепи превращаются в параллельные, магнитный поток превращается в магнитное напряжение, а источник магнитного напряжения превращается в источник тока (рис. 9а).
3. Чтобы убрать витки из источника тока, умножим величины магнитных проводимостей на количество витков N (рис. 9б).
4. Согласно закону Фарадея, в течение времени т всякое изменение магнитного потока ДФ, пронизывающего обмотку N независимо от причины, вызвавшей изменение потока, наведет в обмотке напряжение V [4]:
V = (ШФ)/т. (2)
Таким образом, магнитный поток имеет вольт-секундную зависимость, и его изменение зависит только от количества витков обмотки, намотанной поверх магнитопровода, и времени, в течение которого к этой обмотке приложено напряжение. Материал сердечника, а также его длина и сечение не оказывают какого-либо влияния на изменение потока. Если считать, что нас не ограничивают временные рамки и магнитный поток изменяется от 0, то без всякой потери точности вместо ДФ можно использовать Ф и записать закон Фарадея в несколько упрощенной форме: V = ^Ф. Теперь, чтобы перейти к обмоточному напряжению, каждый элемент схемы и поток Ф умножим на количество витков обмотки N. В результате схема теперь будет представлена обмоточным напряжением и магнитными проводимостями в масштабе № (рис. 9в).
5. Из уравнения (1) мы знаем, что L = №хР. Таким образом, все масштабированные проводимости можно заменить соответствующими индуктивностями (рис. 9г). Таким образом, последовательную магнитную цепь, состоящую из к различных фрагментов, можно моделировать с помощью к отдельных индукторов, каждый из которых соответствует определенному фрагменту магнитной цепи, то есть имеет соответствующую длину, сечение и материал магнитопровода.
Подтвердим полученные выводы с помощью следующей цепочки логических рассуждений.
Так как цепь не разветвленная, то через нее протекает один и тот же поток Ф. Возникает вопрос, как создать одинаковый поток в маг-нитопроводах независимых индукторов? В решении этой задачи нам поможет тот же закон электромагнитной индукции (2), согласно которому мы можем получить равный поток в различных магнитопроводах, если намотаем на них обмотки с равным количеством витков и соединим эти обмотки параллельно.
Теперь, обеспечив равенство потоков, решим проблему корректного моделирования намагничивающего тока.
На каждом из участков различного сечения поток Ф вызывает различную индукцию:
Вк = Ф^к.
В свою очередь, магнитная индукция Бк через кривую намагничивания конкретного фрагмента связана с напряженностью магнитного поля Нквнутри фрагмента магнитопро-вода. Зная напряженность магнитного поля и длину фрагмента, можно найти магнитное напряжение Н^к = Бк, которое прикладывается к каждому из фрагментов. Так как фрагменты магнитопровода включены последовательно, то общее магнитное напряжение будет равно сумме напряжений, требуемых для создания потока в каждом фрагменте:
!Бк = Б.
Ввиду того, что токи намагничивания отдельных фрагментов суммируются, это позволяет моделировать их при помощи той же модели, которая моделирует равные потоки, и в которой отдельные фрагменты моделируются независимыми индукторами с одинаковыми обмотками, включенными параллельно (рис. 10):
1Бк = К^) = !х№
Если на магнитопроводе намотано несколько обмоток, то для получения корректных напряжений можно использовать идеальные трансформаторы напряжения.
Моделирование
разветвленной магнитной цепи
Чтобы преобразовать модель разветвленной магнитной цепи (рис. 11а) в соответствующую ей электрическую схему, необходимо проделать шаги, аналогичные проделанным ранее с неразветвленной магнитной цепью:
1. Исходя из структуры моделируемого электромагнитного устройства, начертить соответствующую схему магнитной цепи (рис. 11б).
Рис. 12. Последовательность получения модели разветвленной магнитной цепи
2. Используя двойственность, конвертировать модель сопротивлений (Я) в модель проводимостей (Р). В этом случае последовательные цепи превращаются в параллельные, магнитный поток превращается в магнитное напряжение, а источник магнитного напряжения превращается в источник тока (рис. 12б).
3. Чтобы убрать витки из источника тока, умножим величины магнитных проводимостей на количество витков N (рис. 12в).
4. Чтобы перейти к обмоточному напряжению, каждый элемент схемы и поток Ф умножим на количество витков обмотки N. В результате схема теперь будет представлена обмоточными напряжениями и магнитными проводимостями в масштабе № (рис. 12г).
5. Из уравнения (1) мы знаем, что L = №хР. Таким образом, все масштабированные проводимости можно заменить соответствующими индуктивностями (рис. 12д). Как и в предыдущем случае, можно рассуждать несколько иначе. Представим себе некую магнитную цепь, состоящую из к (в данном случае к = 2) включенных параллельно фрагментов, имеющих различную длину и сечение (рис. 13).
Так как фрагменты магнитной цепи включены параллельно, то к каждому из них приложено одинаковое магнитное напряжение. Одинаковое магнитное напряжение на каждом фрагменте разветвленной магнитной цепи можно смоделировать, соединив к индукторов, имеющих обмотки с равными витками, последовательно. (Фрагменты реального магнитопровода, не имеющие обмотки, в модели все равно должны быть представлены указанным образом.) Так как при этом через обмотки протекает одинаковый ток I, то к сердечникам прикладывается одинаковое магнитное напряжение !хМ
Моделирование сложной магнитной цепи
Если параллельная магнитная цепь включена последовательно с другой, например, с последовательной цепью (рис. 14), то необходимо обеспечить равенство общего магнитного потока протекающего через них.
Для начала считаем, что Ян на рис. 14 имеет бесконечно большое сопротивление и никак не влияет на картину происходящего. В параллельной магнитной цепи общий поток является суммой потоков, протекающих по каждому параллельному фрагменту к:
Ф = ЕФк-
Параллельная цепь моделируется последовательным включением к дросселей, имеющих одинаковые обмотки. В этом случае элементарные ЭДС Vk, наводящиеся в каждом дросселе, складываются. На основании закона электромагнитной индукции (2), определим суммарный поток для этого случая:
ф = !Фк = ^ктт = ^т)/№
То есть, для обеспечения равенства потоков достаточно для всех дросселей модели указать одинаковое количество витков обмотки и модели отдельных магнитных ветвей включить параллельно. Кроме этого, если фрагмент магнитопровода равного сечения имеет немагнитный зазор, то для его имитации можно использовать один индуктор, что значительно упрощает процесс создания модели. В результате получаем модель (рис. 15).
1_3
Рис. 15. Модель совместной
последовательно-параллельной магнитной цепи
Из модели (рис. 15) видно, что дроссель L3, по сути, моделирует поток рассеяния Ф3 между катушками ю1 и ю2. Поэтому, при подключении нагрузочного резистора Ян к обмотке ю2, напряжение на ней падает.
Выводы
Каждый фрагмент реального магнитопро-вода в модели можно представить элементарной нелинейной индуктивностью (индуктором), имеющей материал магнитопровода, его длину, немагнитный зазор и сечение, соответствующие моделируемому фрагменту. Согласно принципу двойственности, все па-
|?Н
Рис. 14. Пример сложной последовательно-параллельной магнитной цепи
раллельные ветви реального магнитопрово-да в модели представляются соответствующими индукторами, включенными последовательно, и наоборот. На все индукторы модели мотаются обмотки с одинаковым количеством витков, равным количеству витков одной из обмоток прототипа. При необходимости, для согласования модели с остальными элементами схемы, можно использовать дополнительные идеальные согласующие трансформаторы.
Адаптация модели для магнитных цепей, работающих с частичным или полным подмагничиванием
Все предыдущие рассуждения основаны на законе электромагнитной индукции и справедливы для случая, когда к обмоткам магнитного компонента приложено переменное напряжение. Если к обмотке приложено постоянное напряжение или ненулевая постоянная составляющая переменного, то распределение потока в его магнитных цепях определяется уже не законом электромагнитной индукции, а их магнитным сопротивлением и законом Ома для магнитных цепей. Что касается модели, то в ней распределение потоков определяется омическим сопротивлением обмоток индукторов. Следовательно, чтобы модель корректно отражала процессы, протекающие на постоянном токе в реальной магнитной цепи, достаточно последовательно с каждой обмоткой модели включить нелинейное сопротивление, обратно пропорциональное магнитному сопротивлению моделируемого фрагмента магнитной цепи. Так как эта часть модели работает на постоянном токе, то достаточно, чтобы нелинейное со-
У[пОО?|
tran 0 1 0.99 1и
Нс=12 Bs=0.45 Br=0.10 A=1 Lm=1 Lg=0 N=1
Рис. 16. Тестовый узел
Таблица 5. Начальная кривая намагничивания феррита М3000НМС1, полученная программой Иврісе/в^СЛО III
B, Тл -0,465 -0,45 -0,4 -0,3 -0,2 -0,1 0 0,1 0,2 0,3 0,4 0,45 0,465
H, А -10 000 -2500 -320 -85 -35 -15 0 15 35 85 320 2500 10 000
противление отражало начальную кривую намагничивания (без гистерезиса).
Для снятия кривой намагничивания можно использовать тестовый узел (рис. 16). Для примера, проведем тестирование сердечника из феррита М3000НМС1. Сердечник имеет сечение А = 1 м2, среднюю длину магнитной силовой линии Lm = 1 м и на него намотана обмотка в один виток.
Тестовый узел построен на двух источниках тока, G1 и G2, управляемых напряжением, которые используются для измерения и нормирования отображаемых параметров.
Известно, что скорость изменения индукции в сердечнике прямо пропорциональна приложенному напряжению V и обратно пропорциональна площади сечения A и виткам обмотки N:
dB/dt = V/(AN).
Следовательно, подав напряжение V на вход интегратора, на его выходе мы получим сигнал, пропорциональный индукции в сердечнике. В нашем случае роль интегратора возложена на узел, состоящий из источнике G1 и конденсатора С2. Коэффициент передачи управляемого источника тока G1, обеспечивающий выходное напряжение интегратора, равное индукции, можно вычислить по формуле:
KG1 = C2/(AxN) = (1х10-6)/1х1 = = 1х10-6 = 1ц.
Вычисленное значение коэффициента передачи необходимо записать в строке Value меню настройки управляемого источника тока G1.
Коэффициент передачи управляемого источника тока G2, обеспечивающий выходной ток, равный напряженности в сердечнике нелинейного трансформатора, можно вычислить по формуле:
KG2 = N/(R2xLm) = 1/(0,1х1) = 10.
Вычисленное значение коэффициента передачи необходимо записать в строке Value меню настройки управляемого источника тока.
Для устранения постоянной составляющей на выходе интегратора конденсатор С2 зашун-тирован высокоомным резистором R5. В настройках горизонтальной оси, в строчке Quantity Plotted, вместо параметра time впишем параметр I(G2). По вертикали выводим напряжение на выходе интегратора, кликнув по верхнему выводу конденсатора С2 (рис. 15). При этом ток I(G2) равен напряженности магнитного поля, а напряжение на конденсаторе
С2 равно магнитной индукции в магнитопро-воде. В процессе тестирования необходимо снять ряд предельных циклов перемагничи-вания сердечника для различных значений максимальной индукции, вплоть до индукции насыщения. Полученные данные записываются в таблицу 5.
Для создания нелинейного сопротивления используем стандартный источник тока, управляемый напряжением. Для формирования требуемой вольтамперной характеристики (ВАХ) нелинейного сопротивления проще всего использовать табличную форму задания. Для примера ниже приведена строка настройки нелинейного резистора, имитирующего полученную кривую намагничивания магнитопровода из феррита М3000НМС1. В строке настройки указан ряд пар значений, первое значение пары указывает управляющее напряжение, пропорциональное индукции в сердечнике, а второе значение — соответствующую току источника тока, пропорциональному напряженности магнитного поля:
1аЫе=(-0.465 -10000,-0.45 -2500,-0.4 -320,-0.3 -85,-0.2 —35,—0.1 -15,0 0,0.1 15,0.2 35,0.3 85,0.4 320,0.45 2500,0.465 10000)
На рис. 17 показана результирующая ВАХ нелинейного резистора, снятая при помощи программы LTspice/SwCAD III.
Для более детальной аппроксимации достаточно лишь расширить таблицу дополнительными данными.
Так как размеры используемых магнито-проводов отличаются от тех, которые были приняты при тестировании, то перед установкой в модель каждый нелинейный резистор должен быть настроен в соответствии с параметрами моделируемого фрагмента магнито-провода и соответствующей нелинейной индуктивностью, последовательно с которой он будет включен. Величину управляющего напряжения V* для нового сечения А магнито-провода можно определить по формуле:
V* = VхA. (3)
Соответствующую величину тока I* для магнитопровода, имеющего длину Lm, немагнитный зазор Lg и обмотку N можно рассчитать по формуле:
I* = (IхLm+Bх796 000хLg)/N. (4)
Для правильного распределения магнитных потоков важны относительные значения нелинейных резисторов. При выборе абсолютных значений сопротивлений надо стремиться к тому, чтобы общее сопротивление для группового включения нелинейных ре-
Рис. 17. Кривая намагничивания феррита М3000НМС1, имитируемая при помощи нелинейного резистора, построенного на источнике тока, управляемого напряжением
зисторов примерно соответствовало активному сопротивлению соответствующей обмотки моделируемого электромагнитного компонента.
Согласование величин нелинейных сопротивлений с активным сопротивлением обмотки И возможно с помощью некоторого поправочного коэффициента К, на который перемножаются все скорректированные управляющие напряжения всех нелинейных сопротивлений модели. В свою очередь, выбор поправочного коэффициента вызывает определенные трудности, так как после пересчета нелинейных сопротивлений их координаты не совпадают. Для выбора поправочного коэффициента, например, можно определить общее сопротивление всех нелинейных резисторов Инел в режиме без насыщения, и затем поделить на него сопротивление обмотки:
К = И/Инел.
Процесс можно упростить, если в магни-топроводе имеется фрагмент с немагнитным зазором. В этом случае можно выбрать пару У*х1*, соответствующую фрагменту с зазором, которая находится на линейном участке кривой намагничивания. Далее рассчитывается поправочный коэффициент К = (Ях1*)/У*, на который затем умножить все ранее рассчитанные величины У* всех нелинейных сопротивлений модели.
Падение напряжения на нелинейном резисторе пропорционально магнитному потоку, протекающему через соответствующий дроссель. Собственно, при постоянном токе о величине магнитного потока можно судить только по падению напряжения на этом резисторе. Поэтому иногда величину нелинейного резистора стоит выбирать из соображений удобства наблюдения магнитного потока.
Создание модели интегрированного магнитного компонента
Теперь, когда мы определились с принципами моделирования неразветвленных, разветвленных и сложных магнитных цепей, попробуем создать модель трансформатора с магнитопроводом, имеющим сложную структуру.
В настоящее время в источниках вторичного электропитания все чаще используются интегрированные электромагнитные компоненты, совмещающие в себе разделительный трансформатор и дроссель фильтра. Также известны удачные попытки моделирования подобных компонентов средствами 81МРЬОКЕИ [5].
Для примера, создадим модель сложного электромагнитного компонента, изготовленного на ферритовом магнитопроводе Е 42/ 21/20 (рис. 18а). Для удобства керны пронумерованы, крайним присвоены номера 1 и2, а среднему — номер 3. Керн № 3 имеет немаг-
нитный зазор длиной Lg = 0,0005 м. На кернах № 1 и 2 расположено по две обмотки, одна из которых имеет 20 витков, а другая — 6 витков (рис. 18б). Обмотки 1а и 1б имеют сопротивление по 0,02 Ом.
Ранее мы уже моделировали подобную структуру магнитопровода. Нарис. 15 изображена соответствующая схема для переменного тока. На рис. 19 изображена схема модели переменного тока для рассматриваемого примера. Здесь индуктивности L1, L2, L3 соответственно имитируют участки 1, 2, 3 магнитной цепи, изображенной на рис. 18б. Кроме этого, для гальванической развязки обмоток и получения требуемых намоточных данных в модели использовано два идеальных трансформатора К1 и К2.
Так как часть интегрированного компонента используется в качестве дросселя постоянного тока, то необходимо дополнить модель соответствующими нелинейными сопротивлениями, которые сделают ее более реалистичной и достоверной.
Ранее мы определили таблицу настройки нелинейного сопротивления для случая, когда магнитопровод имеет длину 1 м, сечение 1 м2 и на него намотана обмотка в 1 виток:
1аЫ1е=(-0.465 -10000,-0.45 -2500,-0.4 -320,-0.3 -85,-0.2 -35,-0.1 -15,0 0,0.1 15,0.2 35,0.3 85,0.4 320,0.45 2500,0.465 10000)
Скорректируем таблицы нелинейных резисторов с помощью формул (3) и (4). Для кернов 1 и 2:
1аЫ1е=(-0.011625 -34.5,-0.01125 -8.625,-0.01 -1.104,-0.0075 -0.29325, -0.005 -0.12075,-0.0025 -0.05175,0 0,0.0025 0.05175,0.005 0.12075, 0.0075 0.29325,0.01 1.104,0.01125 8.625,0.011625 34.5)
Для керна 3:
1аЫ1е=(-0.022692 -27.2535,-0.02196 -13.455,-0.01952 -8.536,-0.01464 -6.123,-0.00976 -4.043,-0.00488 -2.017,0 0,0.00488 2.017,0.00976 4.043,0.01464 6.123,0.01952 8.536,0.02196 13.455,0.022692 27.2535)
При подгонке величины нелинейных сопротивлений использовался поправочный коэффициент К = 200. Полная модель инте-
Hc=12 Bs=0,45 Br=0,1 A=0,000244 Lm=0,036 Lg=0,0005 N=20 L3
Hc=12 Bs=0,45 Br=0,1 A=0,000125 Lm=0,0069 Lg=0 N=20 Hc=12 Bs=0,45 Br=0,1 A=0,000125 Lm=0,0069 Lg=0 N=20
Hc=12 Bs=0,45 Br=0,1 A=0,000244 Lm=0,036 Lg=0,0005 N=20
Hc=12 Bs=0,45 Br=0,1 A=0,000125 Lm=0,0069 Lg=0 N=20 Hc=12 Bs=0,45 Br=0,1 A=0,000125 Lm=0,0069 Lg=0 N=20
Рис. 19. Модель интегрированного электромагнитного компонента для переменного тока
Рис. 20. Полная модель интегрированного электромагнитного компонента
И
Рис. 21. Графики напряжения, тока и индукции в характерных электрических и магнитных цепях интегрированного электромагнитного компонента
грированного электромагнитного компонента изображена на рис. 20. Отображение таблиц нелинейных сопротивлений G1, G2, G3 отключено с тем, чтобы не загромождать схему лишней информацией.
Как уже говорилось ранее, для контроля индукции в определенном участке магнито-провода удобно контролировать падение напряжения на соответствующем нелинейном резисторе. Крайние значения напряжения управления в настроечных таблицах соответствуют максимальной учтенной индукции в сердечнике. Например, для G1 и G2 падение напряжения, равное 0,011625 В, будет соответствовать индукции 0,465 Тл. Следовательно, чтобы наблюдать сигнал, пропорциональный индукции, достаточно умножить соответствующий график на 0,465/0,011625 = 40.
Соответственно, для G3 максимальное падение напряжения 0,022692, а коэффициент пересчета графика 0,465/0,022692 = 20,5.
Теперь, когда модель интегрированного электромагнитного компонента готова, осталось ее проверить. Для этого создадим простейшую тестовую схему (рис. 21). Напряжение возбуждения для интегрированного компонента формируется с помощью источников V! и V2. Далее это напряжение трансформируется, выпрямляется при помощи диодов D1 и D2, сглаживается и поступает в нагрузку Я1. На рис. 21а изображены графики напряжения возбуждения интегрированного компонента V(n001) и тока в резисторе нагрузки ЦЯ1). На рис. 21б изображены графики индукции (мТл) в крайних (В1 и Б2) и среднем (Б3) кернах интегрированного компонента. ■
Литература
1. IEEE Transactions On Computer-Aided Design, Vol. 10. No. 4, April 1991 (ссылка из HELP-файла программы LTspice/SwCAD III).
2. Бессонов Л. А. Теоретические основы электротехники. Изд. 6-е, перераб. и доп. М.: Высшая школа, 1973.
3. Sandler S. M. SMPS Simulation with Spice 3.
4. Миловзоров В. П. Электромагнитная техника. М.: Высшая школа, 1966.
5. Njiende H., Wetzel H., Froehleke N., Cronje W. A. Models of integrated magnetic components for simulation based design of SMPS with SIMPLORER. http://wwwlea.uni-paderborn.de/fileadmin/ Elektrotechnik/AG-LEA/forschung/ veroeffentlichungen/2003/ 03epe-njiende-wetzel-froehleke-cronje.pdf.
