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

ОБ ИДЕНТИФИКАЦИИ ОБЩИХ ЗАКРЕПЛЕНИЙ ПРЯМОУГОЛЬНОЙ ПЛАСТИНЫ

А. М. Ахтямов ab*, Дж. А. Пардаев c**

a Башкирский государственный университет
Уфа, Россия

b Институт механики им. Р.Р. Мавлютова УФИЦ РАН
Уфа, Россия

c Джизакский государственный педагогический институт
Джизак, Узбекистан

* E-mail: AkhtyamovAM@mail.ru
** E-mail: pardayev.jasurhon@mail.ru

Поступила в редакцию 04.07.2019
После доработки 15.07.2019
Принята к публикации 27.07.2019

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

Аннотация

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

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

1. Введение. Первые систематические исследования по идентификации краевых условий начались в 90-х годах 20 века в работах З.Б. Оганисяна (см., например, [13]). З.Б. Оганесяном исследовались несколько задач идентификации условий закрепления распределенных механических систем: задача идентификации краевых условий круговой пластины [1], задача идентификации краевых условий прямоугольной пластины [2], задача идентификации краевых условий на обоих концах стержня [3]. В работах [46] изучались аналогичные задачи. В них восстанавливались канонические краевые условия. В [7, 8] изучалась идентификация краевых условий, в которых неизвестны все их коэффиценты. Подобные задачи впервые начали изучаться в этих работах и сводятся к идентификации (с точностью до линейных преобразований строк) матрицы из коэффициентов краевых условий по ее минорам. Настоящая работа посвящена идентификации общих видов закреплений прямоугольной пластины на двух противоположных сторонах. Ранее было показано, что 8 неизвестных коэффициентов упругих закреплений для одного из видов канонических краевых условий восстанавливаются с точностью до перестановки закреплений на противоположных концах однозначно по девяти собственным частотам. В работе [8] было показано, что различные виды и параметры общих закреплений могут быть восстановлены по девяти собственным частотам. В частности, в работе было показано, что краевые условия

$\Phi (0) = 0,\quad \Phi {\text{'}}(0) - {\alpha }\Phi {\text{''}}(0) = 0,\quad \Phi (a) = 0,\quad \Phi {\text{'}}(a) - {\beta }\Phi {\text{''}}(a) = 0$
восстанавливаются по двум собственным частотам.

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

2. Постановка задачи. Если толщина $h$ однородной пластины постоянна (цилиндрическая жесткость D = const), то уравнение свободных колебаний прямоугольной пластины имеет следующий вид (см., например, [9] или [10]):

$D{\Delta \Delta }w + {\rho }h\frac{{{{\partial }^{2}}w}}{{\partial {{t}^{2}}}} = 0$

Подстановка $w = {\varphi }({{x}_{1}},{{x}_{2}})\cos ({\omega }t - {\chi })$ приводит к уравнению

${\Delta \Delta \varphi } + {\rho }h{\omega \varphi } = 0$

Если на двух противоположных сторонах пластины реализуются условия свободного опирания

${\varphi } = \frac{{{{\partial }^{2}}{\varphi }}}{{\partial {{t}^{2}}}} = 0\quad {\text{при}}\quad {{x}_{2}} = 0,b,$
то формы колебаний описываются следующими функциями:
${\varphi }({{x}_{1}},{{x}_{2}}) = {{\Phi }_{m}}({{x}_{1}})\sin \frac{{m{\pi }{{x}_{2}}}}{b},\quad m = 1,2,\; \ldots $
где ${{\Phi }_{m}}({{x}_{1}})$ удовлетворяют уравнению
(2.1)
$\Phi _{m}^{{IV}}({{x}_{1}}) - \frac{{2{{m}^{2}}{{{\pi }}^{2}}}}{b}\Phi _{m}^{{''}}({{x}_{1}}) + \left( {\frac{{{{m}^{4}}{{{\pi }}^{4}}}}{{{{b}^{4}}}} - {{{\gamma }}_{m}}} \right){{\Phi }_{m}}({{x}_{1}}) = 0$
с общими видами закреплений
$\sum\limits_{i = 1}^4 {{a}_{{ij}}}{{L}_{i}}{{\Phi }_{m}}({{x}_{1}}),\quad {\text{при}}\;{{x}_{1}} = 0\quad {\text{и}}\quad \sum\limits_{i = 5}^8 {{b}_{{ij}}}{{L}_{i}}{{\Phi }_{m}}({{x}_{1}}),\quad {\text{при}}\quad {{x}_{1}} = a$
где ${{{\gamma }}_{m}} = \frac{{{\rho }h{\omega }{{m}^{2}}}}{D}$,

