Оценка воздействия мощного электромагнитного импульса на антенные системы летательного аппарата (на примере самолета Боинг 737) Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.3

Семёнов А.А., Пивоваров В.В.

Оренбургский государственный университет E-mail: toliksem@rambler.ru

ОЦЕНКА ВОЗДЕЙСТВИЯ МОЩНОГО ЭЛЕКТРОМАГНИТНОГО ИМПУЛЬСА НА АНТЕННЫЕ СИСТЕМЫ ЛЕТАТЕЛЬНОГО АППАРАТА (НА ПРИМЕРЕ САМОЛЕТА БОИНГ 737)

Исследовано влияние мощного электромагнитного импульса ядерного взрыва на радиосвязные и навигационные антенно-фидерные устройства летательного аппарата. Предложены впервые методики оценки воздействия мощных импульсов для разных типов бортовых антенн: вибраторных, апертурных и рамочных. Определены антенные системы, наиболее подверженные влиянию мощных электромагнитных полей.

Ключевые слова: электромагнитный импульс ядерного взрыва, радиосвязные и навигационные антенны, бортовое оборудование.

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

В связи с усовершенствованием и усложнением РЭС, увеличением их быстродействия проблема оценки и обеспечения стойкости к воздействию мощного ЭМИ становится наиболее актуальной [3],[5].

Меры, направленные на повышение стойкости РЭС непосредственно связаны с расчетной оценкой уровней наведенных напряжений и токов на входах антенных систем. Методы оценки стойкости РЭС и антенных систем ЛА, предложенные в [6],[7], основаны на решении интегро-дифференциальных уравнений при использовании численных методов, таких, как метод моментов, метод сингулярных разложений и т. д. Эти методы достаточно сложны, требуют разработки больших программ и значительных затрат машинного времени.

Таким образом, целью статьи является разработка методики оценки воздействия мощных электромагнитных полей по антенно-фидерно-

му тракту. Такая методика должна удовлетворять требованиям универсальности, т. е. возможности её применения для любого типа летательного аппарата, что позволит упростить расчетные соотношения и уменьшить вычислительные затраты.

При этом наибольший интерес представляет решение следующих вопросов:

1. Оценка полной энергии, амплитуды и формы помехи в АФУ в полосе пропускания приемного тракта;

2. Предварительная оценка возможных последствий, вызванных прохождением энергии ЭМИ через АФУ.

1. Источники мощных ЭМИ

Основными источниками мощных электромагнитных излучений естественного и искусственного происхождения являются: грозовые разряды, мощные радиопередающие средства (РПС) и радиолокационные станции (РЛС), высоковольтные линии электропередач (ВЛЭП), контактная сеть железных дорог (КСЖД), а также генераторы сверхширокополосного и микроволнового излучения. Наиболее мощным искусственным источником электромагнитного излучения является ядерный взрыв (ЯВ). Импульсные электромагнитные поля, возникающие при ядерном взрыве, принято называть электромагнитным импульсом ядерного взрыва (ЭМИ ЯВ).

В соответствии со стандартом МЭК серии 61000 [1] выделяют три стадии развития высотного ядерного взрыва: ранняя, промежуточная и поздняя. На рисунке 1 представлен ЭМИ ЯВ на ранней стадии.

В м 50000 45000 40000 35000 30000 25000 20000 15000 10000 5000 О

/ \

' / rsr—? /......... \ ,,,................................ 1

/ \

1....... ■ \ ......

\

г г...........f...............

......f~.............." ^•sj'

1 'l

.......—J

О

10 20 30

40 50

60

SO 90 100 fie-

Рисунок 1. ЭМИ ЯВ на ранней стадии

Для всех стадий поведение электрического поля определяется соотношением [2]

E(t) = Emax k(e-at - е-в), (1)

где Emax - максимальная амплитуда импульса, В/м. Для ранней стадии ЭМИ ЯВ имеет максимальную амплитуду 50 кВ/м; k - нормировочный коэффициент k=1,3(для ранней стадии, k=1 для промежуточной стадии); а, в- показатели экспоненты, с (для ранней стадии а=4-107 с , в=6-108 с ; для промежуточной -а=103 с"1, в=6108 с"1).

