Доклады БГУИР
2013 № 1 (71)
УДК 621.396.677
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ХАРАКТЕРИСТИК ИЗЛУЧЕНИЯ ОДНОЗЕРКАЛЬНОЙ АНТЕННЫ
М.М. ДГАЛИ
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 28 ноября 2012
Приводятся результаты численного моделирования характеристик рассеяния однозеркаль-ной антенны в коммерческой программе FEKO и в специализированной программе. Результаты сравниваются по различным критериям. Делается вывод о том, что при практически одинаковых результатах программа FEKO требует значительно больших компьютерных ресурсов. При сложном облучателе, например в моноимпульсной антенне, моделирование в программе FEKO становится невозможным.
Ключевые слова: зеркальная антенна, диаграмма направленности, коэффициент направленного действия, метод физической оптики, метод моментов.
Введение
При проектировании зеркальных антенн необходимым этапом является численное моделирование, которое проводится с целью определения оптимальных геометрических параметров, при которых удовлетворяются заданные требования к электрическим характеристикам. Эта задача может быть решена с помощью известных коммерческих программ общего назначения, например CST Microwave Studio или FEKO. Если ориентироваться на рекламные материалы, можно сделать вывод о том, что программа FEKO использует максимальное число методов вычислительной электродинамики и, следовательно, обеспечивает наиболее точное моделирование электродинамических структур. Опыт работы с этой программой по моделированию различных типов антенн показывает, что задача моделирования зеркальных антенн требует применения довольно мощных компьютеров. И даже в этом случае время решения задачи при волновом диаметре рефлектора порядка 40-50 составляет несколько часов. Моделирование же, например, моноимпульсной зеркальной антенны с четырехмодовым облучателем на компьютере с 4-х ядерным процессором, тактовой частотой более 3 ГГц и оперативной памятью 32 ГБ требует более суток времени. Это становится неприемлемым при оптимизации антенны методом перебора параметров геометрии. В связи с этим представляет практический интерес вопрос о возможности применения для решения этой задачи специализированной программы, предназначенной для моделирования только зеркальных антенн. В данной статье описываются результаты сравнительного анализа решения задачи моделирования зеркальной антенны с использованием программы FEKO и специализированной программы Ref_ZZ [1]. В программе Ref_ZZ используется комбинация метода физической оптики при определении поля рассеяния рефлектора и метода интеграла Кирхгофа при определении поля рупора облучателя на поверхности рефлектора.
Геометрия задачи
В обеих программах анализируется однозеркальная антенна с рефлектором двух типов - с параболоидом полного профиля и с параболоидом типа «Оффсет». Геометрия задачи поясняется рис. 1, где 1 - параболоид вращения полного профиля; 2 - внеосевая вырезка из па-
раболоида вращения (рефлектор типа «Оффсет» - далее «вырезка»); Д - диаметра параболоида полного профиля; ¥ - фокус параболоида; ¥р - фокусное расстояние параболоида. Д, Д -диаметры эллипса, являющегося контуром раскрыва вырезки; 7ш;п - расстояние от нижней кромки вырезки до фокальной оси параболоида ОТ; ¥1 - точка положения фазового центра рупора облучателя; ДгУ, Дг2 - смещение точки ¥ из фокуса в плоскости У0Ъ (используется для сканирования и фокусировки антенны в ближнюю зону).
Рис. 1. Зеркальная антенна с обозначениями При =0,5 Д внеосевая вырезка превращается в симметричную вырезку относительно фокальной оси с размерами вдоль осей X и У- Дд. и 1)у = 1)р
В качестве облучателя рефлектора использован пирамидальный рупор с волной ТЕ10. Размеры раскрыва рупора: Аг - в плоскости Н; В. - в плоскости Е. Глубина рупора - Яг . Плоскость Е рупора совпадает с плоскостью У ОТ, (рис.1). Плоскость поляризации антенны - плоскость УОЪ (ф = 90 ). В этой плоскости основная составляющая поля - ¥в, составляющая по ортогональной поляризации - ДФ . В плоскости Н (ср = 0) основная составляющая поля - ¥ц,, составляющая по ортогональной поляризации - ¥в.
При моделировании поля рупора на поверхности рефлектора учитывалось нелинейное фазовое распределение на его раскрыве и то, что рефлектор расположен в ближней зоне рупора.
Далее для иллюстрации степени взаимного соответствия результатов, полученных в двух программах, приводятся результаты моделирования зеркальных антенн на частоте 10 ГГц при ¥р = 0.51)^ (как это часто используется). Размеры рупора облучателя выбраны так, чтобы
на краях поверхности рефлектора поле рупора составляло 0,3 от максимума в центре рефлектора. Это примерно соответствует максимуму коэффициента направленного действия (КНД) зеркальной антенны [2, 3].
Результаты численного моделирования
Результаты расчета основных параметров диаграммы направленности (ДН): ширины главного лепестка (2-О 0,5 ) и максимального бокового лепестка (¥Ьт) при двух значениях волнового диаметра рефлектора (Д / X ) для антенны с параболоидом полного профиля представлены на графиках рис. 2 - 3. Результаты показаны в виде зависимостей от ОгУ/Х, где ОгУ -смещение облучателя вдоль оси У (используется для сканирования), X - длина волны.
-4 -б -8 -10
I -20 ^ -22 -24 -26 -28
-1
О
1 2 3 ВгУ/Ъ.
5
^о5 4
3«
Сй ,
Рис.2. Зависимость параметров ДН от смещения облучателя из фокуса: 1.)р / /. =30
о
' -КЕКО '-Ref_ZZ
........;.......г
-2
2 4 ПгУ/Х-
Рис.3. Зависимость параметров ДН от смещения облучателя из фокуса: Ор / к =60
Из рис. 2 и рис. 3, а также из других результатов, не приведенных в статье, следует:
- программы БЕКО и дают одинаковое значение отклонения главного лепестка ДН от фокальной оси при сканировании (изменении величины ВтТ);
- программа дает заниженное значение уровня максимально бокового лепестка, но различие в значениях этого параметра, полученых в двух программах, уменьшается при увеличении волнового диаметра рефлектора и при увеличении угла сканирования. Эти расхождения объясняются тем, что метод физической оптики, использованный в программе
не учитывает краевых эффектов;
- вследствие указанных различий в уровне боковых лепестков программа дает завышенное значение КНД (примерно на 1 дБ).
Рис. 4 иллюстрирует степень различия в ДН, рассчитанных в двух программах.
Принципиальным вопросом сравнения двух программ является сравнение требуемых для работы компьютерных ресурсов. На графиках рис.5 приведены результаты сравнения по времени решения задачи и объему требуемой оперативной памяти.
На графике рис.5, а обозначено отношение Т/1Тге/, где Т/- время решения задачи в программе БЕКО, Тге/- время решения задачи в программе Как видно, программа требует на порядки меньше компьютерных ресурсов.
Аналогичные результаты сравнения получены при моделировании зеркальной антенны типа «Оффсет».
w
ч
-20 -30 -40 -50 -60
Л » i / \ Л •••/• VA 1 li 1
А 11 s
А А /.....j ArJ д У
\ ! £ f ЛА/V, J ■ \ 1 ï il \N\ M U Л V \ /........... ..y
1 У У ï s
-60 -40 -20 0 20 40 Угол наблюдения [град.]
60
80
Thêta | б
Рис. 4. Диаграммы направленности, рассчитанные в программах: а - Ref ZZ; б - FEKO;
Д, / X =60; DrY! X =5,85
а б
Рис.5. Зависимости от волнового диаметра рефлектора (Др/Х):
а - отношения времени счета в двух программах; б - объема оперативной памяти компьютера в двух программах
а
Заключение
Проведенный анализ позволяет сформулировать рекомендации при решении задачи численного моделирования зеркальной антенны. При волновом диаметре рефлектора Д^/Х<20 для моделирования целесообразно использовать программу БЕКО, при большем диаметре целесообразно использовать специализированную программу, построенную с применением метода физической оптики.
NUMERICAL SIMULATION OF SINGLE-DISH ANTENNA RADIATION CHARACTERISTICS
MM. DGHALI Abstract
The numerical simulation results of scattering characteristics of a single-dish antenna obtained by commercial Program FEKO and in the specialized program are given. The results are compared on different criteria. The Conclusion is that in almost identical results FEKO program requires much more computing resources. If the feeder is complex, such as monopulse antenna, has simulation using program FEKO becomes impossible.
Список литературы
1. Юрцев О.А. // Матер. VII межд. НТК «Wyksztalcenie I nauka bez granic-2011». Польша, 7-15 декабря 2011 г. Т. 29. С. 31-37.
2. Шифрин Я.С. Антенны. Харьков, 1978.
3. Stutzman Warren L., Thiele Gary A. Antenna theory and design: third edition. Danvers, 2012.