Об однозначности идентификации параметров упругого закрепления и сосредоточенного инерционного элемента

Автор: Ахтямова Айгуль Азаматовна, Ахтямов Азамат Мухтарович

Журнал: Вычислительная механика сплошных сред @journal-icmm

Статья в выпуске: 1 т.6, 2013 года.

Бесплатный доступ

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

Собственные значения, обратная задача, собственные частоты, балка, сосредоточенный инерционный элемент

Короткий адрес: https://sciup.org/14320658

IDR: 14320658

Текст научной статьи Об однозначности идентификации параметров упругого закрепления и сосредоточенного инерционного элемента

Вопросы вычисления собственных частот распределеных механических систем (в частности, упруго закрепленных стержней и пластин) достаточно хорошо изучены [1–5]. Обратные задачи для таких систем стали решаться сравнительно недавно [6–20]. Так, в работах [6–8] обсуждались задачи идентификации пружинно-массовых систем с конечным числом степеней свободы по собственным частотам их колебаний. В обратных спектральных задачах [6, 9, 10] требуется восстановить коэффициенты дифференциального уравнения и краевых условий. В качестве данных для восстановления краевых условий используется не один спектр или его часть (как в предлагаемой вниманию статье), а несколько спектров или же другие дополнительные спектральные данные (например, спектральная функция, функция Вейля или так называемые весовые числа). К тому же основной целью этих работ является восстановление коэффициентов в уравнении, а не в краевых условиях.

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

Распределенные механические системы (стержни и балки) служат объектами исследования в [16–20], где диагностируются трещины и полости в них по собственным частотам продольных и изгибных колебаний. Но в этих работах определяются не краевые условия, а условия сопряжения.

В [11, 15] и других публикациях (см. [6] и библиографию к ней) решаются задачи идентификации сосредоточенных масс на стержне и балке по собственным частотам продольных и изгибных колебаний. В [12–14] приводятся задачи идентификации видов и параметров закреплений стержней балок. В монографии [12] обобщаются результаты идентификации краевых условий в случае, когда отыскиваются все их коэффиценты. Подобные задачи впервые начали изучаться одним из авторов данной статьи и сводятся к идентификации (с точностью до линейных преобразованиий строк) матрицы из коэффициентов краевых условий по ее минорам. Из упомянутых выше публикаций работы [13–15] наиболее близки к обсуждаемой статье, так как в ней также полагается, что коэффициенты при старших производных y ′′′ (1), y ′′ (1) уже известны. В работах [13, 14] по собственным частотам восстанавливаются

вид и параметры закрепления концов стержня и трубопровода. В [15] по собственным частотам идентифицируются параметры твердого тела (масса, момент инерции, статический момент инерции), прикрепленного к одному из концов балки Тимошенко. Но, в отличие от этих работ, в данной статье идентифицируются как параметры упругого закрепления конца стержня (два относительных коэффициента жесткости пружинок), так и параметры сосредоточенного инерционного элемента (масса и момент инерции), прикрепленного к концу стержня. Ранее подобная задача не рассматривалась. Предложен метод введения дополнительных неизвестных величин, на основе которого доказывается теорема об однозначности восстановления параметров закрепления конца стержня и сосредоточенного инерционного элемента. При доказательстве используется теория некорректных задач [21-27].

Исследуемая проблема часто возникает в технических устройствах, поскольку деталями многих механизмов являются упруго закрепленные стержни с сосредоточенной массой на конце. Такими системами можно считать, в частности, краны и автомобили [8], автоматические записывающие устройства с трубками, из которых вытекает краска, детали некоторых механизмов виброзащиты [14]. Как правило, пружинки (рессоры), служащие составной частью систем, имеют определенную жесткость, влияющую на нагруженность конца. Со временем системы теряют свои массу и жесткость ввиду износа или, наоборот, наращивают их в результате захватывания инородных предметов. Если конец стержня недоступен для визуального осмотра, а разборка механизма представляет собой дорогостоящую процедуру, то для сохранения надежной работы механизма возникает потребность в его ранней и неразрушающей диагностике (например, акустической), то есть возникает задача определения параметров упругого закрепления и нагруженности конца стержня по характеристикам звуковых колебаний, вызванных ударом по стержню. Поэтому становится важной проблема идентификации параметров упругого закрепления и нагруженности конца однородного стержня по собственным частотам его колебаний.

