Научная статья на тему 'К решению обратной статической задачи.'

К решению обратной статической задачи. Текст научной статьи по специальности «Математика»

CC BY
63
17
Читать
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук

Аннотация научной статьи по математике, автор научной работы — Ахтямов А. М.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Предварительный просмотр
i Надоели баннеры? Вы всегда можете отключить рекламу.

On solving of an inverse static problem.

The problem of finding of conditions for fastening the beam endpoints, which are inaccessible to direct observation, from bending of a beam is considered. The problem in question belongs to the class of inverse problems. A theorem on the uniqueness of the solution of this problem is proved and a method for establishing the unknown conditions for fastening the beam endpoints is indicated. The method calculation of unknown boundary conditions is explained with the help of examples.

Текст научной работы на тему «К решению обратной статической задачи.»

К РЕШЕНИЮ ОБРАТНОЙ статической

ЗАДАЧИ

Ахтямов А.М. (AkhtyamovAM@mail.ru) Башкирский государственный университет

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

Суть задачи состоит в следующем: Пусть закрепление стержня на обоих концах неизвестно (концы стержня недоступны для визуального осмотра). Известна длина стержня и интенсивность распределенной нагрузки q. Известны также прогибы стержня в нескольких его точках. Требуется определить вид и параметры закрепления стержня на обоих его концах.

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

2. Математическая постановка обратной задачи в случае стержня постоянного сечения. Дифференциальное уравнение упругой линии стержня постоянного сечения имеет вид

¿4 У{г) = (1)

¿г4 Е ,1Х' ()

где q = qy — интенсивность распределенной нагрузки, Е — модуль упругости Юнга, а ,]Х — момент инерции поперечного сечения относительно оси изгиба.

Краевые условия

ЩУ) = 52 а%з {Ьз У)х=о

3 = 1

О,

1' 2'

и%{У) = X) Ыз {Ьз у)=

3 = 1

О,

для уравнения (??) относятся к значениям: прогиба

Ь1 У

угла поворота

изгибающего момента

и перерезывающей силы

Ь2 У

¿У{г) ¿г

Ьз У = Мх{г ) = ЕЗх

Ь4 У = Яу (г) = Е Зх

¿2У{г) ¿г 2

¿3У{г) ¿г3

3' 4

4

4

Обозначим матрицу, составленную из коэффициентов а^ форм и1(у) и и2(у), через А, а матрицу, составленную из коэффициентов Ь^ форм и3(у) и и4(у), через В:

А :--

ац а12 а1з а 14 а21 а22 а2з а24

В

Ь11 Ь12 Ь1з Ь14 Ь21 Ь22 Ь23 Ь24

В работе [?] приведены различные случаи закрепления стержня. Ниже они перечислены для закрепления стержня на левом конце. Выписаны также их соответствующие матрицы А :

1) жесткое защемление:

2) свободное опирание

3) свободный конец:

4) плавающая заделка:

5) пять различных видов упругого закрепления:

1 0 0 0

0 1 0 0

1 0 0 0

0 0 1 0

0 0 0 1

0 0 1 0

0 0 0 1

0 1 0 0

С1 0 0 1 0 10 0

10 0 0 0 -С2 10

С1 0 0 1 0 0 10

0 0 0 1 0 -С2 10

С1 0 0 1 0 -С2 10

на правом.

Таким образом, имеем 9 видов краевых условий на левом конце и 9 Количество их возможных комбинаций 92 = 81.

Поэтому обратная задача состоит в следующем: правильно распознать одну из 81 комбинации краевых условий и в случае упругого закрепления найти соответствующие коэффициенты с1 и с2.

Заметим, что речь не может идти о восстановлении всех коэффициентов а^ и Ь^, поскольку, например, краевые условия

у(0) = 0, у'(0) = 0

и

4 у(0) = 0, 5 у'(0) = 0

эквивалентны, а их соответствующие коэффициенты а^ различны.

