Известия РАН. Механика твердого тела, 2021, № 4, стр. 18-35

ЗАЩИТА ОБЪЕКТА НА ПОДВИЖНОМ ОСНОВАНИИ С ПОМОЩЬЮ УПРАВЛЕНИЯ С УПРЕЖДЕНИЕМ И ЗАПАЗДЫВАНИЕМ ПРИ ОГРАНИЧЕННЫХ ВОЗМУЩЕНИЯХ

В. А. Корнеев *

Институт проблем механики им. А.Ю. Ишлинского РАН
Москва, Россия

* E-mail: korneev@ipmnet.ru

Поступила в редакцию 03.07.2020
После доработки 13.09.2020
Принята к публикации 24.10.2020

Полный текст (PDF)

Аннотация

Рассматривается задача построения гарантирующего управления для противоударного изолятора, защищающего объект на подвижном основании от ударов, которым подвержено основание. Ударное воздействие на основание заранее неизвестно, задается его ускорением, ограниченной знакопостоянной функцией времени фиксированной длительности, интеграл от которой по времени задан. Управление действует между основанием и защищаемым объектом, ограничено по величине, а абсолютное ускорение основания может превышать эту величину только на одном интервале времени. Минимизируемым критерием качества служит максимальное смещение объекта относительно основания. Разработана методика определения среди возмущений указанного класса наихудших возмущений, доставляющих наибольшее значение критерию качества. Предложен эффективный закон управления с оптимизацией момента начала управления, основанный на комбинации двух методов управления, полученных ранее для прямоугольных возмущений и возмущений с неограниченной амплитудой.

Ключевые слова: противоударная изоляция, оптимальное управление, гарантирующее упреждающее управление, запаздывание, наихудшие возмущения

Введение. Работа посвящена построению способов управления противоударным изолятором при ограниченных наихудших возмущениях. Используется математическая модель системы с одной степенью свободы, которая состоит из основания и защищаемого объекта, находящегося на нем. Рассматривается внешнее воздействие кинематического типа, характеризуемое ускорением основания, описываемым функцией времени. Основание и защищаемый объект движутся вдоль одной прямой, а управляющая сила между основанием и объектом ограничена по величине и создается изолирующим устройством. Впервые задача о минимизации максимума модуля смещения объекта относительно основания при заданном возмущении ставилась в [13]. Дальнейшее развитие теория оптимальной противоударной изоляции получила в [47]. Изоляция объекта, расположенного на подвижном основании, от кратковременных ударных воздействий с помощью активного изолятора с управлением без упреждения изучалась в [8]. Численное решение этой задачи и построение ее аналитических решений для некоторых заданных возмущений с упреждающим управлением рассмотрены в [9]. В [10, 11] аналитически построено гарантирующее упреждающее управление с одним переключением и временем упреждения, позволяющее вычислять минимальную оценку максимальной величины смещения объекта относительно основания для класса возмущений заданной длительности с неограниченной амплитудой. Для этого же класса возмущений в [12, 13] показано, что, используя упреждение и запаздывание для двух простых управлений, предназначенных для конкретных возмущений, и оптимизируя функционал только по моменту начала управления, можно получить значение критерия качества, близкое к оптимальному. Современное описание проблемы противоударной изоляции, включающее в себя публикации по оптимальному управлению противоударными изоляторами, представлено в [14].

В данной работе для возмущений заданной длительности и ограниченной амплитуды доказано, что при заданном управлении наихудшими возмущениями являются непрерывные прямоугольные возмущения наибольшей амплитуды, начинающиеся либо в момент начала возмущения, либо заканчивающиеся в момент окончания допустимого интервала возмущения. Для рассматриваемого класса возмущений используются два управления, полученные в [9, 11] для прямоугольного возмущения и для возмущений заданной длительности с неограниченной амплитудой соответственно. При этом в формулу для управления из [9] подставляется ограничение на амплитуду данной задачи вместо значения амплитуды прямоугольного возмущения. Выбор момента начала действия для этого управления предполагает, что возможно упреждение и запаздывание. Момент начала действия управления выбирается с целью обеспечить минимально возможную величину смещения защищаемого объекта. Проводится сравнение полученных значений критерия качества для управлений из [9, 11] (относительного максимального смещения защищаемого объекта относительно основания). На основе этого сравнения предлагается комбинированное управление, сочетающее достоинства обоих управлений.

1. Механическая система. Пусть механическая система (рис. 1 ) состоит из основания 1 и объекта 2, соединенного с основанием с помощью противоударного изолятора – устройства, генерирующего управляющую силу f между основанием и объектом и созданного для защиты объекта при ударном воздействии на основание. Основание и объект движутся поступательно вдоль одной прямой. Обозначим: z – смещение основания относительно неподвижной (инерциальной) системы отсчета, $x$ – смещение объекта относительно основания, $m$ – масса объекта. Ударное воздействие на основание задается его ускорением $\ddot {z}$, функцией времени, некоторые свойства которой заранее известны.

Рис. 1
Рис. 2
Рис. 3
Рис. 4
Рис. 5

Движение объекта относительно основания описывается уравнением

(1.1)
$\ddot {x} + u = {v}(t),\quad u = - \frac{f}{m},\quad {v} = - \ddot {z}$

Далее полагаем, что сила f удовлетворяет ограничению $\left| f \right| \leqslant {{F}_{0}}$, где ${{F}_{0}}$ – заданная величина, тогда величина $u$ удовлетворяет неравенству

$\begin{gathered} \left| u \right| \leqslant {{u}_{0}},\quad {{u}_{0}} = \frac{{{{F}_{0}}}}{m} \hfill \\ \end{gathered} $

Предполагается, что в начальный момент времени $t = 0$ основание и объект покоятся в положениях, отвечающих нулевым значениям координат $x$ и $z$:

(1.2)
$x(0) = 0,\quad \dot {x}(0) = 0,\quad z(0) = 0,\quad \dot {z}(0) = 0$

В качестве допустимых управлений будем рассматривать кусочно-непрерывные функции u(t), удовлетворяющие ограничению

(1.3)
$\left| {u(t)} \right| \leqslant {{u}_{0}},\quad t \in [0,\infty )$

2. Внешние возмущения. Предполагается, что возмущение ${v}(t)$ имеет вид

(2.1)
${v}(t) = V(t - {{t}_{0}}),\quad {{t}_{0}} \geqslant 0$
где кусочно-непрерывная функция $V(\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$,

(2.2)
$V(\xi ) \leqslant b\quad {\text{для}}\quad 0 \leqslant \xi \leqslant T$
и 4) только на одном интервале, ${{t}_{1}} < \xi < {{t}_{2}}$, величина абсолютного ускорения $V(\xi )$ основания превышает верхнюю границу u0 абсолютного ускорения защищаемого объекта:

$\begin{gathered} V(\xi ) < {{u}_{0}}\quad {\text{для}}\quad 0 \leqslant \xi < {{t}_{1}}\quad {\text{и}}\quad {{t}_{2}} < \xi \leqslant T \\ V(\xi ) > {{u}_{0}}\quad {\text{для}}\quad {{t}_{1}} < \xi < {{t}_{2}} \\ \end{gathered} $

Один или оба из интервалов $0 \leqslant \xi < {{t}_{1}}$ и ${{t}_{2}} < \xi \leqslant T$ могут быть пустыми, если $V(0) > {{u}_{0}}$ или $V(T) > {{u}_{0}}$. В случае, когда $V(\xi )$, $0 \leqslant \xi \leqslant T$, предполагается известной функцией в каждый момент $\xi $ и удовлетворяющей неравенству $V(\xi ) \leqslant {{u}_{0}}$, оптимальное управление определяется тождеством $u(t) \equiv V(t - {{t}_{0}})$ и обеспечивает тождественно равное нулю смещение объекта по отношению к основанию. Далее предполагаем выполнение неравенства

