Форум » Моделирование устройств силовой электроники » Реализация модели Джайлса-Атертона в Microcap 9 » Ответить

Реализация модели Джайлса-Атертона в Microcap 9

Finarfin: К сожалению не могу достать оригинальную статью авторов, но во всех публикациях (для определенности возьмем эту) по теме модель имеет пять параметров, в то время как в microcap их только 4. Обычно во всех статьях: 1) Ms - намагниченность насыщения 2) a - масштабный коэффициент (растягивание-сжатие вдоль оси абсцисс) Ms и a - параметры безгистерезисной кривой 3) k - коэффициент безвозвратной деформации (pinning) стенок, определяет площадь петли гистерезиса (при k=0 гистерезиса нет) 4) alpha - характеризует взаимодействие (coupling, interaction) между доменами 5) c - отношение дифференциальной проницаемости кривой начального намагничивания в начале координат к дифференциальной проницаемости при рассмотрении только безгистерезисной составляющей (в других источниках - коэффициент обратимости перемещаний стенок доменов) Всего пять коэффициентов, подбираемых алгоритмом оптимизации до наилучшего соответствия экспериментальным точкам (пример реализации) В Microcap же видим только Ms, a, с и k. Но где же alpha??? Она полагается равной какому-то постоянному значению? Каково оно? Боюсь что это существенно снижает точность моделирования. По хорошему даже параметр k полагается переменным и переходят от пяти параметров к семи-восьми.

Ответов - 216, стр: 1 2 3 4 5 6 7 8 9 10 11 All

qaki: Кажется нашел. Ток намагничивания первички и, само собой, ее индуктивная реакция считается через DDT(i). Но нужно правильно задавать направление ЭДС самоиндукции.

neon: правильная ссылка: http://www.eecs.berkeley.edu/Pubs/TechRpts/1975/ERL-520.pdf

qaki: Таки прорезалась модель Чана на Микрокапе. Правда пока на костылях и с подпорками, но первый результат выдала.[url=]click here[/url] Это предельный и частный циклы перемагничивания синусоидальным током феррита 2000НМ1 в интертрепации модели Чана. Удобств в модели пока никаких. Но кое-что может быть удастся сделать.

locik: вот что получилось пока у меня. http://zalil.ru/34087763

qaki: locik пишет: вот что получилось пока у меня. Посмотреть не получается. У меня 9.0.7 из Интернета, а у Вас видимо 10. В моей проге даже шейпов таких нет. За 10 охотился, но закончилось тем, что мой мобильник обули на 900 р. и предлагали доплатить еще по 600 в месяц за какие-то непонятные услуги. Еле отбился через службу безопасности. У меня все гораздо проще и недоделаней. Пока только рисуются кривые индукции при синусоидальном намагничивании. Все остальное вручную через задание параметров и убирание ненужных хвостов. Для перехода на намагничивание от генератора напряжения произольной формы нужно, чтобы Микрокап решал нелинейное дифференциальное уравнение. Пока не пробовал. Подобрал всего лишь затычку попроще для имитации эффекта изменения МЮ при намагничивании. Может быть и потянет.

Aml: qaki, см. https://sites.google.com/site/kafeimt/ucebnye-kursy/kaseu/micro-cap-9 Только первый абзац не читай :)

qaki: Aml пишет: qaki, см. Большое спасибо, но дружбы с Google-диском у меня не получается. Просит пароль, которого у меня нет.

Aml: Большое спасибо, но дружбы с Google-диском у меня не получается. Просит пароль, которого у меня нет. Хром, наверное? Черт бы подрал этот Хром... Не нужен там никакой пароль (на самом деле). И никакой броузер (кроме Хрома) его не требует. Просто поражает, что единственный броузер, который глючно работает с диском Google - это собственный броузер Google. Это одна из причин, по которой я его категорически не использую. В общем, под Хромом надо щелкать по ссылке правой клавишей и там в контекстном меню выбирать что-то типа "входить в режиме инкогнито". (точную формулировку не помню) Во всех нормальных броузерах - просто перейти по ссылке.

