CC BY
51
И. Е. Плещинская, Н. Б. Плещинский, И. В. Сабиров
ДИФРАКЦИЯ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ НА ЧАСТИЧНО ЭКРАНИРОВАННОМ СЛОЕ
И ОПТИМИЗАЦИЯ ПАРАМЕТРОВ УСТРОЙСТВА ВВОДА ИЗЛУЧЕНИЯ
Ключевые слова: задача дифракции, периодическая решетка, метод интегрально-сумматорных тождеств, ввод излучения,
частично экранированный слой.
Исследуется процесс возбуждения электромагнитных колебаний в частично экранированной диэлектрической пластине и анализируется возможность оптимизации параметров решетчатого узла ввода излучения в слой. Приводится график зависимости потока энергии через сечение частично экранированного слоя от угла падения плоской волны и ширины щели периодической решетки.
Keywords: diffraction problem, periodical grid, integral-summatorial identity method, radiation input, partially screened layer.
Process of initiation of electromagnetic oscillations in partially screened dielectric plate is investigated and possibility of optimization ofparameters of grating knot of radiation input into a layer is analyzed. The graphic of dependence of energy stream on the angle of a plane wave and width of a crack of periodic grating through the section of partially screened layer is provided.
Введение
Слоистые композитные материалы имеют уникальные физико-механические свойства, что обуславливает их широкое применение в различных конструкциях [1], [2]. В ряде случаев в слоистых пластинах могут распространяться акустические, электромагнитные или упругие волны.
Для возбуждения периодических колебаний в многослойных пластинах внешним полем можно использовать периодические решетчатые структуры. В [3] изучался процесс преобразования электромагнитного излучения слоистым
композитом с тонкими проводящими периодическими включениями. В работе [4] рассматривалась возможность использовать в качестве сканирующего экрана слоистый композит, армированный тонкими проводящими пластинами.
В данной работе исследуется процесс возбуждения электромагнитных колебаний в частично экранированной диэлектрической пластине и обсуждается возможность оптимизации параметров решетчатого узла ввода излучения в слой.
1. Задача дифракции электромагнитной волны на частично экранированном слое
Пусть бесконечный диэлектрический слой О < z < /7 (в декартовой системе координат)
экранирован снизу идеально проводящей пластиной z = О. На верхней поверхности слоя размещена периодическая система из идеально проводящих бесконечно тонких экранов. На слой сверху набегает плоская электромагнитная волна. В результате ее дифракции появляется электромагнитное поле над решеткой (уходящие на бесконечность волны) и поле в слое диэлектрика. Электромагнитные волны в слое в общем случае переносят энергию вдоль него.
Исследуем зависимость потока энергии через поперечное сечение пластины от параметров волноводной структуры и от характеристик возбуждающей волны.
Рассмотрим двумерный (плоский) случай, когда система экранов на поверхности слоя представляет собой £ -периодическую решетку из тонких проводящих лент, параллельных оси .
Выберем гармоническую зависимость поля от времени в виде ехр(- /О)^ .
Для поля ТЕ-поляризации ненулевые
компоненты выражаются через потенциальную функцию и г ) - решение уравнения
Гельмгольца:
,, - 1 ди 1 ди
Еу=и, Нх=-
iou/Уо/У dz
Hz =
/CO/JQ/J дх
д^и
д^и
2 + к и <$,z )= 0, к = ш ц^це^е.
дх дz
Достаточно исследовать процесс дифракции только в плоскости у = 0, а точнее, в
полуплоскости ^ > /7 (область 1) и в полосе
О < z < /7 (область 2). Будем помечать индексами 1 и 2 параметры сред, заполняющих эти области, а также потенциальные функции искомых электромагнитных волн.
Пусть на периодическую структуру падает сверху плоская ТЕ-волна с потенциальной функцией цО ^ 2 $\пд°х-1'к^ соэв0 £-Л)
Нужно найти решения уравнения Гельмгольца в области 1 при к = к\ и в области 2 при к = к2, удовлетворяющие следующим граничным условиям и условиям сопряжения при ^ = 0 и при z = /7:
и (*, 0 + 0 )=0, и (V, /7 - 0 )= 0, и (V, /7 + 0 )= 0 на лентах, и (*, /? - 0 )= и (*, /? + 0 I
ди . , _ . ди . , „
- (Г, /7 — 0 )=- (Г, /7 + 0
дд
на щелях.
Решение уравнения Гельмгольца в области
1 должно также удовлетворять условию излучения,
то есть соответствовать уходящей на бесконечность электромагнитной волне.
2. Парное сумматорное функциональное уравнение
Так как система экранов на поверхности пластины является периодической, то будем искать решения уравнения Гельмгольца в каждой из областей как квазипериодические функции (волны Флоке) вида
+<»
и Xип
п = -ю
к >е
,2тт /-----пх
здесь а - параметр Флоке.
Легко видеть, что коэффициенты Флоке ип £ ) должны удовлетворять дифференциальным уравнениям
и„ +
2тт
/ І I Обозначим
п + а
ип £ )= 0, /7 = 0, + 1,...
V_//7
II
*,2-
277
т
-п + а
и условимся, что или >0 , или ||П)/уя >0 .
В фундаментальной системе решений дифференциального уравнения содержатся две
функции: е ^'1Г,г и 1Г,г. При сделанных
предположениях (ех р(-/о^)) первая из них определяет или приходящие с бесконечности волны, или неограниченно возрастающие при z ^ +да . Вторая функция -- наоборот, или уходящие на бесконечность волны, или затухающие при Z ^ +да .
Условие излучения для области 1 сводится к тому, что в общем решении уравнений должны остаться только слагаемые с экспонентами второго типа. Поэтому будем искать
+го
2>„2>„.
П=-ю
В области 2 должно быть выполнено граничное условие и (*", 0 + 0 )= 0 . Следовательно, произвольные постоянные в общем решении ип£УРпе-'У2 П +Яп е2 должны удовлетворять условию Рп + ЦП = 0 . Тогда потенциальную функцию можно записать в виде
,2
б/1 ^,ЕУЄІ0Х
+ГО
и
Тьпвту2пЕ-еп(!< ) , 0 < z < /7.
П=-ю
Из граничных условий и условий сопряжения при z = /7 следует, что при всех X
и° ^,/7 (М )=и2 )
или
о +^
е1в х+еіах ХапЄп*>
/7=-х
= е/ах ХпугпЬ-еп* ) .
П=-ж
Отсюда
а = ^зіп#0
1 + зо = Аз ап=Ьп елпу2пЬ при /7^0.
Исключим неизвестные е. Для
определения неизвестных е остаются два равенства:
+ГО
^апеп (г )= -1 на лентах,
/7=-х
X а/7 [К1/7 + /'К2пйд К2е] -еп К >
/7=-ю
= К10 - 7К2Ос47 К20л на щелях.
3. Бесконечная система линейных алгебраических уравнений
Так как оба равенства в парном уравнении -равенства периодических функций, то достаточно их рассматривать на частях М и N отрезка [0, £], которые относятся к лентам и щелям соответственно.
Удобнее, когда одно из равенств в парном сумматорном функциональном уравнении является однородным. Введем новые искомые коэффициенты
„ „ Кю -/'К20 с(ду2о
°0 = а0---------:----:-----.
+
сП = а,7 = Ьпъту^^ при /7^0.
Тогда
+Ю
/7=-х
'Х^СпЄп ^ )= ^ на
э
+го
X спПпеп $ )= 0 на /V,
/7=-х
где
^ = 2/К20ФТ20 _ ^ ^^ + ,у2/7С^ у2пь
$0
п = 0, ±1,...
Методом интегрально-сумматорных
тождеств [5] парное сумматорное функциональное уравнение преобразуется в бесконечную систему линейных алгебраических уравнений (БСЛАУ)
,| ■+» +ю ^
-£С/< + — ^ !п-т^т-к = %'!-к,
п=-<х т=-ю
т
к = 0, ±1,...
где
'п=$Єп(1іЯ. ^п = ІЄп(№ п = °.±1— М N
и
Более сложная БСЛАУ в задаче дифракции электромагнитной волны на частично экранированном слое была получена в работе [6].
Приближенное решение бесконечной системы уравнений может быть найдено методом усечения.
4. Поток энергии через сечение слоя
Рассмотрим поперечное сечение плоскостью X = 0 слоя 0 < Z < /7 . Вычислим поток энергии через это сечение (нормированный на единицу длины оси ) электромагнитного поля, возбуждаемого падающей сверху плоской волной. Так как нормальная компонента среднего значения вектора Пойнтинга
пх =
1
(
2ед//
-Re
IU-
д
* '\
д
то поток энергии через сечение л +х>
Fx =
Re
'/71 Х
^ =-а> я2 =-с
/?2
а + ^ГП2 \'Sn''n2'
5. Оптимизация параметров узла ввода излучения
Исследовать на экстремумы функцию Гх с/ ] слишком сложно. Поэтому для определения значений параметров и ,
доставляющих ей глобальный экстремум, можно использовать следующий алгоритм.
На отрезке [[тяглах] выберем
Л/ значений к = 0... N -1. Получим методом проекции градиента решения одномерных задач X , )~> еХ, с/ е [0ггпп’°тах1 к = 0...Л/-1.
Среди найденных пар (#^, с//) к = 0... N -1
можно выбрать наилучшие и уточнить значения методом проекции градиента для целевой функции двух переменных.
Работа выполнена при поддержке РФФИ и АНРТ, грант 12-01-97012-р-поволжье-а.
где
’пьп2
jsinj/21 -SinK2tf2Z 'dz
о
(* - операция комплексного сопряжения).
Вычислительный эксперимент показал, что зависимость потока энергии через сечение частично экранированного слоя от угла падения плоской волны от внешнего источника и ширины щели периодической решетки является
многоэкстремальной функцией = Рх с/ ) (см. рис. 1).
конструкций из Машиностроение,
Литература
[1] В.В. Васильев, Механика
композиционных материалов.
Москва, 1988. 272 с.
[2] В.И. Королев, Слоистые анизотропные пластинки и
оболочки из армированных пластмасс.
Машиностроение, Москва, 1965. 272 с.
[3] И.Е. Плещинская, Н.Б. Плещинский, Вестник Казан. технол. ун-та, №11, 147-153 (2010).
[4] И.Л. Александрова, И.Е. Плещинская, Н.Б. Плещинский, Вестник Казан. технол. ун-та, №7, 37-39 (2012).
[5] I.L. Aleksandrova, E.A. Osipov, N.B. Pleshchinskii, P.A. Rogozhin, Proceedings of PIeRs 2012 in Moscow, 435-439 (2012).
[6] Е.Г. Лапина, Тр. Матем. центра им. Н.И. Лобачевского, 2. Казанск. матем. об-во, Казань, 1999. С. 229-234.
Рис. 1 - Зависимость потока энергии от
параметров и
© И. Е. Плещинская - канд. физ.-мат. наук, доц. каф. информатики и прикладной математики КНИТУ, [email protected]; Н. Б. Плещинский - д-р физ.-мат. наук, проф., зав. каф. прикладной математики К(П)ФУ, [email protected]; И. В. Сабиров - студ. К(П)ФУ.
