Известия РАН. Механика твердого тела, 2022, № 4, стр. 51-69
ЗАЩИТА ОБЪЕКТА НА ПОДВИЖНОМ ОСНОВАНИИ С ПОМОЩЬЮ ПОСТОЯННОГО УПРАВЛЕНИЯ ПРИ ОГРАНИЧЕННЫХ ВОЗМУЩЕНИЯХ
a Институт проблем механики им. А.Ю. Ишлинского РАН
Москва, Россия
* E-mail: korneev@ipmnet.ru
Поступила в редакцию 30.03.2021
После доработки 11.05.2021
Принята к публикации 06.09.2021
- EDN: ENLVUQ
- DOI: 10.31857/S0572329922030096
Аннотация
Рассматривается задача построения гарантирующего управления противоударным изолятором, защищающим объект на подвижном основании от ударных воздействий, которым подвержено основание. Ударное воздействие на основание предполагается заранее неизвестным и определяется его ускорением, ограниченной знакопостоянной функцией времени заданной длительности, интеграл от которой по времени задан. Управляющая сила действует между основанием и защищаемым объектом, ограничена по величине, а абсолютное ускорение основания может превысить максимально допустимую величину абсолютного ускорения объекта только на одном интервале времени. Минимизируемым критерием качества служит наибольшее смещение объекта относительно основания. В качестве ограниченного по величине управления, действующего между основанием и защищаемым объектом, было взято постоянное управление заданной длительности, являющееся оптимальным для мгновенного удара в задаче без упреждения. Получен оптимальный момент начала действия этого управления. Приведены сравнительные оценки по значению критерия качества предлагаемого управления с другими способами управления.
Введение. В работе рассматривается построение способов управления противоударным изолятором при ограниченных наихудших возмущениях. Используемая математическая модель системы состоит из основания и защищаемого объекта, находящегося на нем, и обладает одной степенью свободы. Предполагается, что имеет место внешнее воздействие кинематического типа, задаваемое ускорением основания, определенным как функция времени. Основание и защищаемый объект могут двигаться поступательно вдоль одной заданной прямой. Управляющая сила между основанием и объектом создается изолирующим устройством и ограничена по величине. Впервые задача о минимизации максимума модуля смещения объекта относительно основания при заданном заранее возмущении была сформулирована в [1–3]. В дальнейшем теория оптимальной противоударной изоляции получила развитие в [4–7]. В [8] изучалась задача защиты объекта, расположенного на подвижном основании, от ударных воздействий малой длительности с помощью изолятора с управлением без упреждения. Численное решение задачи из [8] при заданном возмущении и построение ее аналитических решений для некоторых классов возмущений с упреждающим управлением рассмотрены в [9]. Аналитическое гарантирующее упреждающее управление с одним переключением и оптимальным временем упреждения и отвечающее им минимальное значение наибольшей величины смещения объекта относительно основания для класса возмущений заданной длительности с неограниченной амплитудой получены в [10, 11]. Для этого же класса возмущений в [12, 13] показано, что, используя упреждение и запаздывание для двух простых управлений, предназначенных для конкретных возмущений, и оптимизируя функционал только по моменту начала управления, можно получить значение критерия качества близкое к оптимальному. В [14] рассмотрена задача защиты объекта на подвижном основании для класса возмущений заданной длительности и ограниченной амплитуды. Доказано, что при заданном управлении наихудшими возмущениями являются непрерывные прямоугольные возмущения наибольшей амплитуды, начинающиеся в момент начала возмущения или заканчивающиеся в момент окончания допустимого интервала возмущения. Кроме того, в [14] для рассматриваемого класса возмущений изучалась возможность применения двух управлений, полученных в [9, 11] для прямоугольного возмущения и для возмущений заданной длительности с неограниченной амплитудой соответственно. Оптимизация критерия качества проводилась по моменту начала действия этих управлений. Предложено использовать комбинированное управление, сочетающее достоинства обоих управлений.
Обзор современного состояния проблемы противоударной изоляции, включающий в себя публикации по оптимальному управлению противоударными изоляторами, представлен в [15].
В данной работе, так же как и в [14], рассматривается задача защиты объекта на подвижном основании для класса возмущений заданной длительности и ограниченной амплитуды. В классе постоянных управлений ограниченной длительности проводится минимизация возможной величины смещения защищаемого объекта по начальному моменту действия управления. Варьирование момента начала действия управления из рассматриваемого класса предполагает, что возможно упреждение и запаздывание. Приведены сравнительные оценки по значению критерия качества предлагаемого управления со способами управления, рассмотренными в [14]. Выбор для изучения простых законов управления, зависящих от амплитуды возмущения и не зависящих от его длительности, обусловлен тем, что такие способы управления легче реализовать на практике, когда закон управления задан и выбирается только момент начала его действия в зависимости от длительности и амплитуды возмущения.
1. Механическая система. Рассматриваемая механическая система (рис. 1) содержит основание и объект, соединенный с основанием с помощью противоударного изолятора – устройства, которое генерирует управляющую силу f между основанием и объектом и предназначено для защиты объекта при ударном воздействии на основание. Основание и защищаемый объект движутся поступательно вдоль заданной прямой. Пусть $z$ – смещение основания относительно инерциальной (неподвижной) системы отсчета, $x$ – смещение объекта относительно основания, $m$ – масса объекта. Ударное воздействие на основание определяется его ускорением $\ddot {z}$, функцией времени, некоторые свойства которой известны до начала воздействия.
Движение объекта относительно основания выражается уравнением
в котором сила f удовлетворяет ограничению $|{\kern 1pt} {\kern 1pt} f{\kern 1pt} {\kern 1pt} |\, \leqslant {{F}_{0}}$, где F0 – заданная величина. Величина u удовлетворяет неравенствуДалее полагаем, что в начальный момент времени t = 0 основание и объект имеют нулевые значения координат и скоростей:
Допустимыми управлениями будем считать кусочно-непрерывные функции $u(t),$ удовлетворяющие ограничению
2. Внешние возмущения. Будем рассматривать возмущения ${v}(t)$ вида
где $V(\xi )$ – кусочно-непрерывная функция, заданная для всех вещественных $\xi $, причем $V(\xi ) \equiv 0$ для $\xi < 0$, а ${{t}_{0}} \geqslant 0$ – некоторый момент времени, который может быть известен заранее или подлежать определению. Возмущение V начинает действовать на основание только в момент ${{t}_{0}}$ после включения системы противоударной изоляции (управление с упреждением).Будем предполагать, что 1) возмущение $V(\xi )$ действует только в одном заданном направлении и не меняет знака ($V(\xi ) \geqslant 0$), 2) $V(\xi )$ имеет заданную длительность T ($V(\xi ) \equiv 0$, если $\xi > T$), 3) $V(\xi )$ ограничена максимальной величиной воздействия, обозначаемой параметром b,
и 4) только на одном интервале, ${{t}_{1}} < \xi < {{t}_{2}}$, величина абсолютного ускорения $V(\xi )$ основания превышает верхнюю величину u0 абсолютного ускорения защищаемого объекта:Если $V(0) > {{u}_{0}}$ или $V(T) > {{u}_{0}}$, то один или оба из интервалов $0 \leqslant \xi < {{t}_{1}}$, ${{t}_{2}} < \xi \leqslant T$ могут быть пустыми. Если функция $V(\xi ),\;0 \leqslant \xi \leqslant T,$ предполагается известной и удовлетворяет в каждый момент $\xi $ ограничению $V(\xi ) \leqslant {{u}_{0}}$, оптимальное управление задается тождеством $u(t) \equiv V(t - {{t}_{0}})$ и обеспечивает нулевое смещение объекта по отношению к основанию. Далее предполагаем выполнение неравенства
Полагаем также, что заранее известны и заданы слдующие характеристики ударного воздействия: его длительность T, максимально возможная величина воздействия b и параметр ${{{v}}_{0}}$
характеризующий значение скорости, приобретенной (или потерянной) основанием в результате удара. Класс возмущений без ограничений на максимальную величину воздействия $V(\xi )$ с известными и заданными параметрами T, ${{{v}}_{0}}$ обозначим ${{V}_{T}}.$ Обозначим ${{V}_{\infty }}$ класс возмущений с заданной величиной ${{{v}}_{0}}$ из (2.4) без ограничений на максимальную величину воздействия $V(\xi )$ и на максимальную длительность T. Классам воздействий ${{V}_{T}}$, ${{V}_{\infty }}$ соответствует значение $b = + \infty .$ Классы возмущений ${{V}_{b}}$ и ${{V}_{{bT}}}$ определим параметрами b, ${{{v}}_{0}}$ и T, b, ${{{v}}_{0}}$ соответственно. Из неравенств $0 \leqslant V(\xi ) \leqslant b$ и условия (2.4) следует выполнение соотношений т.е. неравенство $bT \geqslant {{{v}}_{0}}$ есть необходимое условие существования возмущений из класса ${{V}_{{bT}}}$. Случай $bT = {{{v}}_{0}}$ соответствует прямоугольным возмущениям и исследован в [9], поэтому далее предполагаем выполнение неравенстваВведенные классы возмущений удовлетворяют включениям ${{V}_{{bT}}} \subset {{V}_{T}} \subset {{V}_{\infty }}$, ${{V}_{{bT}}} \subset {{V}_{b}} \subset {{V}_{\infty }}.$
3. Критерий качества. Далее полагаем, что известно множество $\Omega \subset {{V}_{\infty }}$, которому принадлежат возможные возмущения $V(\xi )$, но само возмущение может быть неизвестно. При заданных управлении $u(t)$ и времени упреждения ${{t}_{0}}$ качество изоляции будем оценивать функционалом J, характеризующим наибольшую величину смещения объекта относительно основания при наихудшем возмущении:
(3.1)
$J(u,{\kern 1pt} \;{{t}_{0}}) = \mathop {max}\limits_{V \in \Omega } \,\mathop {max}\limits_{t \in {{R}_{{0 + }}}} |{\kern 1pt} x(t;\,\,u,V,\,\,{{t}_{0}}){\kern 1pt} |,\quad {{R}_{{0 + }}} = [0,{\kern 1pt} + \infty )$(3.2)
${{J}_{b}}(u,{\kern 1pt} \,{{t}_{0}}) = \mathop {\,max}\limits_{V \in {{V}_{b}}} \mathop {\,max}\limits_{t \in {{R}_{{0 + }}}} {\text{|}}x(t;u,V,\;{{t}_{0}}){\text{|}}$(3.3)
${{J}_{{bT}}}(u,{\kern 1pt} \,{{t}_{0}}) = {\kern 1pt} \mathop {\,max}\limits_{V \in {{V}_{{bT}}}} \mathop {\;max}\limits_{t \in {{R}_{{0 + }}}} |{\kern 1pt} x(t;u,V,\,\,{{t}_{0}}){\kern 1pt} |$(3.4)
${{J}_{T}}(u,\;{\kern 1pt} {{t}_{0}}) = \mathop {\;max}\limits_{V \in {{V}_{T}}} \mathop {\;max}\limits_{t \in {{R}_{{0 + }}}} |{\kern 1pt} x(t;u,V,{{t}_{0}}){\kern 1pt} |$4. Задачи оптимизации. Учитывая различную возможную информированность управляющей стороны о внешнем возмущении и различие допустимых законов управления U, поставим следующие задачи.
Задача 1. Для системы (1.1) с начальными условиями (1.2) построить кусочно-непрерывное управление $u_{{}}^{*}$, удовлетворяющее ограничению (1.3), и время упреждения $t_{0}^{*}$, доставляющие минимум величине (3.1):
Это задача о гарантирующем управления u* с моментом упреждения ${\kern 1pt} t_{0}^{*}$, которые позволяют защищать объект от ударов из множества $\Omega $. Задача 1 обобщает задачу из [1–3] с заданным возмущением и отсутствием упреждения.
Задача 2. Для системы (1.1) при начальных условиях (1.2) и заданном возмущении (2.1) построить кусочно-непрерывное управление u*, удовлетворяющее ограничению (1.3), и время упреждения $t_{0}^{*}$, минимизирующие наибольшую величину смещения объекта относительно основания:
Задачу 2 можно рассматривать как задачу 1, в которой множество допустимых возмущений $\Omega $ содержит лишь один элемент.
В качестве множества допустимых законов управления может быть рассмотрено параметрическое семейство управлений ${{u}_{s}}(t) \in {{U}_{S}}$, зависящих от параметра $s,s \in S.$
Задача 3. Для заданного класса допустимых управлений ${{u}_{s}}(t) \in {{U}_{S}}$ найти время упреждения $t_{0}^{*}$ и значение параметра s*, минимизирующие величину (3.1):
Решение задачи 3 позволяет улучшать качество противоударной защиты с помощью изменения времени упреждения ${{t}_{0}}$ и параметра s при заданном семействе законов управления ${{u}_{s}}(t)$, которое, например, может быть построено на основе оптимального управления для некоторого заданного возмущения $V \in \Omega $ и представлено аналитическими выражениями. При заданном управлении $u(t)$ с условием $u(t) \equiv 0,\;t < 0$ параметрическое семейство управлений может состоять из функций ${{u}_{c}}(t) = u(t - c),\;c \geqslant 0$, и тогда параметр $c$ играет роль запаздывания для управления ${{u}_{c}}(t)$. Заметим, что, если множество $S$ состоит из одного элемента, задача 3 решается для фиксированного управления.
5. Определение наихудшего возмущения. Вычисление функционала (3.1) предполагает процедуру построения наихудшего возмущения $V \in \Omega \subset {{V}_{T}}$, максимизирующего максимум модуля отклонения защищаемого объекта относительно основания (${\text{ma}}{{{\text{x}}}_{t}}{\text{|}}{\kern 1pt} x(t;u,V,{{t}_{0}}){\kern 1pt} |$) при заданных $u(t)$ и ${{t}_{0}}$.
Лемма 1 (см. [9]). При заданных управлении $u(t)$ и времени упреждения ${{t}_{0}}$ среди возмущений $V \in {{V}_{T}}$ наихудшее возмущение представляет собой мгновенный удар интенсивности ${{{v}}_{0}}$, подаваемый в начальный момент ${{t}_{0}}$ или в конечный момент ${{t}_{0}} + T$ допустимого интервала возмущения. Таким образом, наихудшее возмущение есть либо ${{V}_{\delta }}(\xi ) = {{{v}}_{0}}{\kern 1pt} \delta {\kern 1pt} {\kern 1pt} (\xi )$ либо ${{V}_{{\delta T}}}(\xi ) = {{{v}}_{0}}\delta {\kern 1pt} {\kern 1pt} (\xi - T)$.
Следуя лемме 1, для нахождения наихудшего возмущения в классе ${{V}_{T}}$ при заданных $u(t)$ и ${{t}_{0}}$ необходимо решить дифференциальное уравнение (1.1) с начальными условиями (1.2) при ${v}(t) = {{{v}}_{0}}{\kern 1pt} \delta {\kern 1pt} {\kern 1pt} (t - {{t}_{0}})$ и ${v}(t) = {{{v}}_{0}}\delta {\kern 1pt} {\kern 1pt} (t - {{t}_{0}} - T)$, затем для каждого из решений вычислить $ma{{x}_{t}}{\kern 1pt} |{\kern 1pt} x(t;u,{{V}_{\delta }},{{t}_{0}}){\kern 1pt} |$, $ma{{x}_{t}}{\kern 1pt} |{\kern 1pt} x(t;u,{{V}_{{\delta {\kern 1pt} T}}},{{t}_{0}}){\kern 1pt} |$ и выбрать возмущение, отвечающее большему значению абсолютного отклонения.
Лемма 2 (см. [14]). При заданных управлении $u(t)$ и времени упреждения ${{t}_{0}}$ среди возмущений $V \in {{V}_{{bT}}}$ наихудшее возмущение есть либо
Это означает, что при любом заданном управлении наихудшими возмущениями являются прямоугольные возмущения наибольшей амплитуды, либо начинающиеся в момент начала возмущения ${{t}_{0}}$, либо заканчивающиеся в момент окончания допустимого интервала возмущения ${{t}_{0}} + T$.
Примечание. Справедливо равенство ${{V}_{1}}(\xi ) = {{V}_{0}}(\xi - T + {{{v}}_{0}}{\text{/}}b).$
Следуя лемме 2, для нахождения наихудшего возмущения в классе ${{V}_{{bT}}}$ при заданных $u(t)$ и ${{t}_{0}}$ необходимо решить дифференциальное уравнение (1.1) с начальными условиями (1.2) при ${v}(t) = {{V}_{0}}(t - {{t}_{0}})$ и ${v}(t) = {{V}_{1}}(t - {{t}_{0}})$, затем для каждого из решений вычислить $ma{{x}_{t}}{\kern 1pt} {\kern 1pt} |{\kern 1pt} x(t;u,{{V}_{0}},{{t}_{0}}){\text{|}}$, $ma{{x}_{t}}{\kern 1pt} {\kern 1pt} |{\kern 1pt} x(t;u,{{V}_{1}},{{t}_{0}}){\text{|}}$, сравнить эти величины и выбрать возмущение, отвечающее большему значению. Таким образом, критерий качества J при заданном управлении $u$ для наихудших возмущений определяется выражением
(5.1)
${{J}_{u}} = \max \left\{ {\mathop {max}\limits_{t \in {{R}_{{0 + }}}} {\text{|}}x(t;u,{{V}_{0}},{{t}_{0}}){\text{|}},\mathop {\max }\limits_{t \in {{R}_{{0 + }}}} {\text{|}}x\left( {t;u,{{V}_{1}},{{t}_{0}}} \right){\kern 1pt} |} \right\}$Определим безразмерные переменные
Далее при использовании безразмерных переменных опускаем штрихи. В безразмерных единицах имеем ${{u}_{0}} = 1$ и ${{{v}}_{0}} = 1$, а неравенства (2.3), (2.5) приобретают следующий вид
6. Решение задачи 3 для класса ограниченных возмущений при использовании управления для прямоугольного возмущения с дополнительным запаздыванием. Возмущение в виде прямоугольного импульса в размерных переменных задается кусочно-постоянной функцией $V(t)$
(6.1)
$V = V(t) = \left\{ {\begin{array}{*{20}{l}} {b,\quad 0 \leqslant t \leqslant \frac{{{{{v}}_{0}}}}{b}} \\ {0,\quad t > \frac{{{{{v}}_{0}}}}{b}} \end{array}} \right.$(6.2)
$V = V(t) = \left\{ {\begin{array}{*{20}{l}} {b,\quad 0 \leqslant t \leqslant \frac{1}{b}} \\ {0,\quad t > \frac{1}{b}} \end{array}} \right.$Для возмущения (6.2) задача 2 решена в [9] с использованием графо-аналитического метода [1, 2]. Оптимальные значения управления $u(t)$, времени упреждения ${{t}_{0}}$ и критерия качества J задаются соотношениями
(6.3)
${{u}_{r}}(t) = \left\{ \begin{gathered} - 1,\quad 0 \leqslant t \leqslant t_{1}^{r} = \frac{1}{4}\sqrt {1 - {{b}^{{ - 1}}}} ,\quad b > 1 \hfill \\ 1,\quad t_{1}^{r} < t \leqslant t_{2}^{r} = \frac{1}{2}\sqrt {1 - {{b}^{{ - 1}}}} {\kern 1pt} + 1,\quad b > 1 \hfill \\ 0,\quad t > t_{2}^{r},\quad b > 1 \hfill \\ V(t),\quad t \in [0, + \infty ),\quad b \leqslant 1 \hfill \\ \end{gathered} \right.$(6.4)
$t_{0}^{b} = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{2}[(1 - {{b}^{{ - 1}}}) + \sqrt {1 - {{b}^{{ - 1}}}} ],\quad b > 1} \\ {0,\quad b \leqslant 1} \end{array}} \right.$(6.5)
${{J}_{r}} = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{{16}}(1 - {{b}^{{ - 1}}}),\quad b > 1} \\ {0,\quad b \leqslant 1} \end{array}} \right.$Управление из (6.3) назовем r-управлением. Зададим класс управлений ${{U}_{{rC}}} = \{ {{u}_{{rc}}}\} $ параметрическим семейством допустимых управлений ${{u}_{{rc}}}(t)$, где
(6.6)
$\begin{array}{*{20}{c}} {{{u}_{{rc}}}(t) = \left\{ \begin{gathered} 0,\quad t < c \hfill \\ - 1,\quad c \leqslant t \leqslant c + t_{1}^{r} = c + \frac{1}{4}\sqrt {1 - {{b}^{{ - 1}}}} \hfill \\ 1,\quad c + t_{1}^{r} < t \leqslant c + t_{2}^{r} = c + \frac{1}{2}\sqrt {1 - {{b}^{{ - 1}}}} + 1 \hfill \\ 0,\quad t > c + t_{2}^{r} \hfill \\ \end{gathered} \right.} \\ {c \in C,\quad C = \left[ {0,{\kern 1pt} + \infty } \right)} \end{array}$Равенство ${{u}_{{rc}}}(t) = {{u}_{r}}(t - c)$ означает, что управление ${{u}_{{rc}}}(t)$ представляет собой r-управление с запаздыванием c.
В [14] на основе леммы 2 в классах возмущений ${{V}_{{bT}}}$ и управлений ${{U}_{{rC}}}$ получено решение задачи 3:
(6.7)
$\begin{array}{*{20}{c}} {{{t}_{c}} = \left\{ \begin{gathered} c - 1 - T + \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2} + \frac{2}{b} + \frac{{\sqrt {9{{b}^{2}}{\kern 1pt} + 8{{b}^{2}}T - 26b - 8Tb + 17} }}{{2b}},\quad \frac{1}{b} < T \leqslant \frac{7}{8} + \frac{1}{{8b}} \hfill \\ c + \frac{9}{{16}} - \frac{1}{{16b}} + \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2} - \frac{T}{2},\quad \frac{7}{8} + \frac{1}{{8b}} < T \hfill \\ \end{gathered} \right.} \\ {{{c}^{{}}} \geqslant max\left[ {0,\frac{T}{2} - \frac{9}{{16}} + \frac{1}{{16{\kern 1pt} b}} - \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2}} \right]} \end{array}$(6.8)
$\begin{array}{*{20}{c}} {{{J}_{R}}(b,T) = \left\{ {\begin{array}{*{20}{l}} {\frac{{25}}{{16}} - \frac{{41}}{{16{\kern 1pt} b}} + T - \frac{{\sqrt {9{{b}^{2}} + 8{{b}^{2}}T{\kern 1pt} - 26{\kern 1pt} b - 8Tb + 17} }}{{2b}},}&{\frac{1}{b} < T \leqslant \frac{7}{8} + \frac{1}{{8b}}} \\ {\frac{T}{2} - \frac{1}{{2{\kern 1pt} b}},}&{\frac{7}{8} + \frac{1}{{8{\kern 1pt} b}} < T} \end{array}} \right.} \end{array}$Здесь ${{J}_{R}}(b,T) \equiv {{J}_{{bT}}}({{u}_{{rc}}},{{t}_{c}})$, т.е. искомое значение функционала ${{J}_{{bT}}}({{u}_{{rc}}},{{t}_{c}})$ задается функцией ${{J}_{R}}(b,T)$, у которой аргументы $b,T$ указывают на зависимость значения функционала ${{J}_{{bT}}}({{u}_{{rc}}},{{t}_{c}})$ от этих параметров. Оптимальное управление в классах управлений ${{U}_{{rC}}} = \{ {{u}_{{rc}}}\} $ и возмущений ${{V}_{{bT}}}$ для задачи 3 определяется формулой (6.6) со значением параметра c, удовлетворяющему неравенству из (6.7). Значения ${{t}_{c}},c$ из (6.7), доставляющие значение функционалу ${{J}_{R}}(b,T)$ из (6.8), не единственны. Оптимальные минимально возможные момент упреждения $t_{c}^{*}$ и величина запаздывания c* определяются как
(6.9)
$\begin{gathered} t_{c}^{*} = \left\{ \begin{gathered} \frac{{\sqrt {9{\kern 1pt} {\kern 1pt} {{b}^{2}}{\kern 1pt} + 8{{b}^{2}}T{\kern 1pt} - 26b - 8{\kern 1pt} Tb{\kern 1pt} + 17} }}{{2b}} + \frac{2}{b} + \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2} - 1 - {\kern 1pt} \,T,\quad \frac{1}{b} < \,T \leqslant \frac{7}{8} + \frac{1}{{8b}} \hfill \\ \frac{9}{{16}} - \frac{1}{{16b}} + \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2} - \frac{T}{2},\quad \frac{7}{8} + \frac{1}{{8{\kern 1pt} b}} \leqslant T < \frac{9}{8} - \frac{1}{{8b}} + \sqrt {1 - {{b}^{{ - 1}}}} \hfill \\ 0,\quad T \geqslant \frac{9}{8} - \frac{1}{{8b}} + \sqrt {1 - {{b}^{{ - 1}}}} \hfill \\ \end{gathered} \right. \hfill \\ c_{{}}^{*} = \left\{ \begin{gathered} 0,\quad \frac{1}{b} < T \leqslant \frac{9}{8} - \frac{1}{{8b}} + \sqrt {1 - {{b}^{{ - 1}}}} \hfill \\ \frac{T}{2} - \frac{9}{{16}} + \frac{1}{{16b}} - \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2},\quad \frac{9}{8} - \frac{1}{{8b}} + \sqrt {1 - {{b}^{{ - 1}}}} < T \hfill \\ \end{gathered} \right. \hfill \\ \end{gathered} $При значениях $t_{c}^{{}} = t_{c}^{*}$, c = c* из (6.9) оптимальное управление определяется формулой $u_{{rc}}^{*}(t) = {{u}_{r}}(t - c{\text{*}})$ и справедливо равенство $\min \{ t_{c}^{*},c{\text{*}}\} = 0$, т.е. имеет место либо упреждение, либо запаздывание. При не минимально выбранных моменте упреждения $t_{c}^{{}}$ и величине запаздывания c имеем неравенство $\min \{ t_{c}^{{}},c\} > 0$. В этом случае на интервале времени $0 \leqslant t < \min \{ t_{c}^{{}},c\} $ рассматриваемая система бездействует и значения управления и возмущения имеют нулевые значения.
7. Решение задачи 3 для класса ограниченных возмущений при использовании управления, предназначенного для возмущений ограниченной длительности. В [10, 11] на основе леммы 1 получено решение $\{ {{u}_{\tau }},{{t}_{\tau }}\} $ задачи 3 в классе неограниченных возмущений ${{V}_{T}}$ для управлений из класса ${{U}_{\Sigma }} = \{ {{u}_{\tau }}\} $ – параметрического семейства допустимых релейных управлений ${{u}_{\tau }}(t)$ с переключением cо значения $ - {{u}_{0}}$ на значение ${{u}_{0}}$ в момент времени $\tau $ и cо значения ${{u}_{0}}$ на значение 0 в момент времени ${{{v}}_{0}}{\kern 1pt} /{\kern 1pt} {{u}_{0}} + 2\tau $, где
(7.1)
$\begin{array}{*{20}{c}} {{{u}_{\tau }}(t) = \left\{ {\begin{array}{*{20}{l}} { - {{u}_{0}},0 \leqslant t < \tau } \\ { + {{u}_{0}},\tau \leqslant t \leqslant {{T}_{\tau }}} \\ {0,t > {{T}_{\tau }}} \end{array}} \right.} \\ {{{T}_{\tau }} = \frac{{{{{v}}_{0}}}}{{{{u}_{0}}}} + 2\tau ,\quad \tau \in \Sigma ,\quad \Sigma = [0, + \infty )} \end{array}$Временные отрезки со значениями $ \pm {{u}_{0}}$ выбраны таким образом, чтобы обеспечить выполнение равенства
В безразмерных переменных управление ${{u}_{\tau }}(t)$ из (7.1) имеет вид
(7.2)
$\begin{array}{*{20}{c}} {{{u}_{\tau }}(t) = \left\{ {\begin{array}{*{20}{l}} { - 1,\quad 0 \leqslant t < \tau } \\ { + 1,\quad \tau \leqslant t \leqslant {{T}_{\tau }}} \\ {{\kern 1pt} 0,\quad t > {{T}_{\tau }}} \end{array}} \right.} \\ {{{T}_{\tau }} = 1 + 2\tau ,\quad \tau \in \Sigma ,\quad \Sigma = [0, + \infty )} \end{array}$Оптимальное управление в классах управлений ${{U}_{\Sigma }} = \{ {{u}_{\tau }}\} $ и возмущений VT для задачи 3 определяется формулой (7.2) со значением параметра $\tau $ ([10, 11]):
(7.3)
$\tau {\text{*}} = \left\{ {\begin{array}{*{20}{l}} {T{\text{/}}2 + 1{\text{/}}4,}&{T \leqslant 1{\text{/}}2} \\ {1{\text{/}}2,}&{1{\text{/}}2 < T \leqslant 7{\text{/}}2} \\ {\sqrt {(1 + T){\kern 1pt} {\text{/}}2} - 1,}&{T > 7{\text{/}}2} \end{array}} \right.$(7.4)
$\begin{gathered} {{J}_{T}} = \left\{ {\begin{array}{*{20}{l}} {{{{\left( {T{\text{/}}2 + 1{\kern 1pt} {\text{/}}4} \right)}}^{2}},}&{\;T \leqslant 1{\text{/}}2} \\ {T{\text{/}}2,}&{\;1{\text{/}}2 < T} \end{array}} \right. \hfill \\ t_{\tau }^{*} = \left\{ {\begin{array}{*{20}{l}} {T + 1,}&{\;T \leqslant 1{\text{/}}2} \\ {7{\text{/}}4 - T{\text{/}}2,}&{\;1{\text{/}}2 < T \leqslant 7{\text{/}}2} \\ {0,}&{\;T > 7{\text{/}}{\kern 1pt} 2} \end{array}} \right. \hfill \\ \end{gathered} $Решение (7.2)–(7.4) задачи 3 при $T > 1/2$ не единственно. Значение функционала ${{J}_{T}} = T{\kern 1pt} /{\kern 1pt} 2$ при $T > 1/2$ обеспечивается значениями $\tau $ и ${{t}_{\tau }}$, удовлетворяющими соотношениям
Обозначим через $u_{\tau }^{*}(t)$ функцию из (7.2) со значением $\tau = \tau _{{}}^{*}$ из (7.3) и используем функцию $u_{\tau }^{*}(t)$ и значение $t_{\tau }^{*}$ из (7.4) для вычисления значения функционала ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$ из (3.3) в классе ограниченных возмущений ${{V}_{{bT}}}$. Определим функцию ${{J}_{{b\tau }}}\left( {t{\text{*}}} \right)$ соотношениями
Отсюда получаем решение задачи 3 в классе ограниченных возмущений ${{V}_{{bT}}}$ для случая, когда класс допустимых управлений состоит из одного элемента $u_{\tau }^{*}$ и множество выбора моментов упреждения ${{t}_{0}}$ также состоит из одного элемента $t_{\tau }^{*}$.
Введем функцию ${{J}_{\Sigma }}(b,T) \equiv {{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$, у которой аргументы $b,T$ указывают на зависимость значения функционала ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$ от этих параметров.
Значение ${{J}_{\Sigma }}(b,T)$ для всех $1 < b \leqslant 2$ определяется выражением
(7.5)
$\begin{array}{*{20}{c}} {{{J}_{\Sigma }}(b,T) = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{4},\quad }&{\frac{1}{2} \leqslant \frac{1}{b} < T \leqslant \frac{1}{b} + \frac{1}{2}} \\ {\frac{T}{2} - \frac{1}{{2b}},\quad }&{\frac{1}{b} + \frac{1}{2} < T} \end{array}} \right.} \end{array}$Значение ${{J}_{\Sigma }}(b,T)$ для всех $b > 2$ определяется более сложным выражением
(7.6)
$\begin{array}{*{20}{c}} {{{J}_{\Sigma }}(b,T) = \left\{ {\begin{array}{*{20}{l}} {{{{\left( {\frac{T}{2} + \frac{1}{4}} \right)}}^{2}},\quad }&{\frac{1}{b} < T \leqslant \frac{1}{2}} \\ {\frac{1}{4},\quad }&{\frac{1}{2} < T \leqslant \frac{1}{b} + \frac{1}{2}} \\ {\frac{T}{2} - \frac{1}{{2b}},\quad }&{\frac{1}{b} + \frac{1}{2} < T} \end{array}} \right.} \end{array}$8. Решение задачи 3 для возмущений ограниченной длительности из класса VT и ограниченных возмущений из класса ${{V}_{{bT}}}$ при использовании постоянных управлений. Зададим класс управлений ${{U}_{D}} = \{ {{u}_{d}}\} $ – параметрическим семейством допустимых управлений ${{u}_{d}}(t)$, где ${{u}_{d}}(t)$ в размерных и безразмерных переменных задается соответственно соотношениями (8.1) и (8.2)
(8.1)
$\begin{array}{*{20}{c}} {{{u}_{d}}(t) = \left\{ {\begin{array}{*{20}{l}} {0,\;}&{t < d} \\ {{{u}_{0}},\;}&{d\; \leqslant \;t\; \leqslant \;\frac{{{{u}_{0}}}}{{{{{v}}_{0}}}} + d} \\ {0,\;}&{t > \frac{{{{u}_{0}}}}{{{{{v}}_{0}}}} + d} \end{array}} \right.} \\ {d\; \geqslant \;0} \end{array}$(8.2)
$\begin{array}{*{20}{c}} {{{u}_{d}}(t) = \left\{ {\begin{array}{*{20}{l}} {0,\;}&{t < d} \\ {1,\;}&{d \leqslant t \leqslant 1 + d} \\ {0,\;}&{t > 1 + d} \end{array}} \right.} \\ {d \geqslant 0} \end{array}$Управление ${{u}_{d}}(t)$ при d = 0 обозначим ${{u}_{{0d}}}(t)$. Управление ${{u}_{{0d}}}(t)$ с моментом упреждения ${{t}_{0}} = 0$ представляет собой оптимальное управление для задачи 2 без упреждения при возмущении (см. [6]). Управление ${{u}_{d}}(t)$ есть ${{u}_{{0d}}}$-управление с запаздыванием d, так как ${{u}_{d}}(t) = {{u}_{{0d}}}(t - d)$.
Решение задачи 3 в классе неограниченных возмущений ${{V}_{T}}$ для постоянных управлений класса UD. В [13] на основе леммы 1 получено решение $\{ {{u}_{d}},{{t}_{d}}\} $ задачи 3 в классе неограниченных возмущений ${{V}_{T}}$ для управлений класса UD. Это решение в указанных классах для задачи 3 определяется формулой для управления (8.2) и выражениями для момента упреждения ${{t}_{d}}$ и значения функционала ${{J}_{d}}$:
(8.3)
${{t}_{d}} = \left\{ {\begin{array}{*{20}{l}} {d + \sqrt {2T + 2} - T - 1,}&{T \leqslant 1,} \\ {d + 1{\text{/}}2 - T{\text{/}}2,}&{1 < T,} \end{array}} \right.\quad d \geqslant max\left[ {0,\frac{{T - 1}}{2}} \right]$(8.4)
$J({{u}_{d}},{{t}_{d}}) = \left\{ {\begin{array}{*{20}{l}} {3{\kern 1pt} {\text{/}}2 - \sqrt {2T + 2} + T,}&{T \leqslant 1} \\ {T{\text{/}}2,}&{1 < T} \end{array}} \right.$Значение функционала ${{J}_{d}}(T)$ из (8.4) обеспечивается переменными ${{t}_{d}}$, d из (8.3) не единственным образом. При этом профиль управления ${{u}_{d}}(t)$ из (8.2) остается одинаковым, меняется лишь начало его действия. Оптимальные минимально возможные момент упреждения $t_{d}^{*}$ и величина запаздывания d* определяются соотношениями
(8.5)
$t_{d}^{*} = \left\{ {\begin{array}{*{20}{l}} {\sqrt {2T + 2} - T - 1,}&{T < 1,} \\ 0&{T \geqslant 1,} \end{array}} \right.\quad d_{{}}^{*} = \left\{ {\begin{array}{*{20}{l}} {0,}&{T \leqslant 1} \\ {\frac{{T - 1}}{2},}&{1 < T} \end{array}} \right.$При значениях $t_{d}^{{}} = t_{d}^{*}$, $d = d_{{}}^{*}$ из (8.5) оптимальное управление определяется формулой $u_{d}^{*}(t) = {{u}_{{0d}}}\left( {t - d{\text{*}}} \right)$ и справедливо равенство $\min \{ t_{d}^{*},d{\text{*}}\} = 0$, т.е. имеет место либо упреждение, либо запаздывание. При не минимально выбранных моменте упреждения $t_{d}^{{}}$ и величине запаздывания d имеем неравенство $\min \{ {{t}_{d}},d\} > 0$. В этом случае на интервале времени $0 \leqslant t < \min \{ {{t}_{d}},\,d\} $ рассматриваемая система бездействует и значения управления и возмущения имеют нулевые значения.
Решение задачи 3 в классе ограниченных возмущений ${{V}_{{bT}}}$ для постоянных управлений класса UD. При описании этого решения будем использовать обозначения: ${{u}_{{bd}}}$ для элементов UD, ${{t}_{{bd}}}$ для момента упреждения, ${{d}_{b}}$ для величины запаздывания, поскольку при описании решения задачи 3 эти величины и их области изменения могут зависеть от параметра b. Допустимые управления ${{u}_{{b0}}}$, ${{u}_{{bd}}}$ в безразмерных переменных задаются соотношениями (8.6)
(8.6)
$\begin{array}{*{20}{c}} {{{u}_{{b0}}}(t) = \left\{ {\begin{array}{*{20}{l}} {0,\;}&{t < 0,} \\ {1,\;}&{0 \leqslant t \leqslant 1,} \\ {0,\;}&{t > 1,} \end{array}} \right.} \\ {{{d}_{b}} = 0,} \end{array}\quad \quad \begin{array}{*{20}{c}} {{{u}_{{bd}}}(t) = \left\{ {\begin{array}{*{20}{l}} {0,\;}&{t < {{d}_{b}}} \\ {1,\;}&{{{d}_{b}} \leqslant t \leqslant 1 + {{d}_{b}}} \\ {0,\;}&{t > 1 + {{d}_{b}}} \end{array}} \right.} \\ {{{d}_{b}} \geqslant 0} \end{array}$Управление ${{u}_{{bd}}}(t)$ есть ${{u}_{{b0}}}$-управление с запаздыванием db, так как ubd(t) = ${{u}_{{b0}}}(t - {{d}_{b}})$.
Определим функции ${{J}_{{b0}}}\left( {t{\text{*}}} \right),{{J}_{{bd}}}\left( {t{\text{*}}} \right)$ выражениями
Эти функции определяются следующими соотношениями
(8.7)
$\begin{array}{*{20}{c}} {{{J}_{{bd}}}\left( {t{\text{*}}} \right) = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{2}(1 - {{b}^{{ - 1}}}) - t_{{}}^{*} + {{d}_{b}},\quad 0 \leqslant t_{{}}^{*} < (\sqrt 2 - 1)(1 - {{b}^{{ - 1}}}) + {{d}_{b}}} \\ {\frac{{b{{{\left( {t{\text{*}} - {{d}_{b}}} \right)}}^{2}}}}{{2(b - 1)}},\quad (\sqrt 2 - 1)(1 - {{b}^{{ - 1}}}) + {{d}_{b}} \leqslant t_{{}}^{*} < (1 - {{b}^{{ - 1}}}) + {{d}_{b}}} \\ {t_{{}}^{*} - {{d}_{b}} - \frac{1}{2}(1 - {{b}^{{ - 1}}}),\quad t_{{}}^{*} \geqslant (1 - {{b}^{{ - 1}}}) + {{d}_{b}}} \end{array}} \right.} \end{array}$Функции ${{J}_{{b0}}}\left( {t{\text{*}}} \right),{{J}_{{bd}}}\left( {t{\text{*}}} \right)$ отвечают наибольшему относительному смещению защищаемого объекта при использовании соответственно управлений ${{u}_{{b0}}}(t),\;{{u}_{{bd}}}(t)$ и применении возмущения ${v}(t) = {{V}_{0}}\left( {t - t{\text{*}}} \right)$ – прямоугольного возмущения с амплитудой b и началом действия в момент времени $t_{{}}^{*}$. Очевидное равенство ${{J}_{{bd}}}\left( {t{\text{*}}} \right) = {{J}_{{b0}}}\left( {t{\text{*}} - {{d}_{b}}} \right)$ при $t_{{}}^{*} \geqslant {{d}_{b}}$ означает, что график функции ${{J}_{{bd}}}\left( {t{\text{*}}} \right)$ получается смещением графика ${{J}_{{b0}}}\left( {t{\text{*}}} \right)$ на величину ${{d}_{b}}$ вдоль оси абсцисс с доопределением на интервале $[0,{{d}_{b}})$ согласно формуле для ${{J}_{{bd}}}\left( {t{\text{*}}} \right).$ При $b \to + \infty $ функция ${{J}_{{b0}}}\left( {t{\text{*}}} \right) \to {{J}_{{d0}}}\left( {t{\text{*}}} \right)$, где функция ${{J}_{{d0}}}\left( {t{\text{*}}} \right)$ предназначена для дельта возмущения и имеет тот же смысл, что и функция ${{J}_{{b0}}}\left( {t{\text{*}}} \right)$ для ограниченного возмущения. Функция ${{J}_{{d0}}}\left( {t{\text{*}}} \right)$ задается формулой
Функция ${{J}_{{bd}}}\left( {t{\text{*}}} \right)$ линейно убывает при $t{\text{*}} \in [0,(\sqrt 2 - 1)(1 - {{b}^{{ - 1}}}) + {{d}_{b}}]$ от значения $(1 - {{b}^{{ - 1}}}){\text{/}}2$ до значения
(8.8)
${{J}_{{bd}}}((\sqrt 2 - 1)(1 - {{b}^{{ - 1}}}) + {{d}_{b}}) = \left( {\frac{3}{2} - \sqrt 2 } \right)(1 - {{b}^{{ - 1}}}) \approx 0.0858(1 - {{b}^{{ - 1}}})$Согласно лемме 2 величина ${{J}_{{bT}}}({{u}_{{bd}}},\,{{t}_{{bd}}})$ определяется значением выражения
Минимизируя выражение JbT(ubd, tbd) по переменной tbd получаем решение задачи 3 в классе ограниченных возмущений ${{V}_{{bT}}}$
(8.10)
$\begin{array}{*{20}{c}} {{{t}_{{bd}}} = \left\{ {\begin{array}{*{20}{l}} {{{d}_{b}} + \frac{{\sqrt {2(b - 1)(Tb + b - 2)} - Tb - b + 2}}{b},}&{\frac{1}{b} < T \leqslant 1} \\ {{{d}_{b}} + \frac{{1 - T}}{2},}&{1 < T} \end{array}} \right.} \\ {{{d}_{b}} \geqslant max\left[ {0,{\kern 1pt} \;\frac{{T - 1}}{2}} \right]} \end{array}$(8.11)
$\begin{array}{*{20}{c}} {{{J}_{D}}(b,T) = \left\{ {\begin{array}{*{20}{l}} {\frac{{3{\kern 1pt} b + 2Tb - 5 - 2{\kern 1pt} \sqrt {2(b - 1)(Tb + b - 2)} }}{b},}&{\frac{1}{b} < T \leqslant 1} \\ {\frac{T}{2} - \frac{1}{{2{\kern 1pt} b}},{\kern 1pt} }&{1 < T} \end{array}} \right.} \end{array}$Здесь ${{J}_{D}}(b,T) \equiv {{J}_{{bT}}}({{u}_{{bd}}},{{t}_{{bd}}})$, т.е. искомое значение функционала ${{J}_{{bT}}}({{u}_{{bd}}},{{t}_{{bd}}})$ задается функцией ${{J}_{D}}(b,T)$, у которой аргументы b, T указывают на зависимость значения функционала ${{J}_{{bT}}}({{u}_{{bd}}},{{t}_{{bd}}})$ от этих параметров.
Значение функционала ${{J}_{D}}(b,T)$ из (8.11) обеспечивается переменными ${{t}_{{bd}}}$, ${{d}_{b}}$ из (8.10) не единственным образом. При этом профиль управления ${{u}_{{bd}}}(t)$ из (8.6) остается одинаковым, меняется лишь начало его действия. Оптимальные минимально возможные момент упреждения $t_{{bd}}^{*}$ и величина запаздывания $d_{b}^{*}$ определяются соотношениями
(8.12)
$\begin{array}{*{20}{c}} {t_{{bd}}^{*} = \left\{ {\begin{array}{*{20}{l}} {\frac{{\sqrt {2(b - 1)(Tb + b - 2)} - Tb - b + 2}}{b},}&{\frac{1}{b} < T \leqslant 1} \\ {0,}&{1 < T} \end{array}} \right.} \\ {\quad d_{{}}^{*} = \left\{ {\begin{array}{*{20}{l}} {0,}&{T \leqslant 1} \\ {\frac{{T - 1}}{2},}&{1 < T} \end{array}} \right.} \end{array}$При значениях ${{t}_{{bd}}} = t_{{bd}}^{*}$, ${{d}_{b}} = d_{b}^{*}$ из (8.12) оптимальное управление определяется формулой $u_{{bd}}^{*}(t) = {{u}_{{b0}}}(t - d_{b}^{*})$ и справедливо равенство $min\{ t_{{bd}}^{*},d_{b}^{*}\} = 0$, т.е. имеет место либо упреждение, либо запаздывание. При не минимально выбранных моменте упреждения ${{t}_{{bd}}}$ и величине запаздывания db имеем неравенство $min\{ {{t}_{{bd}}},{{d}_{b}}\} > 0.$ В этом случае на интервале времени $0 \leqslant t < min\{ {{t}_{{bd}}},{{d}_{b}}\} $ рассматриваемая система бездействует и значения управления и возмущения имеют нулевые значения.
Решение задачи 3 в классе ограниченных возмущений ${{V}_{{bT}}}$ при использовании постоянного управления без упреждения и запаздывания. В частном случае решения задачи 3 при заданных значениях параметров ${{t}_{{bd}}} = 0$, ${{d}_{b}} = 0$ и при использовании управления ${{u}_{{b0}}}$ из (8.6) получаем значение функционала
(8.13)
$\begin{array}{*{20}{c}} {{{J}_{0}}(b,T) = {{J}_{{bT}}}({{u}_{{b0}}},{\mkern 1mu} 0) = \;\max [{{J}_{{b0}}}(0),\;{{J}_{{b0}}}(T - {{b}^{{ - 1}}})] = } \\ { = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{2}\left( {1 - \frac{1}{b}} \right),\quad \frac{1}{b} < T \leqslant 1} \\ {T - \frac{1}{2} - \frac{1}{{2b}},\quad T > 1} \end{array}} \right.} \end{array}$При этом начало действия наихудшего возмущения определяется соотношениями
9. Сравнение по критерию качества различных способов управления. Значения функционалов ${{J}_{R}}(b,T),{{J}_{\Sigma }}(b,T),{{J}_{D}}(b,T),{{J}_{0}}(b,T)$ из (6.8), (7.5), (7.6), (8.11), (8.13) представляют собой функции от параметров $b,T$, заданные в области
При этом выполнены следующие тождества
Соответственно сравнение функций ${{J}_{R}}(b,T)$, ${{J}_{\Sigma }}(b,T)$, ${{J}_{D}}(b,T)$ можно проводить лишь при
Сравнение функций ${{J}_{R}}(b,T),\;{{J}_{\Sigma }}(b,T)$. Сравнение функций ${{J}_{R}}(b,T)$, ${{J}_{\Sigma }}(b,T)$ проведено в работе [14]. Вычисления показывают, что равенство ${{J}_{R}}(b,T) = {{J}_{\Sigma }}(b,T)$ достигается при $1{\kern 1pt} /{\kern 1pt} {\kern 1pt} b < T < 1$ на кривой $l$, описываемой функцией $b = f(T),\;0 < T \leqslant 1$, где
При этом кривая l стремится к плюс бесконечности при $T \to + 0$, лежит выше кривой $b = 1{\kern 1pt} /{\kern 1pt} T$ при $0 < T < 1$ и пересекается с ней в точке $T = 1,\;b = 1.$ Кривые l и $b = 1{\kern 1pt} /{\kern 1pt} T$ изображены сплошными линиями на рис. 3. Кривая l делит область ${{G}_{{bT}}}$ на две подобласти ${{G}_{R}},{{G}_{\Sigma }},$ располагающиеся ниже и выше кривой l соответственно. При этом в нижней области ${{G}_{R}}$ выполнено неравенство ${{J}_{R}}(b,T) < {{J}_{\Sigma }}(b,T)$, а в верхней области ${{G}_{\Sigma }}$ выполнено неравенство ${{J}_{R}}(b,T) > {{J}_{\Sigma }}(b,T).$ Зададим функцию ${{J}_{{comb}}}(b,T)$ соотношением
(9.2)
${{J}_{{comb}}}(b,T) = min\left\{ {{{J}_{R}}(b,T),{{J}_{\Sigma }}(b,T)} \right\},\quad (b,T) \in {{G}_{{bT}}}$Линия l1 имеет общую точку с кривой $l$:
а кривая l2 имеет общую точку с кривой $l$:В каждой из четырех областей, помеченных на рис. 3 значение ${{J}_{{comb}}}(b,T)$ из (9.2) вычисляется по своей формуле с соответствующим индексом. То есть области со значением индекса $i\;{\kern 1pt} (i = 1,2,3,4)$ соответствует своя функция ${{J}_{i}}(b,T)$. Функции ${{J}_{i}}(b,T)$ имеют следующий вид
(9.3)
$\begin{gathered} {{J}_{1}} = \frac{{25}}{{16}} - \frac{{41}}{{16b}} + T - \frac{{\sqrt {9{{b}^{2}} + 8{{b}^{2}}T - 26b - 8Tb + 17} }}{{2b}} \\ {{J}_{2}} = {{\left( {\frac{T}{2} + \frac{1}{4}} \right)}^{2}} \\ {{J}_{3}} = \frac{1}{4} \\ {{J}_{4}} = \frac{T}{2} - \frac{1}{{2b}} \\ \end{gathered} $Функция ${{J}_{{comb}}}(b,T)$ представляет собой значение функционала (3.3) при использовании в качестве управления и момента упреждения величин $u_{{rc}}^{*},\;t_{c}^{*}$ в области ${{G}_{R}}$ и при использовании в качестве управления и момента упреждения величин $u_{\tau }^{*},t_{\tau }^{*}$ в области ${{G}_{\Sigma }}$. Такой способ управления с моментом упреждения назовем комбинированным управлением.
На рис. 4 изображено разбиение области ${{G}_{{bT}}}$ на две подобласти, которые описывает формула (6.8). Цифры 1 и 4, которыми помечены эти области, означают что функционал ${{J}_{R}}(b,T)$ имеет значения ${{J}_{1}}$ и ${{J}_{4}}$ из (9.3) в соответствующих подобластях.
На рис. 5 изображено разбиение области ${{G}_{{bT}}}$ на три подобласти, которые фигурируют в формулах (7.5), (7.6). Цифры 2–4, которыми помечены эти подобласти, означают что функционал ${{J}_{\Sigma }}(b,T)$ имеет значения ${{J}_{2}}$, ${{J}_{3}}$, ${{J}_{4}}$ из (9.3) в соответствующих подобластях.
В [14] показано, что применение комбинированного управления приводит к улучшению значения функционала, доходящему до 10% по сравнению со значением функционала JR, полученным для управления $u_{{rc}}^{*}$ и момента упреждения $t_{c}^{*}$. Применение комбинированного управления приводит к улучшению значения функционала, доходящему до 700% по сравнению со значением функционала ${{J}_{\Sigma }}(b,T)$, полученным для управления $u_{\tau }^{*}$ и момента упреждения $t_{\tau }^{*}$.
Сравнение функций ${{J}_{R}}(b,T),\;{{J}_{D}}(b,T)$, ${{J}_{{comb}}}(b,T)$. Заметим, что все три функции ${{J}_{R}}(b,T)$, ${{J}_{D}}(b,T)$, ${{J}_{\Sigma }}(b,T)$ при фиксированном $b$ суть возрастающие функции T, а неотрицательная функция ${{J}_{D}}(b,T) - {{J}_{R}}(b,T)$ при $b > 36{\kern 1pt} /{\kern 1pt} 19$ убывает до нуля по T на полуинтервале $(1{\text{/}}b,1]$. Поэтому относительная разность $\gamma (b,T)$ значений ${{J}_{R}}(b,T),\;{{J}_{D}}(b,T)$ при $b > 36{\kern 1pt} /{\kern 1pt} 19$ удовлетворяет соотношениям
Численные расчеты максимальной величины $\gamma (b,T)$ в области $1{\kern 1pt} /{\kern 1pt} T < b \leqslant 36{\kern 1pt} /{\kern 1pt} 19$ также приводят к значению 0.372583, которое тоже достигается при $T = 1{\text{/}}b$ и составляет 37%.
Относительная разность $\mu (b,T)$ значений ${{J}_{{comb}}}(b,T)$, ${{J}_{D}}(b,T)$, задаваемая равенством
удовлетворяет предельным соотношениямЧисленный анализ показал, что при фиксированном b, удовлетворяющем неравенствам $1 < b \leqslant 7/3$, функция $\mu (b,T)$ на полуинтервале $\left( {1{\text{/}}b,1} \right]$ удовлетворяет соотношениям
Было показано, что при b, удовлетворяющем неравенствам $7{\kern 1pt} /{\kern 1pt} 3 < b$, функция $\mu (b,T)$ убывает по T на полуинтервале $(1{\kern 1pt} {\kern 1pt} /{\kern 1pt} {\kern 1pt} b,\;1{\kern 1pt} /{\kern 1pt} 2]$, затем возрастает по T на $(1{\text{/}}2,1{\text{/}}2 + 1{\text{/}}b]$, и затем на на полуинтервале $\left( {1{\text{/}}2 + 1{\text{/}}b,1} \right]$ убывает до нуля. При этом локальный максимум функции $\mu (b,T)$ при $T = 1{\kern 1pt} /{\kern 1pt} 2 + 1{\kern 1pt} /{\kern 1pt} b$ равен
Сравнение функций ${{J}_{0}}(b,T)$, ${{J}_{D}}(b,T)$. Относительная разность $\nu (b,T)$ функций ${{J}_{0}}(b,T)$, ${{J}_{D}}(b,T)$ удовлетворяет соотношениям
Заключение. Рассмотрена задача изоляции объекта на подвижном основании, которое подвержено однонаправленным возмущениям ограниченной длительности и амплитуды с заданным интегралом от возмущения. Для этой модели изучалась возможность применения наиболее простого, постоянного управления заданной длительности, полученного ранее для мгновенного удара для задачи без упреждения. Предложена оптимизация критерия качества по моменту начала действия предложенного управления. Проведено сравнение по значению критерия качества предлагаемой оптимизации с наиболее эффективными управлениями, предложенными автором в предыдущих работах.
Работа выполнена по теме государственного задания АААА-А20-120011690138-6.
Список литературы
Гурецкий В.В. Об одной задаче оптимального управления // Изв. АН СССР. Механика. 1965. № 1. С. 159–162.
Гурецкий В.В. О задаче минимизации максимального смещения // Труды ЛПИ. Механика и процессы управления. 1969. № 307. С. 11–21.
Sevin E., Pilkey W. Optimum Shock and Vibration Isolation. Washington DC: Shock and Vibration Information Analysis Center, 1971. 162 c.
Коловский М.З. Автоматическое управление виброзащитными системами. М.: Наука, 1976. 320 с.
Болотник Н.Н. Оптимизация амортизационных систем. М.: Наука, 1983. 256 с.
Balandin D.V., Bolotnik N.N., Pilkey W.D. Optimal Protection from Impact, Shock, and Vibration. Amsterdam: Gordon and Breach Science, 2001.
Pilkey W.D., Balandin D.V., Bolotnik N.N., Crandal J.R., Purtsezov S.V. Injury Biomechanics and Control: Optimal Protection from Impact. Hoboken, NJ: Wiley and Sons, 2010.
Болотник Н.Н., Корнеев В.А. Анализ предельных возможностей противоударной изоляции при кратковременных внешних воздействиях // Нелинейная динамика. 2015. Т. 11. № 1. С. 147–168. https://doi.org/10.20537/nd1501008
Болотник Н.Н., Корнеев В.А. Противоударная изоляция с упреждающим управлением для внешних возмущений различной формы // Изв. РАН. ТиСУ. 2018. № 3. С. 48–63. https://doi.org/10.7868/S0002338818030058
Болотник Н.Н., Корнеев В.А. Гарантирующее упреждающее управление в задаче противоударной изоляции // Доклады РАН. 2018. Т. 445. № 4. С. 381–385. https://doi.org/10.31857/S086956520001708-1
Корнеев В.А. Защита объекта на подвижном основании с помощью упреждающего управления при наихудших возмущениях // Изв. РАН. ТиСУ. 2019. № 1. С. 89–97. https://doi.org/10.1134/S0002338819010104
Корнеев В.А. Оптимизация управления с упреждением и запаздыванием в задаче противоударной защиты объекта на подвижном основании // Изв. РАН. ТиСУ. 2020. № 3. С. 30–38. https://doi.org/10.31857/S0002338820030075
Корнеев В.А. Использование постоянного управления с упреждением и запаздыванием в задаче защиты объекта от удара на подвижном основании // Изв. РАН. МТТ. 2020. № 2. С. 147–158. https://doi.org/10.31857/S0572329920020087
Корнеев В.А. Защита объекта на подвижном основании с помощью управления с упреждением и запаздыванием при ограниченных возмущениях // Изв. РАН. МТТ. 2021. № 4. С. 18–35. https://doi.org/10.31857/S0572329921030077
Ledezma-Ramirez D.F., Tapia-Gonzalez P.E., Ferguson N., Brennan M., Tang B. Recent advances in shock vibration isolation: an overview and future possibilities // Appl. Mech. Rev. 2019. V. 71. № 6. P. 060802 https://doi.org/10.1115/1.4044190
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Механика твердого тела