Известия РАН. Механика твердого тела, 2023, № 2, стр. 125-135

ЗАДАЧА ГРАНИЧНОГО УПРАВЛЕНИЯ СМЕЩЕНИЕМ НА ДВУХ КОНЦАХ ПРОЦЕССОМ КОЛЕБАНИЯ СТЕРЖНЯ, СОСТОЯЩЕГО ИЗ ДВУХ УЧАСТКОВ РАЗНОЙ ПЛОТНОСТИ И УПРУГОСТИ

В. Р. Барсегян ab*

a Институт механики НАН Армении
Ереван, Армения

b Ереванский государственный университет
Ереван, Армения

* E-mail: barseghyan@sci.am

Поступила в редакцию 29.11.2021
После доработки 14.03.2022
Принята к публикации 20.04.2022

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

Аннотация

Рассмотрена задача граничного управления для одномерного волнового уравнения с кусочно-постоянными характеристиками. При этом полагается, что время прохождения волны через каждый однородный участок одинаково. Управление осуществляется смещением на двух концах. Предложен конструктивный подход построения управляющего воздействия с заданными начальным и конечным условиями. Схема построения заключается в следующем: исходная задача сводится к задаче управления распределенными воздействиями с нулевыми граничными условиями. Далее используется метод разделения переменных и методы теории управления конечномерными системами. Полученные результаты иллюстрируются на конкретном примере.

Ключевые слова: управление колебаниями, продольные колебания кусочно-однородного стержня, поперечные колебания кусочно-однородной струны, граничное управление, разделение переменных

Введение. Задачам граничного управления и оптималного управления колебательными процессами посвящены многие исследования, в частности, работы [112].

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

Необходимость моделирования и управления распределенных колебательных процессов составных систем с кусочно-постоянными (неоднородными) характеристиками возникает во многих теоретических и прикладных областях науки и техники. Однако, научное направление по управлению неоднородными упругими колебаниями пока еще недостаточно исследовано, находится в стадии становления и по нему имеются лишь некоторые результаты. К исследованию решений задач управления подобных разнородных распределенных составных систем посвящены, в частности, работы [612]. В работах [6, 7] (и других работах этого же автора и его учеников) изучены подобные задачи граничного управления неоднородными колебательными процессами. При исследовании этих задач граничного управления использованы метод Даламбера и выведены формулы типа Даламбера. К краевым задачам для уравнения, описывающего процесс продольных колебаний стержня с кусочно-постоянными характеристиками (состоящими из не менее двух участков) со свободным либо закрепленным правым концом, посвящены, в частности, работы [1318]. В этих работах исследования проведены в классе обобщенного решения. В [19] рассматриваются колебания механической системы, представляющей собой два куска струны одинаковой длины, соединенных между собой пружиной. С помощью формулы Даламбера исследована задача граничного управления колебаниями сложносочлененной системы с особенностями. Получен аналог формулы Даламбера.

Цель данной статьи состоит в разработке аналитического подхода построения функции граничного управления одномерными колебательными неоднородными процессами, переводящего колебания за заданный промежуток времени из начального состояния в конечное состояние.

1. Постановка задачи. Рассматриваются продольные колебания кусочно-однородного стержня, расположенного вдоль отрезка $ - {{l}_{1}} \leqslant x \leqslant l$ и состоящего из двух участков: участок $ - {{l}_{1}} \leqslant x \leqslant 0$ с линейной плотностью ${{\rho }_{1}} = {\text{const}}$, с модулем Юнга ${{k}_{1}} = {\text{const}}$ и скоростью прохождения по нему волны, равной ${{a}_{1}} = \sqrt {\frac{{{{k}_{1}}}}{{{{\rho }_{1}}}}} $, и участка $0 \leqslant x \leqslant l$ с линейной плотностью ${{\rho }_{2}} = {\text{const}}$, с модулем Юнга ${{k}_{2}} = {\text{const}}$ и скоростью прохождения по нему волны, равной ${{a}_{2}} = \sqrt {\frac{{{{k}_{2}}}}{{{{\rho }_{2}}}}} $. Как и в работе [6], предполагается, что длины l1 и l участков стержневых систем выбраны так, что время прохождения волны по участку $ - {{l}_{1}} \leqslant x \leqslant 0$ совпадает со временем прохождения волны по участку $0 \leqslant x \leqslant l$, т.е.

(1.1)
$\frac{{{{l}_{1}}}}{{{{a}_{1}}}} = \frac{l}{{{{a}_{2}}}}$

Пусть состояние (продольные колебания) стержня (или поперечные колебания струны), описывается функцией Q(x, t), $ - {{l}_{1}} \leqslant x \leqslant l$, $0 \leqslant t \leqslant T$, а отклонения от состояния равновесия, подчиняются следующему волновому уравнению