(2.3)
$b > {{u}_{0}}$

Основными характеристиками ударного воздействия являются его длительность T, максимально возможная величина воздействия b и параметр ${{{v}}_{0}}$

(2.4)
${{{v}}_{0}} = \int\limits_0^T V(\xi )d\xi $
характеризующий значение скорости, приобретенной (или потерянной) основанием в результате удара. В [10, 11] рассматривался класс возмущений без ограничений на максимальную величину воздействия $V(\xi )$, а параметры $T$, ${{{v}}_{0}}$ предполагались известными и заданными. Такой класс возмущений обозначим ${{V}_{T}}$. Класс возмущений с заданной величиной ${{{v}}_{0}}$ из (2.4) без ограничений на максимальную величину воздействия $V(t)$ и на максимальную длительность T обозначим ${{V}_{\infty }}$. Этим классам воздействий соответствует значение $b = + \infty $. Введем также классы возмущений ${{V}_{b}}$ и ${{V}_{{bT}}}$ с заданными параметрами $b$, ${{{v}}_{0}}$ и $T$, $b$, ${{{v}}_{0}}$ соответственно. Необходимым условием существования возмущений из класса ${{V}_{{bT}}}$ служит условие $bT \geqslant {{{v}}_{0}},$ поскольку ненулевые значения $V(\xi )$ должны располагаться в прямоугольнике с основанием T и высотой b. Поскольку случай $bT = {{{v}}_{0}}$ соответствует прямоугольным возмущениям и разобран в [9], будем полагать выполненным неравенство

(2.5)
$bT > {{{v}}_{0}}$

Очевидно, что ${{V}_{{bT}}} \subset {{V}_{T}} \subset {{V}_{\infty }}$, ${{V}_{{bT}}} \subset {{V}_{b}} \subset {{V}_{\infty }}$.

3. Критерий качества. Будем считать, что возмущение $V(\xi )$ неизвестно, но известно множество $\Omega \subset {{V}_{\infty }}$, которому могут принадлежать возможные возмущения. Качество изоляции при заданных управлении $u(t)$ и времени упреждения ${{t}_{0}}$ будем оценивать функционалом J, характеризующим максимальную величину смещения объекта относительно основания при наихудшем возмущении:

(3.1)
$J(u,{{t}_{0}}) = \mathop {max}\limits_{V \in \Omega } \mathop {max}\limits_{t \in {{R}_{{0 + }}}} \left| {x(t;u,V,{{t}_{0}})} \right|,\quad {{R}_{{0 + }}} = \left[ {0, + \infty } \right)$
где $x(t;u,V,{{t}_{0}})$ – решение уравнения (1.1) с начальными условиями (1.2) для заданных $u(t)$, $V(\xi )$ и ${{t}_{0}}$. Величину J желательно минимизировать выбором оптимального закона управления и времени упреждения. Для случаев когда $\Omega = {{V}_{b}}$, $\Omega = {{V}_{{bT}}}$, $\Omega = {{V}_{T}}$, $\Omega = {{V}_{\infty }}$ для значения функционала из (3.1) введем соответствующие обозначения ${{J}_{b}}(u$, t0), ${{J}_{{bT}}}(u,{{t}_{0}})$, ${{J}_{T}}(u,{{t}_{0}})$, ${{J}_{\infty }}(u,{{t}_{0}})$, т.е.

(3.2)
${{J}_{b}}(u,{{t}_{0}}) = \mathop {max}\limits_{V \in {{V}_{b}}} \mathop {\;max}\limits_{t \in {{R}_{{0 + }}}} \left| {x(t;u,V,{{t}_{0}})} \right|$
(3.3)
${{J}_{{bT}}}(u,{{t}_{0}}) = \mathop {max}\limits_{V \in {{V}_{{bT}}}} \mathop {\;max}\limits_{t \in {{R}_{{0 + }}}} \left| {x(t;u,V,{{t}_{0}})} \right|$
(3.4)
${{J}_{T}}(u,{{t}_{0}}) = \mathop {max}\limits_{V \in {{V}_{T}}} \mathop {\;max}\limits_{t \in {{R}_{{0 + }}}} \left| {x(t;u,V,{{t}_{0}})} \right|$
(3.5)
${{J}_{\infty }}(u,{{t}_{0}}) = \mathop {max}\limits_{V \in {{V}_{\infty }}} \mathop {\;max}\limits_{t \in {{R}_{{0 + }}}} \left| {x(t;u,V,{{t}_{0}})} \right|$

4. Задачи оптимизации. Поскольку условия информированности управляющей стороны о внешнем возмущении и множества допустимых законов управления $U$ могут быть различными, то и задачи оптимального управления формулируются по-разному.

Задача 1. Для системы (1.1) с начальными условиями (1.2) найти кусочно-непрерывное управление u*, удовлетворяющее ограничению (1.3), и время упреждения $t_{0}^{*}$, которые минимизируют величину (3.1):

$J(u{\text{*}},t_{0}^{*}) = \mathop {min}\limits_{u \in U,{{t}_{0}}} J(u,{{t}_{0}})$
где $U$ – множество законов управления $u(t)$, среди которых ищется оптимум.

Это задача о гарантирующем оптимальном упреждающем управлении противоударным изолятором, защищающим объект от ударных воздействий из множества $\Omega $. Она обобщает задачу, рассмотренную в [13] для заданного возмущения в отсутствие упреждения управления.

Задача 2. Для системы (1.1) при начальных условиях (1.2) и заданном возмущении (2.1) найти кусочно-непрерывное управление u*, удовлетворяющее ограничению (1.3), и время упреждения $t_{0}^{*}$, которые минимизируют максимальную величину смещения объекта относительно основания (функционал $J$):

$J(u*,V,t_{0}^{*}) = \mathop {min}\limits_{u,{{t}_{0}}} \mathop {max}\limits_{t \in [0,\infty )} \left| {x(t;u,V,{{t}_{0}})} \right|$

Задачу 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):

$J(u_{s}^{*},t_{0}^{*}) = \mathop {min}\limits_{{{t}_{0}},s} J({{u}_{s}},{{t}_{0}}),\quad u_{s}^{*} = {{\left. {{{u}_{s}}} \right|}_{{s = s*}}}$

Решение задачи 3 дает возможность улучшать качество противоударной защиты путем изменения времени упреждения и параметра при заданном семействе законов управления ${{u}_{s}}(t)$, которое, например, может быть построено на основе оптимального управления для некоторого возмущения $V \in \Omega $ и представлено аналитическими выражениями. Заметим, что когда множество $S$ состоит из одного элемента задача 3 решается для фиксированного управления, а при заданном управлении $u(t)$ с условием $u(t) \equiv 0$, $t < 0$ в параметрическом семействе управлений ${{u}_{c}}(t) = u(t - c)$, $c \geqslant 0$ параметр $c$ может играть роль запаздывания для управления ${{u}_{c}}(t)$.

В [10, 11] получено решение $\{ u_{\tau }^{*},t_{\tau }^{*}\} $ задачи 3 в классе возмущений ${{V}_{T}}$ для управлений из класса ${{U}_{\Sigma }}\,{\kern 1pt} = \,{\kern 1pt} \{ {{u}_{\tau }}\} $ – параметрического семейства допустимых релейных управлений ${{u}_{\tau }}(t)$ с переключением с $ - {{u}_{0}}$ на ${{u}_{0}}$ в момент времени $\tau $ и с ${{u}_{0}}$ на 0 в момент времени ${{{v}}_{0}}{\text{/}}{{u}_{0}} + 2\tau $, где