Для однозначного решения обратной задачи общее представление краевых условий в общем виде (??)-(??) является неэффективным. Более удобным оказывается уточненное представление.

Заметим, что во всех девяти случаях закрепления стержня на левом конце матрица А имеет вид

а1 0 0 а4

0 а2 аз 0

А

где |а 11 + |а41 = 0 и |а2| + |аз| = 0.

Случаи различных видов закрепления стержня на правом конце такие же как и на левом конце. Единственное различие состоит в том, что знаки перед коэффициентами

сх и с2 в матрицах А и В имеют противоположные знаки. Поэтому матрица В может быть записана в виде

" Ъх 0 0 &4 0 Ъ2 Ъз 0

B

где |bi| + \h\ = 0 и \Ь2\ + \Ьз\ = 0.

Эти представления позволяют записать краевые условия (??)-(??) в более простой уточненной форме:

Ui{y) = {ai Li y + a4 L4 y)z=0 = 0, (4)

U2(y) = (a2 L2 У + a3 L3 y)z=0 = 0, (5)

Us(y) = (bi Li y + b4 L4 y)z=i = 0, (6)

U4(y) = (b2 L2 y + Ьз L3 y)z=i = 0, (7)

где \ai\ + \a4\ = 0, \a2\ + \аз\ = 0, \bi\ + ^ = 0 и \b2\ + \Ьз\ = 0.

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

Коэффициент с большим индексом в паре коэффициентов любого из уравнений (??)-(??) будем называть вторым, а коэффициент с меньшим индексом — первым. Так, в паре коэффициентов ai, a4 из уравнения (??) первым коэффициентом является ai, а вторым — a4. В уравнении (??) первым коэффициентом является b2, а вторым —

Ьз.

Каноническими краевыми условиями будем называть такие краевые условия (??)-(??), которые удовлетворяют следующим двум требованиям: 1) второй коэффициент каждого из краевых условий равен единице или нулю; 2) если второй коэффициент равен нулю, то первый коэффициент равен 1.

Канонические краевые условия представляют собой некоторую нормировку краевых условий (??)-(??). Все эквивалентные краевые условия имеют единственное каноническое представление.

Коэффициенты канонических краевых условий будем называть каноническими коэффициентами. Заметим, что все перечисленные выше реальные физические краевые условия являются каноническими.

Обратная статическая задача сводится к распознанию именно таких канонических краевых условий и канонических коэффициентов.

В этих новых терминах поставленная обратная задача формулируется следующим образом: требуется восстановить канонические коэффициенты ai и bi (i = 1, 2, 3, 4) краевых условий (??)-(??) задачи (??), (??)-(??), если известны значения прогибов y(zk) в четырех точках zk (k =1, 2, 3, 4).

3. Метод нахождения канонических краевых условий. Далее для простоты изложения будем считать, что Jx и qy являются постоянными.

Решение обратной задачи будет проведено в два этапа. На первом этапе по прогибам стержня в четырех точках будет найдена функция прогиба. На втором — по функции прогиба будут восстановлены неизвестные канонические краевые условия.

1. Найдем функцию прогиба y(z).

Функция прогиба стержня у(г) должна удовлетворять уравнению

у(г) _ V (8)

Общим решением неоднородного уравнения (??) является функция

2 з 4

у(г) = С\ + С2 г + Сз у + С4 3 + ^ • 4, (9)

где Сг (г =1, 2, 3, 4) — произвольные константы.

Если известны прогибы у(гк) в четырех различных точках

0 < г1 < г2 < гз < г4 < 1,

то функция прогиба у (г) однозначно определяется из (??). Действительно, из (??) следует

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2 з 4

С1 + С2 ^ + Сз | + С4 | = у(гк) - Еу • 4, к = 1, 2, 3, 4. (10)

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