2.    Постановка обратной задачи

Рассмотрим однородную балку Эйлера-Бернулли длиной L , плотностью р и площадью поперечного сечения F , левый конец которой заделан; на правом конце сосредоточен груз, обладающий массой m и моментом инерции I . Груз упруго закреплен на пружинках с жесткостями c и c 2, препятствующими вертикальному смещению балки ( с , ) и повороту ( c 2). Требуется найти массу m , момент инерции I , коэффициенты жесткостей пружинок с , и c 2 по собственным частотам колебаний балки. Отклонения точек оси балки от начального состояния (при t = 0) при поперечных колебаниях однозначно определяются одной функцией двух переменных — осевой координаты X и времени t : U = U ( X , t ) .

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

EI 54 U ( X , t ) F 92 U ( X , t ) =0, д x 4 д t2

где U = U (X, t) — прогиб оси стержня, EI — изгибная жесткость, р — плотность материала, F — площадь поперечного сечения стержня.

При t = 0 должны выполняться начальные условия dU

U (X ,0) = f (X),    д (X ,0) = g (X), д t где f (X), g(X) — функции, определяющие начальное положение оси стержня.

Если левый конец заделан, а на правом конце имеется упруго закрепленный сосредоточенный инерционный элемент, то краевые условия записываются в следующем виде:

U = 0,

д1иг        д1 U

EI —г + c U = - m —7-, д X 3    1 д t 2

  • ^U- = 0 (при X = 0 ), д X

...ди ди . Uu   . Y

EI--7 + c---= -1----7 (при X = L ), дx 2   2 дx     дX д t2

Вводя обозначения x = X)L , и = U/L , запишем уравнение свободных изгибных колебаний однородного стержня и краевые условия следующим образом:

д4и(x, t) рFL4 д2и (x, t) Л Л ди ..

. 4   ■    ,--,     = 0, и = 0,   — = 0 (при x = 0), дx      EI д tд д3 и      з          Зд2 и         д2 и.    Т8 иТт     Т Т д3 и,

EI г + c.L и = - mL— ,   EI—77 + cL—U = - IL---т- (при x = 1).

д x 3    1               д t 2           д x 2    2 д x         1 д x д t 2

Тогда, при замене и ( x , t ) = у ( x )cos( to t ), поставленная выше задача сводится (см., например, [2]) к следующей спектральной задаче [3]:

у (4) = X 4 у ,     U i = у (0) = 0,      U 2 = у '(0) = 0;                             (1)

U з ( у ) = у "(1) + ( a i - a 2 Х 4 ) у (1) = 0,      U 4( у ) = у "(1) + ( a 3 - a 4 X 4 ) у (1) = 0.                 (2)

Здесь a 1 = ( c 1 L 3)/( EI ), a 2 = ( mL 3)/( EI ), a 3 = ( c 2 L )/( EI ), a 4 = ( 1 1 L )/( EI ), X 4 = р FL 4 ю 2/ ( EI ).

Таким образом, имеем краевую задачу (1), (2) со спектральным параметром X и неизвестными коэффициентами a i ( i = 1,2,3,4). Требуется по ее собственным значениям X k найти неизвестные коэффициенты a i ( i = 1,2,3,4), то есть поставленную задачу определения массы, момента инерции груза, сосредоточенного на одном из концов стержня, и коэффициентов жесткости удерживающей его связи, свели к обратной. Покажем, что эту обратную задачу можно переформулировать в терминах характеристического определителя.

