Электродинамическое моделирование СВЧ-структур и антенн – важный этап проектирования сложных радиотехнических систем. Для решения этих задач разработчики используют мощные и эффективные средства САПР, в которых реализован широкий набор методов электромагнитного моделирования. Одним из таких инструментов является программный комплекс FEKO от компании Altair. FEKO – многофункциональная программная среда для численного электромагнитного моделирования, основанная на современных вычислительных технологиях. Рассмотрим, как с помощью FEKO решается задача обеспечения электромагнитной совместимости антенных систем, установленных на летательных аппаратах.
Требования к вычислительным ресурсам компьютеров, используемых для электромагнитного моделирования, возрастают по мере повышения сложности и увеличения физических размеров антенных систем. В определенный момент времени ресурсов, имеющихся для расчета электрически больших объектов, может оказаться недостаточно. Для решения этой проблемы программный комплекс FEKO предоставляет следующие возможности:
Анализ рассеяния электромагнитного поля на электрически больших объектах всегда был сложной задачей. Для достижения высокой точности расчета в FEKO обычно используются методы электродинамического анализа MoM или MLFMM. Для электрически больших структур применяют асимптотические методы, такие как метод PО. Однако применение электродинамических методов на крупных объектах ограничивается вычислительными ресурсами компьютера, в то время как метод PO может использовать очень плотную сетку разбиения, что резко увеличивает время расчета.
Для решения этих проблем в FEKO был реализован метод геометрической оптики с лучевым распространением (Ray Launching Geometrical Optics – RL-GO), который позволяет использовать возбуждение плоской волны. Метод RL-GO основан на моделировании лучей, падающих и отраженных от модели объектов, в соответствии с теорией распространения, отражения и преломления оптических лучей. Воздействие лучей на металлические и диэлектрические структуры моделируется при помощи источников Гюйгенса, помещаемых в каждую точку контакта лучей с материальными границами. Процессом распространения лучей (ray launching) легко управлять на основе углового (для локализованных источников) или поперечного (для источников плоской волны) интервала лучей и количества многократных разрешенных интерференций.
Метод RL-GO обеспечивает высокую степень совпадения результатов с методом MLFMM при меньших требованиях к вычислительным ресурсам компьютера. Пример моделирования рассеяния электромагнитного поля на летательном аппарате (рис.1) показывает, что метод RL-GO требует всего 4,6 Мбайт памяти, в то время как для метода MLFMM нужно 441 Мбайт при существенно большем времени расчета. Столь радикальное снижение требований к вычислительным ресурсам – решающий фактор при решении сложных задач.
Преимущества метода RL-GO наглядно проявляются также при анализе зеркальных антенн и линз. В этом случае оптимальный вариант – сочетание методов RL-GO и MoM. Антенны излучения, освещающие линзу или зеркало, могут быть смоделированы по методу MoM, а при работе с большой структурой очень эффективен метод RL-GO, использующий связь по дальнему или ближнему полю.
Одна из важных задач, стоящих перед проектировщиками радиоканалов летательных аппаратов (ЛА), – обеспечение электромагнитной совместимости (ЭМС) радиоканалов. Необходимость ее решения обусловлена, во-первых, большим количеством разнообразных радиоэлектронных средств (РЭС), в том числе с перекрывающимися частотными диапазонами, а во-вторых, быстро меняющимся взаиморасположением ЛА в пространстве и в группировке.
Для обеспечения электромагнитной совместимости (ЭМС) функционирование РЭС в составе комплекса (включение, выключение, переключение режимов работы) должно регламентироваться по пространственному (диаграммы направленности антенн), частотному и временному разносу. Это можно обеспечить путем создания системы управления параметрами электромагнитного излучения (частотным, мощностным и временным режимами, а также направлением излучения антенн) всех взаимовлияющих РЭС [1].
При решении данной задачи требуется определить в режиме реального времени внутрисистемную электромагнитную обстановку, включая вычисление коэффициентов электромагнитной связи попарно всех антенн комплекса РЭС (каждая с каждой). В то же время, сложность геометрии ЛА, влияющая на характеристики излучения антенн, обуславливает необходимость применения электродинамических методов расчета, требующих значительных вычислительных ресурсов [2].
На предварительном этапе в программной среде МАТLAB формируется ограниченное дискретное множество сценариев взаимного расположения ЛА в пространстве. С помощью ЦЭВМ, обладающей необходимыми вычислительными ресурсами, решается серия соответствующих задач на электродинамическом уровне. По результатам электродинамического расчета коэффициентов электромагнитных связей попарно всех антенн комплекса РЭС для множества сценариев формируется база данных (БД). После необходимой компиляции БД помещается в бортовые вычислители всех летательных аппаратов для последующего использования данных в режиме реального времени.
В рамках рассматриваемой концепции на электродинамическом уровне можно выделить три класса задач:
В качестве антенн в моделях ЛА реализованы четвертьволновые вибраторы на металлическом корпусе ЛА. Проводные порты размещены у оснований вибраторов. Пространственное расположение каждого ЛА (рис.2) определяется координатами начала связанной системы координат (СК) каждого ЛА относительно земной СК и ориентацией ее осей, описываемой тремя углами Эйлера (углом курса, углом тангажа и углом крена).
Таким образом, мгновенное положение каждого ЛА характеризуется шестью обобщенными координатами. Для описания сценариев взаимного расположения ЛА в пространстве используется понятие нечетких обобщенных координат (НОК) каждого ЛА группировки (три декартовы координаты начала связанной СК и три угла Эйлера) [1].
Нечеткие обобщенные координаты – это нечеткое подмножество С универсального (базового) множества Q обобщенных координат q! ∈ Q, имеющее нормальную и выпуклую многомерную функцию принадлежности μc(q!), то есть такую, что, во-первых, существует значение носителя, в котором функция принадлежности равна единице, а во-вторых, при отступлении от своего максимума в любом направлении функция принадлежности не возрастает.
Для рассматриваемой задачи в [1] предложено использовать унимодальные НОК, когда функция принадлежности нечеткой координаты С имеет максимальное значение только в единственной точке (точнее, имеет бесконечно малую область толерантности). На каждом носителе НОК для описания физической величины определяем аппроксимирующую функцию φc(q!). Каждая из этих функций может быть построена различными способами: при помощи разных аналитических функций, либо путем обучения искусственных нейронных сетей (можно получить наиболее точную аппроксимацию). Способ построения аппроксимирующих функций для каждой задачи выбирается исходя из двух показателей качества: точности аппроксимации и минимизации вычислительных ресурсов.
Для решения задачи в условиях случайного распределения положения ЛА и, как следствие, характеристик связи и наведения используется следующий алгоритм [3]:
Рассмотрим монопольную антенну, расположенную на металлической плоскости конечного размера в пространстве (рис.3). У такой антенны излучение вдоль вертикальной оси отсутствует. Можно сказать, что это направление, вдоль которого наблюдается подавление сигнала.
Изменение положения диаграммы направленности (ДН) антенны позволяет отслеживать направление появления помехи и подавлять ее, увеличивая таким образом отношение сигнал / шум на входе приемника. Однако при установке такой антенны на корпус сложного летательного аппарата вид ДН заметно меняется (рис.4). Чтобы менять положение приемной антенны на корпусе, можно ДН, рассчитанную как отдельную задачу, использовать в качестве приемной антенны (receiver antenna) и изменять положение только этого объекта для решения задачи связи и рассеяния поля.
Выполним расчет ЛА длиной 20 м на частоте 200 МГц. Время расчета по методу MLFMM на компьютере с ОЗУ объемом 16 Гбайт и тактовой частотой процессора 1 ГГц составляет 10 мин на частотную точку.
Кроме расчета ДН этот метод позволяет определить связь между антеннами, находящимися на корпусе ЛА, в диапазоне частот, а также оптимизировать положение антенн. Трехмерную ДН, записанную в файл, можно использовать как для описания возбуждения такой антенны в виде точечного источника, так и в качестве приемной антенны с такой же ДН. Таким образом, задача взаимодействия ЛА с учетом НОК упрощается и сводится к анализу характеристик радиоканалов. Для оценки полной мощности, принятой антенной, а также рассеяния, наведенного металлическими объектами, в FEKO предусмотрена специальная опция Include Only the Scattered Part ofthe Field (только рассеянная часть поля). Эта опция служит для оценки полной и рассеянной части мощности, принятой отдельными антеннами, установленными на ЛА сложной формы.
Выполним расчет мощности на отдельных приемных антеннах, развернутых в пространстве под углом 0, 10 и 30°, при облучении ЛА плоской волной, падающей под различными углами в диапазоне Θ = 0–90° (рис.5).
Затем поместим элемент рассчитанной ДН антенны на полную модель ЛА и рассчитаем зависимость рассеянной части мощности от угла падения (рис.6). Рассчитав этот параметр во всех направлениях с заданным угловым шагом, получим необходимую информацию для формирования банка данных, используя которые процессор будет корректировать принятые с заданного направления сигналы.
В последних версиях FEKО реализован еще один метод – расчет во временной области. В этом случае падающий радиоимпульс электромагнитного поля, развернутый по временной оси, позволяет выделить временной процесс как падающего, так и отраженного импульса, появившегося на входе антенны из-за отражения от заднего хвостового оперения и от металлических токонесущих частей корпуса ЛА (рис.7).
***
Мощная САПР электродинамического моделирования FEKO позволяет рассчитать диаграммы направленности крупных антенных систем, установленных на летательных аппаратах, в условиях нечетких систем координат. Когда связь между летательными аппаратами, постоянно меняющими свое положение один относительно другого, носит случайный характер, важно обеспечить высокую точность и скорость расчета. В данной статье показано, как с помощью САПР FEKO можно улучшить связь между радиоканалами и свести к минимуму наводимые пространственные помехи.
Авторы статьи: А. Курушин, к.т.н., И. Мюхкеря, С. Подковырин
Статья опубликована номере No7 (00157) 2016 в журнале ЭЛЕКТРОНИКА наука | технология | бизнес.
ЛИТЕРАТУРА
