Об идентификации переменной жесткости при анализе поперечных колебаний балки

УДК 539.32

А.О. ВАТУЛЬЯН, А.Ю. БУРЬЯН, А.В. ОСИПОВ

ОБ ИДЕНТИФИКАЦИИ ПЕРЕМЕННОЙ ЖЕСТКОСТИ ПРИ АНАЛИЗЕ ПОПЕРЕЧНЫХ КОЛЕБАНИЙ БАЛКИ

Представлен метод нахождения амплитудной характеристики для консольно закрепленной балки переменной жесткости на основе анализа интегрального уравнения Фредгольма 2-го рода, предложен метод определения закона изменения жесткости на основе решения задачи Коши.

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

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

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

Постановка задачи. Рассмотрим установившиеся колебания упругой балки длины l переменной жесткости EI с частотой ю . Уравнение имеет вид [2]:

где Е - модуль Юнга; I - момент инерции; р - плотность; F - площадь поперечного сечения.

Считаем, что балка на конце х1 = 0 жестко закреплена, а на конце х1 = I - действует нагрузка. Рассмотрим два варианта нагружения на конце -задача 1: нагружение сосредоточенным моментом; задача 2: нагружение сосредоточенной силой.

Соответствующие граничные условия имеют вид -задача 1:

Приведем краевую задачу (1), (2) к безразмерному виду, введя безразмерную координату

(EI (Xj)w" (Xj))"-praFw( Xj) = 0,

w(0) = w"(0) = 0 ; EI(l) • w"(l) = M; (EI(xj) • w" (Xj))'\x = = 0;

задача 2:

w(0) = w"(0) = 0 ; EI(l) • w" (l) = 0; (EI(xj) • w"(xj)) "|x = = P.

Таким образом, краевая задача (1), (2) в безразмерных переменных имеет следующий

вид:

Поскольку решение задач 1, 2 пропорционально нагрузке, то примем следующие модельные граничные условия -

задача 1: к(0) = к '(0) = 0; В(1) • к "(1) = 1; (В(х) • к "(х))'|^ = 0;

задача 2: к

Основная цель исследования:

1) разработка метода построения функции к(х,k) при различных значениях параметра

k при произвольных законах изменения жесткости В(х) (монотонных, немонотонных, кусочнонепрерывных);

2) разработка метода построения решения обратной задачи, которая состоит в нахождении функции В(х) по информации о функции к(х,k) при некотором значении k0;

3) проведение серии вычислительных экспериментов, позволяющих оценить эффективность предлагаемых методов.

Исследование прямой задачи. Прямая задача состоит в нахождении функции к , удовлетворяющей краевым условиям при заданном законе изменения В( х) в некотором диапазоне изменения k . Поскольку коэффициенты дифференциальных операторов переменны, то в аналитическом виде найти решение нельзя, возможно только численное решение задачи. Чаще для этого используют метод конечных разностей или метод сведения к интегральному уравнению Фредгольма 2-го рода [3].

Предлагая способ сведения краевой задачи к интегральному уравнению Фредгольма 2-го рода, опишем процедуру сведения на примере задачи 2.

Для построения интегрального уравнения введем функцию:

и(х) = В(х)к"(х) .

Далее перейдем от краевой задачи (3), (4) к системе дифференциальных уравнений 2-го порядка относительно функций и(х) и к(х):

и" (х) - k 4 к( х) = 0,

К "

В

Соответствующие граничные условия приобретают вид

к(0) = к ' (0) = 0,

,

и(1) = 0, и' (1) = 1.

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

х

1(х) = k4 Г(х - + А1 х + А2,

х°

и

0 В

^х) = | (х-^) ^ + Бх + Б2.

Постоянные А, А2, В1, В2 находим из граничных условий (6):

1 11 1