${{L}_{1}}{{\Phi }_{m}} = {{\left( {\frac{{m{\pi }}}{b}} \right)}^{4}}{{\left[ {{{\Phi }_{m}}({{x}_{1}})} \right]}_{{{{x}_{1}} = 0}}},\quad {{L}_{2}}{{\Phi }_{m}} = {{\left( {\frac{{m{\pi }}}{b}} \right)}^{2}}{{[\Phi _{m}^{'}({{x}_{1}})]}_{{{{x}_{1}} = 0}}}$
${{L}_{3}}{{{\Phi }}_{m}} = {{\left[ {\Phi _{m}^{{''}}({{x}_{1}}) - {\nu }{{{\left( {\frac{{m{\pi }}}{b}} \right)}}^{2}}{{{\Phi }}_{m}}({{x}_{1}})} \right]}_{{{{x}_{1}} = 0}}}$
${{L}_{4}}{{\Phi }_{m}} = {{\left[ {\Phi _{m}^{{'''}}({{x}_{1}}) - (2 - {\nu }){{{\left( {\frac{{m{\pi }}}{b}} \right)}}^{2}}\Phi _{m}^{'}({{x}_{1}})} \right]}_{{{{x}_{1}} = 0}}}$
${{L}_{5}}{{\Phi }_{m}} = {{\left( {\frac{{m{\pi }}}{b}} \right)}^{4}}{{\left[ {{{\Phi }_{m}}({{x}_{1}})} \right]}_{{{{x}_{1}} = a}}},\quad {{L}_{6}}{{\Phi }_{m}} = {{\left( {\frac{{m{\pi }}}{b}} \right)}^{2}}{{[\Phi _{m}^{'}({{x}_{1}})]}_{{{{x}_{1}} = a}}}$
${{L}_{7}}{{\Phi }_{m}} = {{\left[ {\Phi _{m}^{{''}}({{x}_{1}}) - {\nu }{{{\left( {\frac{{m{\pi }}}{b}} \right)}}^{2}}{{\Phi }_{m}}({{x}_{1}})} \right]}_{{{{x}_{1}} = a}}}$
${{L}_{8}}{{\Phi }_{m}} = {{\left[ {\Phi _{m}^{{'''}}({{x}_{1}}) - (2 - {\nu }){{{\left( {\frac{{m{\pi }}}{b}} \right)}}^{2}}\Phi _{m}^{'}({{x}_{1}})} \right]}_{{{{x}_{1}} = a}}}$

Фундаментальную систему решений уравнения (2.1) при m = 1 образуют следующие функции

$\begin{array}{*{20}{c}} {{{y}_{1}}(x) = \cos \left( {\sqrt {\sqrt {{{{\gamma }}_{{1k}}}} - \frac{{{{k}^{2}}{{{\pi }}^{2}}}}{{{{b}^{2}}}}} } \right),\quad {{y}_{2}}(x) = \sin \left( {\sqrt {\sqrt {{{{\gamma }}_{{1k}}}} - \frac{{{{k}^{2}}{{{\pi }}^{2}}}}{{{{b}^{2}}}}} } \right)} \\ {{{y}_{3}}(x) = {\text{ch}}\left( {\sqrt {\sqrt {{{{\gamma }}_{{1k}}}} - \frac{{{{k}^{2}}{{{\pi }}^{2}}}}{{{{b}^{2}}}}} } \right),\quad {{y}_{4}}(x) = {\text{sh}}\left( {\sqrt {\sqrt {{{{\gamma }}_{{1k}}}} - \frac{{{{k}^{2}}{{{\pi }}^{2}}}}{{{{b}^{2}}}}} } \right)} \end{array}$

Ранее в работе [5] было показано, что различные случаи закрепления прямоугольной пластины (упругие и неупругие) на двух противоположных краях прямоугольной пластины можно однозначно восстановить по девяти собственным значениям. Суть метода состояла в восстановлении матрицы

$C = \left\| {\begin{array}{*{20}{c}} {{{a}_{{11}}}}&{{{a}_{{12}}}}&{{{a}_{{13}}}}&{{{a}_{{14}}}}&0&0&0&0 \\ {{{a}_{{21}}}}&{{{a}_{{22}}}}&{{{a}_{{23}}}}&{{{a}_{{24}}}}&0&0&0&0 \\ 0&0&0&0&{{{b}_{{11}}}}&{{{b}_{{12}}}}&{{{b}_{{13}}}}&{{{b}_{{14}}}} \\ 0&0&0&0&{{{b}_{{21}}}}&{{{b}_{{22}}}}&{{{b}_{{23}}}}&{{{b}_{{24}}}} \end{array}} \right\|$
с точностью до линейных преобразований строк. Такая матрица восстанавливается с точностью до линейных преобразований строк по вектору, определяемому с точностью до постоянного множителя и составленному из 10 миноров четвертого порядка матрицы C. Поэтому в этом методе для восстановления матрицы требовалось 9 собственных частот.

В настоящей статье предлагается другой метод, основанный на представлении краевых условий в следующем каноническом виде:

(2.2)
$\begin{array}{*{20}{c}} {{{c}_{1}}{{L}_{1}}y - {{L}_{4}}y = 0,\quad {{c}_{2}}{{L}_{2}}y - {{L}_{3}}y = 0} \\ {{{c}_{3}}{{L}_{5}}y + {{L}_{8}}y = 0,\quad {{c}_{4}}{{L}_{6}}y + {{L}_{7}}y = 0} \end{array}$

При этом предполагается, что cj ($j = 1,2,3,4$) могут обращаться в нуль или бесконечность. Так, ${{c}_{j}} = \infty $ ($j = 1,2,3,4$) означает, что реализуется закрепление (заделка)–(заделка).

3. Построение решения и примеры. Собственными значениями задачи являются корни характеристических определителей:

${\Delta }({{{\gamma }}_{1}}) = \left| {\begin{array}{*{20}{c}} {{{c}_{1}}{{L}_{1}}{{y}_{1}} - {{L}_{4}}{{y}_{1}}}&{{{c}_{1}}{{L}_{1}}{{y}_{2}} - {{L}_{4}}{{y}_{2}}}&{{{c}_{1}}{{L}_{1}}{{y}_{3}} - {{L}_{4}}{{y}_{3}}}&{{{c}_{1}}{{L}_{1}}{{y}_{4}} - {{L}_{4}}{{y}_{4}}} \\ {{{c}_{2}}{{L}_{2}}{{y}_{1}} - {{L}_{3}}{{y}_{1}}}&{{{c}_{2}}{{L}_{2}}{{y}_{2}} - {{L}_{3}}{{y}_{2}}}&{{{c}_{2}}{{L}_{2}}{{y}_{3}} - {{L}_{3}}{{y}_{3}}}&{{{c}_{2}}{{L}_{2}}{{y}_{4}} - {{L}_{3}}{{y}_{4}}} \\ {{{c}_{3}}{{L}_{5}}{{y}_{1}} + {{L}_{8}}{{y}_{1}}}&{{{c}_{3}}{{L}_{5}}{{y}_{2}} + {{L}_{8}}{{y}_{2}}}&{{{c}_{3}}{{L}_{5}}{{y}_{3}} + {{L}_{8}}{{y}_{3}}}&{{{c}_{3}}{{L}_{5}}{{y}_{4}} + {{L}_{8}}{{y}_{4}}} \\ {{{c}_{4}}{{L}_{6}}{{y}_{1}} + {{L}_{7}}{{y}_{1}}}&{{{c}_{4}}{{L}_{6}}{{y}_{2}} + {{L}_{7}}{{y}_{2}}}&{{{c}_{4}}{{L}_{6}}{{y}_{3}} + {{L}_{7}}{{y}_{3}}}&{{{c}_{4}}{{L}_{6}}{{y}_{4}} + {{L}_{7}}{{y}_{4}}} \end{array}} \right|$

Если расписать подробно этот определитель, получим достаточно громоздкое выражение, которое составляет пять страниц. Оно содержит как слагаемые без множителей cj ($j = 1,2,3,4$), так и слагаемые с множителями cj ($j = 1,2,3,4$), ${{c}_{1}}{{c}_{2}}$, ${{c}_{1}}{{c}_{3}}$, ${{c}_{1}}{{c}_{4}}$, ${{c}_{2}}{{c}_{3}}$, ${{c}_{2}}{{c}_{4}}$, ${{c}_{3}}{{c}_{4}}$, ${{c}_{1}}{{c}_{2}}{{c}_{3}}$, ${{c}_{1}}{{c}_{3}}{{c}_{4}}$, ${{c}_{2}}{{c}_{3}}{{c}_{4}}$, ${{c}_{1}}{{c}_{2}}{{c}_{3}}{{c}_{4}}$.

К сожалению, решение системы пяти уравнений $\Delta ({{\gamma }_{1}}_{m})$ = 0 (m = 1, 2, 3, 4.5) методом подстановки приводит к алгебраическому уравнению выше пятой степени. Поэтому эта система уравнений не может быть решена в радикалах (аналитически). Однако она может быть решена численно на конкретных примерах.

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

Контрпример. Пусть собственными значениями задачи (1)–(2) с ${\nu } = \frac{1}{3}$, m = 1, a = 1 и $b = \pi $ являются следующие значения ${{{\gamma }}_{{11}}} = 47.930$, ${{{\gamma }}_{{12}}} = 901.12$, ${{{\gamma }}_{{13}}} = 44329$, ${{{\gamma }}_{{14}}}$ = 183 458. Решив систему уравнений ${\Delta }({{{\gamma }}_{{1j}}})$, $j = 1,2,3,4$, получим 24 решения, большинство из которых содержит комплексные и отрицательные числа. Однако есть среди них и решения с положительными числами:

(3.1)
$\begin{array}{*{20}{c}} {\{ {{c}_{1}} = 7.3091,\,\,\,{{c}_{2}} = 1.9908,\,\,{{c}_{3}} = 0.22964,\,\,{{c}_{4}} = 4.0100\} } \\ {\{ {{c}_{1}} = 0.22965,\,\,{{c}_{2}} = 4.0100,\,\,{{c}_{3}} = 3.4819,\,\,{{c}_{4}} = 1.9911\} } \\ {\{ {{c}_{1}} = 3.0000,\,\,{{c}_{2}} = 4.0000,\,\,{{c}_{3}} = 1.0000,\,\,{{c}_{4}} = 2.0000\} } \\ {\{ {{c}_{1}} = 1.0000,\,\,{{c}_{2}} = 2.0000,\,\,{{c}_{3}} = 3.0000,\,\,{{c}_{4}} = 4.0000\} } \end{array}$

Таким образом, имеется четыре решения, удовлетворяющие физическому смыслу задачи.

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

Пример 1. Пусть собственными значениями задачи (1)–(2) с ${\nu } = \frac{1}{3}$, m = 1, a = 1 и $b = {\pi }$ являются следующие значения ${{{\gamma }}_{{11}}} = 47.930$, ${{{\gamma }}_{{12}}} = 901.12$, ${{{\gamma }}_{{13}}} = 44329$, ${{{\gamma }}_{{14}}}$ = = 183 458, ${{{\gamma }}_{{15}}} = 95860$. Решив систему уравнений ${\Delta }({{{\lambda }}_{j}})$, $j = 2,3,4,5$, получим еще 24 решения, большинство из которых содержит комплексные и отрицательные числа. Однако есть среди них и решения с положительными числами. Их всего два:

(3.2)
$\begin{array}{*{20}{c}} {\{ {{c}_{1}} = 1.0000,\,\,\,{{c}_{2}} = 2.0000,\,\,\,{{c}_{3}} = 3.0000,\,\,\,{{c}_{4}} = 4.0000\} } \\ {\{ {{c}_{1}} = 3.0000,\,\,\,{{c}_{2}} = 4.0000,\,\,\,{{c}_{3}} = 1.0000,\,\,\,{{c}_{4}} = 2.0000\} } \end{array}$

Пересечением решений (3.1) и (3.2) является единственное решение (3.2). Оно и является решением, отвечающим физическому смыслу задачи. Упругие закрепления с коэффициентами жесткстей 1, 2, 3, 4 восстанавливаются с точностью до перестановок закреплений на противоположных краях.

Пример 2. Пусть собственными значениями задачи (1)–(2) с ${\nu } = \frac{1}{3}$, m = 1, a = 1 и $b = {\pi }$ являются следующие значения ${{{\gamma }}_{{11}}} = 18.873 \times {{10}^{6}}$, ${{{\gamma }}_{{12}}} = 538.44$, ${{{\gamma }}_{{13}}} = 14874$, ${{{\gamma }}_{{14}}}$ = 89 756, ${{{\gamma }}_{{15}}} = 174744$. Решив систему уравнений ${\Delta }({{{\lambda }}_{j}})$, $j = 1,2,3,4$, получим еще 24 решения, большинство из которых содержит комплексные и небольшие отрицательные числа. Однако есть среди них и такое:

$\begin{array}{*{20}{c}} {\{ {{c}_{1}} = - 7.6028 \times {{{10}}^{{ - 23}}},\,\,\,{{c}_{2}} = - 6.8337 \times {{{10}}^{{ - 23}}}} \\ {{{c}_{3}} = - 4.2899 \times {{{10}}^{{25}}},\,\,{{c}_{4}} = 2.9238 \times {{{10}}^{{23}}}\} } \\ {\{ {{c}_{1}} = - 5.5150 \times {{{10}}^{{25}}},\,\,{{c}_{2}} = - 3.7874 \times {{{10}}^{{23}}}} \\ {{{c}_{3}} = - 5.8543 \times {{{10}}^{{ - 23}}},\,\,{{c}_{4}} = - 5.2671 \times {{{10}}^{{ - 23}}}\} } \end{array}$

Решив систему уравнений ${\Delta }({{{\lambda }}_{j}})$, $j = 2,3,4,5$, получим еще 24 решения, большинство из которых содержит комплексные и небольшие отрицательные числа. Однако есть среди них и такие:

$\begin{array}{*{20}{c}} {\{ {{c}_{1}} = 2.4623 \times {{{10}}^{6}},\,\,\,{{c}_{2}} = 0.11968,\,\,\,{{c}_{3}} = 2.5221,\,\,\,{{c}_{4}} = 2736.1\} } \\ {\{ {{c}_{1}} = 2.5221,\,\,\,{{c}_{2}} = 2736.2,\,\,\,{{c}_{3}} = 2.4623 \times {{{10}}^{6}},\,\,\,{{c}_{4}} = 0.11969\} } \\ {\{ {{c}_{1}} = - 3.0620 \times {{{10}}^{{ - 19}}},\,\,\,{{c}_{2}} = - 1.7473 \times {{{10}}^{{ - 20}}}} \\ {{{c}_{3}} = - 2.4929 \times {{{10}}^{{24}}},\,\,\,{{c}_{4}} = - 5.8345 \times {{{10}}^{{21}}}\} } \\ {\{ {{c}_{1}} = - 2.3505 \times {{{10}}^{{24}}},\,\,\,{{c}_{2}} = - 5.3911 \times {{{10}}^{{21}}},} \\ {{{c}_{3}} = - 3.3482 \times {{{10}}^{{ - 19}}},\,\,\,{{c}_{4}} = - 1.9070 \times {{{10}}^{{ - 20}}}\} } \end{array}$

Вычисления проводились с точностью до 30 значащих цифр. Поэтому с большой степенью точности можно считать, что пересечением этих решений являются следующие решения: ${{c}_{j}} = \infty $, $j = 1,2$, ${{c}_{j}} = 0$, $j = 3,4$ и ${{c}_{j}} = 0$, $j = 1,2$, ${{c}_{j}} = \infty $, $j = 3,4$. Следовательно, реализуется закрепление (заделка)–(свободный край) и (свободный край)–(заделка).

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

5. Благодарности. Работа выполнена с использованием средств государственного бюджета по госзаданию на 2019-2022 годы (№ 0246-2019-0088) и при поддержке грантов Российского фонда фундаментальных исследований (гранты № 18-51-06002-Aз_a, 18-01-00250-a).

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

  1. Гнуни В.Ц., Оганисян З.Б. Определение граничных условий круглой кольцевой пластинки по заданным частотам собственных колебаний // Известия НАН РА, серия “Механика”. 1991. Т. 44. № 5. С. 9–16.

  2. Оганисян З.Б. Об одной задаче восстановления граничных условий на краях пластинки при заданном спектре частот собственных поперечных колебаний // Ученые записки ЕГУ. 1991. № 1. С. 45–50.

  3. Оганисян З.Б. Об одной задаче восстановления граничных условий на концах стержня при заданном спектре частот собственных поперечных колебаний. “Вопросы оптимального управления, устойчивости и прочности механических систем” (научные труды конференции). Ереван. 1997. С. 159–162.

  4. Сергиенко И.В. Системный анализ многокомпонентных распределенных систем / И.В. Сергиенко, В.С. Дейнека. НАН Украины, Институт кибернетики им. В.М. Глушкова. Киев: Наукова думка, 2009. 639 с.

  5. Халилов С.А., Минтюк В.Б. Исследование устойчивости отсека крыла методом идентификации краевых условий на основе упрощенной модели // Авiцiйно-космiчна технiка и технологiя. 2003. Вып. 2. С. 6–10.

  6. Ахтямов А.М., Шагиев В.Р. Идентификация неупругих видов закреплений трубопроводов // Вестник Башкирского университета. 2016. Т. 21. № 1. С. 21–26.

  7. Akhtyamov A.M., Mouftakhov A.V. Identification of boundary conditions using natural frequencies // Inverse Probl. Sci. and Eng-ng. 2004. V. 12. № 4. P. 393–408.

  8. Ахтямов А.М., Муфтахов А.В., Тайхер М., Ямилова Л.С. Об одном методе определения по собственным частотам условий закрепления прямоугольной пластины // Известия РАН. МТТ. 2007. № 1. С. 100–113.

  9. Strutt W. (Lord-Rayleigh) The theory of Sound. V. 1. L.: Macmillan 1926. = Стрэтт Дж.В. (Лорд Рэлей). Теория звука. Т. 1. М.; Л.: Гостехиздат, 1940. 500 с.

  10. Вибрации в технике: Справочник. В 6-ти т. Т. 1. Колебания линейных систем / Под ред. В.В. Болотина. М.: Машиностроение, 1978. 352 с.

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