(1.2)
$\frac{{{{\partial }^{2}}Q(x,t)}}{{\partial {{t}^{2}}}} = \left\{ \begin{gathered} a_{1}^{2}\frac{{{{\partial }^{2}}Q(x,t)}}{{\partial {{x}^{2}}}},\quad - {{l}_{1}} \leqslant x \leqslant 0,\quad 0 \leqslant t \leqslant T \hfill \\ a_{2}^{2}\frac{{{{\partial }^{2}}Q(x,t)}}{{\partial {{x}^{2}}}},\quad 0 \leqslant x \leqslant l,\quad 0 \leqslant t \leqslant T \hfill \\ \end{gathered} \right.$
с начальными (при $t = {{t}_{0}} = 0$) условиями
(1.3)
$Q(x,0) = {{\varphi }_{0}}(x),\quad {{\left. {\frac{{\partial Q(x,t)}}{{\partial t}}} \right|}_{{t = 0}}} = {{\psi }_{0}}(x),\quad - {{l}_{1}} \leqslant x \leqslant l$
с граничными условиями
(1.4)
$Q( - {{l}_{1}},t) = {{\mu }}(t),\quad Q(l,t) = {{\nu }}(t),\quad 0 \leqslant t \leqslant T$
и с условиями сопряжения в точке $x = 0$ соединения участков
(1.5)
$\begin{gathered} Q(0 - 0,t) = Q(0 + 0,t) \\ {{\left. {a_{1}^{2}{{\rho }_{1}}\frac{{\partial Q(x,t)}}{{\partial x}}} \right|}_{{x = 0 - 0}}} = {{\left. {a_{2}^{2}{{\rho }_{2}}\frac{{\partial Q(x,t)}}{{\partial x}}} \right|}_{{x = 0 + 0}}} \\ \end{gathered} $
и конечными (при $t = T$) условиями

(1.6)
$Q(x,T) = {{\varphi }_{T}}(x),\quad {{\left. {\frac{{\partial Q}}{{\partial t}}} \right|}_{{t = T}}} = {{\psi }_{T}}(x),\quad - {{l}_{1}} \leqslant x \leqslant l$

В формуле (1.4) функции ${{\mu }}(t)$ и ${{\nu }}(t)$ – управляющие воздействия (граничные управления). Предполагается, что функция $Q(x,t)\, \in \,{{C}^{2}}({{\Omega }_{T}})$, где ${{\Omega }_{T}} = \{ (x,t):\;x \in [ - {{l}_{1}},l],\;t \in [0,T]\} $, а функции ${{\varphi }_{0}}(x),\;{{\varphi }_{T}}(x) \in {{C}^{2}}\left[ { - {{l}_{1}},l} \right]$, ${{\psi }_{0}}(x),\,\,{{\psi }_{T}}(x) \in {{C}^{1}}[ - {{l}_{1}},l]$.

Предполагается также, что все функции такие, что выполняются условия согласования.

(1.7)
$\begin{gathered} {{\mu }}(0) = {{\varphi }_{0}}( - {{l}_{1}}),\quad {{\dot {\mu }}}(0) = {{\psi }_{0}}( - {{l}_{1}}),\quad {{\nu }}(0) = {{\varphi }_{0}}(l),\quad {{\dot {\nu }}}(0) = {{\psi }_{0}}(l) \hfill \\ {{\mu }}(T) = {{\varphi }_{T}}( - {{l}_{1}}),\quad {{\dot {\mu }}}(T) = {{\psi }_{T}}( - {{l}_{1}}),\quad {{\nu }}(T) = {{\varphi }_{T}}(l),\quad {{\dot {\nu }}}(T) = {{\psi }_{T}}(l) \hfill \\ \end{gathered} $

Задача граничного управления. Требуется найти граничные управления ${{\mu }}(t)$ и ${{\nu }}(t)$ $0 \leqslant t \leqslant T$, переводящие колебания системы (1.2), из заданного начального состояния (1.3) в конечное состояние (1.6).

2. Сведение задачи к задаче с нулевыми граничными условиями. Для решения поставленной задачи перейдем к новой переменной [20]

(2.1)
${{\xi }} = \left\{ \begin{gathered} \frac{{{{a}_{2}}}}{{{{a}_{1}}}}x,\quad - {{l}_{1}} \leqslant x \leqslant 0 \hfill \\ x,\quad 0 \leqslant x \leqslant l \hfill \\ \end{gathered} \right.$
что приводит к растяжению или сжатию отрезка $ - {{l}_{1}} \leqslant x \leqslant 0$ относительно точки $x = 0$. При этом с учетом (1.1), будем иметь, что отрезок $ - {{l}_{1}} \leqslant x \leqslant 0$ переходит к отрезку $ - l \leqslant {{\xi }} \leqslant 0$. Для функции $Q\left( {{{\xi }},t} \right)$ получим на отрезках одинаковой длины одинаковое уравнение
$\frac{{{{\partial }^{2}}Q({{\xi }},t)}}{{\partial {{t}^{2}}}} = \left\{ \begin{gathered} a_{2}^{2}\frac{{{{\partial }^{2}}Q({{\xi }},t)}}{{\partial {{{{\xi }}}^{2}}}},\quad - l \leqslant {{\xi }} \leqslant 0,\quad 0 \leqslant t \leqslant T \hfill \\ a_{2}^{2}\frac{{{{\partial }^{2}}Q({{\xi }},t)}}{{\partial {{{{\xi }}}^{2}}}},\quad 0 \leqslant {{\xi }} \leqslant l,\quad 0 \leqslant t \leqslant T \hfill \\ \end{gathered} \right.$
или
(2.2)
$\frac{{{{\partial }^{2}}Q({{\xi }},t)}}{{\partial {{t}^{2}}}} = a_{2}^{2}\frac{{{{\partial }^{2}}Q({{\xi }},t)}}{{\partial {{{{\xi }}}^{2}}}},\quad - l \leqslant {{\xi }} \leqslant l,\quad 0 \leqslant t \leqslant T$
с соответствующими начальными условиями
(2.3)
$Q({{\xi }},0) = {{\varphi }_{0}}({{\xi }}),\quad {{\left. {\frac{{\partial Q({{\xi }},t)}}{{\partial t}}} \right|}_{{t = 0}}} = {{\psi }_{0}}({{\xi }}),\quad - l \leqslant x \leqslant l$
с граничными условиями
(2.4)
$Q( - l,t) = {{\mu }}(t),\quad Q(l,t) = {{\nu }}(t),\quad 0 \leqslant t \leqslant T$
с конечными условиями
(2.5)
$Q({{\xi }},T) = {{\varphi }_{T}}({{\xi }}),\quad {{\left. {\frac{{\partial Q({{\xi }},t)}}{{\partial t}}} \right|}_{{t = T}}} = {{\psi }_{T}}({{\xi }}),\quad - l \leqslant {{\xi }} \leqslant l$
и с условиями сопряжения в точке ${{\xi }} = 0$ соединения участков

(2.6)
$Q(0 - 0,t) = Q(0 + 0,t),\quad {{a}_{1}}{{\rho }_{1}}{{\left. {\frac{{\partial Q({{\xi }},t)}}{{\partial {{\xi }}}}} \right|}_{{_{{{{\xi }} = 0 - 0}}}}} = {{\left. {{{a}_{2}}{{\rho }_{2}}\frac{{\partial Q({{\xi }},t)}}{{\partial {{\xi }}}}} \right|}_{{{{\xi }} = 0 + 0}}}$

Так как граничные условия (2.4) неоднородны, решение уравнения (2.2) построим в виде суммы

(2.7)
$Q({{\xi }},t) = V({{\xi }},t) + W({{\xi }},t)$
где $V({{\xi }},t)$ – функция с однородными граничными условиями
(2.8)
$V( - l,t) = V(l,t) = 0$
требующая определения, а $W({{\xi }},t)$ – решение уравнения (2.2) с неоднородными граничными условиями

(2.9)
$W( - l,t) = {{\mu }}(t),\quad W(l,t) = {{\nu }}(t)$

Функция $W({{\xi }},t)$ имеет вид

(2.10)
$W({{\xi }},t) = \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\mu }}(t) + (l + {{\xi }}){{\nu }}(t)} \right]$