qaki: to locik Собственно что сделано. Удалось отвязаться от жесткой связи между индуктивным компонентом Микрокапа и его вычислительным процессом. Если мы задаем L=1, то в штатном режиме МС интерпретирует это как 1 Гн. Если мы вызовем компонент К, где эта единица интерпретируется как один виток, то мы автоматом попадаем в модель J&A. Чтобы этого не происходило, в простенькой модели[url=Чан.cir]click here[/url] индуктивная реакция создается с помощью функционального генератора NFV,создающего встречную ЭДС самоиндукции. При этом величина L задается уже через витки и геометрию сердечника, а величина ЭДС равна L*dI/dT. Величина Mu пока константа. Для намагничивания синусоидальным током это не имеет значения. По этой причине источник намагничивания представлен в виде генератора напряжения и последовательного резистора большой величины. Модель рисует верхнюю и нижнюю ветви гистерезисной кривой, их полусумму, которая у Чана играет роль кривой начального намагничивания. Задавая сдвиги по вертикали, можно получать частные циклы. Однако при намагничивании от генератора напряжения, что характерно для силовой электроники, необходимо учитывать эффект снижения Mu при заходе в область насыщения. Для начала хочу добавить простенькую корректирующую кривую. Это один из семидести вариантов определения Mu (по Аркадьеву), где Mu задается отношением В/Н. Просто опасаюсь, что Микрокап встанет раком, если ему еще и подсунуть dB/dH модели Чана. to Aml Спасибо. Попробую пробиться

locik: у меня 9.0.7.0.откуда скачал уже не помню.в некоторых архивах дается ссылка на источник.всегда это удалял.сейчас понял что зря.может причина в отсутствии моделей в схеме.т.е. у меня то эти модели есть в библиотечных файлах CoreCoil и MagCore.это надо спросить у Aml.

qaki: to Aml Еще раз большое спаcибо. Скачал. Начинаю upgrade.

Aml: Если честно, я не в восторге от МС10 как раз в плане использования индуктивных компонентов. Он почему-то обладает свойством "вылетать" даже при простом добавлении индуктивности в схеме. Плюс еще что-то с индуктивностью Марина находила (достаточно серьезное), надо разработчикам писать.

qaki: qaki пишет: Чтобы этого не происходило, в простенькой модели[url=Чан.cir]click here[/url] индуктивная реакция создается с помощью функционального генератора NFV,создающего встречную ЭДС самоиндукции. Какая-то ерунда с этой ссылкой. На Яндекс-народ в "Моих файлах" она есть, а при вызове попадаем на несуществующую страницу. Сейчас перезалью. Чан.cir Фантастика какая-то! Нет ни той, ни другой. Минуту назад вроде бы были. Яндекс видно дурит, хотя и вспышек на Солнце нет http://www.tesis.lebedev.ru/sun_flares.html

locik: чан.cir скачал правда не с 1 раза.со 2 получается.запустил.пока не понял смысла.

qaki: locik пишет: пока не понял смысла. Как-то сразу не соображу, что ответить на вопрос: "Куда лошадь запрягать?" Если есть необходимость исследовать цепь, содержащую ферромагнитный сердечник, не прибегая при этом к модели J&A, то напрямую Микрокап такой возможности не дает. Встроенная модель сердечника жестко привязана к модели J&A со всеми ее достоинствами и недостатками. Для того, чтобы использовать иные модели необходимо организовать вход в вычислитель Микрокапа и реализовать в нем соответствующие математические процедуры. Мне таким удобным шлюзом представляется функциональный источник NFV. Он позволяет формировать напряжение в соответствие с математическим выражением, используемым как его атрибут. Включение индуктивности в любую электрическую цепь может быть представлено как включение эквивалентного источника ЭДС, создающего напряжение равное ЭДС самоиндукции включенной индуктивности. Такая замена сразу открывает возможность задания индуктивности в моделируемую схему как в виде ее номинала в единицах Гн, так и виде ее конструктивных параметров, т.е. числа витков, площади витка, длины магнитной линии и магнитных параметров среды. В этом как раз и заключается смысл представленного файла. Почему такое внимание к индуктивности? Да очень просто. Любой трансформатор и даже дроссель имеет первичную обмотку, т.е. индуктивность, создающую магнитный поток в сердечнике. Если наша модель позволяет вычислить этот поток, то напряжение вторички легко вычисляется в соответствие с законом Фарадея. Тем или иным способом можно учесть связь между обмотками, т.е выйти на полную схему замещения для моделирования полной схемы девайса. При этом нашим первым действием как раз и является получение зависимости индукции в сердечнике от входного воздействия. Сам же файл Чан.cir иллюстрирует процесс перехода от входной цепи к индукции в сердечнике и получаемые при намагничивании петли гистерезиса.