В1 = В2 = 0, А1 = 1 - k4 Гw(^)d^ , А2 = k4Гк(^- 1 - k4Г(1 -£)к(^=£41^(^-1.

0 0 0 0

Таким образом, имеем соотношения для определения искомых функций:

и (У) л

w

0 Б(У

1 1

и( х) = k 41 (х - ^^(^) й£, + х - k 4 х | w(^)d^ + k 41 ^(^)й^ -1 = (х -1) + k 41 (х - ^^(^)й^

00

11 - k41(х - = (х -$) - k41(х - ^М^.

0 х

Объединив эти представления, получим:

w

Поменяем порядок интегрирования в двойном интеграле в представлении (8), переписав его в виде интегрального уравнения Фредгольма 2-го рода:

1

w( х) = у

0

причем его ядро непрерывно для широкого класса функций Б(х) и представимо в форме

K(х,5) = <

4х ^ -У) ^, 0 < 5 < х,

Б(У) ^

Ч х -У)

У(х) = |

Б

0 < х < 5,

0 Б

Для задачи 1 соответствующее уравнение получается аналогично и имеет вид (9), но с другой правой частью, а именно:

'(х-УЬе

1(х)=1 -

0 Б

Уравнения вида (9) можно исследовать либо методом последовательных приближений, либо численно. Наиболее распространенным методом является метод коллокаций. Дискретизацией на равномерной сетке h = 1/ N, хк = кк, k = 0,..., N, хм = 1, заменой интеграла квадратурной формулой Симпсона [5] и использованием метода коллокаций уравнение (9) приводится к исследуемой численно системе линейных алгебраических уравнений:

N

Ч х) - к 4 2 АК

}=1

Анализ конечномерного оператора (10) позволяет эффективно строить численное решение и определять как резонансные частоты, так и амплитудно-частотные характеристики для балок переменной жесткости.

Для оценки точности предложенной схемы вычислений сравним численное решение, полученное в результате использования составленной программы, и точное решение, полученное аналитически при Б = 1, которое выражается через функции Крылова [2].

Расчеты показывают, что относительная погрешность между этими решениями не превосходит 0.5% в диапазоне частот до первого резонанса, что свидетельствует о достаточно эффективной вычислительной схеме.

Результаты расчетов амплитудно-частотной характеристики представлены в виде графиков.

На рис.1 результаты решения задачи 2 изображены сплошной линией при k = 5, а ромбиками - при k = 1 для В(х) = 1+ х2.

На рис.2 изображен результат решения задачи 1, сплошной линией показано решение при k = 5 , а ромбиками - при k = 1 для В(х) = 1 + зт2(10лх).

№(*) _

0,4 -

0,3 -

0,2 -

0,1 -

о --0,1 --0,2 -

Рис.2. Решение прямой задачи 1

Обратная задача. Целью решения обратной задачи является восстановление функции в(х) по известным в наборе точек значениям функции ^(х, k0). При этом рассмотрены варианты постановки обратной задачи:

задача 3: Восстановление функции в(х) по значениям функции ^(х, &0) (которая в рамках вычислительного эксперимента взята из решения прямой задачи 1);

задача 4: Восстановление функции в(х) по значениям функции ^(х, &0) (которая в рамках вычислительного эксперимента взята из решения прямой задачи 2).

Заметим, что подобные обратные задачи рассмотрены и в другой постановке [6, 7]. В данном случае дифференциальное уравнение (3) относительно функции в(х) имеет 2-й порядок

и для ее определения необходимо решать задачу Коши. Некорректность обратной задачи в этой постановке состоит в том, что необходимо находить вторые производные от функции, заданной в наборе точек. Поэтому на первом этапе аппроксимируем функцию ^(х) сплайнами 5-й степени [5], что позволяет находить вторую производную с достаточной степенью точности. На втором этапе построим операторное соотношение, связывающее функции в(х) и ^(х). Осуществим это построение на примере задачи 4. Дважды интегрируя уравнение (3) от 1 до х, получаем:

х х х

| (в (х)w"(х)) dЕ, = | (в (х)w"(х)) - 1 = & 41 w (£)/Е, ,

и далее имеем:

в

*

Для задачи 3 аналогичное соотношение запишется в виде

1

1 - k 41(х -£) *(£)/£

в

*

Дискретизация позволяет найти значения искомой функции:

ІЇХ, )=^

*'' '

Для вычисления определенного интеграла I (х* -£)*(г;У£

1

х*-1 - *4 |(х* -£)*(£)й£

х*

1

'х.

заменим подынтеграль-

х*

ную функцию на каждом отрезке разбиения многочленом Лагранжа 1-й степени с узлами 1 и х [5]. Это соответствует замене на каждом отрезке разбиения кривой на секущую. В результате получаем известную формулу трапеций:

1

I(х* -£)*(£У£ = Ь,[(хк - х*+1)*(х*+1 )+ (х* - х*+ 2)*(х*+2 )+ - +

х

+ (х* - хы-1 >(хы-1 ) +1/2(х* - хы )*(х^ )]

Применение этой формулы к (11) дает расчетные зависимости.

B(xk) = „( ) [xk -1 - k4Щх

і )w(xk+1 )+(х

г)w(xk + 2 ) +

I +... +

Путем аналогичных построений получаем расчетные формулы для задачи 3:

в(Ч1 = •1 - k4^[(xk - xk+і +і1+

М

+ ( Хк - Хк+2^(хк + 2 )+ - + (Хк - ХМ-1)W(XN-1)+ 1/2( Хк - ХМ )W(XN )]•

Как видно из представления (11), задача 4 хуже обусловленная, нежели задача 3, так как на конце х = 1 имеет место неопределенность, поэтому в этой точке необходимо либо задавать значение функции в(х), либо экстраполировать функцию из области х < 1. Этому факту есть простое объяснение в рамках рассмотрения задачи Коши. Для задачи 3 граничные условия (4) дают два необходимых условия в задаче Коши, а для задачи 4 - только одно, поскольку из условия В(1) • w"(1) = 0 и положительности В( х) следует w"(1) = 0, и второе условие пропадает.

На основе описанного метода проведены эксперименты по реконструкции различных видов неоднородностей: гладкой монотонно возрастающей функции, гладкой монотонно убывающей функции, гладкой немонотонной функции, функции, имеющей на отрезке [0, 1] несколько стационарных точек, кусочно-разрывной функции.

На графиках представлены вычислительные эксперименты для различных видов неоднородности.

Рис.3 иллюстрирует результат задачи 4 по восстановлению функции В( х) = 1+ х2 (гладкой монотонно возрастающей функции) при к = 1:

Рис.3. Решение обратной задачи 4

В качестве примера на рис.4 изображен результат восстановления функции в(х)= 1 + эт2(10л;х) (функции, имеющей несколько стационарных точек) в задаче 3 с входными

данными, зашумленными аддитивным образом с амплитудой 5 = 10-5.

>Кх)

Рис.4. Решение обратной задачи 3 с зашумленными входными данными

В случае разрывных законов неоднородности проведено разработанный метод усовершенствован: входной массив данных исследуется на предмет точек разрыва (рассматривается изменение первой и второй производной функции w(x)). При обнаружении точек разрыва в их окрестности проводится экстраполяция функции, что улучшает качество восстановления.

На рис.5 изображен результат восстановления функции:

В(х) =

1 3

1 при 0<х< —, х>—, 3 4

13

5 при — < х < —, 34

к = 0,9 .

Рис. 5. Решение обратной задачи 4 для кусочно-разрывной функции

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

Библиографический список

1. Ватульян А.О. Обратные задачи в механике деформируемого твердого тела / А.О. Ва-тульян. - М.: Физматлит, 2007. - 224 с.

2. Филиппов А.П. Колебания деформируемых систем. - Изд. 2-е перераб. / А.П. Филиппов. - М.: Машиностроение, 1970. - 736 с.

3. Краснов М.Л. Интегральные уравнения (введение в теорию) / М.Л. Краснов. - М.: Наука, 1975. - 303 с.

4. Тихонов А.Н. Методы решения некорректных задач / А.Н. Тихонов, В.Я. Арсенин. - М.: Наука, 1979. - 285 с.

5. Березин И.С. Методы вычислений / И.С. Березин, Н.П. Жидков. - М.: Гос. изд-во физ.-математ. лит., 1959. - 620 с.

6. Бочарова О.В. О реконструкции плотности и модуля Юнга для неоднородного стержня / О.В. Бочарова, А.О. Ватульян //Акустический журнал. - 2009. - Т.55. - №3. - С.281-288.

7. Боброва А.Н. Об определении закона изменения модуля Юнга при анализе продольных колебаний стержня / А.Н. Боброва, А.О. Ватульян // Вестник ДГТУ. - 2009. - №4. - С.613-621.

References

1. Vatul'yan A.O. Obratnye zadachi v mehanike deformiruemogo tverdogo tela / A.O. Vatul'yan. - M.: Fizmatlit, 2007. - 224 s. - in Russian.

2. Filippov A.P. Kolebaniya deformiruemyh sistem. - Izd. 2-e pererab. / A.P. Filippov. - M.: Mashinostroenie, 1970. - 736 s. - in Russian.

3. Krasnov M.L. Integral'nye uravneniya (vvedenie v teoriyu) / M.L. Krasnov. - M.: Nauka, 1975. - 303 s. - in Russian.

4. Tihonov A.N. Metody resheniya nekorrektnyh zadach / A.N. Tihonov, V.Ya. Arsenin. - M.: Nauka, 1979. - 285 s. - in Russian.

5. Berezin I.S. Metody vychislenii / I.S. Berezin, N.P. Jidkov. - M.: Gos. izd-vo fiz.-matemat. lit., 1959. - 620 s. - in Russian.

6. Bocharova O.V. O rekonstrukcii plotnosti i modulya Yunga dlya neodnorodnogo sterjnya / O.V. Bocharova, A.O. Vatul'yan //Akusticheskii jurnal. - 2009. - T.55. - №3. - S.281-288. - in Russian.

7. Bobrova A.N. Ob opredelenii zakona izmeneniya modulya Yunga pri analize prodol'nyh kole-banii sterjnya / A.N. Bobrova, A.O. Vatul'yan // Vestnik DGTU. - 2009. - №4. - S.613-621. - in Russian.

Материал поступил в редакцию 15.05.2010.

A.O. VATULYAN, A.Y. BURYAN, A.V. OSIPOV ON IDENTIFICATION OF VARIABLE RIGIDITY UNDER THE ANALYSIS OF THE BEAN LATERAL OSCILLATION

The method of amplitude characteristic estimation for the cantilever bean of the variable rigidity on the basis of the analysis of Fredholm integral equation of the second kind is presented. The determination method of rigidity changing on the basis of Cauchy sequence analysis.

Key words: beam of the variable rigidity; coefficient inverse problem; Fredholm integral equation of the second kind; spline approximation.

ВАТУЛЬЯН Александр Ованесович (р. 1953), доктор физико-математических наук (1993), профессор (1995), заведующий кафедрой теории упругости ЮФУ (2001), член Российского Национального комитета по теоретической и прикладной механике, член трех докторских диссертационных Советов (ЮФУ, ДГТУ). Окончил механико-математический факультет РГУ (1975).

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

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

БУРЬЯН Алиса Юрьевна, выпускница кафедры «Математика» ДГТУ (2010).

Область научных интересов - обратные и некорректные задачи для дифференциальных операторов, численные методы решения некорректных задач.

ОСИПОВ Алексей Владимирович (р. 1988), выпускник кафедры «Математика» ДГТУ (2010). Область научных интересов - обратные и некорректные задачи для дифференциальных операторов, численные методы решения некорректных задач.

kukuvz@yandex.ru

Alexander O. VATULYAN (1953), Head of the Elasticity Theory Department, South Federal University (2001). PhD in Physics and Maths (1993), professor (1995), member of the Russian National Committee for Theoretical and Applied Mechanics. He graduated from the Mechanics and Mathematics Faculty, South Federal University (1975).

Research interests - investigation of wave processes in bodies with complex mechanical properties (heterogeneity, anisotropy, rheology, piezoeffect occurrence), geometric and coefficient inverse problems of mathematical physics, numerical methods.

Author of more than 300 scientific publications, including 1 monograph.

Alisa Y. BURYAN (1989) graduated from the Maths Department, Don State Technical University (2010).

Research interests - inverse and incorrect problems for differential operators, numerical methods of solving incorrect problems.

Alexey V. OSIPOV (1988) graduated from the Maths Department, Don State Technical University (2010).

Research interests - inverse and incorrect problems for differential operators, numerical methods of solving incorrect problems.