Подставив (2.7) в (2.2) и учитывая (2.10), получим следующее уравнение для определения функции $V({{\xi }},t)$

(2.11)
$\frac{{{{\partial }^{2}}V({{\xi }},t)}}{{\partial {{t}^{2}}}} = a_{2}^{2}\frac{{{{\partial }^{2}}V({{\xi }},t)}}{{\partial {{{{\xi }}}^{2}}}} + F({{\xi }},t),\quad - l \leqslant {{\xi }} \leqslant l,\quad 0 \leqslant t \leqslant T$
где

(2.12)
$F({{\xi }},t) = \frac{1}{{2l}}\left[ {({{\xi }} - l){{\ddot {\mu }}}(t) - ({{\xi }} + l){{\ddot {\nu }}}(t)} \right]$

Функция $V({{\xi }},t)$ удовлетворяет соответствующему (2.6) условию сопряжения в точке ξ = 0 соединения участков. Отметим, что согласно (2.1) будем иметь

(2.13)
${{\varphi }_{0}}( - {{l}_{1}}) = {{\varphi }_{0}}( - l),\quad {{\varphi }_{T}}( - {{l}_{1}}) = {{\varphi }_{T}}( - l),\quad {{\psi }_{0}}( - {{l}_{1}}) = {{\psi }_{0}}( - l),\quad {{\psi }_{T}}( - {{l}_{1}}) = {{\psi }_{T}}( - l)$

В силу начальных и конечных условий, соответственно (2.3) и (2.5), функция $V(\xi ,t)$ должна удовлетворять следующим начальным

(2.14)
$\begin{gathered} V({{\xi }},0) = {{\varphi }_{0}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\mu }}(0) + (l + {{\xi }}){{\nu }}(0)} \right] \\ {{\left. {\frac{{\partial V({{\xi }},t)}}{{\partial t}}} \right|}_{{t = 0}}} = {{\psi }_{0}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\dot {\mu }}}(0) + (l + {{\xi }}){{\dot {\nu }}}(0)} \right] \\ \end{gathered} $
и конечным условиям

(2.15)
$\begin{gathered} V({{\xi }},T) = {{\varphi }_{T}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\mu }}(T) + (l + {{\xi }}){{\nu }}(T)} \right] \\ {{\left. {\frac{{\partial V({{\xi }},t)}}{{\partial t}}} \right|}_{{t = T}}} = {{\psi }_{T}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\dot {\mu }}}(T) + (l + {{\xi }}){{\dot {\nu }}}(T)} \right] \\ \end{gathered} $

С учетом условий (1.7) и (2.13), условия (2.14), (2.15) запишутся следующим образом, соответственно:

(2.16)
$\begin{gathered} V({{\xi }},0) = {{\varphi }_{0}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\varphi }_{0}}( - l) + (l + {{\xi }}){{\varphi }_{0}}(l)} \right] \\ {{\left. {\frac{{\partial V({{\xi }},t)}}{{\partial t}}} \right|}_{{t = 0}}} = {{\psi }_{0}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\psi }_{0}}( - l) + (l + {{\xi }}){{\psi }_{0}}(l)} \right] \\ V({{\xi }},T) = {{\varphi }_{T}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\varphi }_{T}}( - l) + (l + {{\xi }}){{\varphi }_{T}}(l)} \right] \\ \end{gathered} $
(2.17)
${{\left. {\frac{{\partial V({{\xi }},t)}}{{\partial t}}} \right|}_{{t = T}}} = {{\psi }_{T}}({{\xi }}) - \frac{1}{{2l}}\left[ {(l - {{\xi }}){{\psi }_{T}}( - l) + (l + {{\xi }}){{\psi }_{T}}(l)} \right]$

Таким образом, решение задачи сведено к задаче управления движением, описываемым уравнением (2.11) с однородными граничными условиями (2.8), которая формулируется следующим образом: требуется найти такие граничные управления ${{\mu }}(t)$ и ${{\nu }}(t)$, $0 \leqslant t \leqslant T$, под воздействием которых колебание, описываемое уравнением (2.11) с граничными условиями (2.8), переходит из заданного начального состояния (2.16) в конечное состояние (2.17).

3. Решение задачи. Учитывая, что граничные условия (2.8) однородны и выполнены условия согласованности, решение уравнения (2.11) ищем в виде

(3.1)
$V({{\xi }},t) = \sum\limits_{k = 1}^\infty {{{V}_{k}}(t){\text{sin}}\frac{{{{\pi }}k{{\xi }}}}{l}} ,\quad {{V}_{k}}(t) = \frac{1}{l}\int\limits_{ - l}^l {V({{\xi }},t){\text{sin}}\frac{{{{\pi }}k{{\xi }}}}{l}d{{\xi }}} $

Представим в виде рядов Фурье, в базисе $\left\{ {{\text{sin}}\frac{{{{\pi }}k{{\xi }}}}{l}} \right\}$ ($k = 1,\,\,2,...\,$), функции $F({{\xi }},t)$, ${{\varphi }_{0}}({{\xi }})$, ${{\varphi }_{T}}({{\xi }})$, ${{\psi }_{0}}({{\xi }})$ и ${{\psi }_{T}}({{\xi }})$ и, подставив их значения вместе с $V({{\xi }},t)$ в уравнения (2.11), (2.12) и в условия (2.16), (2.17), получим