Функции у 1 ( x , X ) = (cos X x + ch X x )/2, у 2( x , X ) = (sin X x + sh X x )/ ( 2 X ) , у 3( x , X ) = ( - cos X x + ch X x )/ ( 2 X 2 ) , у 4( x , X ) = ( - sin X x + sh X x )/ ( 2 X 3 ) являются линейно независимыми решениями уравнения

у (4) ( x , X ) = X 4 у ( x , X ),

удовлетворяющими условиям у (r "(0,X) = (0,при j * r (j,r = 1,2,3,4)                                 (4)

j ,      |1, при j = r ’      ’ ’ ’

(другими словами, решения у] (x, X) (j = 1,2,3,4) образуют фундаментальную систему Коши и выражаются через функции Крылова [3]).

Общее решение уравнения (3) представляется в следующем виде:

у ( x ) = у ( x , X ) = C 1 у , ( x , X ) + C 2 у 2 ( x , X ) + C 3 у 3 ( x , X ) + C 4 у 4 ( x , X ).

Для нахождения констант C 1 , C 2, C 3, C 4 используем краевые условия (1), (2):

U ( у ) = U ( C 1 у 1 + C 2 у 2 + C 3 у 3 + C 4 у 4 ) = C 1 U ( у 1 ) + C 2 U ( у 2 ) + C 3 U ( у 3 ) + C 4 U ( у 4 ) = 0 ( i = 1,4). (5)

Уравнение для определения собственных значений задачи (1), (2) следует из условия существования ненулевого решения системы (5). Ненулевое решение для C i существует тогда и только тогда, когда равняется нулю определитель системы:

U , ( у 1 ) U 1 ( у 2 ) U 1 ( у 3 ) U 1 ( у 4) U 2 ( у 1 ) U 2 ( у 2 ) U 2 ( у 3 ) U 2 ( у 4 ) U 3 ( у 1 ) U 3 ( у 2 ) U 3 ( у 3 ) U 3 ( у 4 ) U 4 ( у 1 ) U 4 ( у 2 ) U 4 ( у 3 ) U 4 ( у 4 )

Выражение (6) называется характеристическим определителем спектральной задачи (1), (2). Его нули совпадают с собственными значениями этой задачи [28]. Учитывая условия (4), из (6) получаем

10  0  0

0     10     0   uU з ( у з ) и з ( у 4)

U з ( У 1 )   U з ( у 2)   U з ( у з )   U з ( у 4 )     U 4 ( у з )   U 4 ( у 4 )

U 4 ( у 1 )   U 4 ( у 2 )   U 4 ( у з )   U 4 ( у 4 )

Отсюда, с учетом (1), (2), имеем

A ( X ) =

у з” (1) + q у з (1) - а 2 X 4 у з (1) у з" (1) + а з у з ' (1) а 4 X 4 у з ' (1)

у 7(1) + а 1 у 4 (1) а 2 X 4 у 4 (1) у 4 "(1) + а з у /(1) а 4 X 4 у 4 (1)

Следовательно,

А ( Х ) = - f , ( X ) + q f ( X ) + а 2 f , ( X ) + а з f , ( X ) + а 4 f , ( X ) + ( q а 4 + а 2 а з ) f , ( X ) + q а з f , ( X ) + а 2 а 4 f , ( X ), (7)

где f ) ( X )= ( 1 + cos X ch X ) /2; f ( X ) = ( - cos X sh X + sin X ch X )/ ( 2 X 3 ) ; f , ( X ) = -X 4 f 1 ( X ); f , ( X ) = - ( sin X ch X + cos X sh X )/( 2 X ) ; f , ( X ) = -X 4 f , ( X ); f 5( X )= ( cos X ch X- 1 )/ 2; f , ( X ) = - f ;( X )/ X 4 ; f , ( X ) = -X 4 f Д).