2. Найдем неизвестные краевые условия по функции прогиба у (г). Подставим найденную функцию у (г) в краевые условия (??)-(??). Получим четыре уравнения. Каждое уравнение представляет собой уравнение от двух неизвестных коэффициентов. Из (??)-(??) следует, что они должны быть связаны соотношениями:

а1/а4 = - (¿4 у/Ь1 у)х=0 , (11)

а2/аз = - (Ьз у/Ь2 у)х=0 , (12)

Ьх/Ь4 = - (¿4 у/Ь1 у)2=1, (13)

Ь2/Ьз = - (Ьз у/¿2 у)2=1. (14)

Приведем алгоритм определения канонических коэффициентов краевых условий по системе (??)-(??).

Поскольку на концах стержня прогиб ¿1 у и перерезывающие силы ¿4 у не обращаются в нуль одновременно, а также не обращаются в нуль одновременно угол поворота Ь2 у и изгибающий момент Ьз у, то оба коэффициента аг или Ьг из каждого уравнения не могут быть одновременно равняться нулю. То есть случай, когда правая часть одного из равенств равна 0/0 невозможен.

Однако, случай когда знаменатель в правой части какого-либо из равенств (??)-(??) обращается в нуль не исключается. В этом случае считаем, что коэффицент в знаменателе левой части соответствующего равенства (второй коэффициент) равен нулю, а коэффициент в числителе (первый коэффициент) следует принять за единицу.

Если правая часть одного из равенств (??)-(??) равна нулю, то считаем, что коэффициент, стоящий в числителе левой части этого равенства (первый коэффициент), равен нулю, а коэффициент в знаменателе (второй коэффициент) равен единице.

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

Из приведенных выше рассуждений следует, что канонические краевые условия и канонические коэффициенты находятся однозначно.

4. Пример решения обратной задачи по точным значениям прогибов. В

рассматриваемом примере будем полагать, что Ел

1.

Пример 1. Пусть прогибы в точках ¿1 = 1/5, ¿2 = 2/5,

равны

¿3

3/5, ¿4 = 4/5

у(1/5) = 157431770000, у(2/5) = 7771590000, у(3/5) = 9761590000, у(4/5) = 314331770000.

:15)

:1б)

(Значения прогибов выбраны соотвествующими упругому защемлению на левом конце 100 у(0)+у'''(0) = 0, у'(0) = 0 и упругой опоре на правом конце —20 у(1)+у'"(1) = 0,

У"(1) = 0.)

Решим обратную задачу для уравнения у(4) = 1 по значениям прогибов (??), предполагая, что краевые условия неизвестны.

Подставим значения (??) и соответствующие прогибы (??) в систему (??). Получим решение:

С = 31/4720, С2 = 0, Сз = 37/472, С4 = —155/1416. Подставив найденные константы (??) в (??), получим функцию прогиба

¿1 У = У^)

31

+

37

¿2

155

1

1416 ¿ +24 ¿^

4720 472

Поскольку ¿2 У = ^(¿) = 37¿/236 — 155 ¿2/472 + ¿3/6,

Ь3 у = ^(¿) = 37/236 — 155 ¿/236 + ¿2/2, Ь4 у = у'''^) = ¿ — 155/236, то система (??)-(??) имеет вид:

:17)

:18)

ах

а4

100,

а2 аз

Отсюда

ах = 100,

а2

,

1,

Ъ1

Ъ4

20,

Ъ2 Ъз

0.

а3 = 0, а4 = 1,

Ъх = —20, Ъ2 = 0, Ъ3 = 1,

Ъ

4=

1.

Искомые краевые условия соответствуют упругому защемлению на левом конце и упругой опоре на правом конце:

100 у(0)+ у'''(0) = 0, у' (0) = 0,

— 20 у(1)+ у'''(1) = 0, у'' (1) = 0.

Полученные краевые условия в точности совпадают с теми, по которым были выбраны соответствующие прогибы.

5. Пример решения обратной задачи по приближенным значениям прогибов. Как и в предыдущем примере считаем, что = 1-