(3.2)
${{\ddot {V}}_{k}}(t) + {{\lambda }}_{k}^{2}{{V}_{k}}(t) = {{F}_{k}}(t),\quad {{\lambda }}_{k}^{2} = {{\left( {\frac{{{{a}_{2}}{{\pi }}k}}{l}} \right)}^{2}}$
(3.3)
${{F}_{k}}(t) = \frac{{{{a}_{2}}}}{{{{\lambda }}{}_{k}l}}[{{\ddot {\nu }}}(t)(2{{( - 1)}^{k}} - 1) - {{\ddot {\mu }}}(t)]$
(3.4)
$\begin{gathered} {{V}_{k}}(0) = \varphi _{k}^{{(0)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{k}}l}}[{{\varphi }_{0}}( - l) - {{\varphi }_{0}}(l)(2{{( - 1)}^{k}} - 1)] \\ {{{\dot {V}}}_{k}}(0) = \psi _{k}^{{(0)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{k}}l}}[{{\psi }_{0}}( - l) - {{\psi }_{0}}(l)(2{{( - 1)}^{k}} - 1)] \\ \end{gathered} $
(3.5)
$\begin{gathered} {{V}_{k}}(T) = \varphi _{k}^{{(T)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{k}}l}}[{{\varphi }_{T}}( - l) - \,{{\varphi }_{T}}(l)(2{{( - 1)}^{k}} - 1)] \\ {{{\dot {V}}}_{k}}(T) = \psi _{k}^{{(T)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{k}}l}}[{{\psi }_{T}}( - l) - {{\psi }_{T}}(l)(2{{( - 1)}^{k}} - 1)] \\ \end{gathered} $

Здесь коэффициенты Фурье функции $F({{\xi }},t)$, ${{\varphi }_{0}}({{\xi }})$, ${{\varphi }_{T}}({{\xi }})$, ${{\psi }_{0}}({{\xi }})$ и ${{\psi }_{T}}({{\xi }})$ обозначены через ${{F}_{k}}(t)$, $\varphi _{k}^{{(0)}}$, $\varphi _{k}^{{(T)}}$, $\psi _{k}^{{(0)}}$ и $\psi _{k}^{{(T)}}$ соответствственно.

Общее решение неоднородного уравнения (3.2) с условиями (3.4) и его производная имеют вид

(3.6)
$\begin{gathered} {{V}_{k}}(t) = {{V}_{k}}(0){\text{cos}}{{{{\lambda }}}_{k}}t + \frac{1}{{{{{{\lambda }}}_{k}}}}{{{\dot {V}}}_{k}}(0){\text{sin}}{{{{\lambda }}}_{k}}t + \frac{1}{{{{{{\lambda }}}_{k}}}}\int\limits_0^t {{{F}_{k}}} (\tau ){\text{sin}}{{{{\lambda }}}_{k}}(t - \tau )d\tau \\ {{{\dot {V}}}_{k}}(t) = - {{{{\lambda }}}_{k}}{{V}_{k}}(0){\text{sin}}{{{{\lambda }}}_{k}}t + {{{\dot {V}}}_{k}}(0){\text{cos}}{{{{\lambda }}}_{k}}t + \int\limits_0^t {{{F}_{k}}} (\tau ){\text{cos}}{{{{\lambda }}}_{k}}(t - \tau )d\tau \\ \end{gathered} $

Далее, для формулы (3.6), при t = T, учитывая конечные (3.5) условия применяем подходы, приведенные в работах [25, 20]. Тогда получим, что функции управления ${{\mu }}(t)$ и ${{\nu }}(t)$ для каждого $k$ должны удовлетворять следующим интегральным соотношениям:

(3.7)
$\begin{gathered} \int\limits_0^T {{{\mu }}(\tau ){\text{sin}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)d\tau } + {{E}_{k}}\int\limits_0^T {{\nu }} (\tau ){\text{sin}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)d\tau = {{C}_{{1k}}} \\ \int\limits_0^T {{{\mu }}(\tau ){\text{cos}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)d\tau } + {{E}_{k}}\int\limits_0^T {{{\nu }}(\tau ){\text{cos}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)d\tau } = {{C}_{{2k}}} \\ \end{gathered} $
где

${{C}_{{1k}}} = \frac{1}{{{{\lambda }}_{k}^{2}}}\left[ {\frac{{{{{{\lambda }}}_{k}}l}}{{{{a}_{2}}}}{{{\tilde {C}}}_{{1k}}} + {{X}_{{1k}}} + {{E}_{k}}{{Y}_{{1k}}}} \right],\quad {{\tilde {C}}_{{1k}}} = {{{{\lambda }}}_{k}}{{V}_{k}}(T) - {{{{\lambda }}}_{k}}{{V}_{k}}(0){\text{cos}}{{{{\lambda }}}_{k}}T - {{\dot {V}}_{k}}(0){\text{sin}}{{{{\lambda }}}_{k}}T$
${{C}_{{2k}}} = \frac{1}{{{{\lambda }}_{k}^{2}}}\left[ {\frac{{{{{{\lambda }}}_{k}}l}}{{{{a}_{2}}}}{{{\tilde {C}}}_{{2k}}} + {{X}_{{2k}}} + {{E}_{k}}{{Y}_{{2k}}}} \right],\quad {{\tilde {C}}_{{2k}}} = {{\dot {V}}_{k}}(T) + {{{{\lambda }}}_{k}}{{V}_{k}}(0){\text{sin}}{{{{\lambda }}}_{k}}T - {{\dot {V}}_{k}}(0){\text{cos}}{{{{\lambda }}}_{k}}T$
(3.8)
${{X}_{{1k}}} = {{{{\lambda }}}_{k}}{{\varphi }_{T}}( - l) - {{\psi }_{0}}( - l){\text{sin}}{{{{\lambda }}}_{k}}T - {{{{\lambda }}}_{k}}{{\varphi }_{0}}( - l){\text{cos}}{{{{\lambda }}}_{k}}T,\quad {{E}_{k}} = 1 - 2{{( - 1)}^{k}}$
${{X}_{{2k}}} = {{\psi }_{T}}( - l) - {{\psi }_{0}}( - l){\text{cos}}{{{{\lambda }}}_{k}}T + {{{{\lambda }}}_{k}}{{\varphi }_{0}}( - l){\text{sin}}{{{{\lambda }}}_{k}}T$
$\begin{gathered} {{Y}_{{1k}}} = {{{{\lambda }}}_{k}}{{\varphi }_{T}}(l) - {{\psi }_{0}}(l){\text{sin}}{{{{\lambda }}}_{k}}T - {{{{\lambda }}}_{k}}{{\varphi }_{0}}(l){\text{cos}}{{{{\lambda }}}_{k}}T \\ {{Y}_{{2k}}} = {{\psi }_{T}}(l) - {{\psi }_{0}}(l){\text{cos}}{{{{\lambda }}}_{k}}T + {{{{\lambda }}}_{k}}{{\varphi }_{0}}(l){\text{sin}}{{{{\lambda }}}_{k}}T \\ \end{gathered} $