Таким образом, задачу идентификации краевых условий по собственным частотам в терминах функции (7) можно сформулировать следующим образом: требуется найти неизвестные коэффициенты u i ( i = 1,2, з,4) по известным корням X k характеристического определителя (7). Решение можно получить следующим способом.

Пусть X k являются собственными значениями краевой задачи (1), (2). Тогда X k — корни уравнения характеристического определителя (7) [28]. Подставляя значения X k в (7), получаем систему уравнений для отыскания неизвестных а , ( i = 1,2, з,4):

  • а1 fl(Xk ) + а 2 f-(Xk ) + а з fl(Xk ) + а 4 ft(Xk ) + (а1 а 4 + а 2 а з) f5(Xk ) + а1 а з f,(Xk ) + а 2 а 4 f-(Xk ) = f)(Xk ). (8)

  • 3.    Метод нескольких систем нелинейных уравнений и его недостатки

Поскольку число собственных значений X k бесконечно, то (8) можно решать с любыми из них, а значит, система (8) может содержать любое число уравнений. Заметим, что функция (7) является четной целой функцией первого порядка. Отсюда следует, что по всему бесконечному набору корней X k характеристического определителя (7) неизвестные коэффициенты а ( i = 1,2, з,4) восстанавливаются однозначно (см. [15]).

Заметим, что для однозначного восстановления неизвестных G i ( i = 1,2, з,4) нет необходимости использовать весь бесконечный набор корней X k . Для ответа на вопрос, сколько же нелинейных уравнений следует взять, воспользуемся опытом авторов работы [15]. Ими установлено, что требуется решить несколько систем нелинейных уравнений с числом уравнений, совпадающим с числом неизвестных, а затем найти пересечение этих решений.

Подставляя значения X k ( k = 1,2, з,4) в (8), решаем систему четырех уравнений. Получается шесть решений. Добавляем к этой системе еще одно уравнение, используя в (8) еще и X 5. Заменяем в предыдущей системе уравнений четвертое уравнение на пятое. Решаем вторую систему четырех нелинейных уравнений с четырьмя значениями X k ( k = 1,2, з, 5). Опять получаем шесть решений. Не все решения второй системы совпадают с решениями первой. Значит, они не являются решениями всех пяти уравнений (8) со значениями X k ( k = 1,2, з, 4,5). Их можно не рассматривать в качестве решений. Допустим, осталось два общих решения. Можно ли удалить еще одно? Для ответа на этот вопрос добавляем еще одно уравнение: подставляем в (8) X 6 . Заменяем в предыдущей системе уравнений пятое уравнение на шестое. Решаем вторую систему четырех нелинейных уравнений с четырьмя значениями X k ( k = 1,2, з, 6). Снова имеем шесть решений. Если среди этих шести решений нет одного из двух общих решений системы пяти уравнений (8) со значениями X k ( k = 1,2, з, 4,5), то тем единственным решением, которое мы ищем, будет решение, пересекающее все шесть уравнений (8) со значениями X k ( k = 1,2, з, 4,5,6).

Однако этот метод имеет явные недостатки. Так, он не дает ответа на вопросы: сколько собственных значений Xk требуется для однозначной идентификации четырех неизвестных ai (i = 1,2,3,4); достаточно ли шести собственных значений; если этого количества уравнений недостаточно, то сколько собственных значений следует взять для однозначной идентификации; если приближенные решения имеют близкие значения, то можно ли считать их общими (пересекающимися), или это различные решения.

Продемонстрируем эти недостатки на примере.

Пример 1 . Пусть X 1 =3,861483, X 2 = 7,031643, X 3 =10,18502, X 4 =13,33266, X 5 =16,47796, X 6 =19,62205, X 7 = 22,76544 — первые семь собственных значений задачи (1), (2). Подставляя первые четыре из них в уравнение (8), получаем систему четырех нелинейных уравнений с четырьмя неизвестными, имеющую шесть наборов решений (эти решения найдены с помощью пакета аналитических вычислений Maple), среди которых два решения представляются следующим образом:

