Математическое моделирование характеристик рассеяния объектов локации с импедансной поверхностью

Размещено на http://www.allbest.ru/

Математическое моделирование характеристик рассеяния объектов локации с импедансной поверхностью

В задачах расчета входных сигналов систем ближней радиолокации [1] возникает проблема анализа электромагнитных полей, отраженных от объектов сложной геометрической формы и электрофизической структуры. В последние годы в мире накоплен большой опыт по решению задач определения входных сигналов радиолокационных систем с использованием теории дифракции электромагнитных волн на объектах локации. Метод физической теории дифракции (ФТД) [2] позволяет учесть такие дифракционные явления, как рассеяние волн на острых кромках цели и двукратные переотражения волн между ее элементами.

Математическое моделирование характеристик рассеяния цели может быть основано на использовании ее полигональной модели, которая представляет поверхность объекта в виде совокупности простейших геометрических примитивов (плоский треугольный элемент и острая кромка). Наиболее передовым и мощным инструментом 3D моделирования сложных объектов на сегодняшний день является Creo 2.0 (Pro/Engineer), который позволяет синтезировать полигональную (фацетную) модель внешней поверхности объекта максимально точно и эффективно (рис. 1).

В работе [1] на основе метода ФТД обосновано использование выражения суммарного поля, рассеянного объектом локации, в виде когерентной суммы трех компонент [1]:

. (1)

Первая компонента в соотношении (1) представляет собой величину результирующей напряженности электрического поля, рассеянного на всех плоских треугольных элементах, образующих полигональную модель объекта. Вторая составляющая является результирующей напряженностью поля, дифрагировавшего на острых кромках полигональной модели. Третья компонента поля определяется в результате суперпозиции полей рассеяния от треугольных элементов, для которых источники возбуждения вторичных волн наводятся за счет двукратных переотражений от других элементов полигональной модели.

Рис. 1. Полигональная модель истребителя F-22

Компонента, обусловленная рассеянием волн на острых кромках цели, может вносить существенный вклад в результирующее рассеянное поле (1) на некоторых ракурсах наблюдения. Структура острых кромок и изломов на гладкой поверхности различных объектов приведена на рис. 2.

В соответствии с полигональным описанием геометрии цели и с использованием соотношений, полученных в [1], выражение для напряженности поля, рассеянного отдельным элементарным отражателем (треугольник или кромка), имеет вид:

, (2)

, (3)

, (4)

, (5)

где - поляризационная матрица рассеяния (ПМР) элемента, м;

() - компонента вектора напряженности рассеянного поля в направлении векторов поляризации, В/м;

- расстояние между элементарным отражателем (ЭО) и фазовым центром приемной антенны, м;

P_i - мощность сигнала, излучаемого передающей системой, Вт;

D_i - КНД передающей антенны в направлении максимума излучения;

Z₀=120? - волновое сопротивление воздуха, Ом;

- расстояние между фазовым центром антенны и началом системы координат ЭО, м;

- круговая частота электромагнитных колебаний, рад/с;

F_im - значение функции направленности передающей антенны в направлении на m-й ЭО.

Рис. 2. Геометрия изломов гладкой поверхности объектов локации

ПМР описывает рассеивающие свойства элемента в зависимости от поляризационных характеристик передающей и приемной антенн. Каждый элемент ПМР (5) представляет собой комплексный коэффициент отражения [1,3], который определяется из следующего соотношения:

, (6)

где - величина напряженности поля рассеяния, которая вычисляется для случая, когда излучение и прием ведутся на линейной поляризации, В/м.

В соответствии с (6) комплексный коэффициент отражения ЭО выражается в [м] и не зависит от амплитуды электромагнитной волны, падающей на ЭО, а также от расстояния до фазового центра передающей антенны [3]. При этом комплексный коэффициент отражения находится в функциональной зависимости от геометрических размеров ЭО, его электродинамических параметров, а также от его ориентации в пространстве относительно передающей и приемной антенн.

Величина напряженности поля, рассеянного m-м элементом, принятая в апертуре приемной антенны, определяется в результате скалярного произведения вектора напряженности рассеянного поля вблизи приемной антенны на вектор поляризации приемной антенны в направлении на данный отражатель [1]:

. (7)

В результате имеем выражение для вычисления напряженности поля, рассеянного m-м отражателем, в апертуре приемной антенны:

, (8)

. (9)

С использованием соотношения (8), для каждого элемента вычисляется напряженность рассеянного поля в апертуре приемной антенны, а затем эти поля суммируются для получения напряженности результирующего дифракционного поля, рассеянного объектом:

. (10)

Как было показано в [1] ЭПР цели представляется в виде суммы по всем комплексным коэффициентам отражения от элементов полигональной модели:

. (11)

Для треугольника расчет осуществляется с учетом его электродинамических параметров (диэлектрическая проницаемость, удельная проводимость, радиопоглощающие покрытия), а соответствующая математическая модель приведена в [1] и основана на методе физической оптики. При этом для острых кромок (клиновидных структур) в [1] получены соотношения комплексного коэффициента отражения только для случая идеальной проводимости кромки. Поэтому представляет интерес определение комплексных коэффициентов отражения для импедансных клиновидных структур (острых кромок). В основе электродинамической модели взаимодействия плоской электромагнитной волны с острой кромкой лежит физическая теория дифракции (ФТД), которая основана на уточнении поверхностного тока, полученного методом физической оптики. Эта теория известна также под названием метода краевых волн (МКВ). Согласно МКВ, полный ток, наведенный падающим полем на поверхности треугольного элемента полигональной модели, представляется в виде суммы тока , полученного в приближении физической оптики и тока , являющегося поправкой к и обусловленного всеми эффектами, которые не учитываются физической оптикой: . Токи и считаются соответственно равномерной и неравномерной частями полного тока. Поле, рассеянное равномерной частью тока, является суммой геометрооптического и дифракционного поля в приближении физической оптики. Неравномерная часть тока обусловлена возмущениями, создаваемыми отклонениями поверхности от бесконечной плоскости, например из-за ее изломов (ребер). Для вычисления напряженности электрического поля, рассеянного острой кромкой, используется схема радиолокации [3], представленная на рис. 3. Выражение для вектора напряженности поля рассеяния острой кромки с произвольными электрическими свойствами имеет вид [3]:

, (12)

где - направляющий вектор отраженной волны из вершины кромки P₁ в фазовый центр приемной антенны, - функция Грина для данной задачи [1].

Эквивалентные токи на поверхности острой кромки определяются на основе соотношений Нотта и Сеньора [3,4]:

, , (13)

где F, G - угловые функции, которые определяют поле, рассеяное неравномерной частью тока (для вертикальной и горизонтальной поляризации), - пространственные углы между направляющими векторами и и вектором .

В соответствии с [1] комплексный коэффициент отражения острой кромки определяется следующим образом:

, (14)

, (15)

где L - длина ребра, м; ( - проекция разностного вектора на ось x).

Таким образом, задача определения эквивалентных токов сводится к вычислению функций F и G, которые определяют поле, рассеянное неравномерными токами вблизи ребра клина. Для идеально проводящего клина данные функции хорошо известны [2] и определяются в системе координат клина, как это показано на рис.3. Функции F и G вычисляются с использованием следующих обозначений: ? - внешний угол острой кромки P₁P₂ (рис. 4), - углы визирования кромки для падающей и отраженной волн.

Рис. 3. Схема радиолокации острой кромки полигональной модели

Рис. 4. Система координат для определения коэффициентов дифракции

Для идеально проводящего клина в случае моностатической локации () функции F и G определяются соотношениями [2]:

, , (16)