Введем следующие обозначения:

${{\bar {H}}_{k}}(\tau ) = \left( {\begin{array}{*{20}{c}} {{\text{sin}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)}&{{{E}_{k}}{\text{sin}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)} \\ {{\text{cos}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)}&{{{E}_{k}}{\text{cos}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)} \end{array}} \right),\quad {{C}_{k}} = \left( {\begin{array}{*{20}{c}} {{{C}_{{1k}}}} \\ {{{C}_{{2k}}}} \end{array}} \right),\quad U\left( \tau \right) = \left( {\begin{array}{*{20}{c}} {{{\mu }}(\tau )} \\ {{{\nu }}(\tau )} \end{array}} \right)$

Тогда соотношение (3.7) запишется следующим образом

(3.9)
$\int\limits_0^T {{{{\bar {H}}}_{k}}(\tau )} U(\tau )d\tau = {{C}_{k}},\quad k = 1,\;2,...$

Следовательно, для нахождения функции $U\left( \tau \right)$, $\tau \in \left[ {0,\,T} \right]$, получаются бесконечные интегральные соотношения (3.9).

На практике задача синтеза управления решается используя методы теории управления конечномерными системами [1, 21, 22]. Для первых n гармоник введем следующие обозначения блочных матриц

(3.10)
${{H}_{n}}(\tau ) = \left( {\begin{array}{*{20}{c}} {{{{\bar {H}}}_{1}}(\tau )} \\ {{{{\bar {H}}}_{2}}(\tau )} \\ \vdots \\ {{{{\bar {H}}}_{n}}(\tau )} \end{array}} \right),\quad {{{{\eta }}}_{n}} = \left( {\begin{array}{*{20}{c}} {{{C}_{1}}} \\ {{{C}_{2}}} \\ \vdots \\ {{{C}_{n}}} \end{array}} \right)$
с размерностями ${{H}_{n}}\left( \tau \right) - \left( {2n \cdot 2} \right)$, ${{{{\eta }}}_{n}} - \left( {2n \cdot 1} \right)$.

Для первых $n$ гармоник, с учетом (3.10), из (3.9) будем иметь

(3.11)
$\int\limits_0^T {{{H}_{n}}(\tau )} \,{{U}_{n}}(\tau )d\tau = {{{{\eta }}}_{n}}\,$
(здесь и далее обозначение в нижнем индексе буквы “n” будет означать – “для первых $n$ гармоник”).

Таким образом, из (3.11) следует, что первые n гармоники системы (3.2) с условиями (3.3)–(3.5) вполне управляемы тогда и только тогда, когда для любого вектора ${{{{\eta }}}_{n}}$ (3.10) можно найти управление ${{U}_{n}}(t)$, $t \in \left[ {0,\,T} \right]$, удовлетворяющее условию (3.11).

Управляющее воздействие $\,{{U}_{n}}(t)$, удовлетворяющее интегральному соотношению (3.11), представим в виде [21, 22]

(3.12)
$\,{{U}_{n}}(t) = H_{n}^{T}(t)S_{n}^{{ - 1}}{{{{\eta }}}_{n}} + {{f}_{n}}(t)$
где $H_{n}^{T}(t)$ – транспонированная матрица, ${{f}_{n}}(t)$ – вектор-функция и такая, что

(3.13)
$\int\limits_0^T {{{H}_{n}}(t)} {{f}_{n}}(t)dt = 0,\quad {{S}_{n}} = \int\limits_0^T {{{H}_{n}}\left( t \right)H_{n}^{T}\left( t \right)dt} $

Здесь, Sn – известная матрица размерностью $\left( {2n \cdot 2n} \right)$, для которой предполагается, что $\det {{S}_{n}} \ne 0$.

Из формулы (3.12) следует, что существует множество управляющих функций, решающих задачу граничного управления.

Обозначим

$S_{n}^{{ - 1}} = \left( {\begin{array}{*{20}{c}} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(11)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(11)}}}&{...}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(1n)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(1n)}}} \\ {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(11)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(11)}}}&{...}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{21}}^{{(11)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{22}}^{{(11)}}} \\ {...}&{}&{...}&{}&{...} \\ {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(n1)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(n1)}}}&{...}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(nn)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(nn)}}} \\ {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{21}}^{{(n1)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{22}}^{{(n1)}}}&{...}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{21}}^{{(nn)}}}&{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{22}}^{{(nn)}}} \end{array}} \right)$

Из формулы (3.10), (3.12) с учетом введенных обозначений функции граничного управления представляются в виде

(3.14)
$\begin{gathered} {{{{\mu }}}_{n}}(t) = \sum\limits_{k = 1}^n {\left[ {S_{1}^{{(k)}}{\text{sin}}{{{{\lambda }}}_{k}}\left( {T - t} \right) + S_{2}^{{(k)}}{\text{cos}}{{{{\lambda }}}_{k}}\left( {T - t} \right)} \right]} + {{f}_{n}}(t) \\ {{{{\nu }}}_{n}}(t) = \sum\limits_{k = 1}^n {{{E}_{k}}\left[ {S_{1}^{{(k)}}{\text{sin}}{{{{\lambda }}}_{k}}\left( {T - t} \right) + S_{2}^{{(k)}}{\text{cos}}{{{{\lambda }}}_{k}}\left( {T - t} \right)} \right]} + {{f}_{n}}(t) \\ \end{gathered} $
где