{ а 1 =0,00000, а 2 =2,00000, а 3 =0,00000, а 4 =0,00000},

{ а 1 =13653, а 2 =36,044, а 3 = - 0,47369, а 4 =0,00000}.

Решение системы четырех нелинейных уравнений (8) с X 1 =3,861483, X 2 = 7,031643, X 3 =10,18502 и X 5 = 16,47796, снова дает шесть наборов, среди которых два таковы:

{ а 1 =0,00000, а 2 =2,00000, а 3 =0,00000, а 4 =0,00000}, { а 1 =13656, а 2 =36,362, а 3 = - 0,47400, а 4 =0,00000}.

Другие наборы решений второй системы уравнений существенно отличаются от наборов решений первой системы.

Разрешая (8) с X 1 = 3,861483, X 2 = 7,031643 , X 3 = 10,18502 и X 6 = 19,62205 , получаем шесть наборов, среди которых два имеют вид:

{ а 1 =0,00000, а 2 =2,00000, а 3 =0,00000, а 4 =0,00000}, { а 1 =13657, а 2 =36,395, а 3 = - 0,47403, а 4 =0,00000}.

При X 1 =3,861483, X 2 = 7,031643, X 3 =10,18502 и X 7 =22,76544 из решения (8) следуют шесть наборов, и среди них два следующие:

{ а =0,00000, а 2 =2,00000, а 3 =0,00000, а 4 =0,00000}, { а 1 =13656, а 2 =36,383, а 3 = - 0,47403, а 4 =0,00000}.

Видно, что решение { а 1 = 0,00000, а 2 = 2,00000, а 3 = 0,00000, а 4 = 0,00000} является общим для всех систем уравнений, а вторые решения выделенной пары совпадают с точностью до двух-четырех значащих цифр. Можно ли считать вторые решения общими? На этот вопрос трудно ответить без общих теорем. Поэтому ниже предлагается метод решения системы (8), из которого следует общая теорема, утверждающая, что для выделения единственного решения достаточно семи собственных значений. Эта теорема также дает ответ на вопрос об общности вторых решений. И ответ этот — отрицательный.

4.    Метод дополнительных неизвестных

Пусть X1, X2, .. , X7 — собственные значения задачи (1), (2). Вводим дополнительные неизвестные а 5 = а1 а 4 + а 2 а 3, а 6 = а1 а 3, а 7 = а 2 а 4                                       (9)

и подставляем их вместе с собственными значениями X к ( к = 1,2,...,7). в уравнение (8). В результате получаем систему семи линейных уравнений с семью неизвестными n i ( i = 1,2,...,7):

а 1 f3 ( X к ) + а 2 f 2 ( X к ) + а 3 f , ( X к ) + а 4 f , ( X к ) + а 5 f ; ( X к ) + а 6 f ;( X к ) + а 7 f , ( X к ) = f , ( X к ) ( к = 1,2,...,7).     (10)

Из правила Крамера следует, что если определитель системы уравнений (10)

f (X1) fА) . •   f,(X1) D= f (X2) fД2)  . • f,(X2) (11) fl (X7 ) fД7)  . f,(X7) отличен от нуля, то коэффициенты ai (i = 1,2,...,7) однозначно вычисляются по формулам a1 = D1 /D, a 2 = D 2 /D, a 3 = D3 /D, a 4 = D4 /D,                         (12)

a 5= D 5 / D , a 6 = D 6 / D , a 7 = D 7 / D ,                                (13)

где D i — определитель, получающийся из D при замене элементов f ( X 1 ), f ( X 2), ..., f ( X 7) i -го столбца соответствующими свободными членами f , ( X 1 ), f , ( X 2), . , f 0 ( X 7).

В случае, если собственные значения X к ( к = 1,2,...,7) найдены с погрешностью, то может оказаться, что

