PRIKASPIYSKIY ZHURNAL: Upravlenie i Vysokie Tekhnologii (CASPIAN JOURNAL: Management and High Technologies), 2014, 2 (26) THEORETICAL BASES, DEVELOPMENT AND APPLICATION OF MATHEMATICAL MODELS, ALGORITHMS, COMPUTER PROGRAMS
4. Eder J., Olivotto G. E., Graber W. A Data Warehouse for Event Logs. Proceedings of First International Conference on Engineering and Deployment of Cooperative Information Systems (EDCIS 2002). Lecture Notes in Computer Science. China, Beijing, 2002, vol. 2480, pp. 1-15.
5. La Rosa M., Dumas M., ter Hofstede A., Mendling J. Congurable Multi-Perspective Business Process Models. Information Systems, 2011, no. 36 (2), pp. 313-340.
6. Mamaliga. T. Realizing a Process Cube Allowing for the Comparison of Event Data. Master's thesis. Eindhoven, Eindhoven University of Technology, 2013.
7. Ribeiro J. T. S., Weijters A. J. M. M. Event Cube: Another Perspective on Business Processes. Proceedings of 19th International Conference on Cooperative Information Systems. Lecture Notes in Computer Science. Berlin, Springer-Verlag, 2011, vol. 7044, pp. 274-283.
8. van der Aalst W. M. P. Process Cubes: Slicing, Dicing, Rolling Up and Drilling Down Event Data for Process Mining. Proceedings of Asia Pacific Conference on Business Process Management (AP-BPM 2013). Lecture Notes in Business Information Processing, Berlin, Springer-Verlag, 2013, vol. 159, pp. 1-22.
9. van der Aalst W. M. P. Process Mining: Discovery, Conformance and Enhancement of Business Processes. Springer, 2010. 349 p.
10. XES. Available at: http://www.xes-standard.org (accessed 1 December 2013).
УДК 621.391.677: 519.711.3
ПОСТРОЕНИЕ ГЕОМЕТРИЧЕСКОЙ МОДЕЛИ ИЗЛУЧЕНИЯ ЗЕРКАЛЬНОЙ ПАРАБОЛИЧЕСКОЙ АНТЕННЫ
Ширшов Максим Владимирович, аспирант, Пензенский государственный университет, 440026, Российская Федерация, г. Пенза, ул. Красная, 40, e-mail: admin@r4f.su
Юрков Николай Кондратьевич, доктор технических наук, Пензенский государственный университет, 440026, Российская Федерация, г. Пенза, ул. Красная, 40, e-mail: yurkov nk@mail.ru
Якимов Александр Николаевич, доктор технических наук, Пензенский государственный университет, 440026, Российская Федерация, г. Пенза, ул. Красная, 40, e-mail: admin@r4f.su
Для исследования излучения зеркальной параболической антенны (ЗПА) методом математического моделирования необходимо построить её геометрическую модель, заменяющую излучающую поверхность антенны многогранником. При этом поле излучения модели антенны представляется как суперпозиция полей, формируемых токами, протекающими по граням модели с учетом амплитуд и фаз волн, распространяющихся от этих граней. Отличие многогранника от параболоида приводит к возникновению фазовых ошибок, минимизация которых является важной задачей. Авторами предложен алгоритм построения геометрической модели излучающей поверхности ЗПА, основанный на методе триангуляции Делоне, позволяющем произвести дискретизацию излучающей поверхности антенны в виде совокупности равносторонних плоских треугольных конечных элементов; получить совокупности узловых точек, образующих эти конечные элементы и, соответственно, координаты узлов. Этих данных оказывается достаточно для проведения расчета как полей, формируемых конечными элементами излучающей поверхности антенны, так и для нахождения их суперпозиции с учетом пространственной ориентации этих элементов и векторного характера электромагнитного поля. Дана оценка погрешности представленной модели. Предложенный алгоритм, основанный на методе триангуляции Делоне, обеспечивает выполнение требования к максимально допустимой фазовой ошибке при минимальных вычислительных затратах. Это позволяет использовать модель не только при расчете характеристик излучения ЗПА, но и для оценки влияния внешних воздействий на конструкцию антенны конечно-элементным методом.
Ключевые слова: излучающая поверхность, модель, антенна, триангуляция, конечный элемент, пространственная ориентация, суперпозиция, электромагнитное поле, фазовая ошибка, погрешность
ПРИКАСПИЙСКИЙ ЖУРНАЛ: управление и высокие технологии № 2 (26) 2014 ТЕОРЕТИЧЕСКИЕ ОСНОВЫ, РАЗРАБОТКА И ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, АЛГОРИТМОВ, ПРОГРАММ ДЛЯ ЭВМ
BUILDING GEOMETRICAL MODEL OF RADIATION OF THE MIRROR PARABOLIC ANTENNA
Shirshov Maksim K, post-graduate student, Penza State University, 40 Krasnaya St., Penza, 440026, Russian Federation, e-mail: admin@r4f.su
Yurkov Nikolay K., D.Sc. (Engineering), Penza State University, 40 Krasnaya St., Penza, 440026, Russian Federation, e-mail: yurkov_nk@mail.ru
Yakimov Aleksandr N., D.Sc. (Engineering), Penza State University, 40 Krasnaya St., Penza, 440026, Russian Federation, e-mail: admin@r4f.su
To investigate the radiation antenna parabolic mirror (PAD) mathematical modeling, it is necessary to construct a geometric model, replacing the antenna radiating surface polyhedron. The field radiation pattern of the antenna is represented as a superposition of fields generated by currents flowing along the edges of the model taking into account the amplitudes and phases of these faces propagating waves. Unlike polyhedron from parabolic leads to phase errors, minimization of which is an important task. The authors propose an algorithm for constructing a geometric model of the radiating surface of PAD based on Delaunay triangu-lation method, which made possible the sampling of the radiating surface of the antenna, get together nodal points forming the finite elements of the surface in the form of equilateral triangles and flat, respectively, the coordinates of these points. This data is sufficient to perform the calculation as fields formed end surface of the antenna radiating elements, and to locate them in view of the spatial superposition of these elements and the orientation of the vector nature of the electromagnetic field. The estimation error of the presented model. The proposed algorithm based on Delaunay triangulation method, ensures that the requirements for the maximum allowable phase error with minimal computational cost. This enables not only to use the model for calculating specula radiation characteristics of a PAD, but also to evaluate the effect of external influences on the antenna structure finite element method.
Keywords: radiation surface, model, antenna, triangulation, finite element, spatial orientation, superposition, electromagnetic field, phase error, error
Введение. Рассматриваемая зеркальная параболическая антенна (ЗПА) состоит из параболического зеркала и закрепленного в его фокусе рупорного облучателя. Физическое моделирование такой антенны для достижения наилучших электрических характеристик сопряжено с большими материальными и временными затратами. Поэтому перспективным оказывается использование математического моделирования излучения такой антенны.
Построение геометрической модели зеркала антенны сводится к дискретному представлению излучающей поверхности. В связи с необходимостью обеспечения адекватности математической модели линейные размеры дискретных элементов должны быть сравнимы с длиной излучаемых электромагнитных волн [13].
В связи с этим в сверхвысоком частотном диапазоне волн при формировании геометрической модели зеркала ЗПА требуются десятки тысяч дискретных элементов её излучающей поверхности, для каждого из которых необходимо решать интегральную задачу излучения с учетом амплитуд и фаз возбуждающих источников, а также векторного характера формируемого поля.
Из-за большого объема требуемых при этом вычислений целесообразно представление излучающей поверхности в виде многогранника с гранями в виде плоских элементов первого порядка. Отличие этого многогранника от параболоида приводит к возникновению фазовых сдвигов полей, создаваемых отдельными гранями модели излучающей поверхности антенны. При этом результирующее поле излучения антенны является суперпозицией полей, формируемых отдельными гранями с учетом их амплитуд и фаз [2, 5].
PRIKASPIYSKIY ZHURNAL: Upravlenie i Vysokie Tekhnologii (CASPIAN JOURNAL: Management and High Technologies), 2014, 2 (26) THEORETICAL BASES, DEVELOPMENT AND APPLICATION
OF MATHEMATICAL MODELS, ALGORITHMS, COMPUTER PROGRAMS_
При построении геометрической модели излучающей поверхности параболической антенны перспективным оказывается использование метода триангуляции Делоне, обеспечивающего деление этой поверхности на равносторонние плоские треугольные конечные элементы [1, 8].
Целью данной статьи является исследование возможностей метода триангуляции Делоне для построения геометрической модели зеркала ЗПА в виде многогранника с плоскими треугольными гранями и максимально возможным шагом дискретизации при заданной предельно допустимой фазовой ошибке.
Постановка задачи. Алгоритм построения геометрической модели параболической антенны состоит из трёх основных шагов:
• разбиение параболической поверхности на кольцевые зоны;
• триангуляция кольцевых зон;
• объединения элементов дискретизации кольцевых зон в единую модель. Рассмотрим алгоритм построения геометрической модели ЗПА (рис. 1), излучающая
поверхность которой описывается формулой [8] :
z = (х2 + у2)/4/,
где X,y,Z - пространственные координаты, a f - фокусное расстояние. В результате построения геометрической модели излучающая поверхность разбивается на равносторонние или близкие к равносторонним треугольники. В силу криволинейности параболического профиля, эта равносторонность достигается лишь внутри кольцевых зон, определяемых равномерным шагом разбиения параболического профиля в одном из его главных сечений.
Требования к шагу дискретизации обусловлены тем, что в ЗПА максимальна фазовая ошибка, вызванная отклонениями профиля зеркала от параболоида, не должна превышать А/16 [5]. Для выполнения этого требования введем ряд ограничений:
• линейный размер шага дискретизации не должен превышать длину волны, так как проведенные растровым методом исследования показали, что при большем чем X шаге дискретизации диаграмма направленности модели антенны значительно отклоняется от рассчитанной для гладкого параболоида [13];
• размеры сторон треугольных элементов дискретизации должны быть равны между собой или близки, так как есть риск вырождения решения [10].
ПРИКАСПИЙСКИЙ ЖУРНАЛ: управление и высокие технологии № 2 (26) 2014 ТЕОРЕТИЧЕСКИЕ ОСНОВЫ, РАЗРАБОТКА И ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, АЛГОРИТМОВ, ПРОГРАММ ДЛЯ ЭВМ
При растровом методе излучающая поверхность зеркала, описываемая формулой параболоида, в программных пакетах или с помощью специально разработанных программ разбивается на прямоугольные дискретные элементы взаимно перпендикулярными плоскостями, параллельными осям Ох и Оу декартовой системе координат с шагом
АЬ = 0,5Х . Треугольные элементы дискретизации получаются в результате их последующего деления пополам [13].
Опыт проведения вычислительных экспериментов показал, что при отклонениях размеров сторон треугольников до 41 % решение не вырождается.
Построение модели. Для обеспечения разбиения излучающей поверхности антенны на треугольники выбираем уровни сечения параболоида исходя из рекомендованного для растровой дискретизации шага АЬ, длина которого составляет [2]:
АЬ = 0,5-А,,
где X - длина волны электромагнитных колебаний, формируемых зеркалом антенны.
В векторной интерпретации эта кривая сечения параболоида представляет собой годограф векторной функции Г скалярных аргументов X , у , и 2 . Учитывая, что параболическое зеркало относится к осесимметричным излучателям, целесообразно осуществлять равномерное разбиение относительно центра параболы, совмещенного с центром декартовой системы координат (рис. 2). Интервал равномерного разбиения функции при этом определяется как разность Аг радиусов-векторов узловых точек при равномерной дискретизации этой кривой
АГ:
4-1 =
где к = 2, ... , К ; К - максимальный порядковый номер индекса радиуса-вектора узловой точки сечения параболоида в полупространстве 2(Эх с положительной координатой X . При этом индекс к — 1 = 0 соответствует координатам X = 0, 2 = 0.
Рис. 2. Дискретизация параболического профиля в сечении х()г
Получим равномерную линейную дискретизацию годографа векторной функции, описывающей параболический профиль (см. рис. 2), в декартовой системе координат. Тогда шаг дискретизации Л.Т/ будет равным [5]
PRIKASPIYSKIY ZHURNAL: Upravlenie i Vysokie Tekhnologii (CASPIAN JOURNAL: Management and High Technologies), 2014, 2 (26) THEORETICAL BASES, DEVELOPMENT AND APPLICATION
OF MATHEMATICAL MODELS, ALGORITHMS, COMPUTER PROGRAMS_
AL = |Ar| = - ? + (zk - zk_xf . (l)
Так как в плоскости xOz сечение параболоида описывается уравнением:
z = x2/(4-f), (2)
то Zk=x\l{\.f)-,Zk_x=xl_J{\.f).
Используем полученное значение Х^ в плоскости координат хОу как радиус окружности при построении многогранника [3]:
х = хк- cos(a), у = хк • sin (а), (з)
где ОС - угол относительно оси Ох, образованный вращением радиуса-вектора Г вокруг оси Oz, необходимый для определения координат узловых точек граней многогранника:
a = 360°/(4-2¿). (4)
В результате получаем массив координат узловых точек, определяющих геометрию параболоида и составляющих основу сетки с узловыми точками, размещенными на его поверхности. Далее, используя алгоритм, соединяющий треугольники в единую поверхность, производим «сшивание» этой поверхности (рис. 3).
Х,М -,
-0,6
Рис. 3. Результат дискретизации излучающей поверхности антенны в декартовой системе координат
ПРИКАСПИЙСКИЙ ЖУРНАЛ: управление и высокие технологии № 2 (26) 2014 ТЕОРЕТИЧЕСКИЕ ОСНОВЫ, РАЗРАБОТКА И ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, АЛГОРИТМОВ, ПРОГРАММ ДЛЯ ЭВМ
Таким образом, предложенный в данной работе алгоритм позволяет произвести дискретизацию излучающей поверхности параболической антенны, получить совокупности узловых точек, образующих конечные элементы этой поверхности в виде равносторонних треугольников и, соответственно, координаты этих точек. Этих данных оказывается достаточно для проведения расчета как полей, формируемых поверхностными токами конечных элементов излучающей поверхности антенны, так и их суперпозиции с учетом пространственной ориентации этих элементов и векторного характера электромагнитного поля.
Анализ точности модели. Оценим точность построения геометрической модели в рамках теории геометрической оптики по отклонению треугольных граней от параболоида в их центроидах по ходу парциального луча излучения ЗПА.
Центроид м , являющийся центром излучающего элемента дискретизации антенны, определим по пересечению медиан АА' и ВВ' (рис. 4) [7].
С
Точки А ж В определяются полусуммами координат вершин треугольника ABC:
хл = (хс + хв )/2, уА. = (ус + ув )/2, zÁ = (zc + zs )/2, = (Хс + Ха)/2 , yg = (ус + уА)!2, zff = (zc + za)/2 .
Координаты центроида м для каждого дискретного элемента определяются путем совместного решения уравнений, описывающих прямые АА' и ВВ', исходя из координат точек А, А' и В, В .
(Хм - Хл )/(хА - XА) = (ум - уА)/(уА, -уА) = (zM - zA)l(zA, -zA),
СХм - Хв )КХв' ~Хв) = (Ум- У в ЖУв' - У в) = (Zм~ Zb )KZb> ~Zb)-
В результате решения системы уравнений (5) получим координаты ум , Хм , ZM точки м.
у _ (хв ~ха)'(Уа' ~Уа)'(Ув' -Ув) + Уа -(ХА' ~ха)-(Ув' -Ув) (xAf - ХА ) • (ув, - Ув) - (хв, - хв) • (уА, — у А) __У в ' (ХВ' ~Хв)' (у А' ~Уа)_.
(ХА' ~ ХА ) ' (у В' ~ Ув)~ (ХВ' ~ ХВ ) ' (у А' ~У А)
PRIKASPIYSKIY ZHURNAL: Upravlenie i Vysokie Tekhnologii (CASPIAN JOURNAL: Management and High Technologies), 2014, 2 (26) THEORETICAL BASES, DEVELOPMENT AND APPLICATION
OF MATHEMATICAL MODELS, ALGORITHMS, COMPUTER PROGRAMS_
=(УМ-УА)- (XA' - ХА)/(УА' ~УА) + ХА>
ZM =(XM~Xa)- iZA' ~Z A) I (XA' - X A)+ Z A-Погрешность модели определяется разностью хода парциального луча при отражении от дискретного элемента модели и параболического зеркала. Так как отраженные лучи направлены параллельно оси Oz [3], то погрешность определяется суммой отрезков MS и ES (рис. 5).
Точка S на поверхности параболического зеркала является результатом пересечения прямой, проходящей через точки F и М, с сечением параболоида в плоскости xOz . Координаты этой точки могут быть определены из системы уравнений:
<Х - )КХМ - XF ) = (Л - Ук ЖУМ - УР ) = (ZS - ZF )KZM - ZF );
Z, = /(4 • ) '
где XF, yF , Zp - координаты точки F (фокусное расстояние).
В результате решения системы уравнений (6) получаем точку S отражения парциального луча от поверхности параболического зеркала антенны с координатами Xs, ys , zs с использованием следующих формул:
^ = 0,5 -(zM-F))/zlM-,
xs = (д/хм -F + n2 + (Q/z1M )) / zlM;
У8=Ум-(х81хм)> zs=xlK4-F).
ПРИКАСПИЙСКИЙ ЖУРНАЛ: управление и высокие технологии N2 2 (26) 2014 ТЕОРЕТИЧЕСКИЕ ОСНОВЫ, РАЗРАБОТКА И ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, АЛГОРИТМОВ, ПРОГРАММ ДЛЯ ЭВМ
Так как парциальный луч отражается от параболического зеркала параллельно оси ()~, то координаты точки Е (см. рис. 5) можно определить следующим образом. Точка Е лежит в одной плоскости с точкой параллельной плоскости х()~ , а также лежит в одной плоскости с точкой М параллельной плоскости хОу, следовательно, координаты точки Е являются проекциями точек Е, А4 и составляют:
ХЕ = ЛЛ ; У'Е = У Б 11 2Е = гм ■ Дополнительный фазовый набег электромагнитной волны Ар в центроидах элементов дискретизации определяется суммой отрезков А/Л' и Ж:
АР = ДРА£? + АРЖ > (?)
где Арж, Ар
Ш ~ отРезки прямых, описывающие дополнительный фазовый набег
АРд^ = >/(% ~хбУ + (Ум - У б У + (-М - )2 ; АРж =^(хЕ +(УЕ-УБУ + (2Е •
Из теории антенн [3] следует, что фазовая ошибка не должна превышать
АрИ1ЯЧ = VI6
Модельные расчёты погрешности, связанной с дискретизацией задачи, были произведены для ЗПА с диаметром 1 метр, фокусным расстоянием 0,3125 метра и электромагнитной волны длиной 3 см.
Результаты расчета дополнительного фазового набега электромагнитной волны Ар в центроидах элементов дискретизации при заданных условиях представлены на рис. 6.
V
/ / / \ 4 / \ * \
1 I , 1 / 3 / \ N V
\: 1 у *ч ' ■ ц , , * - п „ * \
■ ( и К — —^ " » - - Л \
0,05 0,1 0,1 5 0,2 0,25 0,3 0,35 0,4 0,45 Я, м
Рис. 6. Зависимость дополнительного фазового набега Ар в центроидах элементов дискретизации от расстояния /(
PRIKASPIYSKIY ZHURNAL: Upravlenie i Vysokie Tekhnologii (CASPIAN JOURNAL: Management and High Technologies), 2014, 2 (26) THEORETICAL BASES, DEVELOPMENT AND APPLICATION OF MATHEMATICAL MODELS, ALGORITHMS, COMPUTER PROGRAMS
/ 2 2
Здесь R = J Хт + ут - расстояние от оси Oz до центроидов элементов дискретизации в плоскости хОу, Ар определено согласно уравнению (7).
На рис. 6 кривые 1, 2, 3 и 4 показывают зависимости дополнительного фазового набега Ар в центроидах элементов дискретизации от расстояния R при шаге дискретизации AL , равном 0,5А,, X, 1,5А, и 2Х соответственно.
Для заданной длины волны X = 3 см максимально допустимая фазовая ошибка Apmax = Х/\6 ~ 0,0019 м. Из приведенных зависимостей (см. рис. 6) данному требованию
не удовлетворяет лишь кривая 4, соответствующая шагу дискретизации AI, = 2.Х .
Заключение. Таким образом, приемлемым шагом дискретизации можно считать AL = 1,5А, . Это менее жесткое требование к шагу дискретизации, чем при растровом методе [2]. Предложенный алгоритм построения геометрической модели излучающей поверхности зеркальной параболической антенны, основанный на методе триангуляции Делоне, позволяет использовать модель не только при расчете характеристик излучения зеркальной параболической антенны [2, 5], но и для оценки влияния внешних воздействий на конструкцию антенны конечно-элементным методом, если элементы дискретизации излучающей поверхности антенны являются гранями объемных элементов дискретизации конструкции зеркала антенны [4, 6, 9, 11, 12, 14, 15].
Список литературы
1.Балуков О. Н. Построение геометрической модели антенны с использованием триангуляции Делоне / О. Н. Балуков, А. Н. Якимов // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2009. -№ 1. - С. 109-117.
2. Воскресенский Д. И. Антенны и устройства СВЧ : учеб. для вузов / Д. И. Воскресенский, В. И. Гостюхин, В. М. Максимов, JT. И. Пономарев ; под ред. Д. И. Воскресенского. - Москва : Радиотехника, 2006. - 376 с.
3. Вуд П. Анализ и проектирование зеркальных антенн : пер. с англ. / П. Вуд. - Москва : Радио и связь, 1984. - 208 с.
4. Горбалысов М. С. Моделирование влияния тепловых воздействий на характеристики излучения зеркальной антенны / М. С. Горбалысов // Надежность и качество - 2013 : труды Международного симпозиума : в 2 т. / под ред. Н. К. Юркова. - Пенза : Изд-во ПГУ, 2013. - Т. 1 - С. 92-95.
5. Драбкин А. JT. Антенно-фидерные устройства / А. JT. Драбкин, В. JT. Зузенко, А. Г. Кислов. -Москва : Сов. радио, 1974. - 536 с.
6. Лунин JT. С. Особенности численного моделирования датчиков на основе поверхностных акустических волн (ПАВ) методом конечных элементов / JT. С. Лунин, В. Ф. Катаев, И. В. Крейс // Прикаспийский журнал: управление и высокие технологии. - 2011. - № 3. - С. 39-47.
7. Понарин я. П. Элементарная геометрия : в 2 т. / П. я. Понарин. - Москва : МЦНМО, 2004. -Т. 1: Планиметрия, преобразования плоскости. - 312 с.
8. Скворцов А. В. Триангуляция Делоне и её применение / А. В. Скворцов. - Томск : Изд-во Томск, ун-та, 2002. - 128 с.
9. Камаев В. А. Автоматизированный расчет вариантов электрической составляющей электромеханического привода зеркальных антенн / В. А. Камаев, А. Д. Чивилёв // Прикаспийский журнал: управление и высокие технологии. - 2014. - № 1. - С. 139-147.
10. Сабоннадьер Ж. К. Метод конечных элементов и САПР : пер. с фр. / Ж. К. Сабоннадьер, Ж. Л. Кулон. - Москва : Мир, 1989. - 190 с.
11. Старченко И. Б. Адаптивное управление гидроакустическими системами / И. Б. Старченко, И. А. Кириченко, В. Л. Сахаров, П. С. Голосов // Прикаспийский журнал: управление и высокие технологии. -2013. -№3. - С. 26-33.
ПРИКАСПИЙСКИЙ ЖУРНАЛ: управление и высокие технологии № 2 (26) 2014 ТЕОРЕТИЧЕСКИЕ ОСНОВЫ, РАЗРАБОТКА И ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, АЛГОРИТМОВ, ПРОГРАММ ДЛЯ ЭВМ
12. Шишулин Д. Н. Исследование влияния вибраций на параболическую антенну в ANSYS // Надежность и качество - 2013 : труды Международного симпозиума : в 2 т. / под ред. Н. К. Юркова. -Пенза : Изд-во ПГУ, 2013. - Т. 1 - С. 231-232.
13. Якимов А. Н. Проектирование микроволновых антенн с учетом внешних воздействий : моногр. / А. Н. Якимов. - Пенза : Изд-во Пенз. гос. ун-та, 2004. - 260 с.
14. Maksimov Е. Yu. A finite-element model of the thermal influences on a microstrip antenna / E. Yu. Maksimov, N. K. Yurkov, A. N. Yakimov // Measurement Techniques. - 2011. - Vol. 54, № 2, May. - P. 207-212.
15. Shishulin D. N. Modeling the Radiation of a Mirror Antenna taking Vibration Deformations into Account / D. N. Shishulin, N. K. Yurkov, A. N. Yakimov // Measurement Techniques. - 2014. -Vol. 56, № 11, February. - P. 1280-1284.
References
1. Balukov O. N., Yakimov A. N. Postroenie geometricheskoy modeli antenny s ispolzovaniem tri-angulyatsii Delone [Building a geometric model of the antenna using Delaunay triangulation], Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Tekhnicheskie nauki [News of Higher Educational Institutions. The Volga region. Engineering Sciences], 2009, no. 1, pp. 109-117.
2. Voskresenskiy D. I., Gostyukhin V. I., Maksimov V. M., Ponomarev L. I. Antenny i ustroystva SVCh [Antennas and microwave devices]. Moscow, Radiotekhnika, 2006. 376 p.
3. Vud P. Analiz i proektirovanie zerkalnykh antenn [Analysis and design of reflector antennas]. Moscow, Radio i svyaz, 1984. 208 p.
4. Gorbalysov M. S. Modelirovanie vliyaniya teplovykh vozdeystviy na kharakteristiki izlucheniya zerkalnoy antenny [Modeling the influence of thermal effects on the radiation characteristics of a reflector antenna]. Nadezhnost i kachestvo - 2013: trudy Mezhdunarodnogo simpoziuma [Reliability and quality -2013: Proceedings of the International Symposium], in 2 vol. Penza, Penza State Univ. Publ., 2013, vol. 1, pp. 92-95.
5. Drabkin A. L., Zuzenko V. L., Kislov A. G. Antenno-fidernye ustroystva [Antenna-feeder devices]. Moscow, Sovetskoe radio, 1974. 536 p.
6. Lunin L. S., Kataev V. F., Kreys I. V. Osobennosti chislennogo modelirovaniya datchikov na os-nove poverkhnostnykh akusticheskikh voln (PAV) metodom konechnykh elementov [Peculiarities of numerical modeling of sensors based on surface acoustic waves (SAW) finite element method]. Prikaspiyskiy zhurnal: upravlenie i vysokie tekhnologii [Caspian Journal: Management and High Technologies], 2011, no. 3, pp. 39-47.
7. Ponarin Ya. P. Elementarnaya geometriya [Elementary geometry], in 2 vol. Moscow, 2004. Т. 1: Planimetriya, preobrazovaniya ploskosti [Vol. 1: Planimetry, transformation of the plane]. 312 p.
8. Skvortsov A. V. Triangulyatsiya Delone i ее primenenie [Delaunay triangulation and its application], Tomsk, Tomsk Univ. Publ.], 2002. 128 p.
9. Kamaev V. A., Chivilev A. D. Avtomatizirovannyy raschet variantov elektricheskoy sostavlya-yushchey elektromekhanicheskogo privoda zerkalnykh antenn [Automated calculation of variations of the electric component of an electromechanical drive of reflector antennas]. Prikaspiyskiy zhurnal: upravlenie i vysokie tekhnologii [Caspian Journal: Management and High Technologies], 2014, no. 1, pp. 139-147.
10. Sabonnader Zh. K., Kulon Zh. L. Metod konechnykh elementov i SAPR [The finite element method and CAD], Moscow, Mir, 1989. 190 p.
11. Starchenko I. В., Kirichenko I. A., Sakharov V. L., Golosov P. S. Adaptivnoe upravlenie gidroakusticheskimi sistemami [Adaptive control of hydroacustic systems]. Prikaspiyskiy zhurnal: upravlenie i vysokie tekhnologii [Caspian Journal: Management and High Technologies], 2013, no. 3, pp. 26-33.
12. Shishulin D. N. Issledovanie vliyaniya vibratsiy na parabolicheskuyu antennu v ANSYS [Research of influence of vibrations on antenna dish in ANSYS]. Nadezhnost i kachestvo - 2013: trudy Mezhdunarodnogo simpoziuma [Reliability and quality - 2013: Proceedings of the International Symposium], in 2 vol. Penza, Penza State Univ. Publ., 2013, vol. 1, pp. 231-232.
13. Yakimov A. N. Proektirovanie mikrovolnovykh antenn s uchetom vneshnikh vozdeystviy [Design of microwave antennas taking into account external influences]. Penza, Penza State Univ. Publ., 2004. 260 p.
PRIKASPIYSKIY ZHURNAL: Upravlenie i Vysokie Tekhnologii (CASPIAN JOURNAL: Management and High Technologies), 2014, 2 (26) THEORETICAL BASES, DEVELOPMENT AND APPLICATION OF MATHEMATICAL MODELS, ALGORITHMS, COMPUTER PROGRAMS
14. Maksimov E. Yu., Yurkov N. K., Yakimov A. N. A finite-element model of the thermal influences on a microstrip antenna. Measurement Techniques, 2011, vol. 54, no. 2, May, pp. 207-212.
15. Shishulin D. N, Yurkov N. K., Yakimov A. N. Modeling the Radiation of a Mirror Antenna taking Vibration Deformations into Account. Measurement Techniques, 2014, vol. 56, no. 11, February, pp. 1280-1284.
УДК 004.94, 519.872
ДЕТЕРМИНИРОВАННАЯ ИМИТАЦИОННАЯ МОДЕЛЬ КЛАСТЕРОВ ГРИД-СИСТЕМЫ ДЛЯ СРАВНЕНИЯ ЭФФЕКТИВНОСТИ ИСПОЛЬЗОВАНИЯ ЭВРИСТИК РАСПРЕДЕЛЕНИЯ ЗАДАНИЙ
Гаевой Сергей Владимирович, аспирант, Волгоградский государственный технический университет, 400005, Российская Федерация, г. Волгоград, пр. Ленина, 28, e-mail: gaevserge@mail.ru
Аль-Хадша Фарес Али Хуссейн, аспирант, Волгоградский государственный технический университет, 400005, Российская Федерация, г. Волгоград, пр. Ленина, 28, e-mail: alhadsha@mail.ru
Фоменков Сергей Алексеевич, доктор технических наук, профессор, Волгоградский государственный технический университет, 400005, Российская Федерация, г. Волгоград, пр. Ленина, 28, e-mail: saf@vstu.ru
Лукьянов Виктор Сергеевич ^ дОКТОр технических наук, профессор, Волгоградский
государственный технический университет, 400005, Российская Федерация, г. Волгоград, пр. Ленина, 28
В работе осуществлено прогнозирование выполнения заданий на вычислительной Грид-
системе (Grid-системе), которая состоит из набора кластеров. Для достижения этой цели авторами была разработана детерминированная имитационная модель. В рамках этой модели кластерная система была представлена в виде системы массового обслуживания (СМО). Каждое из них имеет набор следующих параметров: ширина - количество вычислительных машин, необходимое для выполнения работы; длина - время выполнения задания; площадь (сложность) - произведение длины и ширины. Неограниченная по размеру очередь позволяет выполнить все поступающие на вход системы задания.
Такой подход позволяет определить среднюю длину очереди, среднее время ожидания задания и т.д. В настоящее время авторами исследуются методы оптимизации распределения работы между кластерами. Предполагается, что задание выполняется внутри только одного (любого) кластера GRID-системы. Каждый кластер имеет свой собственный входящий поток заданий от пользователей, отделов, организаций и т.п. Межкластерные соединения позволяют кластерам выполнять задания друг друга.
Процедура выбора кластера для исполнения задания требует введения некоторых критериев. В общем случае не существует ответа на вопрос «Какой критерий выбора кластера следует использовать?», нет и аналитического решения. Поэтому целесообразно промоделировать реальные кластеры и испытать различные способы распределения заданий. Из этих способов нужно выбрать оптимальные решения, удовлетворяющие заданным требованиям.
Авторами построена детерминированная имитационная модель для такой вычислительной системы, предложены эвристики для распределения заданий между кластерами Грид-системы и произведено сравнение эффективности этих эвристик для конкретного случая.
Ключевые слова: Грид-системы, брокер заданий, имитационное моделирование, детерминированная имитационная модель, распределение заданий, эвристики распределения заданий, системы массового обслуживания, вычислительный кластер, parallel workloads