$S_{1}^{{(k)}} = \mathop \sum \limits_{j = 1}^n \left[ {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(k\,j)}}{{C}_{{1j}}} + \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(k\,j)}}{{C}_{{2j}}}} \right],\quad S_{2}^{{(k)}} = \mathop \sum \limits_{j = 1}^n \left[ {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{21}}^{{(k\,j)}}{{C}_{{1j}}} + \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{22}}^{{(k\,j)}}{{C}_{{2j}}}} \right]$

Подставляя построенные выражения для функции ${{{{\mu }}}_{n}}(t)$ и ${{{{\nu }}}_{n}}(t)$ в (3.3), а найденное для ${{F}_{k}}(t)$ выражение – в (3.6), получим функцию ${{V}_{k}}\left( t \right)$, $t \in \left[ {0,\,T} \right]$. Далее, из формулы (3.1) будем иметь

(3.15)
${{V}_{n}}({{\xi }},t) = \sum\limits_{k = 1}^n {{{V}_{k}}} (t){\text{sin}}\frac{{{{\pi }}k}}{l}{{\xi }}$
а с помощью (2.7) и (2.10) функция колебания ${{Q}_{n}}\,({{\xi }},t)$ для первых $n$ гармоник запишется в виде
(3.16)
${{Q}_{n}}({{\xi }},t) = {{V}_{n}}({{\xi }},t) + {{W}_{n}}({{\xi }},t)$

(3.17)
${{W}_{n}}({{\xi }},t) = \frac{1}{{2l}}\left[ {(l - {{\xi }}){{{{\mu }}}_{n}}(t) + (l + {{\xi }}){{{{\nu }}}_{n}}(t)} \right]$

Учитывая обозначения (2.1) функция ${{Q}_{n}}\,(x,t)$ при $ - {{l}_{1}} \leqslant x \leqslant l$ представляется в виде:

(3.18)
${{Q}_{n}}(x,t) = \left\{ \begin{gathered} \sum\limits_{k = 1}^n {{{V}_{k}}} (t){\text{sin}}\frac{{{{\pi }}k}}{{{{l}_{1}}}}x + \frac{1}{2}\left[ {\left( {1 - \frac{x}{{{{l}_{1}}}}} \right){{{{\mu }}}_{n}}(t) + \left( {1 + \frac{x}{{{{l}_{1}}}}} \right){{{{\nu }}}_{n}}(t)} \right] \hfill \\ - {{l}_{1}} \leqslant x \leqslant 0,\quad 0 \leqslant t \leqslant T \hfill \\ \sum\limits_{k = 1}^n {{{V}_{k}}} (t){\text{sin}}\frac{{{{\pi }}k}}{l}x + \frac{1}{2}\left[ {\left( {1 - \frac{x}{l}} \right){{{{\mu }}}_{n}}(t) + \left( {1 + \frac{x}{l}} \right){{{{\nu }}}_{n}}(t)} \right] \hfill \\ 0 \leqslant x \leqslant l,\quad 0 \leqslant t \leqslant T \hfill \\ \end{gathered} \right.$

4. Пример. Для иллюстрации вышеизложенного построения предположим, что в граничных условиях (1.4) $Q(l,t) = 0$, $0 \leqslant t \leqslant T$ (т.е. ${{\nu }}(t) = 0$).

В этом случае из формулы (3.3) следует ${{F}_{k}}(t) = - \frac{{2{{a}_{2}}}}{{{{\lambda }}{}_{k}l}}{{\mu }}{\kern 1pt} '{\kern 1pt} '(t)$, а согласно формулам (3.7) будем иметь следующие интегральные соотношения

$\int\limits_0^T {{{\mu }}(\tau ){\text{sin}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)d\tau } = {{C}_{{1k}}}$
$\int\limits_0^T {{{\mu }}(\tau ){\text{cos}}{{{{\lambda }}}_{k}}\left( {T - \tau } \right)d\tau } = {{C}_{{2k}}}$

${{C}_{{1k}}} = \frac{1}{{{{\lambda }}_{k}^{2}}}\left[ {\frac{{{{{{\lambda }}}_{k}}l}}{{{{a}_{2}}}}{{{\tilde {C}}}_{{1k}}} + {{X}_{{1k}}}} \right],\quad {{C}_{{2k}}} = \frac{1}{{{{\lambda }}_{k}^{2}}}\left[ {\frac{{{{{{\lambda }}}_{k}}l}}{{{{a}_{2}}}}{{{\tilde {C}}}_{{2k}}} + {{X}_{{2k}}}} \right]$

Постоянные ${{\tilde {C}}_{{1k}}}$, ${{\tilde {C}}_{{2k}}}$, ${{X}_{{1k}}}$ и ${{X}_{{2k}}}$определяются из формулы (3.8). Следовательно

${{\bar {H}}_{k}}(\tau ) = \left( \begin{gathered} {\text{sin}}{{{{\lambda }}}_{k}}(T - \tau ) \hfill \\ {\text{cos}}{{{{\lambda }}}_{k}}(T - \tau ) \hfill \\ \end{gathered} \right),\quad {{C}_{k}} = \left( \begin{gathered} {{C}_{{1k}}} \hfill \\ {{C}_{{2k}}} \hfill \\ \end{gathered} \right),\quad k = 1,\;2,...$

Для простоты изложения, согласно формуле (3.12) (или (3.14)), (3.13), построим функцию граничного управления ${{{{\mu }}}_{n}}(t)$ при $n = 1$ (следовательно $k = 1$). Тогда, согласно (3.10), будем иметь ${{H}_{1}}(\tau ) = {{\bar {H}}_{1}}(\tau )$, ${{{{\eta }}}_{1}} = {{C}_{1}}$, а из (3.13), элементы матрицы

${{S}_{1}} = \left( {\begin{array}{*{20}{c}} {{{s}_{{11}}}}&{{{s}_{{12}}}} \\ {{{s}_{{21}}}}&{{{s}_{{22}}}} \end{array}} \right)$
получим в виде:
${{s}_{{11}}} = \frac{T}{2} - \frac{1}{{4{{{{\lambda }}}_{1}}}}{\text{sin}}2{{{{\lambda }}}_{1}}T,\quad {{s}_{{12}}} = {{s}_{{21}}} = \frac{1}{{2{{{{\lambda }}}_{1}}}}{\text{si}}{{{\text{n}}}^{2}}{{{{\lambda }}}_{1}}T,\quad {{s}_{{22}}} = \frac{T}{2} + \frac{1}{{4{{{{\lambda }}}_{1}}}}{\text{sin}}2{{{{\lambda }}}_{1}}T$
при этом $\Delta = \det {{S}_{1}} \ne 0$.

Из формулы (3.12) (или (3.14)) следует, что ${{{{\mu }}}_{1}}(\tau ) = H_{1}^{T}(\tau )S_{1}^{{ - 1}}{{{{\eta }}}_{1}} + {{f}_{1}}(\tau )$. Предполагая, что ${{f}_{1}}(\tau ) = 0$, получим, при $\tau \in [0,\,\,T]$

${{{{\mu }}}_{1}}(\tau ) = {\text{sin}}{{{{\lambda }}}_{1}}\left( {T - \tau } \right)[\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{11}}^{{(11)}}{{C}_{{11}}} + \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(11)}}{{C}_{{21}}}] + {\text{cos}}{{{{\lambda }}}_{1}}\left( {T - \tau } \right)[\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{21}}^{{(11)}}{{C}_{{11}}} + \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{22}}^{{(11)}}{{C}_{{21}}}]$