D 5 / D * ( D i / D ) ( D 4 / D ) + ( D 2 / D ) ( D 3 / D ) ,

(D б/ D И Di /D )■( D3/ D), (D 7 /D )*( D 2/ D )■( D4 /D), и тогда система уравнений для определения четырех неизвестных ai (i = 1,2,3,4) становится переопределенной. То есть задача отыскания неизвестных ai (i = 1,2,3,4) по семи значениям Xк (к = 1,2,...,7) может не иметь решения, а поэтому является некорректной по Адамару. Однако она корректна по А.Н. Тихонову.

Действительно, для подхода А.Н. Тихонова к вопросу корректности характерно, что рассматривается некоторое множество M с V , существенно более узкое, чем все пространство V . Пусть образ множества M в пространстве Z при отображении с помощью оператора R есть множество Л , то есть Л = RM .

Задача Rv = z , где v из пространства V , а z из пространства Z , называется корректной по А.Н. Тихонову (условно корректной), если [22, 23, 24]:

  • 1)    априори известно, что решение задачи существует и принадлежит некоторому множеству M пространства V ;

  • 2)    решение единственно на множестве M ;

  • 3)    для любого б > 0 существует такое 5 > 0, что для любых z , z е Л = RM , и таких, что || z - z || z 5 , выполняется неравенство || v - v || <  б .

В обсуждаемом случае под оператором R можно понимать отображение, задаваемое системой уравнений (10), преобразующее семерку чисел a i ( i = 1,2,...,7) в семерку значений X к ( к = 1,2,...,7), при которых матрица системы имеет определитель D * 0 . Множеством корректности M будем называть такое множество семичисельных наборов v = ( a 1 , a 2,..., a 7), для каждого из которых выполняются условия (9). С помощью введенного множества корректности нетрудно показать корректность по А.Н. Тихонову задачи отыскания a i ( i = 1,2,3,4) по значениям X к ( к = 1,2,...,7), при которых система уравнений (10) имеет определитель, отличный от нуля.

Пусть V — это пространство R 7 элементов v = ( v 1 , v 2, . , v 7) с нормой || v || = max(| v 1 |,| v 2 |, . ,| v 7|); Z — это пространство R 7 элементов z = ( z 1 , z 2, . , z 7) с нормой || z || = max(| z 1 |,| z 2 |, . ,| z 7 |) ; образ множества M в пространстве Z при отображении с помощью оператора R есть множество Л , то есть Л = RM . Тогда задача Rv = z будет корректной по А.Н. Тихонову, так как все перечисленные выше условия выполнены (третье условие вытекает из аналитичности f ( X ) ( i = 0,1,...,7) по X ).

Наиболее известны два метода решения корректных по А.Н. Тихонову задач — метод квазирешения и метод подбора. В настоящей статье предлагается метод решения задачи идентификации краевых условий, который по сути представляет собой метод подбора.

Априори известно, что искомые числа ai (i = 1,2,3,4) существуют, так как идентифицируется реально существующая механическая система. Если действительные собственные значения X1, X2, ., X7 и определители D * 0, Di (i = 1,...,7) найдены точно, то условия (9) выполнены, причем числа ai (i = 1,2,3,4) однозначно вычисляются по формулам (12). Все условия корректности по А.Н. Тихонову реализованы, в том числе и третье. Действительно, для любого ε >0 существует такое δ >0, что для любых z = (X1, X 2,..., X7), z = (X 1, X 2,..., X 7) eA = RM , и таких, что ||Z - zl l»’ <δ, соблюдается неравенство ||(D1 /D,D2/D,...,D7/D)-(Ё)1 /D,Ё)2/D,...,Ё)7/D)|| i <ε. Последнее вытекает из аналитичности функций fi (λ) (i = 1,..., 7) по параметру λ .

