(1) |
где ; .
В приведенных формулах , — углы сферической системы координат, определяющие направление фазирования, (, ) — координаты центров излучателей в декартовой системе координат, ось ОХ которой ориентирована параллельно продольной оси вибраторов, а ось OZ — перпендикулярна плоскости решетки.
Возбуждение (1) позволяет рассматривать излучение элемента решетки как возбуждение электромагнитного поля в пространственном волноводе — канале Флоке [1]. Для анализа поля в пространственном волноводе используются известные соотношения
(2) |
где , и u — векторный и скалярный потенциалы.
Интегральное уравнение задачи получено из граничного условия для касательной составляющей электрического поля на поверхности идеально проводящего излучателя
(3) |
где — стороннее электрическое поле, равное нулю всюду, за исключением области зазора между плечами вибратора, где оно принимает известное значение , равное отношению возбуждающего напряжения к ширине зазора. Интегрирование граничного условия (3) с учетом (2) дает
(4) |
где С — постоянная интегрирования, a t = x, y.
Электродинамические потенциалы представляются в виде разложения по плоским волнам [2], а учет граничных условий для нормальных составляющих тока на кромках вибратора и уравнения непрерывности позволяет перейти от поверхностной плотности тока к поверхностной плотности заряда. Рассматривая волну потенциала в виде суперпозиции волн, связанных с H- и E-волнами пространственного волновода, учитываются границы раздела между диэлектрическими слоями. При этом выражение (4) сводится к интегральному уравнению вида
которое получено в предположении равномерного распределения продольной составляющей заряда , вдоль оси у. Как показали численные эксперименты, учет особенности указанной составляющей заряда на кромке узкого вибратора не приводит к заметным изменениям получаемых результатов.
Для решения полученного уравнения был использован метод моментов, преобразующий интегральное уравнение в систему линейных алгебраических уравнений. Для этого искомая функция представлена в виде
где — базисная функция; — неизвестный коэффициент. В качестве базисных функций используются собственные функции области, занимаемой металлической пластинкой (для узкого прямоугольного вибратора — набор тригонометрических функций одной координаты).
Используя базисные функции в качестве весовых функций, получаем следующую систему линейных алгебраических уравнений:
где
Решая полученную систему уравнений, можно найти распределение заряда по металлическому вибратору, через которое находятся все интегральные характеристики излучателя в решетке: парциальная диаграмма направленности, входное сопротивление и другие.
Численные эксперименты проводились при следующих параметрах решетки и излучателя (здесь и далее размеры нормированы к длине волны): шаги решетки = = 0,6; ширина вибратора b = 0,035; = 0,175; = 2,5; = = 0°. Предполагалось наличие диэлектрического покрытия толщиной с диэлектрической проницаемостью .
Рис.1
На рис. 1,а приведены зависимости активной и реактивной составляющих входного сопротивления печатного вибратора в бесконечной решетке от его длины. Параметром является диэлектрическая проницаемость подложки , равная: 1,5 (кривая 1); 2,0 (кривая 2) и 2,5 (кривая 3). При увеличении происходит уменьшение резонансных размеров вибратора, однако при этом сужается рабочая полоса частот. При уменьшении ширины печатного вибратора b наблюдается увеличение активной составляющей входного сопротивления в области второго резонанса и заметное сужение полосы частот. На рис. 1,б b имеет значения: 0,035 (кривая 1); 0,055 (кривая 2) и 0,075 (кривая 3). Можно также отметить некоторое смещение резонансной длины в сторону больших значений длины вибратора. Аналогичное поведение резонансной длины наблюдается при изменении высоты расположения вибратора над экраном. На рис. 1,в имеет значения: 0,150 (кривая 1); 0,175 (кривая 2) и 0,2 (кривая 3).
Входное сопротивление вибраторного излучателя в решетке существенно зависит от шага периодической структуры, что является следствием взаимодействия излучателей в ФАР. На рис. 2,а и б представлены зависимости входного сопротивления печатного вибратора шириной b = 0,075 при изменении шага решетки и соответственно в Е- и H-плоскостях. Кривые на рис. 2, а соответствуют значениям: 1 — = 0,7; 2 — = 0,8 и 3 — = 0,9, а на рис. 2,б = 0,6 (кривая 1); 0,65 (кривая 2) и 0,7 (кривая 3).
Рис.2