2.Методика расчета воздействия ЭМИ

на вибраторные антенны

В тех случаях, когда каналом проникновения ЭМИ в аппаратуру являются вибраторные антенны, могут быть использованы различные модели и методы расчета процессов под влиянием импульсных электромагнитных полей. Наиболее общей исходной моделью являются уравнения Максвелла с соответствующими граничными условиями в виде интегро-дифференциального уравнения [4]. Однако, этот метод достаточно сложен и требует разработки трудоемких программ и значительных затрат машинного времени. Полезным приближением являются модели антенн в виде линий передачи - теория длинных линий [5]. Согласно данному методу через значение действующей длины антенны определяется ЭДС, наводимая на входе АФУ. Достоинство этого подхода заключается в сравнительной простоте и наглядности, возможности получения приближенных аналитических решений.

На вход радиоприемника действует та часть электромагнитного поля (ЭМП), которая проходит через вибраторную антенну. Поэтому наиболее полным будет рассмотрение процессов, возникающих под действием ЭМП в системе «антенна-радиоприемник».

В теории антенн ЭДС на входе вибраторной антенны определяется как [3]

е = ЕтЛ • Р в() (2)

где 1д - действующая длина антенны, м;

Р(в,р) - нормированная диаграмма направленности антенны в направлении, определяемом углами в и р. Будем рассматривать предельный случай, когда Р (в,р) = 1.

Так, действующая длина симметричного полуволнового вибратора 1Д=Х/к, а несимметричного четвертьволнового вибратора 1Д=^/2п. Отметим, что для директорных антенн поле от каждого излучателя складывается синфазно.

3. Методика расчета воздействия ЭМИ

на апертурные антенны

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

V = ^ } е

0

(3)

где Е - напряженность электрического поля ЭМИ, а начало отсчета времени совпадает с моментом прихода ЭМИ;

20 - характеристическое сопротивление свободного пространства, 2=311 Ом; 5 - площадь раскрыва антенны.

Вместо 5 необходимо подставить 5эф - эффективную площадь раскрыва антенны, которая связана с геометрической площадью соотношением [4]

^ ^, (4)

где V- коэффициент использования площади антенны, зависящий от неравномерности амплитудно-фазового распределения поля в апертуре антенны и всегда меньше единицы. При равно-амплитудном синфазном распределении v=1 Принимается V = 0,5-0,6, так как мощность, которая отдается приемной антенной в согласованную нагрузку уменьшается за счет потерь.

Подставив (4) в (3), получим определение энергии в частотной области

V = | |Е(]ю)|2 йю, (5)

о ю

Найдем выражение для энергии Ц в случае воздействия мощного ЭМИ, описываемым дву-хэкспоненциальным выражением. В данном случае спектральная плотность ЭМИ выражается как:

е (ю) = еш

1

1

(6)

Квадрат модуля спектральной плотности определяется:

|е (ю)|2 = е 2

в-а

в + а

1

1

а2 + ю2 в2 + Ю

(1)

Подставив (1) в (5), получим

V =

5

эф

ю

в-а

2Хг,

ю

в + а

1

1

а2 + со2 в2 + ю2

1

Етах 5 эф в —а

220 в + а 1{а2 + ю2 в2 + с

л а+

йю =

йю =

= Ешах 5эф в -а

2г0 в+а

[ —-2\~Б2-2йю

J а1 +сС в +Ю

ю а + ю с

Е1шах 5 эф в —а ■ X

2г0 в + а

х

^2 1 11

--[ —— йю

с а2 + со2 в С со2 , 1

ю ^ ю--+ 1

в2

= \и = С;йи = — йю^ [—т1— = ашя(и) 1 = 1 в в }и2 +1 5

^ 1 J а2 +

ЕШах5эф х в-а /\

220, в + а

а +ю

\arctg (и 2) - аг^и^]

= {/ =ю;й/ = -йю^ ¡-Г— = arctg(Г) 1 = а а */ +1

Е1шах5эф в -а

2г0 в + а

а\arctg Г2) - аге1§(Г!)] - в-х

х[аг^ (и2) - аге1§(и!)]} =