Если числа λ k ( k = 1, 2,..., 7) , а значит, и D 0, D 1, D 2,... , D 7 установлены приближенно, то равенство (14) может не существовать, и поэтому формально по определителям D 0, D 1, D 2,... , D 7 (формулам Крамера) неизвестные ai ( i = 1,2,3,4) вычислить невозможно. Однако для предъявления приближенного решения (14) это и не нужно. Приближенным решением будем считать значения ai ( i = 1,2,3,4), найденные по формулам (12), где определители D 0, D 1, D 2, D 3, D 4 рассчитываются по приближенным λ k ( k = 1, 2,..., 7) . Приближенные решения ai ( i = 1,2,3,4) вместе с числами ai ( i = 5,6,7), определяемыми по формулам (9), а не по (13), лежат во множестве корректности, так как для них справедливы условия (9). Эти приближенные решения для ai ( i = 1, 2,..., 7) тем ближе к точным, чем ближе к точным собственным значениям числа λ k ( k = 1, 2,..., 7) , а значит, чем ближе к точным значениям определители D 0, D 1, D 2, ..., D 7. Таким образом, верна теорема:

Теорема. Если λ k ( k = 1, 2,..., 7) являются действительными собственными значениями краевой задачи (1), (2), причем определитель (11) системы уравнений (10) отличен от нуля, то задача отыскания коэффициентов ai ( i = 1, 2,..., 7) по собственным значениям λ k ( k = 1, 2,..., 7) является корректной по А.Н. Тихонову, и ее единственное решение описывается формулами (12, 13).

Пример 2 . Пусть собственные значения задачи такие же, как в первом примере, а именно λ 1 = 3, 861483 , λ 2 =7,031643 , λ 3 =10,18502, λ 4 =13,33266, λ 5 =16, 47796, λ 6 =19, 62205, λ 7 =22, 76544. Подставляя эти значения в (11), (12) и находя D = - 0, 3680993 10 31 , D 1 = - 0,1608490 10 22 , D 2 = - 0, 7361985 10 31 , D 3 = 0,1224316 10 18 , D 4 = - 0,1467282 10 13 , получаем единственное решение: a 1= D 1 D = 0, 000 , a 2= D 2 D =2,000, a 3= D 3 D =0,000, a 4= D 4 D =0,000 .

Пример 3 . Пусть λ 1 =1, 016782, λ 2 =4,696921, λ 3 =7,831946, λ 4 =10,98039, λ 5 =14,12531, λ 6 = 17, 26906 , λ 7 = 20, 41215 — первые семь собственных значений задачи (1), (2). Используя их, из (11) и (12) имеем D =0,4241723 10 19 , D 1 =0,2120861 10 20 , D 2 =0, 2545034 10 20 , D 3 =0,1272516 10 20 , D 4 = 0,1696689 10 20 , а также единственное решение a 1= D 1 D = 5, 000, a 2= D 2 D = 6, 000 , a 3= D 3 D = 3, 000 , a 4= D 4 D = 4, 000 .

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