qaki: В модели Чана по моим предыдущим ссылкам была ошибка в масштабном коэффициенте. Ошибку устранил и добавил в файл для сравнения модель Джилса-Атертона. Новая ссылка Чан.cir

qaki: Причесал входную часть модели.[url=Чан2.cir]click here[/url] Теперь это представлено в виде двух трансформаторов - в одном сердечник по безгистерезисной модели Чана, в другом по модели J&A. В обоих случаях сердечник имеет единичную длину магнитной линии, единичное сечение магнитопровода и по одному витку. Пришлось придумать новые наименования параметров, например Sech (сечение) вместо Area. Иначе модели конфликтуют. Теперь нужно написать блок управления смещением ветвей в модели Чана. Список опубликованных команд Микрокапа уж очень бедный. Возможно есть еще команды, о которых мало кто знает.

locik: прошу прощения за отсутствие.qaki смещение ветвей это самый интересный момент в модели чана.фактически нужно запомнить величину разности Bup-Bdn в верхней и нижней точке перегиба.и потом вычитать из Bup и Bdn.как запомнить эти величины не ясно.пока получается что нужно использовать управляемые источники.т.е. использовать как их запоминающую ячейку. но это как то коряво.уверен что есть способ проще.используя какие либо функции спайс.

qaki: Добрый день, уважаемый locik. Столкнулся следом за Вами с аналогичной трудностью. Пытался преодолеть ее через условный оператор IF, но Микрокап у меня на эту команду похоже не реагирует. Вроде бы роспись netlist делается по RM, но реакции никакой. Вместо того, чтобы перейти в точке экстремума с безгистерезисной кривой на кривую предельного цикла, рабочая точка остается на безгистерезисной, а на выходе изначально сидит кривая предельного. Натолкнулся в модели IR драйвера IR2104 на несколько иную роспись условного оператора .SUBCKT COMP 1 2 3 4 E1 5 4 VALUE={IF((V(1)>V(2)), V(4)+5V, V(4))} R1 5 3 1 C1 3 4 10P .ENDS Попробую оттестироваться на ней. Появилась мыслишка, а что, если попробовать раздуть Джилсовый сигмоид наподобие Чановской модели. Прикалибровать предельный цикл через Bs, Br и Нс, а дальше деформировать его пропорционально входному воздействию.

locik: я использовал как раз такую интерпретацию IF.она работает.вопрос в том что написать после второй запятой чтобы просто запомнить величину,например Bmax.потому что без третьей части она не работает.пишет не хватает данных.или что то подобное. например .param Bmax={if(ddt(Bdn)>=0,Bup-Bdn,что написать здесь чтобы запомнить Bmax если условие не выполняется)}. если написать Bmax то по сути получается рекурсивный вызов самой Bmax и микрокап ругается.мультисим и лтспайс просто вылетают безо всяких объяснений.наверное есть какой то оператор который означает подстановку не самой функции а ее значения. можно попробовать такую интерпретацию: отвлеченный источник Е1 100 0 value={if(ddt(Bdn)>=0,Bup-Bdn,v(100))} но не уверен что микрокап пропустит такое т.к. фактически Е1 как бы равен самому себе.дальше просто вычитаем: .param B={Bup-v(100)}.так же для другой ветки используя например Е2.у меня просто уже иссяк энтузиазм.жду может кто подскажет.может я уже просто смотрю в книгу но ничего не вижу а решение где то рядом.



полная версия страницы