= ЕШах 5 эф в-а

в + а

аг^

а

- агй§

ю в

- / \ / \ - \

ю ю

arctg 2 - агй§ I

- 1в ) _ /

(8)

Таким образом, энергия помехи на входе приемного тракта в полосе его пропускания для апертурной антенны определяется как

V = Ешах 5эф в -а

220 в + а

arctg

юЛ а

- агй§

/ л

ю в

arctg

Щ

-

ю в

(9)

где ю1 - нижняя рабочая частота системы; ю2 - верхняярабочая частота системы.

4. Методика расчета воздействия ЭМИ на рамочные антенны

В отдельный класс выделяются рамочные антенны. Для таких антенн ЭДС, наводимая в витках рамки, определяется скоростью изменения магнитного потока [5]

йФ

(10>

величина которой, в свою очередь, зависит от числа витков п, их площади 5, абсолютной магнитной проницаемости сердечника рамки ца, напряженности магнитного поля приходящей волны Д ориентации рамки относительно угла в в горизонтальной плоскости между нормалью к фронту падающей волны и плоскостью рамки. С учетом вышесказанного ЭДС в рамке будет равна

■ е = - ja>| anSHsin 8,

(11)

Формула может быть представлена в другом виде для свободного пространства

е = - jknSE sin 8,

(12)

где k - волновое число, k =

2п Я '

Наличие мнимой единицы в формуле (12) свидетельствует о том, что ЭДС в рамке отстает по фазе на п/2 от напряженности электрического поля приходящей волны.

5. Вычисление наводок в антенных системах самолета Боинг 737

Авиационные радиопередачи подразделяются на связь, навигацию и опознавание. Решение столь разнообразных задач обеспечивается антеннами различных типов. На рисунке 2 показано размещение антенн на корпусе самолета Боинг 737.

Рассчитано воздействие мощного ЭМИ на входы приемных трактов АФУ согласно представленной методики. За исходные данные расчетов принято, что максимальная напряжен-нос.ть поля Етах. = 50 кВ/м, соответствующая ГОСТу. Параметры и характеристики оборудования, а также вычисленные наводимая ЭДС (для вибраторных антенн) и энергия помехи (для апертурных антенн) представлены в таблице 1.

6. Вычисление воздействия МЭМП на вибраторные антенны с помощью эквивалентных схем

Оценить влияние ЭМИ на бортовые вибраторные антенны можно, используя методы расчета переходных процессов в антеннах под влиянием импульсных электромагнитных полей. Для этого следует определить ток антенны в приемном режиме. В качестве основы

1 - метеолокатор; 2 - антенна системы предупреждения столкновений; 3 - антенна системы управления воздушным движением; 4 - радиовысотомер; 5 - антенна радиодальномерной системы; 6 - антенна командной станции УКВ; 7 - рамочная антенна радиопеленгатора; 8 - антенна маркерной системы; 9 - GPS антенны; 10 - антенна связной станции КВ; 11 - антенна системы VOR; 12 - глиссадная антенна; 13 - курсовая антенна; 14 - ненаправленная антенна радиопеленгатора.

Рисунок 2. Размещение антенн на корпусе самолета Боинг 737

Таблица 1. Параметры и характеристики оборудования, а также вычисленные наводимая ЭДС (для вибраторных антенн) и энергия помехи (для апертурных антенн)

Тип Частоты, МГц Тип антенны Размер антенн Действующая длина, м Поляризация поля ЭДС (на fc) или энергия помехи

Связная станция KB 2 - 30 пазовая (щель в корпусе) /=1,72 м / = 0 23 1возбуст вертик 22,5 кВ

Командная станция УКВ 118 - 137 четвертьволновый вибратор 1=0,41 м 0,187 вертик 18,7 кВ

Маркерное оборудование 75 вибратор /=0,24 м /Д= /а=0,24 м вертик 24 кВ

Система управления воздушным движением ПРМ-1030 ПРД-1090 лопастного типа четверть волн. вибратор /лоп=0,058м 0,046 вертик 2,3 кВ

Система предупрежд. столкновений ПРМ-1090 РД-1030 директорная из4х элементов /=0,16м /д=4/дод=0,18 вертик 8,8 кВ

Радиовысотомер 4200-4400 рупорная ap= 0,048м bp= 0,069м - вертик 5 ■ 10-15 кВт

Оборудован. радиодальномерной системы 962-1215 лопастного типа /лоп=0,058м 0,044 вертик 2,2 кВ

Оборудование угломерной системы посадки 108-118 четвертьволновый вибратор 0,53 м 0,42 горизонт 21 кВ

108-112 четвертьволновый вибратор 0,5 м 0,43 горизонт 21,5 кВ

328,6-335,4 директорн. из 2х элементов 0,4 м 2/дод= 0,37 горизонт 12 кВ

Радиопеленгатор 0,15-1,75 вертикальный штырь 0,2 м 0,1 кругов 5 кВ

Метеорологический радиолокатор 9300-9500 многощелевая АР с волноводным входом S=0,395 м2 - кругов 2,56 ■ 10"8 кВт

Система GPS 1575,42 микрополосковая a=0,12 м b=0,07 м - эллипт. 7,8 ■ 10-10 кВт

анализа применяется метод, базирующийся на представлении антенн эквивалентными схемами Тевенина.

При этом если линейный размер антенны, взаимодействующей с помехонесущим полем, мал, то анализ её взаимодействия с ЭМП значительно упрощается, т. е. эквивалентные схемы используются, если [4]

1<сЬф (13)

где Ьф - длительностьфронта ЭМИ;

с - скорость распространения электромагнитных волн.

В зависимости от того, будет ли линейный размер антенны меньше, соизмерим или больше длины волны внешней помехи, может быть применен при анализе соответственно квазистатический, резонансный или волновой подход.

Рассмотрим применение резонансного подхода. Антенна в общем случае является резонансной системой с большим числом резонансных точек. Основной вклад в резонансные явления, которые наблюдаются в антенне, вносит точка первого резонанса.

Хорошую аппроксимацию низкочастотных характеристик антенн можно получить с помощью эквивалентной схемы, приведенной на рисунке 3, где первая резонансная точка антенны моделируется введением последовательной индуктивности Ь , образующей резонансный контур с емкостью антенны Са на частоте первого резонанса, возбуждаемого эффективной ЭДС [4]

/ = - = -. Зр Я 41

(14)

Сопротивление (рисунок 3) - сопротивление излучения антенны в первой резонансной точке. Оно составляет 35 Ом для несимметричного вибратора и 70 Ом для полуволнового.

Из эквивалентной схемы видно, что наведенный в антенне ток описывается выражением:

Е (ю) • I д

I (Ю):

(15)

(ЮС) + ю + ки+кн

В общем случае напряжение на нагрузке антенны Я при внешнем воздействии ЭМИ в форме экспоненциального импульса согласно [7] ик (()= Еши1дсаюркн (е- е-|3' )sin ^ (16)

Для определения напряжений и токов, наведенных в бортовых антеннах, необходимо знать их электрические параметры, которые соответствуют элементам, входящим в эквива-

Рисунок 3. Низкочастотная резонансная схема замещения несимметричного вибратора

лентную схему замещения, и находятся в прямой зависимости от линейных размеров и конструктивных особенностей антенн. В таблице 2 приведен расчет электрических характеристик бортовых антенн, наиболее подверженных влиянию мощного ЭМИ.

Таким образом, приближенное аналитическое выражение для максимального наведенного напряжения можно представить:

итах = ЕтаДдСа®р-^Н •

(17)

Как следует из (17) реакция бортовых антенн на импульсное воздействие следующая: максимальное наведенное напряжение для командной станции УКВ принимает значение 72,04 кВ; для оборудования угломерной системы посадки - 22,49 кВ; для бортовой системы предупреждения столкновений - 11,75 кВ; для маркерного оборудования - 7,23 кВ.

Таблица 2. Расчет электрических характеристик бортовых антенн, наиболее подверженных влиянию мощного ЭМИ

Тип Антенна Электрические характеристики

Ьр, Гн Са, Ф

Расчетная формула Значение Расчетная формула Значение

Станция УКВ несимметричный вибратор 1 ЮрСа 2,82 • 10"8 2пе01 6,7 • 10-12

1п ' 21 ^ а V У -1

Маркерное оборудование симметричный вибратор 1 ЮрСа 2,1 • 10-8 пе01 3,07 • 10-12

1п ' 2' а V У -1

Оборудование угломерной системы несимметричный вибратор 1 ЮрСа 2,95 • 10-8 2пе01 9,55 • 10-12

1п ' 21 ^ а V У -1

Система предупреждения столкновений несимметричный вибратор 1 ЮрСа 5,6 • 10-5 2пе01 4,99 • 10-12

1п ' 21 ^ а V -1

Выводы

Результаты проведенных расчетов по методике показывают:

1. Значения наведенных напряжений на нагрузке антенны могут быть значительными, что подтверждается проведенными вычислениями. Так при воздействии на антенну командной станции УКВ электромагнитного поля с амплитудой Етх = 50 кВ/м и временными параметрами а=4-107 с-1, Р=6-108 с-1 можно ожидать значения напряжения, наведенного на ее 50-ом-ной нагрузке, равное 72,04 кВ;

2. На входы угломерной системы посадки и бортовой системы предупреждения столкновений наведенные напряжения могут достигать значений 22,49 кВ и 11,75 кВ соответственно;

3. Незначительное влияние будет оказываться на такие устройства, как радиовысотомер, метеорологический радиолокатор и система GPS, поскольку антенны таких систем относятся к апертурным, а результаты расчета показывают, что данный тип в меньшей степени реагирует на ЭМИ ЯВ;

Полученные результаты показывают, что воздействие ЭМИ ЯВ на антенные системы приведет к образованию сильных помех, что может привести к выходу из строя станция УКВ без установки защитного устройства.

Полученные результаты показывают, что образование помех от воздействия ЭМИ ЯВ на антенные системы без защитного устройства может привести к выходу из строя станции УКВ.

1.05.2014

Список литературы:

1.МЭК 61000-1-3. Электромагнитная совместимость (ЭМС). «Устойчивость к электромагнитному импульсу высотного ядерного взрыва (ЭМИ ВЯВ)», 2003.

2. Балюк, Н.В. Мощный электромагнитный импульс: воздействие на электронные средства и методы защиты / Н.В. Балюк, Л.Н. Кечиев, П.В. Степанов. - М.: ООО «Группа ИДТ», 2007. - 478 с.

3. Balanis, C.A. Antenna Theory. Analysis and Design / C.A. Balanis. - Wiley & Sons, 2ndedition, 1997. - 600 p.

4. Вольман, В.И. Техническая электродинамика / В.И. Вольман, Ю.В. Пименов. - М.: Радио и связь, 1995.-298 с.

5. Мырова, Л.О. Обеспечение стойкости аппаратуры связи к ионизирующим и электромагнитным излучениям / Л.О. Мырова, А.З. Чепиженко. - М.: Радио и связь, 1988. - 296 с.

6. Матье, М. Радиорелейные системы передачи / М. Матье ; пер. Л.О. Мыровой; под ред. В.В Маркова. - М.: Радио и связь, 1982. - 280 с.

7. Кинг, Р. Антенны в материальных средах. 2 книги ; пер с англ. / Р. Кинг, Г. Смит. - М.: Мир, 1984. - 416 с.

Сведения об авторах:

Семенов Анатолий Александрович, аспирант кафедры проектирования и технологии радиоэлектронных средств Оренбургского государственного университета

Пивоваров Валерий Викторович, доцент Оренбургского государственного университета,

кандидат технических наук

460000, г. Оренбург, Шарлыкское шоссе 5, ауд. 14402, тел. (3532) 372556, e-mail: toliksem@rambler.ru