(4.1)
$\begin{array}{*{20}{c}} {{{u}_{\tau }}(t) = \left\{ {\begin{array}{*{20}{l}} { - {{u}_{0}},\quad 0 \leqslant t < \tau } \\ { + {{u}_{0}},\quad \tau \leqslant t \leqslant {{T}_{\tau }}} \\ {0,\quad t > {{T}_{\tau }}} \end{array}} \right.} \\ {{{T}_{\tau }} = \frac{{{{{v}}_{0}}}}{{{{u}_{0}}}} + 2\tau ,\quad \tau \in \Sigma ,\quad \Sigma = \left[ {0, + \infty } \right)} \end{array}$

Соотношение временных отрезков со значениями $ \pm {{u}_{0}}$ выбрано таким образом, чтобы обеспечить выполнение соотношения

$\dot {x}(t;{{u}_{\tau }},V,{{t}_{0}}) \equiv 0,\quad t \geqslant max\left( {{{t}_{0}} + T,{{T}_{\tau }}} \right)$
для того чтобы движение объекта относительно основания прекращалось за конечное время и в дальнейшем не возобновлялось при любых возмущениях класса ${{V}_{T}}$ и временах упреждения t0. Решение $\{ u_{\tau }^{*},t_{\tau }^{*}\} $ будем использовать для вычисления величины ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$, используемой для сравнения с величиной ${{J}_{{bT}}}(u_{{bc}}^{*},t_{c}^{*})$, полученной из решения задачи 3 при $u_{{bc}}^{*}(t) = {{u}_{b}}(t - c*)$, $c* \geqslant 0$, где ${{u}_{b}}(t)$ – управление, полученное из решения задачи 2 для прямоугольного возмущения с амплитудой $b$, a c* – оптимальное значение параметра $c$.

5. Определение наихудшего возмущения. Вычисление функционала (3.1) предполагает определение наихудшего возмущения $V \in \Omega \subset {{V}_{T}}$, которое максимизирует максимум модуля отклонения защищаемого объекта относительно основания ($ma{{x}_{t}}{\text{|}}x(t$; u, V, t0)|) при заданных $u(t)$ и t0.

Лемма 1 (см. [9]). Среди возмущений $V \in {{V}_{T}}$ наихудшее возмущение есть либо ${{V}_{\delta }}(\xi ) = {{{v}}_{0}}\delta (\xi )$ либо ${{V}_{{\delta T}}}(\xi ) = {{{v}}_{0}}\delta (\xi - T)$. Иными словами, наихудшее возмущение есть мгновенный удар интенсивности ${{{v}}_{0}}$, подаваемый в начальный или в конечный момент допустимого интервала возмущения.

Согласно лемме 1 для нахождения наихудшего возмущения в классе ${{V}_{T}}$ при заданных $u(t)$ и ${{t}_{0}}$ надо решить дифференциальное уравнение (1.1) с начальными условиями (1.2) при ${v}(t) = {{{v}}_{0}}\delta (t - {{t}_{0}})$ и ${v}(t) = {{{v}}_{0}}\delta (t - {{t}_{0}} - T)$, для каждого из решений вычислить $ma{{x}_{t}}\left| {x(t;u,{{V}_{\delta }},{{t}_{0}})} \right|$, $ma{{x}_{t}}\left| {x(t;u,{{V}_{{\delta T}}},{{t}_{0}})} \right|$ и выбрать возмущение, отвечающее большему значению абсолютного отклонения.

Лемма 2. Среди возмущений $V \in {{V}_{{bT}}}$ наихудшее возмущение есть либо

${{V}_{0}}(\xi ) = \left\{ \begin{gathered} 0,\quad \xi < 0 \hfill \\ b,\quad 0 \leqslant \xi \leqslant {{{v}}_{0}}{\text{/}}b \leqslant T \hfill \\ 0,\quad {{{v}}_{0}}{\text{/}}b < \xi \hfill \\ \end{gathered} \right.$
либо

${{V}_{1}}(\xi ) = \left\{ {\begin{array}{*{20}{l}} {0,\quad \xi < T - {{{v}}_{0}}{\text{/}}b} \\ {b,\quad T - {{{v}}_{0}}{\text{/}}b \leqslant \xi \leqslant T} \\ {0,\quad T < \xi } \end{array}} \right.$

Иными словами, при заданном управлении наихудшими возмущениями являются прямоугольные возмущения наибольшей амплитуды, начинающиеся либо в момент начала возмущения, либо заканчивающиеся в момент окончания допустимого интервала возмущения.

Примечание. Справедливо равенство ${{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}}\left| {x(t;u,{{V}_{0}},{{t}_{0}})} \right|$, $ma{{x}_{t}}\left| {x(t;u,{{V}_{1}},{{t}_{0}})} \right|$, сравнить получившиеся величины и выбрать возмущение, отвечающее большему значению.

Доказательство. Введем вспомогательные функции ${{f}_{ + }}(t)$ и ${{f}_{ - }}(t)$, определенные на полуоси $0 \leqslant t < + \infty $:

(5.1)
$\begin{gathered} \begin{array}{*{20}{l}} {{{f}_{ + }}(t) = \int\limits_0^t (t - \xi ){{V}_{0}}(\xi - {{t}_{0}})d\xi = \left\{ \begin{gathered} 0,\quad 0 \leqslant t < {{t}_{0}} \hfill \\ \frac{{b{{{(t - {{t}_{0}})}}^{2}}}}{2},\quad {{t}_{0}} \leqslant t < {{t}_{0}} + \frac{{{{{v}}_{0}}}}{b} \hfill \\ {{{v}}_{0}}\left( {t - {{t}_{0}} - \frac{{{{{v}}_{0}}}}{{2b}}} \right),\quad {{t}_{0}} + \frac{{{{{v}}_{0}}}}{b} \leqslant t \hfill \\ \end{gathered} \right.} \end{array} \hfill \\ \hfill \\ \end{gathered} $
(5.2)
${{f}_{ - }}(t) = \int\limits_0^t (t - \xi ){{V}_{1}}(\xi - {{t}_{0}})d\xi = \left\{ \begin{gathered} 0,\quad 0 \leqslant t < {{t}_{0}} + T - \frac{{{{{v}}_{0}}}}{b} \hfill \\ \frac{{{{{\left( {t - {{t}_{0}} - T + \frac{{{{{v}}_{0}}}}{b}} \right)}}^{2}}}}{2}b,\quad {{t}_{0}} + T - \frac{{{{{v}}_{0}}}}{b} \leqslant t < {{t}_{0}} + T \hfill \\ {{{v}}_{0}}\left( {t - {{t}_{0}} - T + \frac{{{{{v}}_{0}}}}{{2b}}} \right),\quad t \geqslant {{t}_{0}} + T \hfill \\ \end{gathered} \right.$

Здесь ${{t}_{0}}$ – время начала действия возмущения, а T – его длительность.

Для возмущения ${v}(t) = V(t - {{t}_{0}})$ длительности T, удовлетворяющего соотношению

$\int\limits_{{{t}_{0}}}^{{{t}_{0}} + T} {v}(t)dt = {{{v}}_{0}}$
справедливы неравенства

(5.3)
$\begin{array}{*{20}{l}} {{{f}_{ - }}(t) \leqslant \int\limits_{{{t}_{0}}}^t (t - \xi ){v}(\xi )d\xi \leqslant {{f}_{ + }}(t)} \end{array}$

Докажем правое неравенство в (5.3) для случая ${{t}_{0}} = 0.$ Рассмотрим некоторую функцию $V( \cdot ) \in {{V}_{b}}$, такую что для нее существуют неотрицательные величины ${{t}_{1}}$, ${{t}_{2}}$, h и функция $\Delta {\kern 1pt} V{\kern 1pt} (t)$, удовлетворяющие условиям

(5.4)
$\begin{gathered} 0 \leqslant {{t}_{1}} \leqslant {{t}_{1}} + h \leqslant {{t}_{2}} \leqslant {{t}_{2}} + h \leqslant T \\ 0 \leqslant \Delta V(t)\quad {\text{при}}\quad 0 \leqslant t \leqslant h \\ \Delta V(t) \equiv 0\quad {\text{при}}\quad t \notin [0,h] \\ 0 \leqslant V(t) + \Delta V(t - {{t}_{1}}) \leqslant b\quad {\text{при}}\quad t \in [{{t}_{1}},{{t}_{1}} + h] \\ 0 \leqslant \Delta V(t - {{t}_{2}}) \leqslant V(t)\quad {\text{при}}\quad t \in [{{t}_{2}},{{t}_{2}} + h] \\ \end{gathered} $

Введем функцию ${{V}_{\Delta }}(t)$, $t \in {{R}_{{0 + }}}$ соотношениями

(5.5)
$\begin{gathered} {{V}_{\vartriangle }}(t) \equiv V(t)\quad при\quad t \notin ({{t}_{1}},{{t}_{1}} + h) \cup ({{t}_{2}},{{t}_{2}} + h) \\ {{V}_{\vartriangle }}(t) = V(t) + \Delta V(t - {{t}_{1}})\quad при\quad t \in \left[ {{{t}_{1}},{{t}_{1}} + h} \right] \\ {{V}_{\vartriangle }}(t) = V(t) - \Delta V(t - {{t}_{2}})\quad при\quad t \in \left[ {{{t}_{2}},{{t}_{2}} + h} \right] \\ \end{gathered} $

Значения функции $\Delta {\kern 1pt} V(t)$ по сравнению с $V(t)$ на отрезке $\left[ {{{t}_{1}},{{t}_{1}} + h} \right]$ увеличены на такие же значения, на которые уменьшены ее значения на отрезке $\left[ {{{t}_{2}},{{t}_{2}} + h} \right]$. Согласно соотношениям (5.4), (5.5) имеет место включение ${{V}_{\Delta }}( \cdot ) \in {{V}_{b}}$ и выполнены равенства

(5.6)
$\begin{array}{*{20}{l}} {\int\limits_0^T V(t)dt = \int\limits_0^T {{V}_{\Delta }}(t)dt = {{{v}}_{0}}} \end{array}$

Введем следующие функции

(5.7)
$\begin{array}{*{20}{l}} {I(t) = \int\limits_0^t (t - \xi )V(\xi )d\xi ,\quad {{I}_{\Delta }}(t) = \int\limits_0^t (t - \xi ){{V}_{\Delta }}(\xi )d\xi ,\quad t \in {{R}_{{0 + }}}} \end{array}$

Тогда

(5.8)
$\begin{array}{*{20}{l}} {{{I}_{\Delta }}(t) = \left\{ \begin{gathered} I(t),\quad 0 \leqslant t < {{t}_{1}}, \hfill \\ I(t) + \int\limits_{{{t}_{1}}}^t {(t - \xi )} \Delta V(\xi - {{t}_{1}})d\xi ,\quad {{t}_{1}} \leqslant t < {{t}_{1}} + h \hfill \\ I(t) + \int\limits_{{{t}_{1}}}^{{{t}_{1}} + h} {(t - \xi )} \Delta V(\xi - {{t}_{1}})d\xi ,\quad {{t}_{1}} + h \leqslant t < {{t}_{2}} \hfill \\ I(t) + \int\limits_0^h {(t - {{t}_{1}} - y)} \Delta V(y)dy - \int\limits_{{{t}_{2}}}^t {(t - \xi )} \Delta V(\xi - {{t}_{2}})d\xi ,\quad {{t}_{2}} \leqslant t < {{t}_{2}} + h \hfill \\ I(t) + ({{t}_{2}} - {{t}_{1}})\int\limits_0^h {\Delta V} (y)dy,\quad {{t}_{2}} + h \leqslant t \hfill \\ \end{gathered} \right.} \end{array}$

Покажем, что

(5.9)
$\begin{array}{*{20}{l}} {{{I}_{\Delta }}(t) \geqslant I(t)} \end{array}$

Согласно (5.8) неравенство (5.9) выполняется для всех $t \in [0, + \infty ){{\backslash }}[{{t}_{2}},{{t}_{2}} + h)$, так как к $I(t)$ добавляется неотрицательная величина.

Поскольку при ${{t}_{1}} \leqslant {{t}_{2}} \leqslant t < {{t}_{2}} + h$

$\begin{array}{*{20}{l}} {0 \leqslant \int\limits_0^{t - {{t}_{2}}} {(t - {{t}_{2}} - y)} \Delta V(y)dy \leqslant \int\limits_0^h {(t - {{t}_{2}} - y)} \Delta V(y)dy \leqslant \int\limits_0^h {(t - {{t}_{1}} - y)} \Delta V(y)dy} \end{array}$
то добавляемая величина к $I(t)$ в (5.8) также неотрицательна

$\begin{gathered} \int\limits_0^h {(t - {{t}_{1}} - y)} \Delta V(y)dy - \int\limits_{{{t}_{2}}}^t {(t - \xi )} \Delta V(\xi - {{t}_{2}})d\xi = \\ = \;\int\limits_0^h {(t - {{t}_{1}} - y)} \Delta V(y)dy - \int\limits_0^{t - {{t}_{2}}} {(t - {{t}_{2}} - y)} \Delta V(y)dy \geqslant 0 \\ \end{gathered} $

Неравенство (5.9) доказано. Нетрудно убедиться, что при кусочно-непрерывной $V(t)$, удовлетворяющей условиям (5.4) и такой что , , $t \in {{R}_{{0 + }}}$ выполнено строгое неравенство ${{I}_{\Delta }}(t) > I(t)$ при $t \geqslant {{t}_{1}} + h$. Таким образом, возможность построения функции приводит к увеличению функции $I(t)$, $t \in {{R}_{{0 + }}}$. Невозможность построения ${{V}_{\Delta }}(t)$ для функции ${{V}_{0}}(t)$, $t \in {{R}_{{0 + }}}$ доказывает невозможность увеличения функции $I(t)$, $t \in {{R}_{{0 + }}}$ и справедливость правого неравенства (5.3) для ${{t}_{0}} = 0$. Аналогично доказывается справедливость левого неравенства (5.3) для ${{t}_{0}} = 0$, а также справедливость неравенств (5.3) для ${{t}_{0}} > 0$.

Справедливость правого неравенства (5.3) для ${{t}_{0}} = 0$ также легко доказывается следующим способом. Из соотношений

$\begin{array}{*{20}{l}} {{{I}_{0}}(t) = \int\limits_0^t {(t - \xi )} {{V}_{0}}(\xi )d\xi ,\quad 0 \leqslant {{V}_{0}}(\xi ) \leqslant b,\quad \int\limits_0^T {{{V}_{0}}(\xi )} d\xi = {{{v}}_{0}},\quad {{{\dot {I}}}_{0}} = \int\limits_0^t {{{V}_{0}}} (\xi )d\xi } \end{array}$
следует что

$\begin{array}{*{20}{l}} {{{{\dot {I}}}_{0}}(t) \leqslant \psi (t),\quad {{I}_{0}}(t) \leqslant \int\limits_0^t {\psi (\xi )} d\xi ,\psi (t) = min\{ bt,{{{v}}_{0}}\} ,\quad \int\limits_0^t {\psi (\xi )} d\xi = {{{\left. {{{f}_{ + }}(t)} \right|}}_{{{{t}_{0}} = 0}}}} \end{array}$

Запишем решение уравнения (1.1) с нулевыми начальными условиями в следующем виде:

${{x}_{u}}(t) = \int\limits_0^t {(t - \xi )} {v}(\xi ){\kern 1pt} d\xi - \int\limits_0^t {(t - \xi )} u(\xi )d\xi $

Применив к этому выражению неравенство (5.3), получим

(5.10)
$\begin{gathered} x_{u}^{ - }(t) \leqslant {{x}_{u}}(t) \leqslant x_{u}^{ + }(t) \\ x_{u}^{ - }(t) = {{f}_{ - }}(t) - \int\limits_0^t {(t - \xi )} u(\xi )d\xi \\ x_{u}^{ + }(t) = {{f}_{ + }}(t) - \int\limits_0^t {(t - \xi )} u(\xi )d\xi \\ \end{gathered} $

Из этих неравенств и выражений (5.1), (5.2) для функций ${{f}_{ + }}(t)$ и ${{f}_{ - }}(t)$ вытекает, что при фиксированном $u = u(t)$ в любой момент времени t величина ${{x}_{u}}(t)$ принимает минимальное и максимальное значения при возмущениях ${{{v}}_{ - }}(t) = {{V}_{1}}(t - {{t}_{0}})$ и ${{{v}}_{ + }}(t) = {{V}_{0}}(t$t0) соответственно. Лемма 2 доказана.

Из (5.10) следует, что максимальное значение, которое может принимать критерий качества J при заданном управлении $u$ и различных возмущениях, определяется выражением

(5.11)
$J_{u}^{ + } = max\left\{ {\mathop {max}\limits_{t \in {{R}_{{0 + }}}} \left| {x_{u}^{ - }(t)} \right|,\mathop {max}\limits_{t \in {{R}_{{0 + }}}} \left| {x_{u}^{ + }(t)} \right|} \right\}$

Безразмерные переменные. Введем безразмерные переменные

$\begin{gathered} x{\text{'}} = \frac{{{{u}_{0}}}}{{{v}_{0}^{2}}}x,\quad t{\text{'}} = \frac{{{{u}_{0}}}}{{{{{v}}_{0}}}}t,\quad t_{0}^{'} = \frac{{{{u}_{0}}}}{{{{{v}}_{0}}}}{{t}_{0}},\quad T{\kern 1pt} {\text{'}} = \frac{{{{u}_{0}}}}{{{{{v}}_{0}}}}T,\quad \xi {\text{'}} = \frac{{{{u}_{0}}}}{{{{{v}}_{0}}}}\xi \\ {v}{\text{'}}(t{\text{'}}) = \frac{1}{{{{{v}}_{0}}}}{v}\left( {\frac{{{{{v}}_{0}}}}{{{{u}_{0}}}}t{\text{'}}} \right),\quad u{\text{'}} = \frac{u}{{{{u}_{0}}}},\quad b{\text{'}} = \frac{b}{{{{u}_{0}}}},\quad J{\text{'}} = \frac{{{{u}_{0}}}}{{{v}_{0}^{2}}}J \\ \end{gathered} $

Далее будем использовать безразмерные переменные, опуская штрихи. В безразмерных единицах имеем ${{u}_{0}} = 1$ и ${{{v}}_{0}} = 1$, а неравенства (2.3), (2.5) приобретают соответственно вид

(5.12)
$b > 1$
(5.13)
$bT > 1$

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.$
в безразмерных переменных. Для этого возмущения задача 2 решена в [9] с помощью графо-аналитического метода [1, 2]. Оптимальное управление $u(t)$, оптимальное время упреждения ${{t}_{0}}$ и оптимальное значение критерия качества $J$ определяются соотношениями

(6.3)
${{u}_{b}}(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}}}} + 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.$

Индекс b указывает, что данные соотношения выражены через значение амплитуды возмущения b. Управление (6.3) будем называть b-управлением. Поскольку формулы (6.3)(6.5) были получены в предположении выполнения равенства

(6.6)
$bT = 1$
их можно переписать в следующем виде

(6.7)
${{u}_{T}}(t) = \left\{ \begin{gathered} - 1,\quad 0 \leqslant t \leqslant t_{1}^{r} = \frac{1}{4}\sqrt {1 - T} ,\quad T < 1 \hfill \\ 1,\quad t_{1}^{r} < t \leqslant t_{2}^{r} = \frac{1}{2}\sqrt {1 - T} + 1,\quad T < 1 \hfill \\ 0,\quad t > t_{2}^{r},\quad T < 1 \hfill \\ V(t),\quad t \in [0, + \infty ),\quad T \geqslant 1 \hfill \\ \end{gathered} \right.$
(6.8)
$t_{0}^{T} = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{2}[\left( {1 - T} \right) + \sqrt {1 - T} ],\quad T < 1} \\ {0,\quad T \geqslant 1} \end{array}} \right.$
(6.9)
${{J}_{r}} = \left\{ {\begin{array}{*{20}{l}} {\frac{1}{{16}}\left( {1 - T} \right),\quad b > 1} \\ {0,\quad b \leqslant 1} \end{array}} \right.$

Индекс T указывает, что данные соотношения выражены через значение длительности T. Управление (6.7) будем называть T-управлением. Такое введение различий для управлений ${{u}_{b}}$, ${{u}_{T}}$ и их упреждающих моментов обусловлено тем, что дальнейшее использование этих управлений для задач гарантированного оценивания (например, в задаче 3), в которых выполнены неравенства $bT > 1$, $b > 1$, приводит к различным результатам. В данной работе изучается эффективность управления ${{u}_{b}}$ с запаздыванием для задачи гарантированного оценивания с построением оптимального момента начала управляющего воздействия.

Введем класс управлений ${{U}_{{bC}}} = \{ {{u}_{{bc}}}\} $, описываемый параметрическим семейством допустимых управлений ${{u}_{{bc}}}(t)$, где

(6.10)
$\begin{gathered} {{u}_{{bc}}}(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 = [0,{\mkern 1mu} + \infty ) \\ \end{gathered} $

Равенство ${{u}_{{bc}}}(t) = {{u}_{b}}(t - c)$ означает, что управление ${{u}_{{bc}}}(t)$ представляет собой b‑управление с запаздыванием $c$.

Введем функцию ${{J}_{{bc}}}(t*)$, определяемую как значение функционала J, отвечающее управлению ${{u}_{{bc}}}(t)$ и возмущению ${v}(t) = {{V}_{0}}(t - t*)$, т.е. прямоугольному возмущению с амплитудой b с началом действия в момент времени t*. Эта функция задается выражением

(6.11)
${{J}_{{bc}}}(t{\text{*}})\, = \,\left\{ {\begin{array}{*{20}{l}} {9t_{1}^{2} + 2{{t}_{1}} - t{\text{*}} + c\quad {\text{при}}\quad 0 \leqslant t* < 8t_{1}^{2} + 2{{t}_{1}} + c} \\ {\frac{{2t_{1}^{2}\, + \,2t_{1}^{2}b\, + \,b{{{(t{\text{*}}\, - \,c)}}^{2}} - 4{{t}_{1}}(t{\text{*}}\, - \,c)b}}{{2(b - 1)}}\quad {\text{при}}\quad {\text{ }}8t_{1}^{2} + 2{{t}_{1}} + c \leqslant t* < {{t}_{2}} - {{b}^{{ - 1}}} + c} \\ {t{\text{*}} - c - 9t_{1}^{2} - 2{{t}_{1}}\quad {\text{при}}\quad t* \geqslant {{t}_{2}} - {{b}^{{ - 1}}} + c} \end{array}} \right.$

Функция ${{J}_{{bc}}}(t*)$ линейно убывает при $t* \in [0,\;8t_{1}^{2} + 2{{t}_{1}} + c]$ от значения $9{\kern 1pt} t_{1}^{2} + 2{\kern 1pt} {{t}_{1}} + c$ до значения

(6.12)
$\begin{gathered} {{J}_{{bc}}}(8t_{1}^{2} + 2{{t}_{1}} + c) = t_{1}^{2} = \frac{{b - 1}}{{16b}} \hfill \\ \end{gathered} $
затем при $t* \in [8t_{1}^{2} + 2{{t}_{1}} + c,\; + \infty )$ монотонно возрастает от значения (6.12) до бесконечности. При этом функция ${{J}_{{bc}}}(t*)$ на отрезке $[8t_{1}^{2} + 2{{t}_{1}} + c,\;{{t}_{2}} - {{b}^{{ - 1}}} + c]$ растет квадратично до значения
(6.13)
$\begin{gathered} {{J}_{{bc}}}({{t}_{2}} - {{b}^{{ - 1}}} + c) = 7t_{1}^{2} = \frac{{7(b - 1)}}{{16b}} \hfill \\ \end{gathered} $
а затем линейно до бесконечности. Найдем значения $t* = {{t}_{p}}$, $T = {{T}_{p}}$, соответствующие значению ${{J}_{{bc}}}$ из (6.13):

${{t}_{p}} = 2t_{1}^{2} + 2{{t}_{1}} + c,\quad {{T}_{p}} = 1 - 2t_{1}^{2} = \frac{1}{8}(7 + {{b}^{{ - 1}}}),\quad c < {{t}_{p}} < c + \frac{5}{8}$

Функция ${{J}_{{bc}}}(t*)$ из (6.11) позволяет оценить влияние ошибки в определении момента удара при расчете упреждающего оптимального управления на величину критерия качества противоударной изоляции. Формула (6.11) является обобщением формулы для величины ${{J}_{{b0}}}(t*)$, определяемой как максимум модуля отклонения объекта относительно основания для возмущения ${v}(t) = {{V}_{0}}(t - t*)$ при управлении ${{u}_{b}}(t)$

(6.14)
${{J}_{{b0}}}(t*) = \left\{ {\begin{array}{*{20}{l}} {9t_{1}^{2} + 2{{t}_{1}} - t{\text{*,}}\quad 0 \leqslant t* < 8t_{1}^{2} + 2{{t}_{1}}} \\ {\frac{{2t_{1}^{2} + 2t_{1}^{2}b + b{{t}^{*}}^{2} - 4{{t}_{1}}t{\text{*}}b}}{{2(b - 1)}}{\text{,}}\quad 8t_{1}^{2} + 2{{t}_{1}} \leqslant t* < {{t}_{2}} - {{b}^{{ - 1}}}} \\ {t{\text{*}} - 9t_{1}^{2} - 2{{t}_{1}}{\text{,}}\quad t* \geqslant {{t}_{2}} - {{b}^{{ - 1}}}} \end{array}} \right.$

Равенство ${{J}_{{bc}}}(t*) = {{J}_{{b0}}}(t{\text{*}} - c)$ при $t* \geqslant c$ означает, что график функции ${{J}_{{bc}}}(t*)$ получается смещением графика на величину $c$ вдоль оси абсцисс с доопределением на интервале [0, c) согласно формуле (6.11). При $b \to + \infty $ функция ${{J}_{{b0}}}(t*) \to {{J}_{{d0}}}(t*)$, где функция ${{J}_{{d0}}}(t*)$ предназначена для дельта возмущения и имеет тот же смысл, что и функция ${{J}_{{b0}}}(t*)$ для ограниченного возмущения. Функция ${{J}_{{d0}}}(t*)$ задается формулой

(6.15)
$\begin{gathered} {{J}_{{d0}}}(t*) = \left\{ \begin{gathered} 17{\text{/}}16 - t*,\quad 0 \leqslant t* < 1 \hfill \\ ({{t}^{*}}^{2} - t*){\text{/}}2 + 1{\text{/}}16,\quad 1 \leqslant t* < 3{\text{/}}2 \hfill \\ t{\text{*}} - 17{\text{/}}16,\quad t* \geqslant 3{\text{/}}2 \hfill \\ \end{gathered} \right. \hfill \\ \hfill \\ \end{gathered} $
и приведена на рис. 2 . Функция ${{J}_{{b0}}}(t*)$ имеет аналогичный вид. Для решения задачи 3 надо, с учетом леммы 2 о наихудшем возмущении, для заданного T найти минимум по переменным ${{t}_{0}}$, c величины

(6.16)
${{f}_{c}}({{t}_{0}},T) = max[{{J}_{{bc}}}({{t}_{0}}),\;{{J}_{{bc}}}({{t}_{0}} + T - {{b}^{{ - 1}}})],\quad {{t}_{0}} \geqslant 0$

Минимум величины из (6.16) достигается при выполнении условий

(6.17)
$c \geqslant max[0,T{\text{/}}2 - 9t_{1}^{2} - 2{{t}_{1}}],\quad {{J}_{{bc}}}({{t}_{0}}) = {{J}_{{bc}}}({{t}_{0}} + T - {{b}^{{ - 1}}}),\quad {{t}_{0}} \geqslant 0$
которые приводят к решению задачи 3 для классов возмущений ${{V}_{{bT}}}$ и управлений ${{U}_{{bC}}}$:

(6.18)
$\begin{gathered} \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}} + 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.} \\ {} \end{array} \\ c \geqslant max\left[ {0,\;\frac{T}{2} - \frac{9}{{16}} + \frac{1}{{16b}} - \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2}} \right] \\ \end{gathered} $
(6.19)
$\begin{array}{*{20}{c}} {{{J}_{{bT}}}({{u}_{c}},{{t}_{c}}) = \left\{ \begin{gathered} \frac{{25}}{{16}} - \frac{{41}}{{16b}} + T - \frac{{\sqrt {9{{b}^{2}} + 8{{b}^{2}}T - 26b - 8Tb + 17} }}{{2b}},\quad \frac{1}{b} < T \leqslant \frac{7}{8} + \frac{1}{{8b}} \hfill \\ \frac{T}{2} - \frac{1}{{2b}},\quad \frac{7}{8} + \frac{1}{{8b}} < T \hfill \\ \end{gathered} \right.} \end{array}$

Оптимальное управление в классах управлений ${{U}_{{bC}}} = \{ {{u}_{{bc}}}\} $ и возмущений ${{V}_{{bT}}}$ для задачи 3 определяется формулой (6.10) со значением параметра c, удовлетворяющим неравенству из (6.18). Заметим, что значение функционала ${{J}_{{bT}}}({{u}_{c}},t_{c}^{{}})$ из (6.19) обеспечивается не единственным образом. Оптимальные минимально возможные момент упреждения $t_{c}^{*}$ и величина запаздывания c* определяются соотношениями

$\begin{array}{*{20}{c}} {t_{c}^{*} = \left\{ \begin{gathered} \frac{{\sqrt {9{{b}^{2}} + 8{{b}^{2}}T - 26b - 8Tb + 17} }}{{2b}} + \frac{2}{b} + \frac{{\sqrt {1 - {{b}^{{ - 1}}}} }}{2} - 1 - 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}{{8b}} \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.} \\ {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.} \end{array}$

При минимально возможных моменте упреждения $t_{c}^{*}$ и величине запаздывания $c{\text{*}}$ оптимальное управление определяется формулой $u_{{bc}}^{*}(t) = {{u}_{b}}(t - c*)$ и справедливо равенство $\min \{ t_{c}^{*},c*\} = 0$, т.е. имеет место либо упреждение, либо запаздывание. При не минимально выбранных моменте упреждения $t_{c}^{{}}$ и величине запаздывания $c$ имеем неравенство $\min \{ {{t}_{c}},c\} > 0$. В этом случае при $0 \leqslant t < \min \{ {{t}_{c}},c\} $ рассматриваемая система бездействует, поскольку значения управления и возмущения имеют нулевые значения.

7. Решение задачи 3 для класса ограниченных возмущений при использовании управления, предназначенного для возмущений ограниченной длительности. В безразмерных переменных управление ${{u}_{\tau }}(t)$ из (4.1) приобретает вид

(7.1)
$\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 }}} \\ {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.1) со значением параметра τ [10, 11]:

(7.2)
$\tau * = \left\{ \begin{gathered} T{\text{/}}2 + 1{\text{/}}4\quad {\text{при}}\quad T \leqslant 1{\text{/}}2 \hfill \\ 1{\text{/}}2\quad {\text{при}}\quad 1{\text{/}}2 < T \leqslant 7{\text{/}}2 \hfill \\ \sqrt {(1 + T){\text{/}}2} - 1\quad {\text{при}}\quad T > 7{\text{/}}2 \hfill \\ \end{gathered} \right.$
а значение функционала и момент упреждения находят как

(7.3)
$\begin{gathered} {{J}_{\tau }} = \left\{ \begin{gathered} {{\left( {T{\text{/}}2 + 1{\text{/}}4} \right)}^{2}}\quad {\text{при}}\quad T \leqslant 1{\text{/}}2 \hfill \\ T{\text{/}}2\quad {\text{при}}\quad 1{\text{/}}2 < T \hfill \\ \end{gathered} \right. \\ t_{\tau }^{*} = \left\{ \begin{gathered} T + 1\quad {\text{при}}\quad T \leqslant 1{\text{/}}2 \hfill \\ 7{\text{/}}4 - T{\text{/}}2\quad {\text{при}}\quad 1{\text{/}}2 < T \leqslant 7{\text{/}}2 \hfill \\ 0\quad {\text{при}}\quad T > 7{\text{/}}2 \hfill \\ \end{gathered} \right. \\ \end{gathered} $

Заметим, что приведенное решение (7.1)–(7.3) задачи 3 при $T > 1{\kern 1pt} /{\kern 1pt} 2$ не единственно. Значение функционала $J = T{\text{/}}2$ для этого случая обеспечивается моментом переключения $\tau $ и моментом упреждения ${{t}_{\tau }}$, удовлетворяющим соотношениям $t_{\tau }^{*}$

$\begin{gathered} max\{ \sqrt {(1 + T){\text{/}}2} - 1,0\} \leqslant \tau \leqslant \sqrt {T{\text{/}}2} \\ {{t}_{\tau }} = 1{\text{/}}2 + {{\tau }^{2}} + 2\tau - T{\text{/}}2 \\ \end{gathered} $

Обозначим через $u_{\tau }^{*}(t)$ функцию из (7.1) со значением $\tau = \tau {\text{*}}$ из (7.2). Введем функцию ${{J}_{{b\tau }}}(t*)$, определяемую как значение функционала J, отвечающее управлению $u_{\tau }^{*}(t)$ и возмущению ${v}(t) = {{V}_{0}}(t - t*)$, т.е. прямоугольному возмущению с амплитудой b с началом действия в момент времени t*. Для вычисления величины ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$ необходимо посчитать значение выражения

${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*}) = max[{{J}_{{b\tau }}}(t_{\tau }^{*}),{{J}_{{b\tau }}}(t_{\tau }^{*} + T - {{b}^{{ - 1}}})]$

Значение J для класса возмущений ${{V}_{{bT}}}$ при управлении $u_{\tau }^{*}$ с моментом упреждения $t_{\tau }^{*}$ для всех $1 < b \leqslant 2$ имеет вид

(7.4)
$\begin{array}{*{20}{c}} {{{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*}) = \left\{ \begin{gathered} \frac{1}{4}\quad {\text{при}}\quad \frac{1}{2} \leqslant \frac{1}{b} < T \leqslant \frac{1}{b} + \frac{1}{2} \hfill \\ \frac{T}{2} - \frac{1}{{2b}}\quad {\text{при}}\quad \frac{1}{b} + \frac{1}{2} < T \hfill \\ \end{gathered} \right.} \end{array}$

Значение J для класса возмущений ${{V}_{{bT}}}$ при управлении $u_{\tau }^{*}$ с моментом упреждения $t_{\tau }^{*}$ для всех $b > 2$ имеет более сложный вид

(7.5)
$\begin{array}{*{20}{c}} {{{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*}) = \left\{ \begin{gathered} {{\left( {\frac{T}{2} + \frac{1}{4}} \right)}^{2}}\quad {\text{при}}\quad \frac{1}{b} < T \leqslant \frac{1}{2} \hfill \\ \frac{1}{4}\quad {\text{при}}\quad \frac{1}{2} < T \leqslant \frac{1}{b} + \frac{1}{2} \hfill \\ \frac{T}{2} - \frac{1}{{2b}}\quad {\text{при}}\quad \frac{1}{b} + \frac{1}{2} < T \hfill \\ \end{gathered} \right.} \end{array}$

8. Качество изоляции при комбинированном способе управления. Значения функционалов ${{J}_{{bT}}}(u_{{bc}}^{*},t_{c}^{*})$, ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$ из (6.19), (7.4), (7.5) представляют собой функции от параметров b, T, заданные в области

${{G}_{{bT}}} = \left\{ {b,T:b > 1,\;Tb > 1} \right\}$

Введем для значений ${{J}_{{bT}}}(u_{{bc}}^{*},t_{c}^{*})$, ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$ другие обозначения, отражающие эту зависимость

${{J}_{ - }}(b,T) = {{J}_{{bT}}}(u_{{bc}}^{*},t_{c}^{*}),{{J}_{ + }}(b,T) = {{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$

Вычисления показывают, что равенство ${{J}_{ - }}(b,T) = {{J}_{ + }}(b,T)$ достигается при $1{\text{/}}b$ < T < 1 на кривой l, описываемой функцией $b = f(T)$, $0 < T \leqslant 1$, где

$f(T) = \left\{ \begin{gathered} \frac{{38 + 59T - 41{{T}^{2}} + 2\sqrt {361 + 528T + 536{{T}^{2}} - 320{{T}^{3}} + 272{{T}^{4}}} }}{{4T(T + 1)({{T}^{2}} - 7T + 4)}} \hfill \\ {\text{при}}\quad 0 < T \leqslant \frac{1}{2} \hfill \\ \frac{{29 + 400T + 32\sqrt {79 - 70T + 8{{T}^{2}}} }}{{160T - 135 + 256{{T}^{2}}}}\quad {\text{при}}\quad \frac{1}{2} < T \leqslant \frac{{13}}{{14}} \hfill \\ \frac{1}{{8T - 7}}\quad {\text{при}}\quad \frac{{13}}{{14}} \leqslant T \leqslant 1 \hfill \\ \end{gathered} \right.$

При этом кривая l стремится к плюс бесконечности при $T \to + 0$, лежит выше кривой $b = 1{\text{/}}T$ при $0 < T < 1$ и пересекается с ней в точке T = 1, b = 1. Кривые l и $b = 1{\text{/}}T$ изображены сплошными линиями на рис. 3 . Кривая l делит область ${{G}_{{bT}}}$ на две подобласти ${{G}_{ - }}$, ${{G}_{ + }}$, располагающиеся ниже и выше кривой l соответственно. При этом в нижней области ${{G}_{ - }}$ выполнено неравенство ${{J}_{ - }} < {{J}_{ + }}$, а в верхней области ${{G}_{ + }}$ выполнено неравенство ${{J}_{ - }} > {{J}_{ + }}$. Зададим функцию ${{J}_{{comb}}}(b,T)$ соотношением

$\begin{gathered} {{J}_{{comb}}}(b,T) = min\left\{ {{{J}_{ - }}(b,T),\;{{J}_{ + }}(b,T)} \right\},\quad (b,T) \in {{G}_{{bT}}} \hfill \\ \end{gathered} $
и поясним как эта функция вычисляется. Дополнительно область ${{G}_{ + }}$ на рис. 3 разделена двумя пунктирными линиями ${{l}_{1}}$, ${{l}_{2}}$, задаваемыми формулами

${{l}_{1}}:T = \frac{1}{2},\quad b \geqslant \frac{{229 + 32\sqrt {46} }}{9} \approx 49.55939549$
${{l}_{2}}:b = \frac{2}{{2T - 1}},\quad \frac{1}{2} < T \leqslant \frac{{13}}{{14}}$

Линия ${{l}_{1}}$ имеет общую точку с кривой $l$:

$T = \frac{1}{2},\quad b = \frac{{229 + 32\sqrt {46} }}{9}$
а кривая ${{l}_{2}}$ имеет общую точку с кривой $l$:

$T = \frac{{13}}{{14}},\quad b = \frac{7}{3}$

В каждой из четырех областей, помеченных на рис. 3 , значение ${{J}_{{comb}}}(b,T)$ вычисляется по своей формуле с соответствующим индексом. То есть области со значением индекса i (i = 1, 2, 3, 4) соответствует своя функция ${{J}_{i}}(b,T)$. Функции ${{J}_{i}}(b,T)$ имеют следующий вид

(8.1)
$\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} $

Функция Jcomb(b, T) представляет собой значение функционала (3.3) при использовании в качестве управления и момента упреждения величин $u_{{bc}}^{*}$, $t_{c}^{*}$. в области ${{G}_{ - }}$ и при использовании в качестве управления и момента упреждения величин $u_{\tau }^{*},t_{\tau }^{*}$ в области ${{G}_{ + }}$. Такой способ управления с моментом упреждения назовем комбинированным управлением.

На рис. 4 изображено разбиение области ${{G}_{{bT}}}$ на две подобласти, которые фигурируют в формуле (6.19). Цифры 1 и 4, которыми помечены эти области, означают, что функционал ${{J}_{{bT}}}(u_{{bc}}^{*},t_{c}^{*})$ имеет значения J1 и J4 из (8.1) в соответствующих подобластях.

На рис. 5 изображено разбиение области ${{G}_{{bT}}}$ на три подобласти, которые фигурируют в формулах (7.4), (7.5). Цифры 2–4, которыми помечены эти области, означают что функционал ${{J}_{{bT}}}(u_{\tau }^{*},t_{\tau }^{*})$ имеет значения J2, J3, J4 из (8.1) в соответствующих подобластях.

Вычисления показывают, что применение комбинированного управления приводит к улучшению значения функционала, доходящему до 10% по сравнению со значением функционала ${{J}_{ - }}$, полученным для управления $u_{{bc}}^{*}$ и момента упреждения $t_{c}^{*}$. Применение комбинированного управления приводит к улучшению значения функционала, доходящему до 700% по сравнению со значением функционала J+, полученным для управления $u_{\tau }^{*}$ и момента упреждения $t_{\tau }^{*}$. Действительно, наихудшие значения относительной разности решений демонстрируют следующие пределы

$\begin{gathered} {{\left. {\mathop {lim}\limits_{b \to + \infty } \frac{{{{J}_{ - }} - {{J}_{{comb}}}}}{{{{J}_{{comb}}}}}} \right|}_{{T = 1/5}}} = {{\left. {\mathop {lim}\limits_{b \to + \infty } \frac{{{{J}_{1}} - {{J}_{2}}}}{{{{J}_{2}}}}} \right|}_{{T = 1/5}}} = 0.1 \\ {{\left. {\mathop {lim}\limits_{b \to 2 + 0} \frac{{{{J}_{ + }} - {{J}_{{comb}}}}}{{{{J}_{{comb}}}}}} \right|}_{{T = 1/2}}} = {{\left. {\mathop {lim}\limits_{b \to 2 + 0} \frac{{{{J}_{3}} - {{J}_{1}}}}{{{{J}_{1}}}}} \right|}_{{T = 1/2}}} = 7.0 \\ \end{gathered} $

Заключение. Рассмотрена задача изоляции объекта на подвижном основании, которое подвержено однонаправленным возмущениям ограниченной длительности и амплитуды с заданным интегралом от возмущения. Для этой модели изучалась возможность применения управления, предназначенного для прямоугольного возмущения с оптимизацией начала его действия, а также управления, предназначенного для ударов с неограниченной амплитудой, но с ограниченной длительностью. Было показано, что при больших длительностях возмущения предложенные управления и значения функционалов совпадают, а при относительно небольших длительностях возмущений первое управление предпочтительнее использовать для небольших амплитуд возмущения, а при относительно больших амплитудах возмущения предпочтительнее использовать второе управление. На основе этих двух управлений было построено комбинированное управление, сочетающее преимущества обоих управлений.

Работа выполнена по теме государственного задания АААА-А20-120011690138-6.

Список литературы

  1. Гурецкий В.В. Об одной задаче оптимального управления // Изв. АН СССР. Механика. 1965. № 1. С. 159–162.

  2. Гурецкий В.В. О задаче минимизации максимального смещения // Труды ЛПИ. Механика и процессы управления. 1969. № 307. С. 11–21.

  3. Sevin E. and Pilkey W. Optimum Shock and Vibration Isolation. Washington DC: Shock and Vibration Information Analysis Center, 1971, 162 c.

  4. Коловский М.З. Автоматическое управление виброзащитными системами. М.: Наука, 1976, 320 с.

  5. Болотник Н.Н. Оптимизация амортизационных систем. М.: Наука, 1983, 256 с.

  6. Balandin D.V., Bolotnik N.N., and Pilkey W.D. Optimal Protection from Impact, Shock, and Vibration. Amsterdam: Gordon and Breach Science, 2001, 440 p.

  7. Pilkey W.D., Balandin D.V., Bolotnik N.N., Crandal J.R., Purtsezov S.V. Injury Biomechanics and Control: Optimal Protection from Impact. Wiley and Sons, Inc., Hoboken, NJ, 2010. 286 p.

  8. Болотник Н.Н., Корнеев В.А. Анализ предельных возможностей противоударной изоляции при кратковременных внешних воздействиях // Нелинейная динамика. 2015. Т. 11. № 1. С. 147–168.

  9. Болотник Н.Н., Корнеев В.А. Противоударная изоляция с упреждающим управлением для внешних возмущений различной формы // Известия РАН. Теория и системы управления. 2018. № 3. С. 48–63.

  10. Болотник Н.Н., Корнеев В.А. Гарантирующее упреждающее управление в задаче противоударной изоляции // Доклады РАН. 2018. Т. 445. № 4. С. 381–385.

  11. Корнеев В.А. Защита объекта на подвижном основании с помощью упреждающего управления при наихудших возмущениях // Изв. РАН. ТиСУ. 2019. № 1. С. 89–97.

  12. Корнеев В.А. Оптимизация управления с упреждением и запаздыванием в задаче противоударной защиты объекта на подвижном основании // Изв. РАН. ТиСУ. 2020. № 3. С. 30–38.

  13. Корнеев В.А. Использование постоянного управления с упреждением и запаздыванием в задаче защиты объекта от удара на подвижном основании // Изв. РАН. Механика твердого тела. 2020. № 2. С. 147–158.

  14. 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. https://doi.org/10.1115/1.4044190

Дополнительные материалы отсутствуют.