CC BY
18
Владикавказский математический журнал 2012, Том 14, Выпуск 3, С. 31-44
УДК 519.46
ОБРАТНЫЕ КОЭФФИЦИЕНТНЫЕ ЗАДАЧИ ДЛЯ ДИССИПАТИВНЫХ ОПЕРАТОРОВ И ИДЕНТИФИКАЦИЯ СВОЙСТВ ВЯЗКОУПРУГИХ МАТЕРИАЛОВ1
И. В. Богачев, А. О. Ватульян
Представлена общая постановка обратной задачи идентификации неоднородных характеристик вяз-коупругого тела, сформулировано вариационное уравнение. Рассмотрена задача о восстановлении функций-коэффициентов диссипативных операторов, возникающих при решении ряда задач об идентификации свойств слоистых неоднородных вязкоупругих структур на основе анализа спектральных характеристик. Представлены способы построения итерационного процесса, приведены результаты вычислительных экспериментов по восстановлению функций различных типов.
Ключевые слова: обратные задачи, диссипативные операторы, вязкоупругость, идентификация, итерационный процесс.
Введение
Исследование характеристик материалов со сложными неоднородными свойствами [1], таких как полимеркомпозиты, функционально-градиентные материалы, геологические породы, биологические ткани, в настоящее время является одним из важнейших направлений механики сплошной среды. Вследствие сложности прямых экспериментальных оценок механических свойств таких материалов со сложной реологией важна разработка иных методов идентификации неоднородных характеристик, основанных на различных моделях вязкоупругости. Кроме того, в связи со спецификой самих материалов (например, биологических тканей) интерес представляют неинвазивные методы [2]. Одним из способов воздействия является акустическое зондирование, при специальной обработке результатов которого для упругих тел [3, 4] удается восстанавливать неизвестные функции по информации об амплитудно-частотных характеристиках, измеренных в некоторых точках исследуемого объекта [5]. Отметим, что исследование установившихся колебаний тел в рамках линейной вязкоупругости для модели стандартного вязкоупру-гого тела приводит к краевым задачам с переменными характеристиками для диссипа-тивных операторов. Соответственно, решение обратных задач об определении функций координат, характеризующих неоднородность, приводит к некоторым обобщениям методов, использовавшихся ранее для положительных операторов.
В случаях исследования более сложных моделей, учитывающих неоднородную структуру материалов, ситуация значительно усложняется и требуется более глубокий анализ
© 2012 Богачев И. В., Ватульян А. О.
1 Работа выполнена при частичной поддержке Российского фонда фундаментальных исследований, проект № 10-01-00194-а, и Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг., госконтракт № П596.
проблемы, основанный на исследовании решений общих краевых задач с переменными коэффициентами и изучение зависимости следов решений от коэффициентов [6].
На начальных этапах исследования таких задач в зарубежной литературе была широко распространена постановка, в которой известны физические поля внутри исследуемого объекта [7, 8]. Обратная задача об определении характеристик материалов в таком случае оказывается линейной и сводится к решению системы дифференциальных уравнений первого порядка или к интегральному уравнению Фредгольма первого рода. Проведено исследование и построение численных решений с использованием различных регуляризованных подходов (метод А. Н. Тихонова [10, 11], метод минимизации расширенного лагранжиана [9], резуляризованные методы обращения разностных схем [8]). Если информацию о следах решений можно получить только на границе, то обратная задача существенно нелинейна. В последнее время исследования связаны с более сложной постановкой, в которой известны (измерены) граничные поля в некотором диапазоне изменения спектрального параметра (частоты колебаний зондирующего возмущения). Она сводится к нелинейным операторным уравнениям, которые содержат промежуточные переменные — компоненты полей смещений. Задачи в такой постановке могут быть исследованы лишь на основе некоторых итеративных процедур, основные принципы построения которых опираются на метод линеаризации и слабую постановку [12, 13].
Результаты исследования таких задач на примере определения свойств одномерных упругих объектов (стержней, балок) представлены в статьях [11, 13, 19]. В работах для исследования задач о восстановлении модулей упругости и плотности упругого или вязкоупругого неоднородного стержня при возбуждении продольных и/или изгибных колебаний предложен итерационный алгоритм, основанный на аппарате интегральных уравнений Фредгольма первого и второго рода.
Для областей типа слоя прямые и обратные коэффициентные задачи исследованы, например, в работе [15], где представлен алгоритм восстановления характеристик-коэффициентов, произвольно меняющихся по глубине, построен итерационный процесс, на каждом шаге которого решаются интегральные уравнения Фредгольма первого рода относительно функций-поправок.
Также распространение получили задачи восстановления коэффициентов вязкоупру-гих моделей и неоднородных характеристик тела при возбуждении колебаний в нем [16, 17]. В работе [18] исследованы математические аспекты ультразвукового сканирования биологических тканей и реконструкции их свойств, в частности модуля сдвига, при этом в качестве регуляризующего алгоритма при нахождении градиентов полей использован метод В-сплайнов.
Отметим принципиальную разницу в процедуре идентификации упругих и вязкоуп-ругих моделей. Для вязкоупругих моделей в рамках концепции динамических модулей [1] процесс колебаний описывается комплекснозначными операторами, обладающих свойством диссипативности.
В настоящей работе приведена общая постановка обратной коэффициентной задачи для трехмерного вязкоупругого тела, для которой с помощью процедуры варьирования сформулированы вариационные уравнения, позволяющие строить итерационные процессы. В качестве примера рассмотрена задача о восстановлении свойств изотропного неоднородного по толщине вязкоупругого слоя на основе данных о полях смещений, измеренных на верхней границе слоя. Схема решения основана на предварительном сведении с помощью преобразования Фурье к двум более простым однотипным несвязанным краевым задачам относительно осредненных характеристик, каждая из которых сводится к коэффициентной обратной задаче для диссипативного оператора, что позволяет исполь-
зовать построенные ранее итерационные схемы и регуляризационные алгоритмы для вещественных операторов [11-14, 19]. Представлены некоторые общие свойства решений и результаты вычислительных экспериментов.
1. Общая постановка обратной коэффициентной задачи
Общая постановка прямой задачи об определении поля смещений и напряженно-деформированного состояния при анализе установившихся колебаний с частотой ш вязко-упругого тела, занимающего область V и ограниченного кусочно-гладкой границей Б, в рамках концепции динамических модулей имеет вид [4]:
+ рш2щ = 0, (1)
(?тз = Ст^к1(х1гш)пку1, (2)
иг к =0, Щ к = Рг. (3)
Здесь Сгцк1(х,гш) — компоненты тензора динамических модулей, являющиеся кусочно-непрерывными комплексными функциями координат и частоты, р — плотность, пц — компоненты единичного вектора внешней нормали к Б .
В обратной задаче, используя дополнительную информацию,
иг = ¡г, ш е [Ш1 ,Ш2], (4)
характеризующую общие сведения о граничных полях смещений, которые доступны для непосредственного измерения в некотором частотном диапазоне, требуется определить не только поля смещений, но и законы изменения Сцк1(х,гш) и плотности как функций координат.
Введем множество функций пространственных координат Нт(V), характеризующих законы изменения динамических модулей и плотности, которые могут иметь конечное число разрывов первого рода на некоторых поверхностях внутри области V, которые разбивают ее на конечное число подобластей Vm; внутри каждой из этих подобластей функции бесконечно дифференцируемы. Эти функции могут быть кусочно-постоянными в случае включения или равны нулю в случае полости на некоторых из этих подобластей.
В работе [4] сформулированы и доказаны некоторые утверждения, необходимые для исследования задачи в упругом случае. Если рассматривать в них Съцы(х,гш) как функции не только координат, но и частоты, можно установить некоторые аналогичные свойства и для вязкоупругого случая. Введем понятие возможного поля. Под возможным полем понимается любое непрерывное в V поле смещений, поле компонент тензора динамических модулей и плотности с характеристиками неоднородности из Нт{У), которые удовлетворяют уравнениям (1)-(2) и граничным условиям (3). Используя принцип соответствия, легко установим, что имеет место
Свойство 1. Для любого возможного поля выполняется равенство
J 2Ь(щ,Уг,Сгзы(х,ш),р(х)) (IV + Jры йБ = 0, (5)
где 2Ь(щ,Уг,СцЫ(х,гы), р) = рш2щы - Сгцк1(х1гш)ик^гц.
Проварьировав соотношение (5) по всем переменным, используя формулу Гаусса — Остроградского, можно установить
Свойство 2. Имеет место следующее вариационное уравнение
J 2 Ь(щ ,иг,5Сцы(х,ш),5р(х)) йУ — J Ргбщ йБ = 0. (6)
V 5СТ
Уравнение (6) является линейным относительно вариаций и порождает линейный оператор, действующий из Нт(У) в Сте([шьш2]),
Т(5Сцы,5р) = J (^р(х) ш2ии — 5Сцы(х, гш) и^и^-) йУ.
V
Ядра этого оператора при любых ш (в отличие от упругого случая) представляют собой функции, аналитические по ш и суммируемые по пространственным координатам, поэтому соответствующий оператор является вполне непрерывным. Использование уравнения (6) позволяет сразу строить итерационную процедуру одновременного определения динамических модулей и плотности, заключающуюся в построении последовательности решений и(п), С^^х, гш), р(п). На первом шаге выбирается некоторое начальное
распределение модулей и плотности С(0)г(х, гш) и р(0) (например, в простейшем классе функций — линейных или кусочно-постоянных из условия минимума функционала невязки). Далее, если на п — 1 шаге и(п 1) есть решение краевой задачи (1)-(3) с известными С^ь1) (х,гш) и р(п_1), то элементы последовательности должны удовлетворять следующему линейному операторному уравнению первого рода
12Ь(и(Г~1),иг(п"1) ,С<$(х,«ш),р(п)) йУ — |рг(,П — и^) йБ = 0, ш £ [ш1,Ш2], (7)
V £ст
которое является следствием (6).
Отметим, что одного уравнения (7) недостаточно для однозначного определения всего набора искомых функций С^^х, гш), р(п), в общем случае требуется использовать различные способы нагружения, но в ряде ситуаций этого уравнения вполне достаточно.
Далее рассматривается задача о реконструкции вязкоупругих характеристик в неоднородном по толщине слое, для которой можно значительно упростить процедуру построения решения.
2. Постановка двумерной задачи об идентификации параметров вязкоупругого слоя
Рассмотрим установившиеся колебания с частотой ш вязкоупругого изотропного неоднородного по толщине слоя, занимающего область У = {х1 £ (—го, го), х2 £ (—то, го), хз £ [0, Л,]} в условиях плоской деформации. Нижняя грань слоя Б1 жестко защемлена, на верхней границе приложены нагрузки, определяемые вектором р, где р = (р1,0,рз).
Считая, что у вектора смещений ненулевыми являются компоненты и1 = и1(х1, хз), и3 = и3(х1, хз), получим краевую задачу:
(А + 2^)и1Д1 + Аиз,з1 + 0и(и1,з + из,1)) ,з +рш2и1 = 0,
((А + 2^)из,з);з + (Аи1,1 ),з + (^(и1,з + изд)),1 +рш2из = 0,
и^ =0, г = 1,3, (8)
Аим + (А + 2^)из,з|^2 = Рз,
^(и1,з + и1,з)|^2 = Р1,
где х3,ш) = ^ = и _ неизвестные ком.
плекснозначные функции (аналоги характеристик Ляме) и плотность неоднородного слоя соответственно. Функции »1(х3), А1(х3) и »2(х3), А2(х3) характеризуют мгновенные и длительные модули упругости соответственно, п > 0 — время релаксации. При этом из общих свойств операторов вязкоупругости [1] следуют ограничения 0 < А1(х3) < А2(х3) и 0 < ^1(хз) < ^2(хз).
Обратная коэффициентная задача в рассматриваемом случае сводится к определению четырех функций »1(х3), А1 (х3), »2(х3), А2(х3) по известной информации о полях смещений, заданных на верхней границе слоя при известном законе изменения плотности в некотором частотном диапазоне [ш1,ш2].
3. Формулировка операторных соотношений
Основной проблемой, возникающей при исследовании коэффициентных обратных задач, является построение операторных соотношений, связывающих искомые и заданные функции. На основе свойств 1, 2, сформулированных в общем случае, можно построить итерационные процессы. Для рассматриваемой постановки возможно значительное упрощение ситуации и разделение общей задачи на две подзадачи. Применяя интегральное преобразование Фурье по переменной х1 к задаче (8) и проводя несложные преобразования, получим систему обыкновенных дифференциальных уравнений второго порядка относительно трансформант компонент поля смещений и1(а1,х3,ш), и3(а1 ,х3,ш):
-а? (А + 2/л)и1 - га1Аиз,з + (^1,3 - га.1 Щ)),з +рш2171 = 0, ((А + 2р)и3,з),з - га1 ((АЩ),з + »и13) - а2»Щ + рш2и = 0, Цг|хз=0 = 0, (А + 2р)Щ,з - га1Аи 1Хз=н = рз, р(и1,з - га1Щ)1хз =ь = Р1, г = 1,3.
Полагая в задаче (9) а1 = 0, получаем две однотипные несвязанные краевые задачи относительно усредненных характеристик смещений, в одну из которых входит неизвестная комплекснозначная функция » = »(х3,гш), а во вторую — А(х3,гш) + 2»(хз,гш).
Задача 1. Первая краевая задача имеет вид:
(»и1,з),з +рш21/1 =0,
Щ |хз=0 =0, »Щ1,з1хз =ь = Р1-
Задача 2. Вторая краевая задача имеет вид:
((А + 2»)из,з),з + рш2из = 0,
из|хз=о =0, (А + 2»)из,з|хз=н = Рз, ( )
где для простоты записи обозначим иг(х3,ш) = иг(0,х3,ш), г = 1, 3.
Обратную коэффициентную задачу будем решать на основе информации об интегральных характеристиках полей смещений в случаях задач 1, 2, измеренных на верхней границе слоя.
Восстановление законов изменения комплексных функций по данным о спектральных свойствах решений сводится к последовательному решению двух обратных коэффициентных задач.
Обратная задача 1 (ОЗ 1). Определить пару функций »1 (х3), »2(х3) из (10) по информации
и (Н,ш)= ¡1(ш), Ш £ [Ш1,Ш2]. (12)
Далее, если ^1(хз) и (хз) определены из задачи 1, то формулируется Обратная задача 2 (ОЗ 2). Определить пару функций А1(хз), А2(хз) из (11) по информации
Цз(Л,ш) = /з(ш), ш £ [шз,ш4]. (13)
Отметим, что (10) и (11) представляют собой однотипные краевые задачи для обыкновенных дифференциальных уравнений второго порядка с переменными комплексными коэффициентами относительно усредненных характеристик полей смещений. Для произвольных законов изменения ^1(хз), ^2(хз), А1(хз), А2(хз) решение возможно лишь численно с привлечением аппарата численного исследования интегральных уравнений Фредгольма 1-го и 2-го родов. Подобная обратная задача возникает при описании продольных колебаний неоднородного по длине упругого стержня, которая была исследована ранее в работе [14].
Решение обратной задачи о реконструкции вязкоупругих характеристик приводит к нелинейной системе интегральных уравнений, которая может быть исследована на основе построения итерационного процесса. В то же время необходимо изучить некоторые общие свойства решений, важные при анализе обратных задач при произвольных законах изменения искомых характеристик.
В силу однотипности задач (10), (12) и (11), (13) построим операторные соотношения
для вспомогательной задачи, которая получается после введения безразмерных перемен-
22
ных х = Хз/к, и = Ц~г/к, к2 = — и безразмерных характеристик
(С(х,гк)и' (х, к))' + к2 г(х)и(х, к) = 0, и(0, к) = 0, С(1,гк)и' (1, к) = 1, ( )
при наличии дополнительного условия
и(1,к)= /(к), к £ [к1, к2], (15)
причем диапазон изменения спектрального параметра к выбирается в окрестности одного из максимумов амплитудно-частотной характеристики.
Здесь г (ж) = рор(х), С(х,ж) = гткэ(х)+и(х) ^ где ^^ = ^^ = ддя задачи х
И 9(х) = Л2(^(0)2(Ж)^ ^ = Х1(ТтФ) для задачи 2, г = п^.
Обратная задача сводится к определению пары функций Л(х), #(х), 0 < Л(х) < #(х), (г(х) — известная положительная функция) из краевой задачи (14)-(15).
Проведем исследование рассматриваемой задачи. Отметим некоторые свойства функции С(х,гк) и решения и(х, к).
Отделим вещественную и мнимую части в выражении для функции комплексного модуля
С(х, гк) = И,е(С(х, гк)) + г 1ш(С(х, гк)) = (^(х)г2 (к) + Л(х)г1 (к)) + г(#(х) — Л,(х))гз (к),
где *1(/с) = г2(к) = г3(к) =
Свойство 3. Имеет место И,е(С(х,гк)) > 0, 1ш(С(х,гк)) ^ 0. Введем в рассмотрение следующий оператор
Аи = —(Си')', и(0, к) = 0, и'(1, к) = 0.
Отметим, что введенный оператор обладает свойством диссипативности.
Свойство 4. Оператор А — диссипативен.
Поскольку согласно определению, свойством диссипативности обладает линейный оператор А с областью определения О а, плотной в гильбертовом пространстве Н, такой, что 1ш(Ах,х) ^ 0 при х £ О а, проверим его выполнение. Действительно, находим
1
1т((Аи,и)) = 1т((-(Си')', и)) = -(Си')' и йх
0
1 1 1
, \ / [ гткд + Н \ [ (д - К)тк
= 1т / Си и йх = 1т / -и и йх = / -и йх ^ 0.
\J ) \] 1 + гтк ) } 1 + т2 к2
0 0 0
Построим решение краевой задачи (14) в виде разложения в ряд по степеням к
те
и(х,к) = ^ Щ (х)к3. (16)
3=0
Учитывая вид комплексного модуля С(х, гк), перепишем уравнение задачи (14) в виде
((гткд(х) + Н(х))и'(х, к))' + к2г(х) и(х, к) + гтк3г(х) и(х, к) = 0. (17)
Свойство 5. Имеют место следующие рекуррентные соотношения
ип+1 =
причем
х 1
~Ш^и'п(л) + щ/ (г(0ип-Л0 + гтг(0ип-2(0)й^ о п
х х
[ йп Г тд(п) .
ПО = ]Щу
00 х 1
йп, п = 2,3,
и2
Н(п) Н(п)
о п
йп.
Эти соотношения легко устанавливаются при подстановке (16) в (14), формулировке краевых задач, получающихся при одинаковых степенях к, и нахождении их решений.
4. Формулировка итерационного процесса
Итерационный процесс, подобный описанному в [11-14, 19], может быть построен путем последовательного решения систем интегральных уравнений Фредгольма 1-го и 2-го рода, полученных с использованием метода линеаризации.
Итерационный процесс можно построить следующим образом:
1. Из интегрального уравнения Фредгольма 2-го рода (18) по известным п -1 приближениям неизвестных функций Нп-1 (х) и дп-1 (х) (из которых можно найти Сп-1(х,гк)) определим соответствующую функцию ип-1 (х, к)
1 х
ип-г(х) = ! Кп_1 (£, х)ип-х (£) й£ + j ——^ , (18)
0 0 п- 1
где = г^)^-), £ € [0,1].
Для решения интегрального уравнения (18) использован метод коллокаций, при этом интегральный оператор аппроксимируется конечномерным и задача сводится к решению алгебраической системы относительно узловых неизвестных.
2. Построим следующие приближения неизвестных функций Нп (ж) = Нп-1 (ж) + ¿Н(ж) и дп(ж) = дп-1(ж) + ¿д(ж), где $Н(ж) и $д(ж) — поправки для длительного и мгновенного модуля, входящие в функцию $С(ж,т), поправки комплексного модуля, должны удовлетворять интегральному уравнению Фредгольма 1-го рода (см. свойство 2)
1
/«(*,ж),,п-12(*,к)¿ж = „п-1(1,к) - /(К), К 6 [К1,К2]. (19)
0
Тогда из интегрального уравнения (19) получим систему: 1
0
1
(5д(ж),г2(к) + 5к{х)х\{к))и'п_1{х, к)и'п_1 (ж, к) с1х = 11е(ига_1(1, к) — /(к)), 1
/ (<**) - = 1т(..„-1(1,к) - /(к)).
(20)
Система (20) представляет собой систему интегральных уравнений Фредгольма 1-го рода с суммируемыми ядрами (при любых законах изменения неоднородности), решение которой является некорректной задачей и требует использования регуляризационных методов, например, метода регуляризации А. Н. Тихонова [10].
3. В результате реализации итерационного процесса получим функции Нп (ж) и дп(ж), которые дают приближенное решение обратной коэффициентной задачи (14) при удовлетворении условию выхода, в противном случае позволяют скорректировать ядра интегральных операторов (20) для построения решения на следующей итерации. В качестве условия выхода могут быть использованы следующие:
[0,1] <е, Рд||с[0,1] <е, (21)
К-1(1,к) - /(к)||с[0,1] <е, (22)
п > N (23)
где С[0,1] — пространство непрерывных на [0,1] функций; N — предельное значение максимально допустимого количества итераций.
5. Численная реализация
Следуя изложенной выше схеме, был проведен вычислительный эксперимент по восстановлению неизвестных безразмерных функций А1 = А1(ж3), А2 = А2(ж3), = ^1(ж3), № = (ж3) для слоя толщины Н = 1. В серии расчетов плотность выбиралась постоянной р = 1, параметр, характеризующий время релаксации, полагался т = 0.1.
Итерационный процесс (19)-(23) был применен сначала к задаче (10), (12) для восстановления функций = ^1(ж3), = ^2(ж3) в наборе точек, затем к задаче (11), (13) для определения функций А1 = А1(ж3), А2 = А2(ж3).
Необходимо отметить, что первый этап итерационного процесса требует знания начального приближения. В настоящей работе начальное приближение Со(х,гк) строилось в классе линейных комплексных ограниченных функций, коэффициенты которых находились из условия минимума функционала невязки на построенном из априорной информации компактном множестве:
К2
Ф = ! Щ1,к) - ¡(к)|2йк,
(24)
К1
где и(1, к) — трансформанта функции смещения на поверхности слоя при некотором законе изменения комплексного модуля. Минимизация функционала осуществлялась на равномерной сетке области [0,1] х [^ ,к2].
Далее приведены примеры восстановления функций различного характера.
Пример 1. Случай монотонных законов изменения неоднородных характеристик слоя
»1(х3) = 1 + 0.4х3, »2(х3) = 3 - 1.5х3, А1(х3) = 2 - х3, А2(х3) = 2.1 + х\.
Рис. 1.
Для первого этапа начальные приближения найдены в виде »0(х3) = 0.45 + 0.85х3 и »0(х3) = -1.6х3 + 3.2. В ходе анализа амплитудно-частотной характеристики (АЧХ) вещественной и мнимой частей трансформанты смещения при х3 = Н (рис. 1 (линией обозначена АЧХ, соответствующая точным функциям, точками — начальному приближению)), позволяющего выявить наиболее эффективные с точки зрения идентификации диапазоны изменения волновых чисел (между резонансными числами), для восстановления был выбран отрезок к £ [9.3,10.5]. Выход из итерационного процесса во всех экспериментах производился либо по числу итераций (Ж = 20), либо по условию (22) при е = 10-3. В приведенном случае выход осуществлен в соответствии со вторым условием, для выполнения которого потребовалось 11 итераций (рис. 2). Относительная погрешность для первого этапа была меньше 8%.
На рис. 2 представлены графики точного решения — сплошная линия, начального приближения — пунктир, восстановленной функции — квадратиками.
Рис. 2.
Для 2-го этапа были определены начальные приближения в виде
А0(ж3) + 2^?(ж3) = 2.9 + 0.6 ж3, А2 Ы + = -1.7 ж3 + 8.2.
Идентификация производилась на отрезке к £ [7.5,11.5], потребовалось 20 итераций, на рис. 3 представлены результаты восстановления исходных функции А 1(ж3), А2(ж3). При этом относительная погрешность восстановления второй функции не превышает 10%.
Рис. 3.
Пример 2. В случае немонотонных законов изменения (ж3) = 2—еоз(5ж3), ^2(ж3) = 3 + 8т(3ж3), А1(ж3) = 2.5 — 8т(5ж3), А2(ж3) = 5 — 2зш(5ж3 — 1) определялись начальные приближения, частотные диапазоны, ^°(ж3) = 2, ^°(ж3) = 3.5, к £ [5,8.5] и А^(ж3) + 2^°(ж3) = 6, А°(ж3) + 2^2(ж3) = 12, к £ [6.1,15.74]. Выход из итерационного процесса был произведен, соответственно через 15 и 30 итераций. Погрешность составила не более 12%. Результаты реконструкции приведены на рис. 4-5.
Рис. 4.
Рис. 5.
Рассмотрим влияние зашумления дополнительной информации на погрешность восстановления, что моделирует реальный эксперимент с установкой датчиков, измеряющих смещение. Каждый из таких датчиков снимает информацию с некоторой погрешностью, которая моделируется аддитивным зашумлением. Реальное зашумление будем моделировать с помощью соотношения / (к) = / (к)(1 + 57), где 5 — величина зашумления, 7 — случайная величина с равномерным законом распределения.
Далее рассмотрен эксперимент на зашумление дополнительной информации в 1% (5 = 0.01) в примере 1 (рис. 6-9).
Первый этап. Относительная погрешность для первого этапа составила менее 10%.
Второй этап. Относительная погрешность для второго этапа составила менее 12%.
Рис. 6.
Рис. 7.
Из приведенных результатов можно делать вывод о работоспособности предлагаемого подхода в условиях реального зашумления входной информации.
6. Заключение
Представлен способ исследования обратной коэффициентной задачи о восстановлении функций, характеризующих неоднородные свойства вязкоупругого слоя при анализе плоской деформации. С помощью интегрального преобразования Фурье задача сведена к двум более простым задачам относительно усредненных характеристик слоя. Решение построено с помощью итерационного процесса, сочетающего на каждом шаге решение прямой задачи (сведенной к интегральному уравнению Фредгольма 2-го рода) и определение поправок на основе решения интегрального уравнения Фредгольма 1-го рода с суммируемым (гладким) ядром. Рассмотрены примеры идентификации вязкоупругих
характеристик неоднородных материалов для различных законов изменения комплексного модуля по толщине слоя. В результате вычислительных экспериментов установлено, что восстановление неизвестных законов изменения параметров слоя происходит с погрешностью менее 8% при условно точной входной информации, а при зашумле-нии порядка 1% точность реконструкции составляет 12%, что позволяет заключить, что предложенный в работе [14] подход, примененный к одномерным обратным задачам для упругих стержней, может быть успешно применен и для областей типа слоя, свойства которого меняются по толщине при наличии диссипации. Проведенные вычислительные эксперименты реконструкции неизвестных функций по амплитудно-частотным характеристикам свидетельствуют об эффективности предложенного подхода для монотонных функций.
Замечание. Следует отметить, что предложенный алгоритм восстановления свойств материала существенно опирается на возбуждение колебаний двухкомпонентной нагрузкой. В случае действия лишь нормальной нагрузки при сведении к более простым задачам одна из задач получается однородной и имеет лишь тривиальное решение, из нее восстановить вторую неизвестную функцию не удастся.
Литература
1. Кристенсен Р. Введение в механику композитов.—М.: Мир, 1974.—338 с.
2. Gladwell G. M. L. Inverse problems in vibration.—Waterloo: Kluwer Academic Publishers, 2004.—471 с.
3. Ватульян А. О. Обратные задачи в механике деформируемого твердого тела.—М.: Физматлит, 2007.—223 с.
4. Ватульян А. О. О вариационном подходе при исследовании обратных коэффициентных задач в теории упругости // Владикавк. мат. журн.—2009.—Т. 11, вып. 1.—С. 3-8.
5. Kim S., Kreider K. L. Parameter identification for nonlinear elastic and viscoelastic plates // Applied Numerical Math.—2006.—Vol. 56.—P. 1538-1554.
6. Ватульян А. О. К теории обратных коэффициентных задач в линейной механике деформируемого тела // Прикладная математика и механика.—2010.—№ 6.—С. 911-918.
7. Jadamba B., Khan A. A., Raciti F. On the inverse problem of identifying Lame coefficients in linear elasticity // J. Computers and Math. with Appl.—2008.—Vol. 56.—P. 431-443.
8. McLaughlin J., Yoon J.-R. Unique identifiability of elastic parameters from time-dependent interior displacement measurement // Inverse Problems.—2004.—Vol. 20.—P. 25-45.
9. Zhiming Chen, Zou J. An augmented Lagrangian method for identifying discontinuous parameters in elliptic systems // J. Control and Optimization.—1999.—Vol. 37.—P. 892-910.
10. Тихонов А. Н., Арсенин В. Я. Методы решения некорректных задач.—М.: Наука, 1979.—288 с.
11. Ватульян А. О., Шевцова М. С. Об идентификации свойств неоднородных вязкоупругих материалов // Механика композиционных материалов и конструкций.—2009.—Т. 15, № 4.—С. 475-485.
12. Ватульян А. О. Об итерационных процессах в коэффициентных обратных задачах // Исслед. по мат. анализу, дифференц. уравнениям и их приложениям.—Владикавказ: ВНЦ РАН и РСО-А, 2010.—С. 20-32.—(Итоги науки. Юг России. Мат. форум. Т. 4).
13. Баранов И. В., Ватульян А. О., Соловьев А. И. Об одном численном алгоритме и его применении в обратных задачах идентификации упругих сред // Вычислительные технологии.—2008.—№ 3.— С. 14-26.
14. Бочарова О. В., Ватульян А. О. О реконструкции плотности и модуля Юнга для неоднородного стержня // Акустический журнал.—2009.—Т. 55, № 3.—С. 281-288.
15. Ватульян А. О., Сатуновский П. С. Об определении упругих модулей при анализе колебаний неоднородного слоя // Докл. РАН.—2007.—Т. 414, № 1.—С. 36-38.
16. Ara'ujo A. L., Mota Soares C. M., Mota Soares C. A., Herskovits J. Inverse estimation of elastic, viscoelastic and piezoelectric properties of anisotropic sandwich adaptive structures // 7th Euromech Solid Mechanics Conference.—2009.—Vol. 20.—P. 25-45.
17. Polansky J., Boiron O., Novacek V. Identification of viscoelastic properties of artificial materials simulating vascular wall // J. Computer Methods in Biomechanics and Biomedical Engineering.—2005.— Vol. 8.—P. 10-41.
18. McLaughlin J., Ji L. Recovery of the Lam'e parameter ^ in biological tissues // Inverse Problems.— 2004.—Vol. 20.—P. 1-24.
19. Аникина Т. А., Богачев И. В., Ватульян А. О. Об определении неоднородных реологических свойств балок // Вестн. ДГТУ.—2011.—Т. 10, № 7.—С. 1016-1023.
20. Колмогоров А. Н., Фомин С. В. Элементы теории функций и функционального анализа.—М.: Наука, 1976.—543 с.
Статья поступила 26 апреля 2011 г.
Богачев Иван Викторович
Южный федеральный университет,
аспирант кафедры теории упругости
РОССИЯ, 344090, Ростов-на-Дону, ул. Мильчакова, 8 а
E-mail: bogachev89@yandex.ru
Ватульян Александр Ованесович Южный математический институт ВНЦ РАН и РСО-А, заведующий отделом дифференциальных уравнений РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22; Южный федеральный университет, заведующий кафедрой теории упругости РОССИЯ, 344090, Ростов-на-Дону, ул. Мильчакова, 8 а E-mail: vatulyan@math.rsu.ru
THE INVERSE COEFFICIENT PROBLEM FOR DISSIPATIVE OPERATORS AND IDENTIFICATION OF THE PROPERTIES OF VISCOELASTIC MATERIALS
Bogachev I. V. and Vatulyan A. O.
We give some general formulation and the variational equation of the inverse problem of identifying the inhomogeneous characteristics of three-dimensional viscoelastic body. Under consideration is the problem of reconstruction of the functional coefficients of dissipative operators arising in solving several problems of identification of the properties of layered inhomogeneous viscoelastic structures in the analysis of spectral characteristics. We suggest a method for constructing an iterative process and present the results of recovering functions of different types.
Key words: inverse problems, dissipative operators, viscoelasticity, identification, iterative process.
