УДК 517.977.56
ГРНТИ 27.31.15
УПРАВЛЕНИЕ КОЛЕБАНИЯМИ МАЧТОВЫХ АНТЕННЫХ КОНСТРУКЦИЙ
В.В. ПРОВОТОРОВ, доктор физико-математических наук, доцент
ВУНЦВВС «ВВА имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж)
А.П. ЖАБКО, доктор физико-математических наук, профессор
Санкт-Петербургский государственный университет (г. Санкт-Петербург)
Работа посвящена задаче управления антенной системой типа «мачта-растяжки», а именно, гашению несанкционированных колебаний конструкции (успокоению конструкции), инициированных внешними воздействиями. Формируется математическая модель процесса колебаний и решение поставленной задачи (имеется ввиду численная реализация задачи успокоения антенной конструкции). Осуществляется замена континуальной задачи для системы дифференциальных уравнений конечной («усеченной») системой алгебраических уравнений. Разрешимость такой системы представлена в терминах необходимых и достаточных условий для использования метода моментов. Таким образом, решение «усеченной» системы можно рассматривать как приближенное решение исходной континуальной задачи гашения колебаний с наперед заданной точностью. Представленный подход может использоваться в задачах оптимального управления, определения условий осцилляции и стабилизации и задачи сетевой экономики.
Ключевые слова: эволюционные системы, начально-краевые задачи, гиперболическое уравнение, распределенные параметры на графе, граничное управление, метод моментов.
ANTENNA MAST STRUCTURES OSCILLATION CONTROL
V.V. PROVOTOROV, Doctor Physico-mathematical Sciences, Assistant Professor
MESC AF «N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh)
A.P. ZHABKO, Doctor Physico-mathematical Sciences, Professor
St. Petersburg State University (St. Petersburg)
The work is devoted to the task of controlling the antenna system of the type «mast-stretching», namely, to quench unauthorized vibrations of the structure (to calm the structure) initiated by external influences. A mathematical model of the oscillation process and the solution of the problem (for the numerical implementation of the problem of calming the antenna structure) is being formed. The continual problem for the system of differential equations is replaced by a finite («truncated») system of algebraic equations. The solvability of such a system is presented in terms of the necessary and sufficient conditions for using the method of moments. Thus, the solution of the «truncated» system can be considered as an approximate solution of the initial continual problem of damping oscillations with a predetermined accuracy. The presented approach can be used in problems of optimal control, determination of conditions for oscillation and stabilization, and tasks of the network economy.
Keywords: evolutionary systems, initial-boundary value problems, hyperbolic equation, distributed parameters on a graph, boundary control, method of moments.
Введение. Исследование, представленное в работе, посвящено часто встречающейся в прикладных вопросах теории сплошных сред задаче перевода сложносочлененной промышленной конструкции, состоящей из системы упругих фрагментов, из фиксированного состоя-
ния в состояние покоя - задача успокоения нежелательных колебаний промышленной конструкции. В качестве сложносочлененной конструкции рассматриваются мачтовые антенные конструкции, управляющие воздействия на конструкцию осуществляются в местах закрепления антенных растяжек (граничное управление системой). Представлена математическая модель малых упругих поперечных колебаний конструкции, которая описывается системой волновых уравнений с распределенными параметрами на геометрическом графе. Работа разбита на две части: 1) представлена математическая модель, описывающая колебания антенной конструкции, где тело мачты имеет систему крепежных растяжек, зафиксированных в одной точке мачты - односекционная антенна и 2) математическая модель двухсекционной антенны с несколькими системами крепежных растяжек, зафиксированных в различных точках тела мачты, число таких точек равно Ь > 2. В обоих случаях фрагмент мачтовой конструкции, находящийся выше места закрепления (или всех мест закрепления) растяжек интерпретируется в виде сосредоточенной массы (точечной массы) М; фрагмент мачтовой конструкции, находящийся ниже всех мест закрепления растяжек, испытывает только продольные малые колебания [1]. Используемые для исследования антенные конструкции относятся к числу так называемых конструкций типа «мачта-растяжки». В качестве естественного инструмента математического описания рассматриваемых антенных конструкций используется ограниченный связный граф типа цепочки графов-звезд с фиксированным числом ребер у каждой звезды, в качестве математической модели, описывающей волновой процесс - начально-краевая задача для уравнения с частными производными гиперболического типа с распределенными параметрами на графе с обобщенными условиями Кирхгофа в местах сопряжения звезд. Следует отметить, что для исследуемого процесса становятся важными вопросы устойчивости и стабилизации получаемых решений [2-4].
Актуальность. Рассматривается актуальная в теории упругости задача о гашении колебаний сложносочлененной конструкции, состоящей из упругих линейных фрагментов, к таковым относятся мачтовые антенны различного типа, функционирующие в экстремальных условиях. Колебания конструкции инициируются внешними воздействиями, математический формализм которых описывается классом суммируемых функций с нетривиальным носителем на конечном связном ограниченном графе. Строится математическая модель мачтовой антенны с одним и несколькими местами креплений удерживающих растяжек. Используются формализмы эволюционных уравнений математической физики гиперболического типа с пространственной переменной, изменяющейся на графе [1].
Теоретической базой исследования являются методы теории оптимального управления. Указаны исчерпывающие условия применимости классического метода моментов для решения актуальной прикладной задачи: гашение нежелательных (и даже опасных) колебаний конструкции в результате воздействий сторонних (несанкционированных) внешних сил. Такие колебания гасятся с помощью специально рассчитанных внешних управляющих динамических воздействий, приложенных в наиболее доступных точках конструкции, в терминах математической модели таковыми являются граничные точки упругой системы. Представленный подход допускает анализ другой актуальной задачи, когда влияние на упругую систему осуществляется с временным запаздыванием, при этом используются уравнения с отклоняющимся (запаздывающим) временным аргументом. Это направление в данном исследовании не обсуждается и является темой специального анализа, частично проведенного в работах [2-4].
1. Односекционная антенная конструкция. Введем в рассмотрение геометрический связный граф-звезда Г с одним узлом % и примыкающими к нему т ребрами ук (к = 1, т) . Для класса функций с носителем на графе определим параметризацию и установим ориентацию на каждом ребре Г : параметризация и ориентация на ребрах ук (к = 1, т -1) осуществляется отрезком [О, I / 2] , а на ребре ут - соответственно отрезком [£ /2,£].
Обозначим через и(х,г), х,г еГх[0,Т] функцию, определяющую достаточно малые отклонения малых (точечных) элементов мачтовых конструкций (линейных и узловых фрагментов) от исходного состояния. Функция и(х, г) удовлетворяет системе дифференциальных уравнений гиперболического типа:
^и(х,г\к = и(х,г)Гк -д(х)Г1 и(х,г)^, х,геГ0х(0,Т),
'ч
'Гк'-
(1)
( и(х, г) - сужение функции и(х, г) на ребро ук, к = 1, т; Г0 - совокупность всех ребер графа
Г, не содержащих концевых точек), а также условиями сопряжения (условиями согласования):
(2)
5 ттЛ л V1 5 ттА л л
в узле %, несущем сосредоточенную массу М (фрагмент тела мачты, находящийся выше узла % заменен массой М ). Функции д( х)у = д( х)у =... = д( х)у - рассматривается достаточно часто встречающаяся на практике ситуация идентичности растяжек мачты [1, 5]; упругие характеристики фрагмента мачты, находящиеся выше узла % креплений растяжек мачты определяются функцией д(х)у , символ /у означает сужение функции / на ребро у. При этом остают-
/ т '
ся понятия и обозначения, используемые в работах [1, 5]. Присоединяя к соотношениям (1), (2) начальные:
и (х, 0) = К х)
(3)
и граничные:
и( х, г) = ^(г), х едГ
(4)
условия ( дГ - множество граничных узлов графа Г ). Начально-краевая задача (1) - (4) определяет математическую модель процесса колебаний мачтовой системы.
Управление колебаниями (гашение колебаний) системы (1)-(4) состоит в отыскании управляющих функций вк (г) (управлений) таких, что в момент времени г = Т имело место
д
и(х,Т) = 0, — и(х,Т) = 0, то есть мачтовая конструкция примет состояние динамического покоя.
дг
Как хорошо известно, применяя метод разделяющих переменных (метод Фурье) к (1)-(4) [6], получаем спектральную задачу (задачу Штурма-Лиувилля) на графе-звезде Г . В работе [6] содержится обоснование метода разделяющих переменных для эволюционных уравнений с частными производными 2-го порядка с распределенными параметрами на графе.
Задача Штурма-Лиувилля в пространстве С (Г) С 2[ Г] представляет собой совокупность уравнений:
-у" + д(х)у = Ху, х еук (к = 1, т),
(5)
во внутренней части каждого ребра ук (к = 1, т), при этом функция у(х, Х) удовлетворяет условиям согласования:
-у\£!2), +^У\£/2) +М(д(£/2)г -Л)у(£/2), =0
/ т к / т / т
к=1
во внутреннем узле % и однородным условиям в граничных узлах:
(6)
№Гк=0 (к = 1,т-1), у[ /) 0.
(7)
где Л - спектральный параметр, д(х) е С[Г]. Уравнением на звезде Г, как и выше, называется совокупность уравнений (5), (6), краевая задача на Г определяется соотношениями (5) - (7). Рассмотрим вспомогательное уравнение:
- г" + 0( х) г = Л, х е (0, л / 2) и (л / 2, л),
где искомая функция х) допускает разрыв:
- (г '(л /2 + 0) - г '(л /2 - 0)) + М (0(л / 2) - Л) г (л / 2) = 0,
(8)
(9)
пропорциональный значению г(л / 2) ; здесь:
0(х) = д(х), хе[0,л/2)и(л/2,л] = (^1 иут)\{%}, 0(л/2) = д(л/2)г .
Систему соотношений (8), (9) назовем дифференциальным уравнением на интервале (0,л), а функцию г(х) класса С[0,л]оС:([0,л]\{л/2})оС2((0,л)\{л/2}), удовлетворяющую системе (8), (9), назовем решением системы (8), (9).
Теорема 1. Пусть г0, г'0 - произвольные действительные числа. Система (8), (9) имеет
единственное решение г(х,Л) класса С[0,л]оС:([0, л]\{л/2})оС2((0, л)\{л/2}), удовлетворяющее начальным условиям г(0,Л) = г0, г' (0,Л) = г'0 (либо условиям г(л/2, Л) = г0, г ' ( л / 2 - 0 , Л) = г'). Это решение для каждого х из отрезка [0, л] является целой аналитической функцией параметра Л . Многообразие решений системы (8), (9) имеет размерность два.
Для уравнения (8) введем два линейно независимых решения х,Л) и г/(х, Л) класса
С[0,л]оС:([0, л]\{л/2})оС2((0, л)\{л/2}) ,которые удовлетворяют начальным условиям:
ц(И2-0,Х) = \, ¡Л{£ /2-ОД) = 0, г}(И2-0,Х) = 0, ?/(>/2-0Д) = 1.
Из непрерывности функций /л(х,Л), г/(х, Л) и из соотношения (9) вытекают следующие соотношения:
¡л{£! 2, Л) = /л{£! 2 + ОД) = 1, /л\£ / 2 + 0, А) = -М(Л - ()(£ / 2)),
г}(£/2,Л) = г}(£/2 + 0,А) = 0, г]\ж/2 + 0,Л) = \.
Заметим, что функции:
Э1
и
¿кЛ(х,ЛХ х
(г =1,т-1)
а (,
(к=1,т-1) Ы( х ,Л), х С/т ,
а
,(х,Л) = ц(х,Л), х еГ,
Г1, г = к,
(¿к =\ ) линейно независимы и, очевидно, образуют фундаментальную систему реше-
[0, г * к,
ний уравнения (5), удовлетворяющую условиям согласования (6). Предположим, что функции и(х,Л), V(х,Л) С[0,л]елС1[0, л]\{л/2})елС2((0, л)\{л/2}) с начальными условиями:
и
(0,Л) = 1, и(0,Л) = 0,
у(£,А) = 1, у,(£,Л) = 0,
являются также решениями уравнения (5), удовлетворяющих условиям согласования (6). В соответствии с утверждением теоремы 1 для каждого фиксированного действительного числа Л решения и (х, Л), V(х, Л) являются целыми аналитическими по Л функциями. Рассмотрим систему функций {Ок (х, Л)}т=1 вида:
вк (х,Л) = и\£ /2-0, (х,Л)~ сок (х, Л)), к = 2, т -1,
1 т— 1
ет (х,Л) =--у\£ /2-0, (х,Л) + ^ / 2' ¿К,(*> Л).
т-1 А.=1
1 т—\
вя (хД) =--у\£ /2-0, ЛУ£сок (х, А) + у(£ / 2, ЯК(х, Л).
т-1 к=1
явное представление которых определено соотношениями:
3 (х,Л) = •
и
(х,Л), х е/, г =1, 2,..., т-1,
(т - 1)и'(л/2 - 0, Л)^( х, Л) + и(л/ 2, Л)^ (х,Л), х е /т,
'(л/2-0,Л)ф,Л), хе/,
вк (х, Л) = < , ,
(к=2Щ {-¿к,и (л/2-0,Л)^(х,Л), х е/, г = 2,3,..., т,
(10)
От (х, Л) =
V (л/2 10,Л) х,Л) + V(л/2,Л)^(х,Л), х е/, г = 1,2,..., т-1,
т-1
V (х, Л) , х е /т .
Нетрудно показать, что функции Ок (х, Л) (к = 1, т) являются решениями уравнения (5), (6) и удовлетворяют краевым условиям:
Э1
и
01(0,Л)Л =0, вк(л,Х\к =0 (к = 2,т).
Систему линейно независимых функций {вк (х, Л)}т назовем фундаментальной системой
решений уравнения (5), (6). Множество решений уравнения (5), (6) образует т -мерное линейное многообразие, то есть произвольное решение уравнения (5), (6) является линейной комбинацией функций (10), принадлежащих системе {вк (х, Л)}т . Многообразие {вк (х, Л)}т используется для построения системы собственных функций спектральной задачи (5)-(7).
Теорема 2. Краевая задача (5)-(7) имеет только дискретный спектр, его элементы - вещественные числа. Собственные функции задачи (5)-(7) вещественные. При этом, собственные функции, соответствующие отличным друг от друга собственным значениям, ортогональны в
? (Г).
Доказательство указанных в условии теоремы спектральных свойств собственных значений и собственных функций почти полностью повторяет рассуждения классической теоремы на отрезке.
Обозначим через А следующее выражение:
Д = (т - 1)и' (ж/2 - 0, Л)у(ж /2, Л) - и (ж / 2, Л)у'(ж /2 - 0, Л) = = (т -1 )и' (ж /2 - 0, Л)у(ж / 2, Л) - и (ж / 2, Л)у '(ж /2 + 0, Л) --М (Л - £ (ж / 2))и (ж / 2, Л)т<ж / 2, Л),
и пусть О - спектр краевой задачи (5)-(7), т.е. множество собственных значений, через О', О'', О''' обозначим совокупности вида:
О' = {Л:В(Л) = 0,и(ж / 2, Л) * 0},
О '' = {Л: В(Л) = 0,и (ж/ 2, Л) = 0},
О ' '' = {Л: В (Л) * 0,и (ж/2, Л) = 0}.
Теорема 3. Множества О', О'', О' ' ' собственных значений Л, попарно не пересекаются и О = О' о>О'' о>О''', для любого Л справедливы утверждения:
1) пусть Л = Л'еО', тогда Л' простое,
2) пусть Л = Л'еО'', тогда Л'имеет кратность, равную т-1,
3) пусть Л = Л'''еО''', тогда Л''' имеет кратность, равную т - 2.
Доказательство теоремы 3 представлено в [6].
Замечание. Из утверждений теоремы 3 следует представление собственных функций задачи (5)-(7):
1) пусть ЛеО', тогда ему соответствует собственная функция в1(х, Л') ;
2) пусть Л' еО'', тогда ему соответствуют т-1 собственных функций вк (х, Л") (к = Гт-1):
Х( х, X' ) = %( х, X''),
Х2 (х, X' ) = ^ (х, X'),
к—1
хк (х, X') = —£ - X (X, X') + вк (х, X'), к = 3, т -1;
3) пусть X ' 'ей''', тогда ему соответствуют т — 2 собственные функции %к(х,X''),
к = 2, т—1 (Хк (х, XX '') совпадают с вк (х, X'), если X'n''=X'' (к = 2, т—1)).
Очевидно, что, используя процесс ортогонализации Сонина-Шмидта для функций 0к(х, X'), можно считать все собственные функции взаимно ортогональными.
Нетрудно проверить, что при условии непрерывности на каждом ребре функции д (х) (т.е. условии выполнения д(х) е С[Г]), дифференциальный оператор, порожденный краевой задачей (5) - (7), положительно определен, а значит, положительным является его спектр, т.е. являются положительными все собственные значения и на бесконечности находится их точка сгущения.
Расположим все собственные значения краевой задачи (5)-(7) в порядке возрастания и с учетом их кратностей:
0<Я0 <Я1 <л2 < ... <\ <...
(заметим, что в приведенной цепочке неравенств число Xk присутствует ровно столько раз, какова его кратность - см. утверждение теоремы 3). В связи с этим, собственные функции нумеруются также, как нумеруются им соответствующие собственные значения Хп (х) (для одинаковых собственных значений им соответствующие собственные функции берутся в произвольном порядке (см. замечание к теореме 3)).
Обоснованием использования метода Фурье (метода разделения переменных) для начально-краевой задачи (1)-(4) являются приводимые ниже утверждения (полные доказательства см. в работах [5, 6]).
Теорема 4. Множество {Хп (х)} 0 обладает свойством спектральной полноты в пространстве (Г) , то есть (Г) является сепарабельным пространством с ортонормированным базисом {Хп (х)}п^0 .
Теорема 5. Любая абсолютно непрерывная функция /(х), х еГ разлагается в обобщенный ряд Фурье по системе {Хп (х)}я>0:
а =
/ (х) = ^апХп (х),
п=0
- ¡/ (г) Хп (г) йг, ®п = \\Х] (г) йг,
'^п Г УГ
(11)
сходимость ряда (11) равномерная относительно хеГ .
Ниже представлен анализ математической модели колебательного процесса антенной конструкции, описываемой начально-краевой задачей (1)-(4). Пусть достаточно гладкая функция и (х, г) является решением задачи (1)-(4) и ее разложение в обобщенный ряд Фурье имеет
вид:
1 = 2
ад
и (х, г) = £ии (г) Х (х), ип (г) = \и (х, г) Х (х) дх,
п=0 г
здесь базисные функции Хп (х) определяются системой функций {вк (х, Х)}т (см. представления (10)). Дважды продифференцируем функцию ип (г) и, используя уравнение (1), (2), получаем дифференциальное уравнение 2-го порядка для коэффициентов ип (г) :
<(г)+ р2(г) = ^(г), г>0, Л = рр,
(12)
т-1
где 2 (/) = ~УЦк^)Х'(0) +у(1)Х'(1) . Решение уравнения (12) для произвольного п, учиты-
< ' "к 'т
к=1
вая начальные условия при г = 0, определяются соотношениями:
1 г
и. (г) = рп соб рпг+— sin рпг+— \гп (г) sin рп (г - г)дг,
рп Рп 0
(13)
через рп, —п обозначены обобщенные коэффициенты Фурье начальной и граничной функций
ад ад
р(х), — (х), разложенных по базису {Хп (х)} о: р(х) = ^Р„Хп (х), —(х) = Ё—пХп (х). Значит,
п=0
п=0
ди
решение и(х, г) начально-краевой задачи (1)-(4) вместе с производной -(х, г) имеют пред-
дг
ставление в виде рядов:
и (х, г ) = ¿р соб рпг+—бш рпг]Хп (х)+£ ¡гя (фт рп (г-г)дг, (14)
рп П=0 рп 0
п=0
ди_
дг
ад ад 1
(x, г) = £[-РпРп Sin Лпг + —п С0Б Рпг]Хп (х) + ЁХп (х)\2п (г) С0Б Рп (г - г)дг. (15)
ди
Динамический покой системы (1), (2) (т.е. и (х, г) =-(х, г ) = 0, х, г еГх[0, Т]), определяется
дг
совокупностью уравнений:
£р соб Рпг+—sin Рпг ] Хп (х)+Ё Хп(х) \гп (г) sin Рп (г - г)дг = 0,
рп п=0 рп о
п=0
ад ад 1
Ё[-РпРп Рпг + —п С0Б Рпг]Хп (х) + ЁХп (х)\2п (г) С0Б Рп (г - г)дг = 0,
(16)
п=0
п=0 0
которая определяет совокупность условий на коэффициенты рп, —п - условия гашения колебаний в системе (1), (2) для произвольного временного момента г. Соотношения (16), а значит, и систему (1), (2), удобно рассматривать в фиксированный момент времени г = Т. Сравнивая коэффициенты при функциях Хп (х) системы {Хп (х )}яй 0, получаем совокупность соотношений ( п = 0,1,2,...):
T m—1
-срп cos рпТ - ^ sin рпТ = — (г)х; (0) - v{r)X'n {£)} sin рп (Т - z)dz,
Pn Pn 0 k=1
(17)
T m—1
ap„<p„ sin pT - y/n eos pT = J*[^y/Á (r)X'n (0) - v(t)X'„ (£)] cos pn (T - z)dz.
0 k=1 k
Рассматривая это как систему соотношений относительно выражений:
T m—1
{[!>, (0):. - v(r)X'n (/)] cos pnr)dT 0 k=1
и
T m—1
(r)X,; (0) - v(r)X;; (/)] Sin pnz)dz
0 k=1
( и = 0,1,2,...) и разрешая ее (определитель при этом равен единице!), приходим к новой системе соотношений ( п = 0,1,2,...):
с с
IX М (r)cos KTdT - К eos Ánzdz = -у/п,
k=1 0 0 »1-1 т т sin P„T^jX'„ (0)/t Ja (T) ~ X'„ (£)jv(z) sin pnzdz = pncpn.
(18)
Таким образом, получена удобная для практической реализации система уравнений (18) относительно так называемых моментов функций / (/) (к = 1, т — 1), у(/) (система моментных равенств), являющихся граничными управляющими воздействиями. И задача гашения колебаний системы (1)-(4) свелась к известной «¿'-проблеме моментов» (см. [4, 5]): граничные управляющие воздействия /к (/) (к = 1, т — 1), у(/) определяются решением системы (18) с любой, наперед заданной точностью, зависящей только от числа взятых номеров п = 0,1,2,....
Замечание 1. Колебания антенных фрагментов можно погасить, используя в качестве управляющей функции только одно из граничных управляющих воздействий
/к (х) (к = 1, т — 1), например, / (х) (остальные фиксированы). Моментная система для этой
функции аналогична системе (18).
Замечание 2. Сетеподобные антенные конструкции, представляют собой набор звезд типа Г с общими ребрами (разумеется, масса в узлах крепления растяжек исключается). При этом, возникают так называемые циклы (петли) [1] - подграфы, содержащие замкнутые контуры, состоящие из последовательно связных между собой ребер (начало реберного контура совпадает с его концом). Анализ колебаний такой сетеподобной антенной конструкции аналогичен приведенному выше.
2. Ь -уровневая антенная конструкция. Для упругой Ь -уровневой антенной конструкции ниже рассматривается аналогичная задача гашения колебаний. Математическая модель малых поперечных колебаний использует формализмы краевых задач для эволюционных уравне-
ний с частными производными гиперболического типа, пространственная переменная которых изменяется на множестве, представляющем собой ориентированный граф-цепочку (обозначим его через 3 ). Такой граф структурирован Ь последовательно соединенными между собой гра-
фами-звезда Г, (£ = 1,1,):3 = (^Г/ . Ребра у'к (к = \,т,, £ = \,Ь), а также узлы (£ = \,Ь),
/=1
графа 3 параметризованы и ориентированы по аналогии с параметризацией и ориентацией графа-звезды Г раздела 1.
Обозначим через П(х,г), х,г е3х[0,Т] достаточно гладкую функцию, которая описывает амплитудные изменения малых перемещений достаточно малых фрагментов антенной конструкции при произвольном значении пространственной переменной х и временной переменной г (см. обозначения работы [1]). Математическая модель колебательного процесса в антенной конструкции описывается совокупностью уравнений:
д2 д2
— П(хг) , =— П(хг) , - х) 1:П(хг)
дГ Ч дх~ Ч Ч Ч
(19)
на каждом ребре ук {к = \,т,, I = соотношениями (условиями трансверсальности):
Ж
ж
П(- г), = П(- г ) , , к = 1, т -1,
2 Ч 2 у£
тр
дС1 ,,ж . д£1,.ж , дС1 .. п .
(20)
в узлах £ = 2, Ь и соотношениями:
ж
ж
П(-г)! = П(-г), , к = 1,т -1,
2 у к 2 'тк
дП ж ч т1-1 дП ж ч , , ж. ч , ,д2П ж .
7Т (9 , г V -£т~ (9 , г V , г V = (9 , г V
дх 2 у-1 к=1 дх 2 2 7т1 2 у-1 д7 2 у-1
(21)
в узле ^, несущему сосредоточенную массу М (количественная характеристика оконечности антенны). При этом имеют место начальные условия (начальное состояние антенны):
д -
О (х, 0) = т( х), — О (х, 0) = т( х), дг
(22)
при х е 3, ^ = 0 и граничные условия:
0(^-1)^/2,0, =/а0 (* = \,т-\(. = \,Ь\ ЩЬ + \)тг/2,Г) = КО, (23)
/ т/. / т^
в граничных узлах цепочки 3 (г е[0, Т]). Отметим, что возможны и иные граничные условия, определяемые конструкцией крепежа растяжек.
Ь
Спектральная задача (задача Штурма-Лиувилля) на цепочке 3 для функций класса С (3) С 2[ 3] определяется совокупностью уравнений:
~У / + Ф) /■ У < = Лу , (к = \,т„ £ = 1,1),
гк гк гк гк
(24)
на внутренних частях ребер у'к при выбранной параметризации (Л- спектральный параметр), соотношениями:
»Я,—1
у\£п / 2) + ¿>'(^/2), =У(^/ 2), ,
у" у' у"
в узлах , £ = 2,Ь каждой звезды Г, = 2,Ь), соотношением:
(25)
т
-1
£ /(л/2)1 М (д (л/2)1 -X) у (л/2)1 = / (л/2) 1
£=1 ^ ^^
(26)
в узле ^ звезды Г1 и краевыми условиями:
у((£-\)ж12) г= О, А: = 1,£ = 2,1, +1)^/2)^ = 0.
Гк
(27)
Заметим, что априори задаваемое семейство функций //ДО (к = 1,т, -1, £ = фигурирующих в формализмах (23) граничных условий, представляет множество управлений изучаемой системы (19) - (23), состояние которой является функцией 0(х,г), х,г е3х[0,Т] класса С2 ( 3х[0, Т ]).
Введем следующие ограничения:
тах | //ДО |<С {к = \,т, -1, £ = \,1),
(28)
где С >0 - фиксированная постоянная.
Задача гашения колебаний Ь -уровневневой антенной конструкции, т.е. дифференциальной системы (19) - (23) состоит в определении управляющих воздействий
¡л'к(0 (к = 1,т, — 1, £ = 1, V), удовлетворяющих ограничениям (28) и таких, чтобы условия:
О (х,Т ) = 0, — О (х,Т ) = 0, дг
(29)
(условия успокоения системы (19) - (23)) выполнялись при минимальном значении временной переменной, т.е. минимальном значении Т .
Таким образом, учитывая начальные распределения (22) амплитуд и скоростей изменения параметров колебаний, необходимо через минимальное время Т привести к динамическому покою все фрагменты антенной конструкции. Это означает равенство нулю как распределение амплитуд колебаний, так и скоростей изменения колебаний - эволюционная система (19)-(23) переведена из возбужденного состояния (22) в нулевое (29) или состояние покоя.
В рамках вышеприведенного метода моментных соотношений предлагается универсальный способ вычислительных операций, не зависящий ни от сложности и порядка системы дифференциальных уравнений, ни от числа функций, являющихся управляющими воздействиями. Этот метод, требующий только знание множества собственных функций системы (19)—(23), то есть множества {Xn (x)}n>1 собственных функций задачи Штурма-Лиувилля (24)—(27), является основополагающим в задаче гашения колебаний антенной системы произвольной мачтовой конструкции. При этом точность отыскания решения определяется только числом собственных функций.
Исходя из представления функции Q(х, t), аналогичное функции U(х, t) (раздел 1), и представления правых частей zn (t) для уравнения (12), условия (29) дают следующую совокупность соотношений ( n = 1,2,...):
~Tn C0s PnT--- sln PnT =
Pn
1 T. L '", -1
— {["EE/4(0^-1^/2) f+v{QXn{{L + \)xl2) ]sinPn{T-OdC
~Tn cos PnT--sin pnT =
Pn
1 T. L "Y~l
, +v(C)Xn((L + l)7Tl2) ]sinp„(T-C)dC-
Pn 0 /■=! *=1 Г*
Рассматривая эту систему (определитель ее для каждого фиксированного п = 0,1,2,... равен единице) относительно величин:
Т L п,-1
{[-££^(0^-1)^2) , +v(£)Xn((L + 1)х/2) ]cosрпСЖ
J 77 1-г ri- ""
О '=1 к=1 *
и
Т L п,-1
, +v(OX„((Z + l)T/2b ]sm
о Ml k=\ ,mi
получаем новую систему (n = 0,1,2,...) вида:
L >»i
-1
£ Ё xn ((£ -1 / 2)f \t4 (O cos k£dZ + Xn ((L +1)n / 7) \v(£) cos = -z„,
/=1 k=1 )k А 'т1 A
^ 1 1 sin р/Ё Ё -1 / 2V К (О sin + ((L + l)xt 2)}i jv(£) sin pn£d£ = pnzT.
Ml k = 1 ^ A ™L A
или
ш-1 1 1
Yx'n (0\ \рк (г) COS A„TdT - Х'„ Щу{т) cos Хптйт ■
k=1 0 0
m-\
sin pjJ^K (0\ fa (т) - X'n {l)\v{z) sin pjc/r = p
k=1 о 0
Vn
Тем самым получена удобная (хотя и достаточно громоздкая) для компьютерной реализации аналогичная (17) моментная система уравнений для функции р'к (?) {к = 1, — 1, £ = 1,Ь), ), являющихся граничными управляющими функциями.
Выводы. Дальнейшее движение на пути получения решения поставленной задачи (имея ввиду численную реализацию задачи успокоения антенной конструкции) - замена полученной бесконечной системы алгебраических уравнений конечной («усеченной») системой (возникает конечная проблема моментов). Разрешимость такой системы представлена в терминах необходимых и достаточных условий в работе [7; § 2, гл. III]. Таким образом, решение «усеченной» системы можно рассматривать как приближенное решение исходной задачи гашения колебаний с наперед заданной точностью (см. задачи оптимального управления [8-12], задачи определения условий осцилляции и стабилизации [13-16] и задачи сетевой экономики [17-22]).
СПИСОК ЛИТЕРАТУРЫ
1. Провоторов В.В. Моделирование колебательных процессов «мачта-растяжки» // Системы управления и информационные технологии. 2008. № 1 (31). С. 272-277.
2. Aleksandrov A.Y., Aleksandrova E.A., Zhabko A.P. Asymptotic stability conditions for certain classes of mechanical systems with time delay // WSEAS Transactions on Systems and Control. 2014. vol. 9. P. 388-397.
3. Александров А.Ю., Жабко А.П. Об устойчивости решений одного класса нелинейных систем с запаздыванием // Автоматика и телемеханика. 2006. № 9. С. 1355-1365.
4. Aleksandrov A.Y., Zhabko A.P., Hu G.-D. Delay-independent stability conditions for some classes of nonlinear systems // IEEE Transactions on Automatic Control On stability and dissipativity of some classes of complex systems. 2014. vol. 59. no. 8. P. 2209-2214.
5. Провоторов В.В. Метод моментов в задаче гашения колебаний дифференциальной системы на графе // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2010. № 2. С. 60-69.
6. Провоторов В.В. Собственные функции задачи Штурма-Лиувилля на графе-звезде // Математический сборник. 2008. Т. 199. № 10. С. 105-126.
7. Бутковский А.Г. Теория оптимального управления системами с распределенными параметрами. М.: Наука, 1965. 474 с.
8. Подвальный С.Л., Провоторов В.В. Определение стартовой функции в задаче наблюдения параболической системы с распределенными параметрами на графе // Вестник Воронежского государственного технического университета. 2014. Т. 10. № 6. С. 29-35.
9. Provotorov V.V., Ryazhskikh V.I., Gnilitskaya Yu.A. Unique weak solvability of a nonlinear initial boundary value problem with distributed parameters in a netlike region // Vestnik of Saint Petersburg University. Applied mathematics. Computer sciense. Control processes. 2017. Volume 13. Issue 3. P. 264-277.
10. Podvalny S.L., Provotorov V.V., Podvalny E.S. The controllability of parabolic systems with delay and distributed parameters on the graph // В сборнике: Procedia Computer Sciense 12th. Сер. «12th. International Symposium Intelligent Systems, INTELS 2016». 2017. P. 324-330.
11. Провоторов В.В., Провоторова Е.Н. Синтез оптимального граничного управления параболической системы с запаздыванием и распределенными параметрами на графе // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2017. Т. 13. №. 2. С. 209-224.
12. Карелин В.В. Штрафные функции в задаче управления процессом наблюдения // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2010. № 4. С. 109-114.
13. Kamachkin A.M., Yevstafyeva V.V. Oscillations in a relay control system at an external disturbance // Control Applications of Optimization 2000: Proceedings of the 11th IFAC Workshop. 2:(2000). P. 459-462.
14. Веремей Е.И., Сотникова М.В. Стабилизация плазмы на базе прогноза с устойчивым линейным приближением // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2011. № 1. С. 116-133.
15. Aleksandrov A.Y., Aleksandrova E.A., Zhabko A.P. Asymptotic stability conditions for certain classes of mechanical systems with time delay // WSEAS Transactions on Systems and Control. 2014. Vol. 9. P. 388-397.
16. Aleksandrov A.Y., Aleksandrova E.A., Zhabko A. P. Asymptotic stability conditions of solutions for nonlinear multiconnected time-delay systems // Circuits Systems and Signal Processing. 2016. Vol. 35. no. 10. P. 3531-3554.
17. Krasnov S.V., Sergeev S.M., Mukhanova N.V., Grushkin A.N. Methodical forming business competencies for private label // Reliability, Infocom Technologies and Optimization (Trends and Future Directions) 6th International Conference ICRITO. 2017. P. 569-574.
18. Сергеев С.М. Выбор инновационной маркетинговой стратегии предприятий на основе экономико-математического моделирования // Инновации. 2013. № 3 (173). С. 116-119.
19. Сергеев С.М. Математическое моделирование потоков через POS-терминалы // Вестник Тамбовского университета. Серия: Естественные и технические науки. 2013. т. 18. № 1. С. 227-229.
20. Iliashenko O.A., Krasnov S.V., Sergeev S.M. Calculation of high-rise construction limitations for non-resident housing fund in megacities. E3S Web of Conferences. Volume 33. 6 March 2018. Paper number 030062017 International Scientific Conference on High-Rise Construction. HRC 2017. P. 425-437.
21. Сергеев С.М. Формирование структуры управления сложными коммерческими сетями // Современные методы прикладной математики, теории управления и компьютерных технологий (ПМТУКТ - 2015). Сборник трудов VIII Международной конференции. 2015. С. 324-326.
22. Borisoglebskaya L.N., Sergeev S.M. Model of assessment of the degree of interest in business interaction with the university // Journal of Applied Economic Sciences 2018. Volume XII. Issue 8 (54) Winter 2017. P. 2423-2448.
REFERENCES
1. Provotorov V.V. Modelirovanie kolebatel'nyh processov «machta-rastyazhki» // Sistemy upravleniya i informacionnye tehnologii. 2008. № 1 (31). pp. 272-277.
2. Aleksandrov A.Y., Aleksandrova E.A., Zhabko A.P. Asymptotic stability conditions for certain classes of mechanical systems with time delay // WSEAS Transactions on Systems and Control. 2014. vol. 9. pp. 388-397.
3. Aleksandrov A.Yu., Zhabko A.P. Ob ustojchivosti reshenij odnogo klassa nelinejnyh sistem s zapazdyvaniem // Avtomatika i telemehanika. 2006. № 9. pp. 1355-1365.
4. Aleksandrov A.Y., Zhabko A.P., Hu G.-D. Delay-independent stability conditions for some classes of nonlinear systems // IEEE Transactions on Automatic Control On stability and dissipativity of some classes of complex systems. 2014. vol. 59. no. 8. pp. 2209-2214.
5. Provotorov V.V. Metod momentov v zadache gasheniya kolebanij differencial'noj sistemy na grafe // Vestnik Sankt-Peterburgskogo universiteta. Seriya 10: Prikladnaya matematika. Informatika. Processy upravleniya. 2010. № 2. pp. 60-69.
6. Provotorov V.V. Sobstvennye funkcii zadachi Shturma-Liuvillya na grafe-zvezde // Matematicheskij sbornik. 2008. T. 199. № 10. pp. 105-126.
7. Butkovskij A.G. Teoriya optimal'nogo upravleniya sistemami s raspredelennymi parametrami. M. : Nauka, 1965. 474 p.
8. Podval'nyj S.L., Provotorov V.V. Opredelenie startovoj funkcii v zadache nablyudeniya parabolicheskoj sistemy s raspredelennymi parametrami na grafe // Vestnik Voronezhskogo gosudarstvennogo tehnicheskogo universiteta. 2014. T. 10. № 6. pp. 29-35.
9. Provotorov V.V., Ryazhskikh V.I., Gnilitskaya Yu.A. Unique weak solvability of a nonlinear initial boundary value problem with distributed parameters in a netlike region // Vestnik of Saint Petersburg University. Applied mathematics. Computer sciense. Control processes. 2017. Volume 13. Issue 3. pp. 264-277.
10. Podvalny S.L., Provotorov V.V., Podvalny E.S. The controllability of parabolic systems with delay and distributed parameters on the graph // V sbornike: Procedia Computer Sciense 12th. Ser. «12th. International Symposium Intelligent Systems, INTELS 2016». 2017. pp. 324-330.
11. Provotorov V.V., Provotorova E.N. Sintez optimal'nogo granichnogo upravleniya parabolicheskoj sistemy s zapazdyvaniem i raspredelennymi parametrami na grafe // Vestnik Sankt-Peterburgskogo universiteta. Seriya 10: Prikladnaya matematika. Informatika. Processy upravleniya. 2017. T. 13. №. 2. pp. 209-224.
12. Karelin V.V. Shtrafnye funkcii v zadache upravleniya processom nablyudeniya // Vestnik Sankt-Peterburgskogo universiteta. Seriya 10: Prikladnaya matematika. Informatika. Processy upravleniya. 2010. № 4. pp. 109-114.
13. Kamachkin A.M., Yevstafyeva V.V. Oscillations in a relay control system at an external disturbance // Control Applications of Optimization 2000: Proceedings of the 11th IF AC Workshop. 2:(2000). pp. 459-462.
14. Veremej E.I., Sotnikova M.V. Stabilizaciya plazmy na baze prognoza s ustojchivym linejnym priblizheniem // Vestnik Sankt-Peterburgskogo universiteta. Seriya 10: Prikladnaya matematika. Informatika. Processy upravleniya. 2011. № 1. pp. 116-133.
15. Aleksandrov A.Y., Aleksandrova E.A., Zhabko A.P. Asymptotic stability conditions for certain classes of mechanical systems with time delay // WSEAS Transactions on Systems and Control. 2014. Vol. 9. pp. 388-397.
16. Aleksandrov A.Y., Aleksandrova E.A., Zhabko A. P. Asymptotic stability conditions of solutions for nonlinear multiconnected time-delay systems // Circuits Systems and Signal Processing. 2016. Vol. 35. no. 10. pp. 3531-3554.
17. Krasnov S.V., Sergeev S.M., Mukhanova N.V., Grushkin A.N. Methodical forming business competencies for private label // Reliability, Infocom Technologies and Optimization (Trends and Future Directions) 6th International Conference ICRITO. 2017. pp. 569-574.
18. Sergeev S.M. Vybor innovacionnoj marketingovoj strategii predpriyatij na osnove 'ekonomiko-matematicheskogo modelirovaniya// Innovacii. 2013. № 3 (173). pp. 116-119.
19. Sergeev S.M. Matematicheskoe modelirovanie potokov cherez POS-terminaly // Vestnik Tambovskogo universiteta. Seriya: Estestvennye i tehnicheskie nauki. 2013. t. 18. № 1. pp. 227-229.
20. Iliashenko O. Krasnov S.V., Sergeev S.M. Calculation of high-rise construction limitations for non-resident housing fund in megacities. E3S Web of Conferences https://www.scopus.com/sourceid/21100795900?origin=recordpage. Volume 33. 6 March 2018. Paper number 030062017 International Scientific Conference on High-Rise Construction. HRC 2017. pp. 425-437.
21. Sergeev S.M. Formirovanie struktury upravleniya slozhnymi kommercheskimi setyami //
Sovremennye metody prikladnoj matematiki, teorii upravleniya i komp'yuternyh tehnologij (РМТ1ЖТ-2015). Sbornik trudov VIII Mezhdunarodnoj konferencii. 2015. pp. 324-326.
22. Borisoglebskaya L.N., Sergeev S.M. Model of assessment of the degree of interest in business interaction with the university // Journal of Applied Economic Sciences 2018. Volume XII. Issue 8 (54) Winter 2017. pp. 2423-2448.
© Провоторов B.B., Жабко А.П., 2019
Провоторов Вячеслав Васильевич, доктор физико-математических наук, доцент, старший научный сотрудник 22 отдела научно-исследовательского 2 управления научно-исследовательского научно-исследовательского центра (проблем применения, обеспечения и управления авиацией ВВС), Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, [email protected].
Жабко Алексей Петрович, доктор физико-математических наук, профессор, заведующий кафедрой теории управления, Санкт-Петербургский государственный университет (г. Санкт-Петербург), Россия, 199034, г. Санкт-Петербург, Университетская набережная, 7-9, /habко.арmath. spbu a mail, ru.