Пример 2. Пусть найденные прогибы в точках

z1 = 1/5, ¿2 = 2/5, ¿3 = 3/5, ¿4 = 4/5 (19)

равны:

у(1/5) = 0.00889435, у(2/5) = 0.01317119, ()

у(3/5) = 0.01654407, у(4/5) = 0.01775876. ()

(Приближенные значения совпадают с прогибами из примера 1 с точностью до 10_8.)

Решим обратную задачу для уравнения у(4) = 1 по этим приближенным значениям прогибов (??).

Подставим значения (??) и соответствующие прогибы (??) в систему (??), найдем неизвестные константы С1, С2, С3 и С4. Подставив найденные константы в (??), получим функцию прогиба

Ь1 у = у^) = 0.006567780000 + 0.125 10"6 z

+ .07838958333 ¿2 - 0.1094631250 ¿3 + ¿4/24.

Поскольку

Ь2 у = 0.125 • 10"6 + 0.1567791667 z - 0.3283893750 ¿2 + ¿3/6, Ь3 у = у'\г) = 0.1567791667 - 0.6567787500 z + г2/2, Ь4 у = у'"^) = z - 0.6567787500,

то система (??)-(??) имеет вид:

— = 100.0001142, — = -0.1254233334 • 107,

«4 «3

Ь Ь

-1 = -20.00003788, ^ = 0.8429392740 • 10"4.

Ь4 Ь3

Поэтому

а1 = 100.0001142, а2 = -0.1254233334 • 107,

«3 = 1, «4 = 1,

Ь1 = -20.00003788, Ь2 = 0.8429392740 • 10"4,

Ь3 = 1, Ь4 =1. Искомые краевые условия имеют вид:

100.0001142 • у(0)+ у'''(0) = 0,

-0.1254233334 • 107 • у'(0) + у'''(0) = 0,

-20.00003788 • у(1) + у'''(1) = 0, 0.8429392740 • 10"4 • у'(1) + у''(1) = 0.

(21)

Получившиеся канонические коэффициенты совпадают каноническими коэффициентами из примера 1 с точностью до 10_3.

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

6. Заключение. Рассмотренная задача может быть отнесена к задачам технической диагностики (см., например, [?]). Однако, ранее в такой постановке она не ставилась. Наиболее близки ей по постановке задачи акустической диагностики закреплений (см. [?, ?]. Найденный метод решения можно распространить для решения обратных статических задач для более сложных случаев прогиба стержня, а также для случаев прогибов пластин и оболочек. Следует лишь заменить соответствующие дифференциальные выражения согласно подобающих этим в случаях формулам (см., например, [?]-[?]).

Автор выражает признательность М. А. Ильгамову за постановку задачи и полезные обсуждения. Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и ФЦП "Интеграция".

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

[1] Прочность, устойчивость, колебания. Справочник в 3-х т. / Под ред. И. А. Бир-гера и Я. Г. Пановко. М.: Машиностроение, 1968. Т. 1. 831 с.

[2] Биргер И.А. Техническая диагностика. М.: Машиностроение, 1978. 239 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

[3] Ахтямов А. М. О единственности восстановления краевых условий задачи Штурма-Лиувилля по ее спектру // Математическое моделирование. T. 12. № 3, 2000. С. 6.

[4] Ахатов И. Ш., Ахтямов А. М. Определение вида закрепления стержня по собственным частотам его изгибных колебаний // Прикладная математика и механика. 2001. Т. 65. Вып. 2. C. 290-298.

[5] Тимошенко С. П. Пластины и оболочки. М., Л.: ОГИЗ. 1948. 460 с.

[6] Ильгамов M. A. Введение в нелинейную гидроупругость. — М.: Наука, 1991. 200 с.

[7] Ilgamov M. A. Static Problems of Hydroelasticity. — Moscow: Nauka. Fizmatlit, 1998. 208 p.

i Надоели баннеры? Вы всегда можете отключить рекламу.