Прикладная математика и механика, 2023, T. 87, № 4, стр. 661-669
Обратные задачи для уравнения колебаний консольной балки по отысканию источника
1 Самарский государственный технический университет
Самара, Россия
* E-mail: faoks@yandex.ru
Поступила в редакцию 15.05.2023
После доработки 15.06.2023
Принята к публикации 20.06.2023
- EDN: MMTCAD
- DOI: 10.31857/S0032823523040057
Аннотация
Для уравнения колебания балки изучаются обратные задачи по отысканию правой части, т.е. источника колебаний. Решения задач методами спектрального анализа и интегральных уравнений Вольтерра построены в явном виде как суммы рядов и доказаны соответствующие теоремы единственности и существования. При обосновании существования решения обратной задачи по определению сомножителя правой части, зависящей от пространственной координаты, возникает проблема малых знаменателей. В связи с этим установлены оценки знаменателей, гарантирующие их отделенность от нуля, с указанием соответствующей асимптотики. На основании этих оценок обоснована сходимость рядов в классе регулярных решений уравнения колебаний балки.
1. Введение. В строительной механике, авиастроении, машиностроении и других областях значимую роль играют задачи о колебаниях балок, стержней и пластин. Описание таких колебательных процессов приводит к дифференциальным уравнениям более высокого порядка, чем уравнение колебания струны ([1], с. 143–145), ([2], с. 276–277), ([3], с. 314–315), ([4], с. 34), ([5], с. 75–77).
Рассмотрим однородную балку длины $l$, один конец которой свободен, а другой наглухо заделан. Ее вынужденные изгибные поперечные колебания под действием непрерывной внешней силы $G(x,t)$, при отсутствии вращательного движения, описываются уравнением
где $\rho $ – линейная плотность балки, $S$ – площадь поперечного сечения, $E$ – модуль упругости материала, $J$ – момент инерции сечения относительно своей горизонтальной оси. Это уравнение можно переписать в виде где ${{\alpha }^{2}} = EJ{\text{/}}\rho S$, $F(x,t) = G(x,t){\text{/}}\rho S$.Отметим, что при изучении задач расчета устойчивости вращающихся валов и вибрации кораблей так же приходят к уравнению (1.1) ([6], с. 326).
В данной работе рассматривается уравнение (1.1) в области
где $l$ и $T$ – заданные положительные действительные числа, с граничными условиями, соответствующими консольной балке,(1.2)
$u(0,t) = {{u}_{x}}(0,t) = {{u}_{{xx}}}(l,t) = {{u}_{{xxx}}}(l,t) = 0;\quad 0 \leqslant t \leqslant T,$В данной работе исследуются следующие задачи.
Задача 1. В области $D$ найти решение $u(x,t)$ уравнения (1.1), такое, что
(1.4)
$u(x,t) \in C_{{x,t}}^{{4,2}}\left( D \right) \cap C_{{x,t}}^{{2,1}}\left( {\overline D } \right),$Задача 2. Пусть $F(x,t) = f(x)g(t)$. Найти функции $f(x)$ и $u(x,t)$, удовлетворяющие уравнению (1.1) в области $D$, такие, что $f(x) \in С\left[ {0,l} \right]$, а функция $u(x,t)$ удовлетворяет условиям (1.2)–(1.4) и, кроме того,
где $g(t)$, $\varphi (x)$, $\tau (x)$ и $\nu (x)$ – заданные достаточно гладкие функции, ${{t}_{0}}$ – заданная точка из $\left( {0,\,T} \right]$.Задача 3. Пусть $F(x,t) = f(x)g(t)$. Найти функции $g(t)$ и $u(x,t)$, удовлетворяющие уравнению (1.1) в области $D$, такие, что $g(t) \in С\left[ {0,T} \right]$, а функция $u(x,t)$ удовлетворяет условиям (1.2)–(1.4) и, кроме того,
где $f(x)$, $h(t),$ $\tau (x)$ и $\nu (x)$ – заданные достаточно гладкие функции, ${{x}_{0}}$ – заданная точка из $\left( {0,l} \right)$.Из постановок задач видно, что задача 1 представляет собой прямую начально-граничную задачу для неоднородного уравнения колебаний балки (1.1). Задачи 2 и 3 являются обратными, поэтому условия (1.5) и (1.6) являются дополнительными для определения соответственно сомножителей $f(x)$ и $g(t)$ правой части $F(x,t)$ уравнения (1.1).
Данная работа является продолжением исследований автора [7], посвященных обоснованию корректности постановки начально-граничной задачи 1 для уравнения (1.1). Основное внимание здесь уделяется задачам 2 и 3, постановки которых аналогичны работам [8–12], где похожие задачи изучались для уравнений теплопроводности, колебаний струны и других дифференциальных уравнений второго порядка и высоких порядков.
В данной работе на основе прямой задачи 1, изученной в [7], доказаны теоремы единственности и существования решений обратных задач 2 и 3. При этом решения построены в виде суммы рядов. Отметим, что при обосновании сходимости рядов в задаче 2 возникает проблема малых знаменателей, создающая дополнительные трудности. Установлены оценки, гарантирующие отделенность от нуля знаменателей, с указанием соответствующей асимптотики. На основании этих оценок обоснована сходимость рядов в классе регулярных решений уравнения (1.1).
2. Построение решения прямой задачи. В этом пункте кратко приведены результаты исследования задачи 1, полученные в [7], где методом интегралов энергии доказана единственность решения этой задачи. Методом разделения переменных построено решение задачи (1.2)–(1.4) для уравнения (1.1) в явном виде как суммы ряда по системе собственных функций следующей спектральной задачи:
Собственные значения этой спектральной задачи находятся по формуле ${{\lambda }_{n}} = - d_{n}^{4}$, где ${{d}_{n}}$ – корни уравнения $chdl \cdot \cos dl = - 1$, для которых справедлива асимптотическая формула
(2.1)
${{d}_{n}} = \frac{\pi }{l}\left( {n - \frac{1}{2} + {{{( - 1)}}^{n}}{{\Theta }_{n}}} \right),\quad {{\Theta }_{n}} \in \left( {0,\frac{1}{2}} \right),\quad {{\Theta }_{n}} = O\left( {\frac{1}{{{{n}^{2}}}}} \right)$Система собственных функций имеет вид
Нормируя эту систему, получаем
Отметим, что система функций ${{Y}_{n}}(x)$ полна и образует ортонормированный базис в пространстве ${{L}_{2}}\left[ {0,l} \right]$.
Тогда решение задачи 1 определяется в виде суммы ряда
где ${{u}_{n}}(t)$ определяются по формуле и после вычисления принимают вид(2.3)
${{u}_{n}}(t) = {{\tau }_{n}}\cos \alpha d_{n}^{2}t + \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}t + \frac{1}{{\alpha d_{n}^{2}}}\int\limits_0^t {{{F}_{n}}(s)\sin \left[ {\alpha d_{n}^{2}(t - s)} \right]ds} ,$Справедливо следующее утверждение.
Теорема 1. Если функция $\tau (x)$ принадлежит классу ${{C}^{6}}[0,l]$, $\tau (0) = \tau {\kern 1pt} '(0)$ = $\tau {\kern 1pt} '{\kern 1pt} '(l) = \tau {\kern 1pt} '{\kern 1pt} '{\kern 1pt} '{\kern 1pt} (l)$ = = ${{\tau }^{{(4)}}}(0) = {{\tau }^{{(5)}}}(0)$ = 0, функция $\nu (x)\;$ – классу ${{C}^{4}}[0,l]$, $\nu (0) = \nu {\kern 1pt} '(0)$ = $\nu {\kern 1pt} '{\kern 1pt} '(l) = \nu {\kern 1pt} '{\kern 1pt} '{\kern 1pt} '(l)$ = 0, а функция $F(x,t)$ – классу $C\left( {\bar {D}} \right) \cap C_{x}^{4}\left( {\bar {D}} \right)$ и $F(0,t) = {{F}_{x}}(0,t)$ = ${{F}_{{xx}}}(l,t) = {{F}_{{xxx}}}(l,t)$ = 0 при $0 \leqslant t \leqslant T$, то существует единственное решение задачи 1 и оно определяется рядом (2.2).
Полное доказательство этой теоремы приведено в работе [7].
3. Исследование обратной задачи 2. Пусть $F(x,t) = f(x)g(t)$. В силу теоремы 1 функции $f(x)$ и $g(t)$ должны удовлетворять условиям:
Тогда функции ${{F}_{n}}(t)$ и ${{u}_{n}}(t)$ принимают вид
(3.1)
${{F}_{n}}(t) = g(t){{f}_{n}},\quad {{u}_{n}}(t) = {{\tau }_{n}}\cos \alpha d_{n}^{2}t + \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}t + {{f}_{n}}{{g}_{n}}(t),$(3.2)
${{f}_{n}} = \int\limits_0^l {f(x){{Y}_{n}}(x)dx} ,\quad {{g}_{n}}(t) = \frac{1}{{\alpha d_{n}^{2}}}\int\limits_0^t {g(s)\sin \left[ {\alpha d_{n}^{2}\left( {t - s} \right)} \right]ds} $Удовлетворяя функцию (2.2) условию (1.5), получим уравнение
Тогда с учетом (2.3) найдем
(3.3)
${{f}_{n}} = \frac{1}{{{{g}_{n}}({{t}_{0}})}}\left( {{{\varphi }_{n}} - {{\tau }_{n}}\cos \alpha d_{n}^{2}{{t}_{0}} - \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}{{t}_{0}}} \right),$Подставляя найденные значения ${{f}_{n}}$ в равенство (2.3), построим в явном виде функции
(3.4)
$\begin{gathered} {{u}_{n}}(t) = {{\tau }_{n}}\cos \alpha d_{n}^{2}t + \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}t + \frac{{{{g}_{n}}(t)}}{{{{g}_{n}}({{t}_{0}})}}\left( {{{\varphi }_{n}} - {{\tau }_{n}}\cos \alpha d_{n}^{2}{{t}_{0}} - \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}{{t}_{0}}} \right) = \\ = {{\tau }_{n}}\left( {\cos \alpha d_{n}^{2}t - \frac{{{{g}_{n}}(t)}}{{{{g}_{n}}({{t}_{0}})}}\cos \alpha d_{n}^{2}{{t}_{0}}} \right) + \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\left( {\sin \alpha d_{n}^{2}t - \frac{{{{g}_{n}}(t)}}{{{{g}_{n}}({{t}_{0}})}}\sin \alpha d_{n}^{2}{{t}_{0}}} \right) + {{\varphi }_{n}}\frac{{{{g}_{n}}(t)}}{{{{g}_{n}}({{t}_{0}})}} \\ \end{gathered} $Тогда решение задачи 2 находится как сумма рядов
(3.5)
$u(x,t) = \sum\limits_{n = 1}^\infty {{{u}_{n}}(t){{Y}_{n}}(x)} \quad {\text{и}}\quad f(x) = \sum\limits_{n = 1}^\infty {{{f}_{n}}{{Y}_{n}}(x)} ,$Из однозначного характера построения решения задачи 2 следует его единственность. Действительно, пусть ${{g}_{n}}({{t}_{0}}) \ne 0$ при всех $n \in N$ и $\varphi (x) = \tau (x) = \nu (x) \equiv 0$. Тогда из (3.3) и (3.4) вытекает, что ${{f}_{n}} = 0$ и ${{u}_{n}}(t) \equiv 0$ при всех $n \in N$. В силу чего из формул (2.2) и (3.2) следуют равенства
из которых в силу полноты системы функций ${{Y}_{n}}(x)$ на ${{L}_{2}}\left[ {0,l} \right]$ следует, что $u(x,t) = 0$ и $f(x) = 0$ почти всюду на $\left[ {0,l} \right]$ при любом $t \in \left[ {0,T} \right]$. Тогда из условия (1.4) и непрерывности функции $f(x)$ на $\left[ {0,l} \right]$ следует, что $u(x,t) \equiv 0$ в $\bar {D}$ и $f(x) \equiv 0$ на $\left[ {0,l} \right]$ при любой функции $g(t) \in C\left[ {0,T} \right]$.Если при некотором ${{t}_{0}}$ и $n = m$ выражение ${{g}_{m}}({{t}_{0}}) = 0$, то однородная задача 2 (при $\varphi (x) = \tau (x) = {v}(x) \equiv 0$) при любой непрерывной функции $g(t)$ имеет ненулевое решение
Выражение ${{g}_{n}}({{t}_{0}})$ за счет $\sin \left( {\alpha d_{m}^{2}(t - s)} \right)$ может обращаться в нуль не более чем в счетном числе точек.
Поэтому возникает проблема малых знаменателей и необходимо установить оценки величин ${{g}_{n}}({{t}_{0}})$, гарантирующие их отдаленность от нуля, и указать асимптотику этих оценок при больших $n$.
Пусть $g(s)$ монотонна на отрезке $\left[ {0,{{t}_{0}}} \right]$. Тогда на основании второй теоремы о среднем для некоторой точки $\xi \in (0,{{t}_{0}})$ справедливо:
Пусть $g(s)$ на отрезке $\left[ {0,{{t}_{0}}} \right]$ возрастает и неотрицательна. Тогда $g({{t}_{0}}) = g(0) + \beta $, $\beta \geqslant 0$, и тогда имеем
Таким образом, при $g(s) \equiv \operatorname{const} \ne 0$ равенство ${{g}_{n}}({{t}_{0}}) = 0$ возможно только при $\alpha d_{n}^{2}{{t}_{0}} = 2\pi k$, $k \in N$. В силу (2.1) при больших $n$: ${{d}_{n}} \approx \frac{\pi }{{2l}}(2n - 1)$, тогда при $d = \frac{{\alpha {{t}_{0}}\pi }}{{8{{l}^{2}}}}$ = = $\frac{k}{{{{{(2n - 1)}}^{2}}}}$ нарушается единственность решения задачи 2.
Итак, установлен критерий единственности решения задачи 2.
Теорема 2. Если существует решение задачи 2, то оно единственно тогда и только тогда, когда ${{g}_{n}}({{t}_{0}}) \ne 0$ для всех $n \in N$.
Лемма 1. Пусть $g(t)$ – возрастающая и положительная функция на отрезке $\left[ {0,T} \right]$ и $d$ – алгебраическое число степени $s \geqslant 2$. Тогда при больших $n$ справедливы оценки
Доказательство проводится аналогично работе [12].
Лемма 2. Пусть $d$ – алгебраическое число степени $s \geqslant 2$. Тогда при любом $t \in \left[ {0,T} \right]$ и больших $n$ справедливы оценки
Справедливость этих оценок следует из формул (3.4) и (3.3) на основании леммы 1.
Формально почленным дифференцированием первого ряда из (3.5) составим ряды
В силу леммы 2 эти ряды при любых $(x,t) \in \bar {D}$ мажорируются числовым рядом
(3.6)
${{C}_{5}}\sum\limits_{n = 1}^\infty {{{n}^{{8 + \varepsilon }}}\left( {\left| {{{\varphi }_{n}}} \right| + \left| {{{\tau }_{n}}} \right| + {{n}^{{ - 2}}}\left| {{{\nu }_{n}}} \right|} \right)} ,$(3.7)
$\begin{gathered} \varphi (x) \in {{C}^{{10}}}\left[ {0,l} \right] \\ \varphi (0) = \varphi {\kern 1pt} '(0) = \varphi {\kern 1pt} '{\kern 1pt} '(l) = \varphi {\kern 1pt} '{\kern 1pt} '{\kern 1pt} '(l) = {{\varphi }^{{(4)}}}(0) = {{\varphi }^{{(5)}}}(0) = {{\varphi }^{{(6)}}}(l) = {{\varphi }^{{(7)}}}(l) = {{\varphi }^{{(8)}}}(0) = {{\varphi }^{{(9)}}}(0) = 0 \\ \end{gathered} $При выполнении этих условий ряд (3.6) мажорируется сходящимся рядом
Таким образом, доказана следующая
Теорема 3. Пусть функции $\varphi (x)$, $\tau (x)$ $\nu (x)\;$ удовлетворяют условиям (3.7), кроме этого, непрерывная функция $g(t)$ и число $d$ удовлетворяют условиям леммы 1. Тогда существует единственное решение задачи 2 и оно определяется рядами (3.5).
4. Исследование обратной задачи 3. Удовлетворим функцию (2.2) граничному условию (1.6):
С учетом (3.1) из последнего равенства получаем интегральное уравнение Вольтерра первого рода относительно неизвестной функции $g(t)$:
где(4.2)
$K(t,s) = \sum\limits_{n = 1}^\infty {\frac{{{{f}_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}(t - s){{Y}_{n}}\left( {{{x}_{0}}} \right)} ;\quad 0 \leqslant s \leqslant t \leqslant T$(4.3)
$\tilde {h}(t) = h(t) - \sum\limits_{n = 1}^\infty {\left( {{{\tau }_{n}}\cos \alpha d_{n}^{2}t + \frac{{{{\nu }_{n}}}}{{\alpha d_{n}^{2}}}\sin \alpha d_{n}^{2}t} \right)} {{Y}_{n}}\left( {{{x}_{0}}} \right);\quad 0 \leqslant t \leqslant T$Формально почленным дифференцированием ряда (4.2) составим ряды
(4.4)
$K_{t}^{'}(t,s) = \sum\limits_{n = 1}^\infty {{{f}_{n}}\cos \alpha d_{n}^{2}(t - s){{Y}_{n}}\left( {{{x}_{0}}} \right)} ,\quad K_{{tt}}^{{''}}(t,s) = - \alpha \sum\limits_{n = 1}^\infty {d_{n}^{2}{{f}_{n}}\sin \alpha d_{n}^{2}(t - s){{Y}_{n}}\left( {{{x}_{0}}} \right)} $Аналогично из (4.3) получим
(4.5)
$\begin{gathered} \tilde {h}{\kern 1pt} '(t) = h{\kern 1pt} '(t) + \sum\limits_{n = 1}^\infty {\left( {\alpha d_{n}^{2}{{\tau }_{n}}\sin \alpha d_{n}^{2}t - {{\nu }_{n}}\cos \alpha d_{n}^{2}t} \right){{Y}_{n}}({{x}_{0}})} \\ \tilde {h}{\kern 1pt} '{\kern 1pt} '(t) = h{\kern 1pt} '{\kern 1pt} '(t) + \sum\limits_{n = 1}^\infty {\left( {{{{\left( {\alpha d_{n}^{2}} \right)}}^{2}}{{\tau }_{n}}\cos \alpha d_{n}^{2}t + \alpha d_{n}^{2}{{\nu }_{n}}\sin \alpha d_{n}^{2}t} \right)} {{Y}_{n}}\left( {{{x}_{0}}} \right) \\ \end{gathered} $На основании леммы 2 при любых $(x,t) \in \overline D $ ряды (4.2) и (4.4) мажорируются числовым рядом ${{C}_{7}}\sum\nolimits_{n = 1}^\infty {{{n}^{2}}\left| {{{f}_{n}}} \right|} $, а ряды (4.3) и (4.5) – рядом ${{C}_{7}}\sum\nolimits_{n = 1}^\infty {{{n}^{4}}\left( {\left| {{{\tau }_{n}}} \right| + {{n}^{{ - 2}}}\left| {{{\nu }_{n}}} \right|} \right)} $.
Для сходимости указанных числовых рядов достаточно потребовать выполнения условий теоремы 3. Тогда эти ряды мажорируются сходящимся рядом ${{C}_{8}}\sum\nolimits_{n = 1}^\infty {\frac{1}{{{{n}^{2}}}}} $.
Дифференцируя интегральное уравнение (4.1) по t, будем иметь
Учитывая, что $K(t,t) = 0$, при повторном дифференцировании получим
Поскольку из (4.2) следует, что
(4.6)
$g(t)f\left( {{{x}_{0}}} \right) + \int\limits_0^t {g(s)K_{{tt}}^{{''}}(t,s)ds} = \tilde {h}{\kern 1pt} '{\kern 1pt} '(t);\quad 0 \leqslant t \leqslant T$Тогда, если $f\left( {{{x}_{0}}} \right) \ne 0$, то уравнение (4.6) является интегральным уравнением Вольтерра второго рода с непрерывным ядром и непрерывной правой частью при условии $h(t) \in {{C}^{2}}\left[ {0,T} \right]$. Следовательно, это уравнение имеет единственное решение $g(t)$ в классе $C\left[ {0,T} \right]$.
Теорема 4. Пусть функции $\tau (x)$, $\nu (x)$ удовлетворяют условиям теоремы 1, $f(x) \in {{C}^{4}}\left[ {0,l} \right]$, $f(0) = f{\kern 1pt} '(0) = f{\kern 1pt} '{\kern 1pt} '(l) = f{\kern 1pt} '{\kern 1pt} '{\kern 1pt} '(l)$ = 0, $h(t) \in {{C}^{2}}\left[ {0,T} \right]$ и $h(0) = \tau \left( {{{x}_{0}}} \right)$, h'(0) = = $\nu \left( {{{x}_{0}}} \right)$. Тогда, если $f\left( {{{x}_{0}}} \right) \ne 0$, то задача 3 имеет единственное решение, которое определяется формулой (2.2), а функция $g(t)$ находится из интегрального уравнения (4.6).
Выясним, насколько существенно условие $f\left( {{{x}_{0}}} \right) \ne 0$ в теореме 4. Пусть для некоторых $n = m$ выполняется равенство ${{Y}_{m}}({{x}_{0}}) = 0$. Тогда для функции $f(x) = {{Y}_{m}}(x)$ при любой функции $g(t) \in C\left[ {0,\,T} \right]$ существует ненулевое решение задачи 3 (где τ(x) = = $\nu (x) = h(t)$ = 0)
Список литературы
Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1966. 724 с.
Рэлей Л. Теория звука. Т. 1. М.: Гостехиздат, 1955. 503 с.
Тимошенко С.П. Колебания в инженерном деле. М.: Физматлит, 1967. 444 с.
Филиппов А.П. Колебания деформируемых систем М.: Машиностроение, 1970. 736 с.
Доннел Л.Г. Балки, пластины и оболочки. М.: Наука, 1982. 568 с.
Крылов А.Н. Вибрация судов. М.: Гостехиздат, 2012. 447 с.
Сабитов К.Б., Фадеева О.В. Начально-граничная задача для уравнения вынужденных колебаний консольной балки // Вестн. Самар. Гос. тех. ун-та. Сер.: Физ.-мат. науки. 2021. Т. 25. № 1. С. 51–66.
Романов В.Г. Обратные задачи уравнений математической физики. М.: Наука, 1984. 264 с.
Денисов А.М. Введение в теорию обратных задач. М.: МГУ, 1994. 208 с.
Prilepko A.I., Orlovsky D.G., Vasin I.A. Method for Solving Inverse Problems in Mathematical Physics. New York; Basel: 1999. 709 c.
Кабанихин С.И. Обратные и некорректные задачи. Новосибирск, Сиб. научн. изд-во, 2009. 457 с.
Сабитов К.Б. Обратные задачи для уравнения колебаний балки по определению правой части и начальных условий // Дифф. уравн. 2020. Т. 56. № 6. С. 773–785.
Дополнительные материалы отсутствуют.
Инструменты
Прикладная математика и механика