Отметим, что согласно формулам (3.15)(3.16) будем иметь

${{Q}_{1}}({{\xi }},t) = {{V}_{1}}(t){\text{sin}}\frac{{{\pi }}}{l}{{\xi }} + \frac{1}{{2l}}(l - {{\xi }}){{{{\mu }}}_{1}}(t)$

Предположим, что $T = 2\frac{l}{{{{a}_{2}}}}$. Тогда, с учетом ${{{{\lambda }}}_{1}} = \frac{{{{a}_{2}}{{\pi }}}}{l}$, получим $T{{{{\lambda }}}_{1}} = 2{{\pi }}$. Для матриц ${{S}_{1}}$ и $S_{1}^{{ - 1}}$ будем иметь:

${{S}_{1}} = \left( {\begin{array}{*{20}{c}} {\frac{l}{{{{a}_{2}}}}}&0 \\ 0&{\frac{l}{{{{a}_{2}}}}} \end{array}} \right),\quad S_{1}^{{ - 1}} = \left( {\begin{array}{*{20}{c}} {\frac{{{{a}_{2}}}}{l}}&0 \\ 0&{\frac{{{{a}_{2}}}}{l}} \end{array}} \right),\quad \det {{S}_{1}} = {{\left( {\frac{l}{{{{a}_{2}}}}} \right)}^{2}}$

Учитывая, что $\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{12}}^{{(11)}} = \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{s} _{{21}}^{{(11)}} = 0$, получим

${{C}_{{11}}} = \frac{l}{{{{a}_{2}}}}\left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right),\quad {{C}_{{21}}} = \frac{l}{{{{a}_{2}}{{{{\lambda }}}_{1}}}}\left( {\psi _{1}^{{(T)}} - \psi _{1}^{{(0)}}} \right)$
${{{{\mu }}}_{1}}(\tau ) = \frac{1}{{{{{{\lambda }}}_{1}}}}\left( {\psi _{1}^{{(T)}} - \psi _{1}^{{(0)}}} \right){\text{cos}}{{{{\lambda }}}_{1}}\tau - \left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right){\text{sin}}{{{{\lambda }}}_{1}}\tau $
следовательно

${{F}_{1}}(\tau ) = \frac{{{{a}_{2}}}}{l}\left( {\psi _{1}^{{(T)}} - \psi _{1}^{{(0)}}} \right){\text{cos}}{{{{\lambda }}}_{1}}\tau - \frac{{{{a}_{2}}{{{{\lambda }}}_{1}}}}{l}\left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right){\text{sin}}{{{{\lambda }}}_{1}}\tau $

Подставляя найденное для ${{F}_{1}}(\tau )$ выражение в (3.6), получим функцию ${{V}_{1}}\left( t \right)$, $t \in \left[ {0,\,T} \right]$ в явном виде

$\begin{gathered} {{V}_{1}}(t) = \left[ {\varphi _{1}^{{(0)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{1}}l}}{{\varphi }_{0}}( - l) + \frac{{{{a}_{2}}}}{{2l}}\left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right)t} \right]{\text{cos}}{{{{\lambda }}}_{1}}t + \\ \; + \frac{1}{{{{{{\lambda }}}_{1}}}}\left[ {\psi _{1}^{{(0)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{1}}l}}{{\psi }_{1}}( - l) - \frac{{{{a}_{2}}}}{{2l}}\left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right) + \frac{{{{a}_{2}}}}{{2l}}\left( {\psi _{1}^{{(T)}} - \psi _{1}^{{(0)}}} \right)t} \right]{\text{sin}}{{{{\lambda }}}_{1}}t \\ \end{gathered} $

Следовательно, для функций ${{Q}_{1}}\,({{\xi }},t)$, получим следующее явное выражение

$\begin{gathered} {{Q}_{1}}({{\xi }},t) = \left\{ {\left[ {\varphi _{1}^{{(0)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{1}}l}}{{\varphi }_{0}}( - l) + \frac{{{{a}_{2}}}}{{2l}}\left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right)t} \right]{\text{cos}}{{{{\lambda }}}_{1}}t + } \right. \\ \;\left. { + \frac{1}{{{{{{\lambda }}}_{1}}}}\left[ {\psi _{1}^{{(0)}} - \frac{{{{a}_{2}}}}{{{{{{\lambda }}}_{1}}l}}{{\psi }_{1}}( - l) - \frac{{{{a}_{2}}}}{{2l}}\left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right) + \frac{{{{a}_{2}}}}{{2l}}\left( {\psi _{1}^{{(T)}} - \psi _{1}^{{(0)}}} \right)t} \right]{\text{sin}}{{{{\lambda }}}_{1}}t} \right\}{\text{sin}}\frac{{{\pi }}}{l}{{\xi }} + \\ \; + \frac{1}{{2l}}\left( {l - {{\xi }}} \right)\left[ {\frac{1}{{{{{{\lambda }}}_{1}}}}\left( {\psi _{1}^{{(T)}} - \psi _{1}^{{(0)}}} \right){\text{cos}}{{{{\lambda }}}_{1}}t - \left( {\varphi _{1}^{{(T)}} - \varphi _{1}^{{(0)}}} \right){\text{sin}}{{{{\lambda }}}_{1}}t} \right] \\ \end{gathered} $

