22. Бабич В.М., Булдырев В.С., Молотков И.А. Пространственновременной лучевой метод. Линейные и нелинейные волны. -Л.: Изд-во ЛГУ, 1985. - 272 с.
23. Молотков И.А., Вакуленко С.А. Сосредоточенные нелинейные волны. - Л.: Изд-во ЛГУ, 1988. - 240 с.
24. Вакуленко С.А., Маслов В.П., Молотков И.А., Шафаре-вич А.И. Асимптотические решения уравнения Хартри, сосредоточенные при в малой окрестности кривой // Доклады РАН. - 1995. - Т 345. - № 6. - С. 743-745.
25. Во Хань Фук, Четвериков В.М. Обобщенные солитоны уравнения Шредингера с унитарной нелинейностью // Теоретическая и математическая физика. - 1978. - Т. 36. - № 3. - С. 345-351.
26. Sanchez A., Bishop A.R. Collective coordinates and length-scale competition in spatially inhomogenious soliton-bearing equations // SIAM Rev. - 1998. - V. 40. - № 3. - P. 579-615.
27. Mertens F.G., Schnitzer H.J., Bishop A.R. Hierarchy of equations of motion for nonlinear coherent excitations applied to magnetic vortices // Phys. Rev. B. - 1997. - V. 56. - № 5. - P. 2510-2520.
28. Rice M.J. Physical dynamics of solitons // Phys. Rev. B. - 1983. -V. 28. - № 6. - P. 3587-3589.
29. Rice M.J., Mele E.J. Phenomenological theory of soliton formation in lightly-doped polyacetylene // Solid State Commun. - 1980. -V. 35. - № 6. - P. 487-491.
30. McLaughling D.W., Scott A.C. Perturbation analysis of fluxon dynamics // Phys. Rev. A. - 1978. - V. 18. - № 4. - P. 1652-1680.
31. Malomed B.A. Perturbative analysis ofthe interaction of a phi4 kink with inhomogeneities // J. Phys. A: Math. Gen. - 1992. - V. 25. -№ 4. - P. 755-764.
32. Kivshar Y.S., Fei Z., Vazquez L. Resonant soliton-impurity interactions // Phys. Rev. Lett. - 1991. - V. 67. - № 10. - P. 1177-1180.
33. Fei Z., Kivshar Y.S., Vazquez L. Resonant kink-impurity interactions in the sine-Gordon model // Phys. Rev. A. - 1992. - V. 45. -№ 8. - P. 6019-6030.
34. Белов В.В., Доброхотов С.Ю., Тудоровский Т.Я. Асимптотические решения нерелятивистских уравнений квантовой механики в искривленных нанотрубках: I. Редукция к пространственно одномерным уравнениям // Теоретическая и математическая физика. - 2004. - Т. 141. - № 2. - С. 267-303.
35. Олвер П. Приложения групп Ли к дифференциальным уравнениям. - M.: Наука, 1989. - 639 с.
36. Fushchych W.I., Shtelen W.M., Serov N.I. Symmetry analysis and exact solutions of equations of nonlinear mathematical physics. -Dordrecht: Kluwer, 1993. - 435 p.
37. Fushchich W.I., Nikitin A.G. Symmetries of equations of quantum mechanics. - N.Y.: Allerton Press Inc., 1994. - 480 p.
38. Shvedov O.Yu. Semiclasical symmetries // Ann. Phys. - 2002. -V. 296. - P. 51-89.
Поступила 11.02.2009 г.
УДК 537.874.4
ИССЛЕДОВАНИЕ ВЛИЯНИЯ ТОНКИХ ПРОВОДНИКОВ НА БИСТАТИЧЕСКИЕ СЕЧЕНИЯ РАССЕЯНИЯ ДИЭЛЕКТРИЧЕСКОГО ЭЛЛИПСОИДА
Ю.А. Келлер
Томский государственный университет E-mail: [email protected]
На основе метода вспомогательных источников получено решение задачи рассеяния электромагнитных волн на трехмерном магнитодиэлектрическом теле при наличии вблизи него тонких проводников, расположенных произвольным образом в пространстве относительно тела. Построенный алгоритм реализован в виде компьютерной программы для расчета характеристик рассеяния ряда структур, отличающихся взаимным расположением тел, входящих в них. Исследовано влияние тонких проводников на бистатические сечения рассеяния диэлектрического эллипсоида.
Ключевые слова:
Метод вспомогательных источников, диэлектрическое тело, тонкие проводники, математическое моделирование электромагнитного рассеяния.
Задачи рассеяния электромагнитного поля на структурах, состоящих из диэлектрического тела и тонких проводников, возникают в различных областях науки и техники, например, в антенной технике и радиолокации. Тонкие проводники часто используют в качестве передающих и приёмных антенн. При расположении таких антенн вблизи диэлектрических тел возникает проблема оценки влияния диэлектрических тел на параметры антенны, решение которой требует решения поставленной задачи рассеяния. В радиолокации при оценке радиолокационной заметности сложного объекта часто возникает ситуация, когда часть объекта - это
диэлектрическое тело с расположенными вблизи него тонкими проводниками. Расчет бистатическо-го сечения рассеяния (БСР) такой части объекта также требует решения поставленной задачи. Если расстояние между диэлектрическим телом и проводниками меньше или сравнимо с длиной волны, то корректная постановка исследований подобного рода приводит к необходимости решения граничных задач теории рассеяния с учетом электромагнитного взаимодействия между рассеивателями. Существующие численные методы [1] позволяют решать подобные задачи. Применительно к задачам, рассматриваемым в данной статье, из наиболее
популярных методов предпочтительным является метод интегральных уравнений [2, 3]. Однако соответствующие вычислительные алгоритмы получаются чрезвычайно емкими по затратам компьютерной памяти и времени вычислений.
В последние годы для решения задач электромагнитного рассеяния на структурах, состоящих из взаимодействующих тел, начали использовать метод вспомогательных источников [3, 4]. В частности, в [3] предложен вариант этого метода для численного решения задач электромагнитного рассеяния, составленных из конечного числа трехмерных идеально проводящих тел.
Данная статья является продолжением серии работ, посвященных применению метода вспомогательных источников к анализу рассеяния на структурах из нескольких тел. В ней предложен вариант метода вспомогательных источников для численного решения задач электромагнитного рассеяния на структурах, состоящих из диэлектрического эллипсоида и расположенных вблизи него тонких проводников конечной длины. Дана математическая формулировка метода и краткое описание возможностей реализованной на его основе программы для расчёта компонент рассеянного поля. Приведены примеры численных результатов, характеризующих влияние тонких проводников на БСР диэлектрического эллипсоида.
1. Формулировка задачи и метод её решения
Геометрия задачи показана на рис. 1. Будем рассматривать стационарную^задачу рассеяния электромагнитного поля {Е 0,Н 0} на структуре, состоящей из объемного диэлектрического тела D¡, ограниченного поверхностью S, с диэлектрической и магнитной д проницаемостями и П тонких проводников, ограниченных поверхностями Su (и=1,2,...,П) и расположенных произвольным образом по отношению к телу D¡ (зависимость рассеяния от времени выбрана в виде ехр(-ш/)). Под объемным телом будем понимать тело, максимальный и минимальный поперечный размеры которого сравнимы между собой, а под тонким проводником - идеальный проводник круглого сечения, диаметр которого мал по сравнению с длиной проводника и длиной волны. Эта структура размещена в однородной безграничной среде De с диэлектрической и магнитной проницаемостями ее, д в декартовой системе координат с центром, выбранным внутри диэлектрического тела. Требуется найти рассеянное поле {Ее,Не}.
Математическая постановка задачи имеет следующий вид:
Ух Ее = тдНе Ух Не = -тЕ,Е„
Ух Е = гюдЦ Ух Н = -те,Е
(1)
Рис. 1. Геометрия задачи
п х (Е, - Ее) = п х Е0 п х (Н, - Не) = п х Н0
п., х Ее = -п, х Е0
на и = 1,2,..., и,
(2)
{&Е. -»¡ДН,} х Я / Я + {^ДНе; -у/^Е } =
= 0(Я _1), Я (3)
где Ее,Не, и Е,Ц - поля в областях De и D¡, п - единичный вектор нормали к поверхности S, пи - единичные векторы нормалей к поверхностям тонких проводников; Л=(х2+/+г2)1/2 - расстояние от рассматриваемой точки до центра системы координат, а хЬ - векторное произведение.
Решение сформулированной выше задачи получим следующим образом. Построим модель рассеянного поля. Аналогично тому, как это сделано в работе [5], введем (рис. 1) две вспомогательные поверхности S¡ и Se, подобные поверхности диэлектрического тела 8 в смысле гомотетии с центром в точке О, расположенной внутри тела и являющейся началом системы координат. Если поверхность S является центральной, то центр гомотетии выбираем так, чтобы он совпадал с центром поверхности. Поверхность 8е=Ке8 расположена внутри диэлектрического тела и характеризуется коэффициентом гомотетии (подобия) Ке, меньшим единицы; поверхность S¡=K¡S расположена вне тела и характеризуется коэффициентом подобия К, большим единицы. Если К=К=1, вспомогательные поверхности 8е и S¡ совпадают с S.
Выберем на вспомогательной поверхности Se конечную совокупность точек {М„,е}^= 1, в каждой из которых разместим пару независимых элементарных электрических диполей с моментами Рй=Рй ер 5е и р 2е=р2ер 2е, а на вспомогательной поверхности S¡ - конечную совокупность точек {М„,,)* 1, в каждой из которых разместим пару независимых вспомогательных электрических диполей с моментами р 2¡=р2¡е и р'2¡=р2¡е 2¡. Единичные векторы е6 п1е, е6 п2 ¡ выбраны в плоскости, касательной к Se в точке Ып е, а единичные векторы е ^¡, е 2‘
- в плоскости, касательной к S¡ в точке Мй>¡. Предполагается, что диполи, размещенные на Se, излучают в однородную среду с параметрами ее, д, а диполи, размещенные на S¡, - в однородную среду с параметрами е;, д.
Внутри каждого из тонких проводников на оси разместим неп6рерывно распределенный вспомогательный ток 1и.
р Представим неизвестное рассеянное поле {Ее,Не} в De в виде суммы полей вспомогательных диполей, расположенных на вспомогательной поверхности Se, и вспомогательных токов:
Е.(М) = ф{±Ух(УхЙи,) + ^Ух(Ух 11.)},
ке п=1 и= 1
1 N и
Н, (М ) = —{¿УхЙи е +£УхЦ, },
Де п=1 и=1
(4)
пя,е =т.(м,мпе)рпе, г/ = рп^ + рх, п. = |т, (М,м,и)ЗЛ М е ц,
1и
а поле Е,Н, в D¡ - в виде суммы полей вспомогательных диполей, расположенных на вспомогательной поверхности S¡:
N
Iф
Е(М) = — ¿Ух (УхИи ,),
к, п=1
Н,(М) = -¿УхП,
Д п=1
II и, = т, (М, мп i) г/,
рп/ = p”1^e”lj + рх>гх*г , М е ц.
В представлениях (4), (5)
Т, (М, Мп,е) = ехр(1кеЯММе) / 4пЯммп Т, (М,М,и ) = ехр(1кеЯММы ) / 4пЯмм1 и
Т, (М,Мп,) = ехр(/кЯмм„) / 4пЯм ,
(5)
к,
= аф
Кмм„,е и К-ти - расстояния от точки М„е на вспомогательной поверхности Se и точки М,и на оси проводника с номером и до точки наблюдения М в De; К-ММп,1 - расстояние от точки М^ на вспомогательной поверхности 8 до точки М в D¡; р^, р 2е (2=1,2,...,Ие) и рТ{, р2 (n=1,2,...,N¡) - неизвестные дипольные моменты; Ие и - числа точек размещения д6иполей на вспомогательных поверхностях 8е и 1и (и=1,2,..., и) - неизвестные осевые вспо-
могательные токи; интегрирование в (4) проводится вдоль осей проводников 1и.
Представления (4), (5) удовлетворяют уравнениям Максвелла (1) и условиям излучения (3). Для того, чтобы удовлетворить граничным условиям (2), необходимо соответствующим образом выбрать значения дипольных моментов р2е, р 2е (2=1,2,...,Же) и рТ{,Рр2 (2=1,2,...,Л9 и распределения осевых токов 1и (и=1,2,...,П).
Введем кусочно-постоянную аппроксимац6ию осевых токов. Разобьем линию 1и каждого тока 1и на Ыи малых участков, в пределах каждого из которых ток м6ожно считать постоянным. Тогда выражение для Пи в (4) приближенно можно записать в виде
N.
I, и
п =у 3 я , Г Т (М, м, ы ,
и ¿ш^ и, и, I е у ’ I и 7 ,=1 1,-1 и
(6)
где I - величина тока на ¡-м участке проводника с номером и; - единичный вектор, направление
которого совпадает с направлением касательной в средней точке рассматриваемого участка. При таком подходе нахождение неизвестных распределений осевых токов сводится к нахождению значе-
ний
элементов тока.
и=1
Для определения величин дипольных моментов и элементов тока используем граничные условия (2), удовлетворяя им в соответствии с методом коллокаций. Пусть М (/=1,2,...,Ь) - точки коллокации на поверхности диэлектрического тела 5, а М/(/=1,2,...,Ьи) - точки коллокации на поверхности проводников 5/; Ь - число точек коллокации на 5, а Ьи - число точек коллокации на 5/. В силу предположения о малости диаметра проводника по сравнению с длиной проводника и длиной волны будем считать, что вкладом в рассеянное поле азимутальных составляющих токов на поверхностях тонких проводников можно пренебречь. Тогда для нахождения неизвестных рь р«/ («=1,2,...Д), рЦ, р2 (и=1,2,...Д) и (и=1,2,...Д г=1,2,...,Ии) получим сле-
дующую систему линейных алгебраических уравнений:
п1 х (Ё/ - Ё]е) = П х Ё0, П х (Щ - Щ ) = П х Н0,
1 = 1,2,..., Ь, (7)
Ё^, = -0=,, и = 1.2,..., и, ] =1,2,..., Ь,
где р/' - значение единичного вектора нормали в т очк6е М/ н а п оверхности диэлектрического тела; Е/, Н/ , и Е/, ИЦ, - значения компонент внешнего и внутреннего полей в точке М; Е0!, Н0!, - значения компонент возбуждающего поля в этой же точке; Е/ и>; и Е/и1 - значения составляющих рассеянного и возбуждающего полей вдоль оси проводника с номером и в точках коллокации на его поверхности.
Отметим, что в отличие от обычно используемых систем линейных алгебраических уравнений метода дискретных источников [6], элементы матриц которых выражены в конечном виде,
и \ и
4Ь + ^Ьи |х£N элементов матрицы системы
и=1 / и=1
(7) содержат интегралы, которые возникают при применении в соответствии с (4) необходимых дифференциальных операторов к функции Пи (это элементы матрицы, стоящие перед /„,,-).
Решение системы (7) определяется путем минимизации функционала
ф = ¿{| п х Ё -ЁЦ) -п х Ё112 +
]=1
+-^ | п1 х (Н1 -Н1) -п хЩ |2} +
и Ьи
+1111 + Ё^Т- (8)
и=1 1=1
После решения задачи минимизации (определения неизвестных дипольных моментов рТ, р 2е («=1,2,...,А/), рТ/, Рй' («=1,2,...,N1) и элементов тока 1^ (и=1,2,...,и, 1=1,2,...,Д,) необходимые характеристики рассеянного поля определяются из (4).
Контроль точности модели (4), (5) осуществляется путем вычисления относительного значения функционала (8) на сетке точек, промежуточных по отношению к точкам коллокации, выбираемых как на поверхности диэлектрического тела 5, так и на поверхностях £/ всех проводников, входящих в структуру:
Д = (Ф/ Ф 0)1/2,
Ф0 = X{| п хЁ012 + ^| п х 110 |2} + ии Ё^, |2,
1=1 и =1
где Ф' - значение функционала (8) на указанной выше совокупности точек; Ф0 - значение соответствующей нормы падающего поля на этой же совокупности точек; Ь - число промежуточных точек на поверхности диэлектрического тела; Ьи' - число промежуточных точек на поверхности проводника с номером и.
2. Численные результаты
На основании изложенной выше модели создана программа для расчета компонент рассеянного поля и контроля точности полученного решения. Диэлектрическое тело может быть трехосным эллипсоидом. Максимальное число тонких проводников в структуре, предусмотренное программой, равно трем. Предполагается, что все проводники являются прямолинейными. Положение проводников относительно диэлектрического тела, а также их длина определяются заданием координат начальной и конечной точек проводника. Ниже приведены некоторые результаты, касающиеся влияния тонких проводников на БСР диэлектрических тел.
Методика исследований заключалась в сравнении БСР структур, состоящих из диэлектрического тела и тонких проводников, с БСР одиночного диэлектрического тела. Сравнение проводилось для структур, содержащих один проводник, расположенный со стороны падения возбуждающей волны (рис. 2) и в «зоне тени» эллипсоида (рис. 3).
Рис. 2. Структура, состоящая из эллипсоида и тонкого проводника, расположенного на расстоянии Дz со стороны падающей волны
Ф=сош1 точек размещения диполей, так и в полусе-чениях, проведенных посередине между ними.
Рис. 3. Структура, состоящая из эллипсоида и тонкого проводника, расположенного на расстоянии Az в «зоне тени» эллипсоида
Рис. 4 относится к структуре, представленной на рис. 2, а рис. 5 - к структуре на рис. 3. Геометрические параметры эллипсоида следующие: кеа=кеЬ=3, кее=4; относительная диэлектрическая проницаемость эллипсоида е/ее выбрана равной 4; относительная магнитная проницаемость ц/ре выбрана равной 1. Длина проводника равна 0,91, радиус проводника г равен 0,021. Кривые 1 на этих рисунках - это БСР одиночного эллипсоида, кривые 2 - БСР структуры из эллипсоида и проводника, когда последний расположен на расстоянии Аг=0,011 от эллипсоида, кривые 3 - те же характеристики, когда проводник расположен на расстоянии Аг=0,11 от эллипсоида. БСР приведены в Е-плоскости - плоскости, в которой лежат векторы Е 0 и ке падающей волны.
При получении результатов, представленных на рис. 4, 5, число элементов разбиения осевого тока проводника N выбрано равным 35, число точек коллокации на поверхности проводника - равным 140. Параметры подобия вспомогательных поверхностей выбраны следующими: Ке=0,6, К=4, количества точек размещения диполей на внутренней и внешней вспомогательных поверхностях выбраны одинаковыми: N=N¡=484; распределены эти точки следующим образом. В каждом из 22 полусечений Ф=сош1, отстоящих одно от другого на угловое расстояние Аф=16,36°, равномерно по углу в выбраны 22 точки размещения диполей. Число точек колло-кации Ь на поверхности эллипсоида выбрано равным 968; алгоритм их размещения по углу в выбран таким же, как для точек размещения диполей, но располагаются эти точки как в полусечениях
Рис. 4. БСР в E-плоскости эллипсоида с параметрами кеа=кеЬ=3, kec=4, е/ее=4, ц,/у.е=1 и структуры, состоящей из такого же эллипсоида и расположенного со стороныы падения возбуждающей волныы проводника длиной ¡=0,91. Кривая: 1) одиночный эллипсоид, 2) эллипсоид и проводник, А1=0,011,3) эллипсоид и проводник, Az=0/11
Рис. 5. БСР в Е-плоскости эллипсоида с параметрами кеа=кеЬ=3, кес=4, е/ее=4, ц/у.е=1 и структурыы, состоящей из такого же эллипсоида и расположенного в «зоне тени» проводника длиной ¡=0,91. Кривая: 1) одиночный эллипсоид, 2) эллипсоид и проводник, А1=0,011,3) эллипсоид и проводник, А1=0,11
Показано, что для эллипсоида с диэлектрической проницаемостью е/е=4 присутствие тонкого проводника в наименьшей степени влияет на сечение рассеяния в направлениях 70°<в<100°, прилегающих к направлению бокового рассеяния в=90°. Проводник, расположенный в «зоне тени» эллипсоида, влияет на БСР диэлектрического тела в меньшей степени, чем такой же проводник, расположенный на таком же расстоянии со стороны падения возбуждающей волны. Это объясняется тем, что в последнем случае токи, наводимые на проводнике падающей волной, больше по величине, чем токи, наводимые на проводнике ближним рассеянным полем эллипсоида, когда проводник находится в «зоне тени». Как показывает сравнение кривых 1-3 на рис. 4, 5 между собой, величина сечения рассеяния существенно зависит от расстояния между по-
верхностью эллипсоида и проводником. Например (см. рис. 4), при расположении проводника на расстоянии 0,011 от поверхности эллипсоида сечение рассеяния в направлении в=140° равно -9 дБ, а при расположении этого же проводника на расстоянии 0,11равно -19 дБ (уменьшилось на 10 дБ).
Выводы
На основе метода вспомогательных источников получено решение задачи рассеяния электромаг-
нитных волн на трехмерном магнитодиэлектрическом теле при наличии вблизи него тонкого проводника, расположенного в «зоне тени» и со стороны падения возбуждающего поля. Построенный алгоритм реализован в виде компьютерной программы для расчета характеристик рассеяния структур, отличающихся взаимным расположением входящих в них тел. Исследовано и показано существенное влияние тонких проводников на би-статические сечения рассеяния диэлектрического эллипсоида.
СПИСОК ЛИТЕРАТУРЫ
1. Kahnert F.M. Numerical methods in electromagnetic scattering theory // Journal of Quantitative Spectroscopy & Radiative Transfer. -2003. - V. 79-80. - P. 775-824.
2. Вычислительные методы в электродинамике // Под ред. Р. Митры. Пер. с англ. - М.: Мир, 1977. - 485 с.
3. Хижняк Н.А. Интегральные уравнения макроскопической электродинамики. - Киев: Наукова думка, 1986. - 278 с.
4. Еремин Ю.А., Орлов Н.В. Анализ рассеяния волн на нескольких магнитодиэлектрических телах методом дискретных источников // Радиотехника и электроника. - 1994. - Т. 39. -№ 5. - С. 740-748.
5. Дмитренко А.Г., Колчин В.А. Численный метод анализа электромагнитного рассеяния структурами из конечного числа трехмерных идеально проводящих тел // Радиотехника и электроника. - 2001. - Т. 46. - № 3. - С. 277-282.
6. Еремин Ю.А., Свешников А.Г. Метод дискретных источников в задачах электромагнитной дифракции. - М.: Изд-во МГУ, 1992. - 181 с.
Поступила 16.01.2009 г.
УДК 53.072;531.43;539.62
ФАКТОРЫ, ОПРЕДЕЛЯЮЩИЕ ГЕНЕРАЦИЮ УПРУГИХ ВОЛН ПРИ ТРЕНИИ. МОДЕЛИРОВАНИЕ НА ОСНОВЕ ДИСКРЕТНО-КОНТИНУАЛЬНОГО ПОДХОДА
А.Ю. Смолин1-3, С.А. Добрынин1, С.Г. Псахье1-3
’Институт физики прочности и материаловедения СО РАН, г. Томск 2Томский политехнический университет, г. Томск 3Томский государственный университет, г. Томск E-mail: [email protected]
На основе дискретно-континуального подхода рассмотрена модель механического взаимодействия в пятне контакта при трении скольжения. Проанализирован частотный спектр возникающих упругих волн, показано наличие частот, зависящих от шероховатости профиля взаимодействующих поверхностей. Анализ, выполненный с помощью Фурье- и вейвлет-преобразований, позволил выявить сложную структуру возникающих при трении упругих волн. Сделан вывод о том, что закономерности процесса изнашивания могут быть изучены на основе анализа соответствующих акустических спектров.
Ключевые слова:
Трение, моделирование, метод частиц, упругие волны, Фурье-анализ, вейвлет-анализ, сигнал.
Введение
Одна из наиболее острых современных проблем
- износ машин и механизмов - связана с трением. Трудности в изучении процессов трения и изнашивания во многом обусловлены тем, что они представляют собой сложное существенно динамическое явление. Так, непрерывное измерение коэффициента трения скольжения показывает, что эта характеристика пары не является постоянной величиной даже при установившемся режиме, а периодически изменяется, достигая максимумов и
минимумов [1, 2]. Проявлением динамической сущности трения является также возникновение акустических колебаний, которые играют важную роль в поведении взаимодействующих тел, в том числе в их изнашивании. Срок службы деталей узлов трения можно существенно повысить только за счет гашения этих колебаний демпферами, нанесением поверхностных слоев или за счет подбора геометрии деталей таким образом, чтобы возникающие при трении колебания опасной частоты быстро затухали [3].