5.    Заключение

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

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

  • Стрэтт Дж.В. (лорд Рэлей) Теория звука. Т. 1. -М.: Гостехиздат, 1955. -504 с.
  • Коллатц Л. Задачи на собственные значения (с техническими приложениями). -М.: Наука, 1968. -504 c.
  • Вибрации в технике: справочник в 6 томах. Т. 1. Колебания нелинейных систем/Под ред. В.В. Болотина. -М.: Машиностроение, 1978. -352 с.
  • Гонткевич В.С. Собственные колебания пластинок и оболочек. -Киев: Наукова думка, 1964. -288 с.
  • Акуленко Л.Д., Нестеров С.В. Частотно-параметрический анализ собственных колебаний неоднородных стержней//ПММ. -2003. -Т. 67, № 4. -С. 588-602.
  • Гладвелл Г.М.Л. Обратные задачи теории колебаний. -М.-Ижевск: НИЦ «Регулярная и хаотическая динамика», Институт компьютерных исследований, 2008. -608 с.
  • Movahhedy M., Ismail F. and Gladwell G.M.L. Reconstruction of mass-spring system from spectral data. II: Experiment//Inverse Probl. Eng. -1995. -V. 1, N. 4. -P. 315-327.
  • Иванов В.Н., Домбровский И.В., Шевелев Н.А. Численная идентификация параметров динамического поведения элементов машиностроительных конструкций//Вычисл. мех. сплош. сред. -2011. -Т. 4, № 3. -С. 58-67.
  • Юрко В.А. Введение в теорию обратных спектральных задач. -М: Физматлит, 2007. -384 c.
  • Садовничий В.А., Султанаев Я.Т., Ахтямов А.М. Обратные задачи Штурма-Луивилля с нераспадающимися краевыми условиями. -М.: Изд-во МГУ, 2009. -184 с.
  • Morassi A., Dilena M. On point mass identification in rods and beams from minimal frequency measurements//Inverse Probl. Eng. -2002. -V. 10, N. 3. -P. 183-201.
  • Ахтямов А.М. Теория идентификации краевых условий и ее приложения. -М.: Физматлит, 2009. -272 c.
  • Ахтямов А.М., Муфтахов А.В., Ямилова Л.С. Идентификация вида и параметров закрепления стержня по собственным частотам его колебаний//Акустический журнал. -2008. -Т. 54, № 2. -C. 181-188.
  • Ахтямов А.М., Сафина Г.Ф. Определение виброзащитного закрепления трубопровода//ПМТФ. -2008. -Т. 49, № 1. -C. 139-147.
  • Ахтямов А.М., Урманчеев С.Ф. Определение параметров твердого тела, прикрепленного к одному из концов балки, по собственным частотам колебаний//Сибирский журнал индустриальной математики. -2008. -Т. XI, № 4. -C. 19-24.
  • Ватульян А.О. Обратные задачи в механике деформируемого твердого тела. -М.: Физматлит, 2007. -224 с.
  • Ермолов И.Н., Алешин Н.П., Потапов А.И. Неразрушающий контроль: В 5 книгах. Кн. 2. Акустические методы контроля: Практ. пособие/Под ред. В.В. Сухорукова. -М.: Высшая школа, 1991. -283 с.
  • Вибродиагностика качества механизмов приборов: Сб. статей/Под ред. К.Н. Явленского. -Л.: ЛИАП, 1987. -144 с.
  • Ильгамов М.А., Хакимов А.Г. Диагностика закрепления и повреждений балки на упругих опорах//Контроль. Диагностика. -2010. -№ 9. -С. 57-63.
  • Ильгамов М.А. Диагностика повреждений вертикальной штанги//Тр. ин-та механики УНЦ РАН. -Уфа: Гилем, 2007. -Вып. 5. -С. 201-211.
  • Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. -М.: Наука, 1979. -284 с.
  • Иванов В.К., Васин В.В., Танана В.П. Теория линейных некорректных задач и ее приложения. -М.: Наука, 1978. -206 с.
  • Лаврентьев М.М., Романов В.Г., Шишатский С.П. Некорректные задачи математической физики и анализа. -М.: Наука, 1980. -288 с.
  • Тихонов А.Н., Гончарский А.В., Степанов В.В., Ягола А.Г. Численные методы решения некорректных задач. -М.: Наука, 1990. -232 с.
  • Тихонов А.Н., Леонов А.С., Ягола А.Г. Нелинейные некорректные задачи. -М.: Наука, 1995. -312 с.
  • Лаврентьев М.М., Резницкая К.Х., Яхно В.Г. Одномерные обратные задачи математической физики. -Новосибирск: Наука, 1982. -88 с.
  • Лаврентьев М.М. Теория операторов и некорректные задачи. -Новосибирск: Изд-во Ин-та математики, 1999. -702 с.
  • Наймарк М.А. Линейные дифференциальные операторы. -М.: Наука, 1969. -526 c.
Еще
Статья научная