Имея явные виды функции ${{{{\mu }}}_{1}}(t)$, ${{V}_{1}}(t)$ или ${{Q}_{1}}\,({{\xi }},t)$ с помощью формул (3.15)–(3.18) или формулы замены переменных (2.1) можно записать явное выражение функции ${{Q}_{1}}\,(x,t)$, $ - {{l}_{1}} \leqslant x \leqslant l$, $0 \leqslant t \leqslant T$.

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

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

  1. Бутковский А.Г. Методы управления системами с распределенными параметрами. М.: Наука, 1975. 568 с.

  2. Barseghyan V.R. The control problem for stepwise changing linear systems of loaded differential equations with unseparated multipoint intermediate conditions // Mech. Solids. 2018. V. 53. № 6. P. 615–622. https://doi.org/10.21538/0134-4889-2019-25-3-24-33

  3. Barseghyan V.R. The problem of optimal control of string vibrations // Int. Appl. Mech. 2020. V. 56. № 4. P. 471–480. https://doi.org/10.1007/s10778-020-01030-w

  4. Barseghyan V.R. The problem of optimal control of vibrations of a string with non-separated conditions into state functions at given intermediate time instants // Autom. Remote Control. 2020. V. 81. № 2. P. 226–235. https://doi.org/10.31857/S0005231020020038

  5. Barseghyan V., Solodusha S. Optimal boundary control of string vibrations with given shape of deflection at a certain moment of time // Mathematical Optimization Theory and Operations Research. MOTOR 2021. Lecture Notes in Computer Science. V. 12755 / Ed. by P. Pardalos, M. Khachay, A. Kazakov. Cham: Springer, 2021. P. 299–313. https://doi.org/10.1007/978-3-030-77876-7_20

  6. Ильин В.А. Оптимизация граничного управления колебаниями стержня, состоящего из двух разнородных участков // Доклады РАН. 2011. Т. 440. № 2. С. 159–163.

  7. Ильин В.А. О приведении в произвольно заданое состояние колебаний первоначально покоящегося стержня, состоящего из двух разнородных участков // Доклады РАН. 2010. Т. 435. № 6. С. 732–735.

  8. Егоров А.И., Знаменская Л.Н. Об управляемости упругих колебаний последовательно соединенных объектов с распределенными параметрами // Тр. ИММ УрО РАН. 2011. Т. 17. № 1. С. 85–92.

  9. Провоторов В.В. Построение граничных управлений в задаче о гашении колебаний системы струн // Вестн. Санкт-Петерб. ун-та. Сер. 1. С. 62–71.Вып. 1.

  10. Amara J. Ben, Bouzidi H. Null boundary controllability of a one-dimensional heat equation with an internal point mass and variable coefficients // J. Math. Phys. 2018. V. 59. № 1. P. 011512. https://doi.org/10.1063/1.5021947

  11. Amara J. Ben, Beldi E. Boundary controllability of two vibrating strings connected by a point mass with variable coefficients // SIAM J. Control Optim. 2019. V. 57. № 5. P. 3360–3387. https://doi.org/10.1137/16M1100496

  12. Mercier D., Régnier V. Boundary controllability of a chain of serially connected Euler-Bernoulli beams with interior masses // Collect. Math. 2009. V. 60. № 3. P. 307–334. https://doi.org/10.1007/BF03191374

  13. Кулешов А.А. Смешанные задачи для уравнения продольных колебаний неоднородного стержня и уравнения поперечных колебаний неоднородной струны, состоящих из двух участков разной плотности и упругости // Доклады РАН. 2012. Т. 442. № 5. С. 594–597.

  14. Кулешов A.A. Смешанные задачи для уравнения продольных колебаний неоднородного стержня со свободным либо закрепленным правым концом, состоящего из двух участков разной плотности и упругости // Доклады РАН. 2012. Т. 442. № 4. С. 451–454.

  15. Рогожников А.М. Исследование смешанной задачи, описывающей процесс колебаний стержня, состоящего из нескольких участков, при условии совпадения времени прохождения волны по каждому из этих участков // Доклады РАН. 2011. Т. 441. № 4. С. 449–451.

  16. Рогожников А.М. Исследование смешанной задачи, описывающей процесс колебаний стержня, состоящего из нескольких участков с произвольными длинами // Доклады РАН. 2012. Т. 444. С. 488–491.

  17. Аниконов Д.С., Коновалова Д.С. Прямая и обратная задачи для волнового уравнения с разрывными коэффициентами // Науч.-тех. ведом. СПбГПУ. Физ.-мат. науки. 2018. Т. 11. № 2. С. 61–72.

  18. Смирнов И.Н. Смешанные задачи для телеграфного уравнения, в случае системы, состоящей из двух участков, имеющих разные плотности и разные упругости, но одинаковые импедансы // Доклады РАН. 2010. Т. 435. № 2. С. 172–177.

  19. Зверева М.Б., Найдюк Ф.О., Залукаева Ж.О. Моделирование колебаний сингулярной струны // Вестн. Воронеж. гос. ун-та. Сер.: Физ., мат. 2014. № 2. С. 111–119.

  20. Холодовский С.Е., Чухрий П.А. Задача о движении неограниченной кусочно-однородной струны // Уч. зап. Забайкальск. гос. ун-та. Сер. Физ. мат. техн. технол. 2018. Т. 13. № 4. С. 42–50. https://doi.org/10.21209/2308-8761-2018-13-4-42-50

  21. Барсегян В.Р. Управление составных динамических систем и систем с многоточечными промежуточными условиями. М.: Наука, 2016. 230 с.

  22. Зубов В.И. Лекции по теории управления. М.: Наука, 1975. 496 с.

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