Автоматизированное проектирование высокодобротной колебательной системы транзисторного генератора. Численные методы расчета электромагнитного поля
В коротковолновой части сантиметрового диапазона длин волн уровень собственных шумов генератора растет из-за падения добротности колебательной системы ТрАГ. Высокую стабильность частоты генератора можно достичь применением фазовой автоподстройки частоты (ФАПЧ). Однако для реализации такого технического решения требуются дорогостоящие высокостабильные микросхемы. Кроме того, ФАПЧ – крупногабаритная система, в то время как ТрАГ без такой системы выгодно отличаются миниатюрностью и дешевизной. Так, ТрАГ на основе ДР [3] на частоту 16 ГГц размером 2,5x2,5x0,8 см обеспечивает выходную мощность не менее 10 дБм при токе потребления 30 мА и напряжении питания 3 В (рис.1). Однако в отсутствие ФАПЧ при отстройке от несущей на 10 кГц уровень частотных шумов такого ТрАГ составляет 80 дБ/Ц––Гц , что для большинства применений явно недостаточно.
ПОСТАНОВКА ЗАДАЧИ
Недостаточная кратковременная стабильность частоты ТрАГ в коротковолновой части сантиметрового диапазона длин волн связана с уменьшением добротности традиционных колебательных систем (КС) генераторов вследствие роста потерь в диэлектрике при повышении частоты, а также ухудшением добротности микрополосковых резонаторов (МПР), входящих в состав КС. Необходимая добротность колебательной системы ТрАГ достигается за счет высокой собственной добротности ДР (в 3-см диапазоне длин волн равной ~5000 [2]) при оптимальном включении ДР в систему отрезков микрополосковых линий (МПЛ), формирующих МПР. Значение добротности МПР (определяемой потерями в проводниках и на излучение) не превышает 100–200 [4] и в первом приближении не зависит от свойств материала подложки. Поскольку из-за частотной дисперсии в коротковолновой части сантиметрового диапазона значения тангенса угла потерь современной СВЧ-керамики с высокой относительной диэлектрической проницаемостью сравнительно высокие (например, при er = 41 на частоте 15 ГГц tgd » 3·10-4 [2]), добротность колебательной системы с ее применением снижается.
Возможный путь повышения кратковременной стабильности без использования ФАПЧ – замена микрополосковых резонаторов более добротными полосковыми резонаторами (ПР) на линиях с подвешенной подложкой. Добротность ПР достигает 2000, что значительно выше чем у МПР [4]. Однако проектирование подобных колебательных систем затруднено из-за сложного спектра собственных частот системы, не допускающей применения приближенных расчетных оценок. Достоверный расчет систем с диэлектрическими и полосковыми резонаторами возможен с применением численных методов расчета электромагнитных полей. Это и обуславливает актуальность разработки метода автоматизированного проектирования колебательных систем высокостабильных ТрАГ с применением современного программного обеспечения, реализующего численные методы расчета электромагнитных полей.
Без снижения общности ограничим рассмотрение метода проектирования примером расчета колебательной системы ТрАГ на основе последовательной обратной связи с варакторной перестройкой частоты (рис.2). В такой системе ДР электромагнитно связан с размещенными на подложке отрезками МПЛ 1и 2. МПЛ 2 соединена с затвором полевого СВЧ-транзистора, а линия 3 служит для вывода мощности. Отрезок МПЛ 1 соединен с варикапом. На противоположных концах МПЛ 1 и 2 разомкнуты. Перестройка частоты генератора осуществляется с помощью варикапа, создающего дополнительную удлиняющую емкость и тем самым смещающего резонансные частоты полосковой линии.
Для повышения добротности в такой колебательной системе ТрАГ'а целесообразно применять несимметричные линии с подвешенной подложкой (рис.3), поскольку они более технологичны, чем симметричные, и их проще стыковать с микрополосковыми линиями в трактах передачи сигнала. Образец такого генератора, разработанного в СССР в середине 80-х годов прошлого века [5], обеспечивал электронную перестройку частоты на 40 МГц относительно несущей частоты ~15 ГГц при выходной мощности 150 мВт. Его долговременная относительная нестабильность несущей частоты Df/f составляла ~ 2·10-4, а спектральная частота – -90 дБ при отстройке от несущей на 10 кГц. Габаритные размеры ТрАГ с вентильной развязкой – 50x30x15 мм. На основе известных размеров топологии этого ТрАГ и примененного в нем ДР можно экспериментально проверить достоверность метода проектирования.
ОБОСНОВАНИЕ ВЫБОРА ПРОГРАММНОГО ПРОДУКТА
Как отмечалось, расчет колебательной системы, содержащей ДР и ПР на линиях с подвешенной подложкой, требует применения численных методов расчета электромагнитных полей для анализа спектра собственных частот. При выборе численного метода и его программной реализации следует руководствоваться следующими критериями:
· достоверностью анализа электромагнитного поля диэлектрического резонатора;
· достоверностью анализа поля полосковой линии;
· возможностью расчета собственных частот электродинамической системы;
· возможностью оценки добротности системы;
· возможностью управления точностью/трудоемкостью решения;
· возможностью автоматизации многократного анализа.
Сегодня на рынке предлагается множество коммерческих программ, реализующих различные численные методы. Наиболее популярны программные продукты, основанные на методе моментов (Method of Moments, MoM) и методе конечных разностей во временной области (Finite Differences Time Domain Method, FDTDM) (табл.1). Достоинства таких программ – дружественный интерфейс, простота работы, сравнительно малые аппаратные требования. Менее популярны "тяжелые" конечно-элементные пакеты (Finite Elements Method, FEM), для которых требуются значительно большие вычислительные мощности. К тому же они более сложны в освоении, хотя и обладают широкими возможностями.
Всем перечисленным требованиям в наибольшей мере удовлетворяет программа ANSYS, реализующая метод конечных элементов. В ней предусмотрены возможности управления густотой конечно-элементной сетки, автоматизации многократного анализа, вычисления собственных частот и добротностей электродинамической системы, что подтверждает опыт ее применения авторами при расчете полей ДР. Корректность расчета поля ДР подтверждена экспериментальными исследованиями [11].
АЛГОРИТМ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ
Исходные данные для численного моделирования колебательной системы автогенератора – требуемая частота (в предположении ее соответствия несущей частоте генератора, т.е. без учета вносимых реактивностей со стороны активного прибора, который может и должен быть выполнен на последующем этапе проектирования собственно автогенератора) в сочетании с максимальной добротностью. В случае синтеза перестраиваемой колебательной системы добавляется еще и требование обеспечения необходимого диапазона перестройки частоты. Несмотря на достоинства программы ANSYS, при непосредственном расчете электродинамической системы она не позволяет получить ясного результата и автоматизировать вычисления. Эффективный анализ и синтез колебательной системы высокостабильного транзисторного автогенератора с помощью программы ANSYS можно провести только с применением специального алгоритма. Этот алгоритм состоит из следующих этапов:
1) создание и расчет исходной модели:
· создание геометрической модели с некоторыми (начальными) конструктивными параметрами;
· создание конечно-элементной сетки для модели, проверка ее качества;
· исследование численной сходимости при определении спектра собственных частот – многократное вычисление первых нескольких резонансных частот при различных настройках конечно-элементного разбиения, выбор оптимальных по соотношению точность/трудоемкость настроек;
2) анализ спектра собственных частот исходной модели:
· вычисление (с выбранными настройками) спектра собственных частот колебательной системы;
· вычисление добротности для каждой собственной частоты;
· анализ картины электромагнитного поля для каждой собственной частоты;
· составление по результатам расчета значений добротности и анализа картин поля вывода о типе каждого резонанса и о его влиянии на технические характеристики ТрАГ;
3) автоматизированное исследование колебательной системы:
· составление программы на внутреннем алгоритмическом языке программирования ANSYS (с использованием вспомогательных файлов, созданных программой для исходной модели), позволяющей автоматически генерировать геометрическую модель, конечно-элементную сетку и граничные условия;
· тестовый прогон созданной программы – вычисление спектра собственных частот исходной модели;
· вычисление с помощью разработанной программы требуемых характеристик системы (например, зависимости значений собственных частот системы от параметров материала диэлектрического резонатора);
· расчет АЧХ системы при нагруженных полосковых линиях, вычисление коэффициентов связи между ПР и ДР;
4) анализ полученных данных, выбор оптимальных настроек.
При реализации алгоритма целесообразно руководствоваться следующими рекомендациями:
· ограничить число собственных частот, анализируемых при проектировании колебательной системы и не рассматривать частоты, резко отличающиеся (в основном вверх по частоте) от предполагаемой несущей частоты генератора;
· при вычислении добротности рассматривать в качестве единственного источника потерь токи в проводниках, поскольку программа ANSYS не позволяет вычислять собственные частоты электродинамической системы, содержащей диэлектрики с потерями. Такое допущение возможно, так как добротность колебательной системы определяется прежде всего добротностью полосковых линий, потери в которых обусловлены относительно высокой плотностью токов на краях проводников;
· делать вывод о природе резонанса для каждой из найденных собственных частот системы, сравнивая добротности при резонансах на данной частоте и на остальных собственных частотах и анализируя структуру поля, соответствующую конкретной частоте. Резонансам, контролируемым ДР, соответствуют большие значения добротности, а резонансам, определяемым только полосковыми линиями, – меньшие ее значения;
· для рассматриваемой колебательной системы наблюдаются парные (для четного и нечетного типа возбуждения) резонансы полосковых линий, их расстройка по частоте определяется значением коэффициента связи между полосковыми резонаторами;
· моделировать перестройку частоты, осуществляемую варикапом, меняя длину соответствующего полоскового резонатора.
С учетом сказанного блок-схему системы автоматизированного проектирования можно представить в виде, показанном на рис.4.
Предлагаемый алгоритм, основанный на программе ANSYS, поддерживается программным обеспечением CADGen, работающим под управлением Windows 95 и позволяющим оперативно проектировать колебательную систему ТрАГ. Для вычисления s-параметров системы при нагруженных полосковых линиях предусмотрена возможность присоединения 50-Ом коаксиальных портов. Разработанная база знаний позволяет автоматически создавать высококачественное конечно-элементное разбиение, задавать граничные условия, настраивать программу ANSYS. Номограммы и табличные данные, составляющие базу данных, обеспечивают оперативное получение начальной топологии и ее оптимизацию. Заметим, что с помощью предлагаемого алгоритма возможен и синтез колебательных систем на основе микрополосковых резонаторов.
ПРИМЕР ПРАКТИЧЕСКОГО ПРИМЕНЕНИЯ АЛГОРИТМА ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ
Рассмотрим этапы проектирования колебательной системы ТрАГ на частоту 14,7 ГГц с помощью программы CADGen и комплекса ANSYS (рис.5). При численном моделировании в качестве нулевого приближения приняты следующие электрофизические параметры материалов и исходные конструктивные размеры:
· er подложки ПР равно 9, материала резонатора – 35, диэлектрической подставки – 4; s металлического покрытия проводников – 5,8·107 См/м, mr = 1;
· толщина подложки – 0,5 мм; ширина выступа, образующего линию с подвешенной подложкой, – 1,75 мм; длина и ширина полосков – 2 см и 0,75 мм, соответственно; диаметр ДР – 3,9 мм, высота – 2 мм; высота диэлектрической подставки – 0,5 мм; зазоры между ДР и полосками – 0,3 мм; размер модели (с учетом воздушного заполнения) 10x20x5,5 мм.
Конечно-элементное разбиение модели создавалось с учетом малого зазора между ДР и ПР. В диэлектрическом резонаторе и в зазоре между ним и полосками сетка имеет наибольшую плотность размещения конечных элементов. Поскольку топология геометрической модели иррегулярна, для построения сетки применялись тетраэдральные конечные элементы и устанавливался нерегулярный ("свободный" – free) способ разбиения. Изучались различные варианты разбиения. Оптимальным по соотношению точности и трудоемкости оказалось разбиение, представленное на рис.5.
Поскольку анализ проводился для колебательной системы ТрАГ'а с рабочей частотой ~14,7 ГГц, целесообразно ограничиться рассмотрением первых десяти собственных колебаний (табл.2).
Выводы о природе резонансов, обнаруженных в системе, можно сделать на основе анализа значений добротности и картин поля, соответствующих нескольким характерным резонансам системы, представленным на рис.6–9 в виде изоповерхностей напряженности электрического поля.
С помощью программы CADGen исследовалось влияние различных конструктивных факторов на значения собственной частоты в системе (табл. 3).
При оценке чувствительности топологии колебательной системы к точности исполнения, а также для построения электронной (варакторной) перестройки или подстройке частоты ТрАГ можно использовать зависимость резонансных частот колебательной системы от длины ПР (рис.10).
Создание эффективной колебательной системы ТрАГ предполагает отбор наиболее добротных резонансов. Для рассматриваемой системы – это собственные колебания под номерами 6, 7 и 8 (см. табл.2). А поскольку заданная для проектирования ТрАГ рабочая частота составляет 14,7 ГГц, то реализовывать колебательную систему целесообразно на основании собственного колебания № 7. Следует отметить, что результаты расчета привели к синтезу колебательной системы ТрАГ, идентичной системе, примененной в созданном ранее ТрАГ с варакторной перестройкой частоты [5].
Таким образом, разработана методика автоматизированного синтеза колебательной системы высокодобротного ТрАГ на основе программы ANSYS, реализующей метод конечных элементов. Предложенный алгоритм позволяет проектировать колебательные системы, включающие в себя диэлектрические резонаторы и резонаторы на отрезках линий с подвешенной подложкой. Возможен синтез систем, содержащих резонаторы иного типа. Полученные с помощью численных расчетов характеристики колебательных систем хорошо согласуются с данными для использующихся на практике генераторов. Следовательно, можно говорить о решении задачи оперативного проектирования колебательных систем для миниатюрных высокодобротных транзисторных автогенераторов.
Литература
1. Абраменков, И.И. Бродуленко, В.М. Геворкян, Д.А. Ковтунов А.И. Состояние и перспективы применения миниатюрных диэлектрических резонаторов в перестраиваемых полупроводниковых генераторах. – Обзоры по электронной технике. Сер. Электроника СВЧ/Под ред. В.М. Геворкяна – М.: ЦНИИ "Электроника", 1988, вып. 5(1359). – 70 с.
2. Высокочастотные керамические материалы и микроволновые элементы: Каталог. – СПб.: ООО "Керамика", 2000.
3. New Products. Phase Locked Oscillators for LMDS, VSAT and POINT TO POINT Wireless Applications. – Microwave-Magazine, 2000, Dec., v.1, N4, p.104–105.
4. Малорацкий Л.Г. Микроминиатюризация элементов и устройств СВЧ. – М.: Сов. радио, 1976.
5. Геворкян В.М., Ковтунов Д.А. Генератор СВЧ. Авт. св. № 1788664, 05.07.90 г.; МКИ 5 Н03В 5/18: – Бюл. № 2 от 15.01. 1993.
6. ANSYS Electromagnetic field analysis guide. – ANSYS Inc, 1998.
7. Microwave Office – EMSight Reference. – Applied Wave Research, 2000.
8. IE3D User's manual. – Zeland Software, 2000.
9. Fidelity User's manual// Zeland Software, 1999.
10. New features in CONCERTO V3. – Vector, 2002, v.17, N1, p.6.
11. Казанцев Ю.А., Бунин А.В., Геворкян В.М., Вишняков С.В., Капустин Ю.Г., Полукаров В.И. СВЧ фильтры на диэлектрических резонаторах. – Конференция "Излучение и рассеяние электромагнитных волн", 2001. – Таганрог: изд-во ТГРУ, с.191–193.