Научный отдел
МЕХАНИКА
Известия Саратовского университета. Новая серия. Серия: Математика. Механика. Информатика. 2022. Т. 22, вып. 2. С. 180-195 Izvestiya of Saratov University. Mathematics. Mechanics. Informatics, 2022, vol. 22, iss. 2, pp. 180-195 https://mmi.sgu.ru
https://doi.org/10.18500/1816-9791-2022-22-2-180-195
Научная статья УДК 539.3
Решение обратной задачи об идентификации двух термомеханических характеристик функционально-градиентного стержня
.20
А. О. Ватульян1,2, С. А. Нестеров
1 Южный федеральный университет, Институт математики, механики и компьютерных наук им. И. И. Воровича, Россия, 344090, г. Ростов-на-Дону, ул. Мильчакова, д. 8-А
2Южный математический институт - филиал Владикавказского научного центра Российской академии наук, Россия, 362025, г. Владикавказ, ул. Ватутина, д. 53
Ватульян Александр Ованесович, доктор физико-математических наук, 1 заведующий кафедрой теории упругости, 2заведующий отделом дифференциальных уравнений, [email protected], ИПрБ: //orcid.org/0000-0003-0444-4496, АиШогГО: 3469 Нестеров Сергей Анатольевич, кандидат физико-математических наук, старший научный сотрудник отдела дифференциальных уравнений, [email protected], https://orcid.org/0000-0003-3780-5104, АиШогГО: 808934
Аннотация. Представлен подход к решению обратной задачи об одновременной идентификации двух термомеханических характеристик функционально-градиентного стержня. Рассмотрены две задачи термоупругости с различной тепловой нагрузкой на торцах стержня. В качестве входной информации выступают данные измерений температуры на торце стержня на конечном временном интервале. Прямые задачи после применения преобразования Лапласа решены на основе метода пристрелки и обращения трансформант на основе разложения оригинала в ряд по смещенным многочленам Лежандра. Проведен анализ влияния законов изменения переменных характеристик на значения снимаемой в эксперименте входной информации. Решение нелинейной обратной задачи построено на основе итерационного процесса. Начальное приближение для итерационного процесса находится в классе линейных функций, коэффициенты которых определяют из условия минимального значения функционала невязки. Для нахождения поправок к законам
изменения термомеханических характеристик на основе слабой постановки каждой прямой задачи и метода линеаризации получена система интегральных уравнений Фредгольма 1-го рода. Регуляризация системы интегральных уравнений осуществлена на основе метода А. Н. Тихонова. Проведены вычислительные эксперименты по одновременной реконструкции двух теплофизических характеристик при известных законах изменения остальных характеристик. Восстанавливались пары как монотонно возрастающих, так и монотонно убывающих функций. Ключевые слова: функционально-градиентный стержень, коэффициент теплопроводности, удельная теплоемкость, коэффициент температурного напряжения, метод пристрелки, идентификация, коэффициентная обратная задача термоупругости, система интегральных уравнений Фредгольма 1-го рода
Благодарности: Работа выполнена при поддержке Южного математического института — филиала Владикавказского научного центра РАН.
Для цитирования: Ватульян А. О., Нестеров С. А. Решение обратной задачи об идентификации двух термомеханических характеристик функционально-градиентного стержня // Известия Саратовского университета. Новая серия. Серия: Математика. Механика. Информатика. 2022. Т. 22, вып. 2. С. 180-195. https://doi.org/10.18500/1816-9791-2022-22-2-180-195 Статья опубликована на условиях лицензии Creative Commons Attribution 4.0 International (CC-BY 4.0)
Article
Solution of the inverse problem of two thermomechanical characteristics identification of a functionally graded rod
A. O. Vatulyan1, S. A. Nesterov20
1Southern Federal University, Institute of Mathematics, Mechanics and Computer Sciences named after I. I. Vorovich, 8-A Milchakova St., Rostov-on-Don 344090, Russia
2Southern Mathematical Institute, Vladikavkaz Scientific Center of Russian Academy of Sciences, 53 Vatutina St., Vladikavkaz 362025, Russia
Alexander O. Vatulyan, [email protected], https://orcid.org/0000-0003-0444-4496, AuthorID: 3469 Sergey A. Nesterov, [email protected], https://orcid.org/0000-0003-3780-5104, AuthorID: 808934
Abstract. An approach to solving the inverse problem of the simultaneous identification of two thermomechanical characteristics of a functionally graded rod is presented. Two problems of thermoelasticity with different heat loads at the ends of the rod are considered. The input information is the temperature measurement data at the end of the rod over a finite time interval. Direct problems after applying the Laplace transform are solved on the basis of the shooting method and inversion of transformants based on the expanding the original in a series in the shifted Legendre polynomials. The analysis of the influence of the variable characteristics change laws on the values of the input information taken in the experiment is carried out. The solution to the nonlinear inverse problem is based on an iterative process. The initial approximation for the iterative process is in the class of linear functions, the coefficients of which are determined from the condition of the minimum value of the residual functional. To find corrections to the laws of change in thermomechanical characteristics on the basis of a weak statement of each direct problem and the linearization method, a system of Fredholm integral equations of the 1st kind is obtained. The system of integral equations is regularized based on the method of A. N. Tikhonov. Computational experiments on the simultaneous reconstruction of two thermophysical characteristics with known laws of change in other characteristics are carried out. Pairs of both monotonically increasing and monotonically decreasing functions were reconstructed.
Keywords: functionally graded rod, heat conduction coefficient, specific heat capacity, thermal stress coefficient, shooting method, identification, coefficient inverse problem of thermoelasticity, system of Fredholm integral equations of the 1st kind
Acknowledgements: This work was supported by the Southern Mathematical Institute, a branch
of the Vladikavkaz Scientific Center of the Russian Academy of Sciences.
For citation: Vatulyan A. O., Nesterov S. A. Solution of the inverse problem of two thermo-
mechanical characteristics identification of a functionally graded rod. Izvestiya of Saratov
University. Mathematics. Mechanics. Informatics, 2022, vol. 22, iss. 2, pp. 180-195 (in Russian).
https://doi.org/10.18500/1816-9791-2022-22-2-180-195
This is an open access article distributed under the terms of Creative Commons Attribution 4.0 International License (CC-BY 4.0)
Введение
Обычно расчеты, связанные с нахождением напряженно-деформированного состояния стержневых конструкций, находящихся в условиях высокотемпературного нагружения, проводят для однородных материалов. Однако в настоящее время в различных областях техники все шире стали использовать функционально-градиентные материалы (FGM) — неоднородные материалы с заданными свойствами [1]. В этих материалах, в отличие от слоистых композитов, удается избежать скачка термомеханических характеристик через поверхность раздела. Из-за сложной и многоступенчатой технологии изготовления FGM в конечном изделии могут присутствовать отклонения от заданных законов неоднородности. Поэтому требуется идентификация реальных свойств изделия после изготовления. Однако из-за того, что термомеханические характеристики являются функциями координат, они могут быть определены только неразрушающими методами контроля, теоретически опирающимися на аппарат коэффициентных обратных задач (КОЗ) термомеханики.
КОЗ термоупругости [2] — это задачи нахождения термомеханических характеристик (модулей упругости, плотности, коэффициента теплопроводности, удельной теплоемкости, коэффициента температурного напряжения) как коэффициентов дифференциальных операторов уравнений термоупругости по дополнительной информации о полях перемещений или температуры, измеренной на части границы объекта, на некотором временном отрезке.
Из-за некорректности и нелинейности КОЗ в настоящее время актуальной задачей является построение экономных и устойчивых алгоритмов их решения. Наиболее распространенным методом решения КОЗ является построение функционала невязки и его минимизация каким-либо из градиентных методов [3, 4] или генетическим алгоритмом [5]. На основе такого подхода выполнено довольно большое количество исследований КОЗ теплопроводности [3-8] и теории упругости [9-11]. Однако для ряда материалов необходимо учитывать связанность полей и решать КОЗ термоупругости [12]. При решении КОЗ термоупругости возникает проблема построения операторных уравнений, связывающих искомые и измеряемые в эксперименте функции.
Ранее авторами был разработан и апробирован численный метод решения одномерных нелинейных КОЗ теории упругости [13] и термоупругости [12, 14, 15], основанный на построении итерационного процесса. Для этого на основе слабой постановки прямой задачи и линеаризации была получена система линеаризованных операторных уравнений в трансформантах Лапласа [12]. Ограничившись восстановлением только одной из термомеханических характеристик при известных других, после перехода
от трансформант Лапласа к оригиналам была проведена идентификация термомеханических характеристик стержня [14] и цилиндра [15]. Начальное приближение определялось в классе линейных функций, коэффициенты которых определяют из условия минимального значения функционала невязки. Для нахождения поправок термомеханических характеристик на каждом шаге итерационного процесса решалось интегральное уравнение Фредгольма 1-го рода.
Однако обычно на практике неизвестными являются несколько термомеханических характеристик. Поэтому обратная задача о реконструкции нескольких термомеханических характеристик неоднородных тел является актуальной.
В данной работе рассматривается КОЗ об одновременной идентификации двух термомеханических характеристик функционально-градиентного стержня. Рассмотрены две постановки задач термоупругости с различной тепловой нагрузкой на торцах стержня. Прямые задачи после применения преобразования Лапласа решены с использованием метода пристрелки и обращения трансформант на основе метода разложения оригинала в ряд по смещенным многочленам Лежандра. Решение обратных задач построено на основе итерационного процесса, на каждом шаге которого численно решалась система интегральных уравнений Фредгольма 1-го рода. В ходе вычислительных экспериментов восстанавливались пары как монотонно возрастающих, так и монотонно убывающих функций.
1. Постановка обратной задачи термоупругости для функционально-градиентного стержня
Рассмотрим КОЗ об определении законов изменения двух теплофизических характеристик функционально-градиентного стержня длины I. Торец стержня х = 0 жестко закреплен и на нем поддерживается нулевая температура, на другом торце х = I, свободном от напряжений, действует тепловой поток. В рамках тепловой диагностики идентификация переменных свойств может быть осуществлена по данным о температуре, измеренной на торце стержня, на конечном временном интервале. Поскольку неизвестными являются две теплофизические характеристики, то рассмотрим две задачи термоупругости с разной тепловой нагрузкой на торце х = /.В первой задаче на торце действует тепловой поток д01 ево второй задаче — тепловой поток д0Н(¿), где Н(¿) — функция Хевисайда.
Постановка задачи 1 имеет вид
дох . . д2и г ,, г
ах = р(х)д^, х е [0,^ ге [0, (x), д и
ах = Е(х) — - ^(х)в, д д д д2 и
дх(к(х)дх) =С£(х) т +То1 (х)х е [0'1], и(0, г) = 0, ах(1, г) = 0, д
9(0, Ь) = 0, -к(1)^(/) = qotе-1, х
ди
6(х, 0) = и(х, 0) = — (х, 0) = 0.
г е [0, то),
(1) (2)
(3)
(4)
(5)
(6)
Постановка задачи 2 совпадает с постановкой задачи 1 за исключением тепловых граничных условий (5), вместо которых ставятся условия
0(0, ¿) = 0,
- к(1) £ (1) = доН (г).
(7)
у(г) = ^,
I \ / 70 >
Г =
£о = Ц •
Перейдем в (1)-(7) к безразмерным параметрам и переменным, обозначая [14]:
* = ?, ВД = , Ф) = ^, Р(*) = ^, Е(г) = Щ?) ^ - ^ ь = ^ г = £, Ш = ^, и = |, П = ^ 8о = Ц, ^ = Йй,
Здесь Е0, р0, к0, с0, 70 — некоторые характерные величины.
В дальнейшем в обезразмеренных формулах обозначим функции с индексами «I» и «II» как функции, отвечающие первой и второй постановке задач соответственно. Тогда задача 1 в безразмерной форме имеет вид
дг
д ш
д П/ 2 д2П/ = ?2Р
д 2
ге [0,1], те [0, то),
дП
П/ = Е(г) длт - 7(¿)Ш1,
дШ/ д2П/
ФО^- + 6о7(?)
д
П/(0, т) = 0,
д
П/(1, т) дШ/
ге [0,1], те [0,то),
= 0,
Ш/(0, т) = 0, -к(1)-^-(1, т)=(3те1
Ш/(г, 0) = П/(г, 0)
дП/
( , 0) = 0.
(8) (9)
(10) (11) (12) (13)
Обезразмеренная постановка задачи 2 отличается от задачи 1 только тепловыми граничными условиями. Вместо граничных условий (12) имеем
дШ/
Шп(0, т) = 0, -к(1)^(1, т)=РН(т).
д
(14)
Множитель е0, стоящий в правой части уравнения движения (8), очень мал (10-16..10-12), поэтому в дальнейшем для упрощения расчетов будем рассматривать квазистатические постановки прямых задач термоупругости, не учитывающих инерционные члены, т.е. ^ = 0, ^^ = 0.
Прямая задача термоупругости состоит в нахождении функций Ш/(х,т), П/(х,т) из (8)—( 13) по известным термомеханическим характеристикам (Е(г), р(г), с(г), к(г),
7(*)).
Особенности идентификации механических характеристик неоднородных тел исследованы в [13]. Мы же подробно рассмотрим особенности идентификации тепло-физических характеристик.
В обратной задаче требуется одновременно восстановить две теплофизические характеристики при известных остальных по дополнительной информации о торцевой температуре /( ), / ( ), измеренной в ходе проведения двух экспериментов с тепловой нагрузкой, моделируемой граничными условиями (12) и (14) соответственно:
Ш/ (1, т) = // (т), те [а 1, 61], Ш// (1, т) = /// (г), те [«2, 62].
(15)
(16)
* т)
р). (19)
2. Решение прямой задачи термоупругости
При произвольных законах изменения термомеханических характеристик прямые задачи вида (8)—( 13) могут быть решены лишь численно, например методом пристрелки. Применим его к решению поставленных прямых задач. Для этого применим преобразование Лапласа к уравнениям равновесия (8), теплопроводности (10), граничным условиям (11), (12) с учетом начальных условий (13) и того, что £1(г,р) = 0, после некоторых преобразований перейдем к канонической системе обыкновенных дифференциальных уравнений 1-го порядка с переменными коэффициентами:
^ = Ш«' (*. Р). (17)
аг к(г)
чг=" (к«+* Ш>и/' , ^ (|8)
аи1 = кМ_
Ах Е(г)
Далее для 1-й задачи методом Рунге-Кутта решается вспомогательная задача Коши, состоящая из системы дифференциальных уравнений (17)-(19) со следующими условиями на торце г = 0: £Л(0) = 0, Й^(0) = 0, (1(0) = 1. Составляются выражения III = а11!1, Ш = а1У°1, (/ = а-^Я1. Неизвестная константа а1 определяется путем удовлетворения условию на торце г = 1: Я(1, р) = .
Для 2-й задачи решается вспомогательная задача Коши, состоящая из системы дифференциальных уравнений (17)-(19) со следующими условиями на торце г = 0: и2(0) = 0, 1°2(0) = 1, (2(0) = 0. Составляются выражения ип = а2и2, = <^2Уй2, Яп = а2(°2. Неизвестная константа а2 определяется путем удовлетворения условию на торце г = 1: Ян(1, р) = |.
После нахождения решений прямых задач в трансформантах для каждого значения параметра р = 1, 2.... для нахождения оригиналов функций в работе применялось обращение преобразования Лапласа на основе метода разложения оригинала в ряд по смещенным многочленам Лежандра. На примере однородного стержня проведено сравнение результатов, полученных предложенным методом решения прямой задачи с аналитическим решением. Выяснено, что для того чтобы погрешность вычисления оригиналов температуры для времени т ^ 0.01 не превосходила 1%, достаточно ограничиться 50 членами ряда в разложении по смещенным многочленам Лежандра.
Важным этапом при решении обратных задач является анализ влияния законов изменения переменных характеристик на значения снимаемой в эксперименте входной информации.
Проведен анализ влияния законов изменения теплофизических характеристик и параметра связанности 80 на температуру, измеренную на торце стержня г = 1, которая выступает при решении обратных задач в качестве входной информации. Выяснено, что законы изменения коэффициента теплопроводности и удельной теплоемкости для рассматриваемых задач оказывают сильное влияние на торцевую температуру при любом параметре связанности 80 = 0.01...0.4. Влияние же коэффициента температурных напряжений на торцевую температуру сильно зависит от величины параметра связанности и оказывается существенным только при 80 ^ 0.2.
На рис. 1-3 показаны графики Ш(1, г) и Шц(1, г) для различных наборов законов изменения свойств материала при @ =1.
Так, на рис. 1, для того чтобы показать влияние коэффициента теплопроводности к (г) на торцевую температуру, были построены графики Ш/ (1, т) и Ш// (1, т) при 80 = 0.03 и двух наборах законов изменения свойств материала: 1) к(г) = 1, с(г) = 7(г) = Е(г) = 1 + 3г (сплошная линия); 2) с(г) = 7(г) = Е(г) = к(г) = 1 + 3г (пунктир).
а / a б / b
Рис. 1. Графики торцевой температуры, построенные при различных законах изменения коэффициента теплопроводности: а — WI(1 ,т); б — WII(1 ,т) Fig. 1. End temperature graphs plotted for various laws of change in the heat conduction
coefficient: a — WI(1,r); b — WII(1,r)
Для оценки влияния законов изменения удельной теплоемкости с(г) аналогично на рис. 2 построены графики торцевой температуры при 0 = 0.03 и двух наборах законов изменения свойств материала: 1) с(г) = 1, к(г) = 7(%) = Е(г) = 1 + 2е-2г (сплошная линия); 2) с(г) = к (г) = 7(%) = Е (г) = 1 + 2е-2г (пунктир).
а / a б / b
Рис. 2. Графики торцевой температуры, построенные при различных законах изменения удельной теплоемкости: а — WI(1,т); б — WII(1,т) Fig. 2. End temperature graphs plotted for various laws of change in the heat capacity
change: a — Wi(1,r); b — Wii(1,r)
На рис. 3 показано влияние коэффициента температурного напряжения у(г) на торцевую температуру неоднородного стержня при большом параметре связанности 50 = 0.4. Рассмотрены следующие два набора законов изменения свойств материала: 1) у(г) = 1, к (г) = с(г) = Е (г) = 2 + ео8(3г) (сплошная линия); 2) к(г) = с(г) = Е (г) = у (г) = 2 + еов(3г) (пунктир).
а / a б / b
Рис. 3. Графики торцевой температуры, построенные при различных законах изменения коэффициента температурного напряжения: а — Wi(1 ,т); б — Wu(1 ,т) Fig. 3. End temperature graphs plotted for various laws of variation of the thermal stress
coefficient: a — WI(1,т); b — WII(1,r)
Графики торцевой температуры на рис. 2 для различных наборов законов удельной теплоемкости и на рис. 3 для различных наборов законов коэффициента температурных напряжений отличаются друг от друга незначительно, что делает невозможным одновременную реконструкцию удельной теплоемкости и коэффициента температурных напряжений.
3. Решение обратной задачи термоупругости
КОЗ термоупругости является нелинейной задачей. Построим ее решение на основе итерационного процесса, воспользовавшись техникой, описанной в работе [12]. Записав для задач 1 и 2 слабые постановки в трансформантах Лапласа и используя метод линеаризации для нахождения поправок теплофизических характеристик ôa(n~1) (ôk(n~1), 5с(п~1), 5y(n~1)) на (п — 1)-й итерации, была получена система двух операторных уравнений 1-го рода:
1 i (п 1) \ 2 1
(p+1)2J Sk(n- 1М j dz+p(p+1)2 J ôc(n~ 1)( w(n~ 1))2dz+
П v / П
1 пи (п~ 1) +5пр(р + 1)2 / Sy(n~1)^-™(n"^ - «'^ ™(n"1)'
-Wf" 1)dz = 3 (fi (p) — W)n" 1)(1,р)),
(20)
р J s k(n~ ^f <Wr~) dz + p2 j
ПП
+Snp2
dz
ЙТ~Т(n"1)
5y(n"1) ^^W,
dz+ p2 Sc(n~ 1)( Win~ 1])2dz+
(n-1)
ii
dz = 3(fii (p) — Win" 1)(1, p)).
(21)
1
Исходя из анализа влияния законов изменения теплофизических характеристик на торцевую температуру, рассмотрено два различных случая одновременной идентификации двух теплофизических характеристик.
1. Пусть 7(г) известна. Необходимо определить к (г) и с(г). В этом случае, полагая в (20), (21) б7(п-1) = 0, получим
1 ( (п 1) \ 2 1
(р + 1)2 ! 5к(п-1)( ) йг + р(р + 1)2 ! 5с(п-1){\¥(1п-1))2йг =
00
= Р(I/(Р) -Шп-1)(1,Р)), (22)
1 // 2 1 ( п- 1) 2
Р У $к{п-1)(+ р2 у 6с(-п-1)(\¥/п-1))2йг = Р(¡//(р) -Ш%-1)(1,р)). (23) 00
Обращая операторные уравнения в трансформантах (22), (23), получим систему операторных уравнений в оригиналах, которая является системой интегральных уравнений Фредгольма (ИУФ) 1-го рода:
1
(дк(п-1)Кп(г, т) + д с(п-1)Я12(г, т))йг = ¡/ (т) -Ш/п-1)(1, т), те К, Ь], (24)
j (5к(п-1)П21(г, т) + 6¿п-1)П22(г, т))йг = ¡П(т) - ш/п-1)(1, т), т е [а.2, Ь]. (25) 0
Здесь ядра (24), (25) имеют вид
0
1
К (7 )=1 [ (дШ(п-1)(г, п) + 2д2Ш/п~1\г, п) + д 3Ш[п-1){г, п)\
К11(", Т) = ^ [ дх дхдт + дхдт2 Х
0
дШ}п-1)(г1 т - п) д
х ——--^^-,
, „=1 / , Ч)++х
0
дШ(гп-1)(г, г-п) ,
Х---^-¿Т1,
д
К =1 [ д2Ш/п-1)(г,п) дШ/п-1 (г- п)
К21(Х, Т) = ^ дЖт дг ¿Ти
0
К (гт)=1 / дШ/Г%, П) дШ/Т%,Т - ¿т 0
2. Пусть с(г) известна. Необходимо определить к (г) и у(г). В этом случае, полагая в (20), (21) ёс(п-1) = 0, после обращения трансформант получим систему операторных уравнений в оригиналах для нахождения поправок
8к(п-1), 67(га-1):
(бк(п-1) Мп(г, т) + 8у(п-1)М12(г, т))йг = Ц (т) -Ш(п-1)(1, т), те [а 1, Ь], (26)
(5к(п-1)М21(г, т) + 5у(п-1)М22(г, т))<1г = ¡п(т) - Ш%-1)(1, т), т е [а2, Ь2]. (27)
Здесь ядра (26), (27) имеют вид
Мц(г, т) = Пц(г, т), М21 (г, т) = Д21 , т),
М,2,= | / (шТ"(,, П) + 2+ • ^
о ^
)
х
х
д2и(п-1)(г,т - п)
дхдт
т(п 1)
¿Тл
5о ГдШ(/;-1)(г, Т1)д2и(п-1)(г ,т - п)
М22(г, т) =
дт
дхдт
1.
На основе системы операторных уравнений вида (24)-(27) можно организовать итерационный процесс для определения искомых характеристик.
Итерационный процесс стартует с выбора начальных приближений двух функций, которые определяются в классе линейных функций а10) = т1 + т2г, а20) = т3 + т4г на основе минимизации функционала невязки
3 = (¡1 (г) -Ш.I (1, т))2^ + (¡П (т) -Шп (1, т))2^,
(28)
«1 «2 на двух множествах А1 и А2, построенных по априорной информации об ограниченности искомых функций 0 < а- ^ а^г) ^ а+, 0 < а- ^ а2 (г) ^ а+, г е [0,1]. Здесь Ш. (1,т) и Шц(1,т) — решения задач 1 и 2 в точке г =1 при заданных законах а!0), а20), где параметры т1, т2, т3, т4 выбраны из соответствующих множеств А1 и А2. Численная реализация этого подхода осуществляется путем замены интегральных операторов в правой части (28) на конечные суммы, вычислением их значений для каждого набора точек на сетках разбиений множеств А1 и А2 и выбором точек, на котором достигается наименьшее значение функционала 3.
Далее осуществляется корректировка текущих законов изменения термомеханиче-
г— (п— 1) г-(п-1)
скои характеристики путем нахождения двух поправок оа\ , оа'2 из решения соответствующих систем ИУФ 1-го рода. После нахождения поправок строятся новые приближения а!_п) = аЩ-1) + баЩ-1), а2п) = о^-1) +6 о^-1). Далее процесс замыкается и повторяется до тех пор, пока не будет достигнуто условие остановки: достижения некоторого порогового значения функционала невязки (28), равного 10-4, или достижения предельного количества итераций, равного 20.
1
1
т
При этом следует отметить, что основная сложность реализации этого процесса состоит в построении численного решения систем ИУФ 1-го рода. Поскольку решение системы ИУФ 1-го рода является некорректной задачей, то для его регуляризации применяется метод А. Н. Тихонова [16]. Для этого вводится равномерное разбиение отрезка [0,1] на n отрезков точками Zj = Az (j — 1), j = 1...n, где Az = nh — шаг разбиения по безразмерной координате и равномерное разбиение безразмерных временных отрезков [a1, Ь1] и [а2, Ь2] на m отрезков точками щ = AT\(i — 1),r2i = Ат2(% — 1), г = 1...m, где Ат\ = , Ат2 = ^zf — шаги разбиения по безразмерному времени.
Согласно методу Тихонова решение обратной задачи сводится к решению регуля-ризованной системы уравнений:
Вау = Fa, (29)
где Ва = В + аС, В = АТА, Fa = ATb, Ат — транспонированная матрица А, С = Е + Со, Е — единичная матрица, а — параметр регуляризации, который выбирается по обобщенной невязке согласно [16].
При этом компоненты векторов у, b и матрицы А имеют вид
j
а1( ^), а2(Zj—n) ,
1 ... п,
п + 1...2п,
bi
{
h (n),
fii (Ti-m).,
1 ... m,
m + 1...2m,
А
Rll( Zj, n),
Rn(Zj-n, n), R21(Zj , Ti-m), R22(Zj-n, Ti—m), г = m +1...2m, j = n + 1...2n,
= 1... m, = 1 ... п, i = 1...m, j = n + 1...2n, i = m + 1...2m, j = 1...n,
С0 — матрица размером 2n x 2n, которая имеет вид
( __L_
' Л ^2
Az^
0.
0
__^ __L_ 0
Az2 Л -2 л -2..........W
Az 2
V
Az^ AZ2 'AZ2
A ,
Az2 /
4. Результаты вычислительных экспериментов
В данном разделе представлены результаты вычислительных экспериментов по одновременной реконструкции двух теплофизических характеристик в классах степенных, логарифмических, тригонометрических и экспоненциальных функций. В вычислениях принято @ = 1, п = 31.
Найдены наиболее информативные интервалы для измерения входной информации. В случае тепловой торцевой нагрузки, соответствующей 1-й постановке задачи, наиболее информативным является отрезок [а-\^, Ь^ = [0.04,0.72], а при 2-й постановке задачи — отрезок [а2, Ь2] = [0.03,0.48]. Съем дополнительной информации происходил в 11 точках (параметр т =11) внутри выбранных информативных отрезков.
На рис. 4-7 представлены результаты реконструкции пар теплофизических характеристик (с(г) и к(г)), ( к (г) и 7(гУ). При этом сплошной линией изображен точный закон неоднородности, точками — восстановленный.
В первой серии экспериментов восстанавливалась пара функций ( ), ( ) при известных 7(г) = Е(г) = 1 + 2г и 50 = 0.03.
Сначала были проведены эксперименты по одновременной реконструкции двух характеристик, имеющих монотонно возрастающие законы неоднородности: к(г) = = 1 - 0.5(г - 1)2, с(г) = 0.5(1 + г2) (рис. 4). Начальные приближения коэффициентов:
l
к0(г) = 0.45) + 0.5)г, со(г) = 0.5) + 0.45)г. Потребовалось 16 итераций; максимальная погрешность восстановления двух функций не превысила 6%.
а / a б / b
Рис. 4. Результаты реконструкции возрастающих функций: а — k(z) = 1 — 0.5(2 — 1)2;
б — ф) = 0.5(1 + z2) Fig. 4. Results of reconstruction of increasing functions: a — k(z) = 1 — 0.5(z — 1)2;
b — c(z) = 0.5(1 + z2)
Затем были проведены эксперименты по одновременной реконструкции двух характеристик, имеющих монотонно убывающие законы неоднородности: k(z) = 1 + + 0.5cos(3z), c(z) = 1.2cos(z) (рис. 5). Начальные приближения коэффициентов: k0(z) = 1.45 — 0.9z, с0(z) = 1.25 — 0.6z. Потребовалось 14 итераций; максимальная погрешность восстановления двух функций не превысила 5%.
о Ш 04 С Е- ВД 2 3 ...... 0.2 0.4 0 6 0.8 z
а / a б / b
Рис. 5. Результаты реконструкции убывающих функций: а — k(z) = 1 + 0.5 cos(3z);
б — c(z) = cos(z)
Fig. 5. The results of the reconstruction of decreasing functions: a — fc(z)=1+0.5cos(3z);
b — c(z) = cos(z)
Во второй серии экспериментов восстанавливалась пара функций к(г) и у(г) при известных с(г) = Е(г) = 1 + 2г. При этом одновременная реконструкция к(г) и у(г) оказалась возможной только при параметре связанности 80 ^ 0.2.
На рис. 6 показаны результаты вычислительных экспериментов по одновременной реконструкции двух теплофизических характеристик, имеющих монотонно возрастающие законы неоднородности: ВД = 1 + 0.2е2^, 7(г) = 1 + >/2г при 8о = 0.4. Начальные приближения коэффициентов: /с0(^) = 1.2 + 1.25г, 70(2) = 1 + 1.45г. Выход из итерационного процесса произошел по предельному количеству итераций. При этом на последней итерации максимальная погрешность восстановления двух функций не превысила 8%.
а / a б / b
Рис. 6. Результаты реконструкции возрастающих функций: а — k(z) = 1 + 0.2e2z;
б — 7(2) = 1 + V2z
Fig. 6. Results of reconstruction of increasing functions: a — k(z) = 1 + 0.2e2^;
b — k(*) = 1 + V2z
На рис. 7 показаны результаты вычислительных экспериментов по одновременной реконструкции двух теплофизических характеристик, имеющих монотонно убывающие законы неоднородности: 7(г) = 1 + е-4г, 7(2) = 2 — 1п(1 + 5г) при 50 = 0.4. Начальные приближения коэффициентов: 7(0)(г) = 2 — 0.9г, 70(2) = 2.1 — 1.85г.
о 0.2 Щ аз; о 0.2 о о в Ш z
а / a б / b
Рис. 7. Результаты реконструкции убывающих функций: а — k(z) = 1 + e-4z;
б — k(^) = 2 - ln(1 + 5z) Fig. 7. The results of the reconstruction of decreasing functions: a — k(z) = 1 + e-4z;
b — kCz) =2 - ln(1 + 5z)
Максимальная погрешность восстановления двух функций на 20-й итерации не превысила 7%.
Заключение
Исследована обратная задача об одновременной идентификации двух термомеханических характеристик функционально-градиентного стержня по дополнительной информации, полученной в результате проведения двух вычислительных экспериментов с различным типом торцевой нагрузки. Прямые задачи после обезразмеривания и преобразования Лапласа решались с использованием метода пристрелки и обращения трансформант на основе разложения оригинала в ряд по смещенным многочленам Лежандра. Для нахождения двух поправок термомеханических характеристик в итерационном процессе получена система интегральных уравнений Фредгольма 1-го рода. В ходе вычислительных экспериментов выяснено, что: 1) на результат реконструкции пары функций (c(z) и k(z)) величина параметра связанности 50 оказывает небольшое влияние, в то же время успешная реконструкция пары ( k(z) и j(z)) возможна только при 50 ^ 0.2; 2) пары функций c(z) и k(z) восстанавливаются с меньшей погрешностью (не превышающей 6%), чем пары k(z) и j(z).
Список литературы
1. Birman V., Byrd L. W. Modeling and analysis of functionally graded materials and structures // Applied Mechanics Reviews. 2007. Vol. 60, iss. 5. P. 195-216. https: //doi.org/10.1115/1.2777164
2. Ломазов В. А. Задачи диагностики неоднородных термоупругих сред. Орел : Изд-во ОрелГТУ, 2002. 168 с.
3. Алифанов О. М., Артюхин Е. А., Румянцев С. В. Экстремальные методы решения некорректных задач. Москва : Наука, 1988. 288 с.
4. Razzaghi H., Kowsary F., Ashjaee M. Derivation and application of the adjoint method for estimation of both spatially and temporally varying convective heat transfer coefficient // Applied Thermal Engineering. 2019. Vol. 154. P. 63-75. https://doi.org/10.1016/jj.applther maleng.2019.03.068
5. Raudensky M., Woodbary K. A., Kral J. Genetic algorithm in solution of inverse heat conduction problems // Numerical Heat Transfer, Part B : Fundamentals. 1995. Vol. 28, iss. 3. P. 293-306. https://doi.org/10.1080/10407799508928835
6. Dulikravich G. S., Reddy S. R., Pasqualette M. A., Colaco M. J., Orlande H. R., Coverston J. Inverse determination of spatially varying material coefficients in solid objects // Journal of Inverse and Ill-posed Problems. 2016. Vol. 24. P. 181-194. https: //doi.org/10.1515/jiip-2015-0057
7. Cao K., Lesnic D. Determination of space-dependent coefficients from temperature measurements using the conjugate gradient method // Numerical Methods for Partial Differential Equations. 2018. Vol. 34, iss. 4. P. 1370-1400. https://doi.org/10.1002/num. 22262
8. Helmig T., Al-Sibai F., Kneer R. Estimating sensor number and spacing for inverse calculation of thermal boundary conditions using the conjugate gradient method // International Journal of Heat and Mass Transfer. 2020. Vol. 153. Art. 119638. https: //doi.org/10.1016/j.ijheatmasstransfer.2020.119638
9. Geymonat G., Pagano S. Identification of mechanical properties by displacement field measurement: A variational approach // Meccanica. 2003. Vol. 38. P. 535-545. https: //doi.org/10.1023/A:1024766911435
10. Grediac M., Hild F., Pineau A. Full-Field Measurements and Identification in Solid Mechanics. Great Britain : Wiley-ISTE, 2013. 485 p. https://doi.org/10.1002/9781118578469
11. Avril S., Pierron F. General framework for the identification of constitutive parameters from full-field measurements in linear elasticity // International Journal of Solids and Structures. 2007. Vol. 44, iss. 14-15. P. 4978-5002. https://doi.org/10.1016/j.ijsolstr.2006.12.018
12. Ватульян А. О., Нестеров С. А. Коэффициентные обратные задачи термомеханики. Ростов-на-Дону ; Таганрог : Изд-во Южного федерального ун-та, 2019. 146 с.
13. Dudarev V. V., Vatulyan A. O, Mnukhin R. M., Nedin R. D. Concerning an approach to identifying the Lame parameters of an elastic functionally graded cylinder // Mathematical Methods in the Applied Sciences. 2020. Vol. 43, iss. 11. P. 6861-6870. https://doi.org/10. 1002/mma.6428
14. Nedin R., Nesterov S., Vatulyan A. On an inverse problem for inhomogeneous thermoelastic rod // International Journal of Solids and Structures. 2014. Vol. 51, iss. 3. P. 767-773. https://doi.org/10.1016/j.ijsolstr.2013.11.003
15. Ватульян А. О., Нестеров С. А. О задаче идентификации термомеханических характеристик конечного функционально-градиентного цилиндра // Известия Саратовского университета. Новая серия. Серия: Математика. Механика. Информатика. 2021. Т. 21, вып. 1. С. 35-47. https://doi.org/10.18500/1816-9791-2021-21-1-35-47
16. Тихонов А. Н., Гончарский А. В., Степанов В. В., Ягола А. Г. Численные методы решения некорректных задач. Москва : Наука, 1990. 230 с.
References
1. Birman V., Byrd L. W. Modeling and analysis of functionally graded materials and structures. Applied Mechanics Reviews, 2007, vol. 60, iss. 5, pp. 195-216. https://doi.org/ 10.1115/1.2777164
2. Lomazov V. A. Zadachi diagnostiki neodnorodnykh termouprugikh sred [Diagnostics Problems for Inhomogeneous Thermoelastic Media]. Orel, OrelGTU Publ., 2002. 168 p. (in Russian).
3. Alifanov O. M., Artyukhin E. A., Rumyantsev S. V. Ekstremal'nye metody resheniya nekorrektnykh zadach [Extreme Methods of Solving Ill-posed Problems]. Moscow, Nauka, 1988. 288 p. (in Russian).
4. Razzaghi H., Kowsary F., Ashjaee M. Derivation and application of the adjoint method for estimation of both spatially and temporally varying convective heat transfer coefficient. Applied Thermal Engineering, 2019, vol. 154, pp. 63-75. https://doi.org/10.1016/_j. applthermaleng.2019.03.068
5. Raudensky M., Woodbary K. A., Kral J. Genetic algorithm in solution of inverse heat conduction problems. Numerical Heat Transfer, Part B: Fundamentals, 1995, vol. 28, iss. 3, pp. 293-306. https://doi.org/10.1080/10407799508928835
6. Dulikravich G. S., Reddy S. R., Pasqualette M. A., Colaco M. J., Orlande H. R., Coverston J. Inverse determination of spatially varying material coefficients in solid objects. Journal of Inverse and Ill-posed Problems, 2016, vol. 24, pp. 181-194. https://doi.org/10.1515/ jiip-2015-0057
7. Cao K., Lesnic D. Determination of space-dependent coefficients from temperature measurements using the conjugate gradient method. Numerical Methods for Partial Differential Equations, 2018, vol. 34, iss. 4, pp. 1370-1400. https://doi.org/10.1002/num.22262
8. Helmig T., Al-Sibai F., Kneer R. Estimating sensor number and spacing for inverse calculation of thermal boundary conditions using the conjugate gradient method. International Journal of Heat and Mass Transfer, 2020, vol. 153, Art. 119638. https: //doi.org/10.1016/j.ijheatmasstransfer.2020.119638
9. Geymonat G., Pagano S. Identification of mechanical properties by displacement field measurement: A variational approach. Meccanica, 2003, vol. 38, pp. 535-545. https: //doi.org/10.1023/A:1024766911435
10. Grediac M., Hild F., Pineau A. Full-Field Measurements and Identification in Solid Me-
chanics. Great Britain, Wiley-ISTE, 2013. 485 p. https://doi.org/10.1002/9781118578469
11. Avril S., Pierron F. General framework for the identification of constitutive parameters from full-field measurements in linear elasticity. International Journal of Solids and Structures, 2007, vol. 44, iss. 14-15, pp. 4978-5002. https://doi.org/10.1016/j-.ijsolstr.2006.12.018
12. Vatulyan A. O., Nesterov S. A. Koeffitsiyentnye obratnye zadachi termomekhaniki [Coefficient Inverse Problems of Thermomechanics]. Rostov-on-Don, Taganrog, Southern Federal University Publ., 2019. 146 p. (in Russian).
13. Dudarev V. V., Vatulyan A. O., Mnukhin R. M., Nedin R. D. Concerning an approach to identifying the Lame parameters of an elastic functionally graded cylinder. Mathematical Methods in the Applied Sciences, 2020, vol. 43, iss. 11, pp. 6861-6870. https://doi.org/10. 1002/mma.6428
14. Nedin R., Nesterov S., Vatulyan A. On an inverse problem for inhomogeneous thermoelastic rod. International Journal of Solids and Structures, 2014, vol. 51, iss. 3, pp. 767-773. https://doi.org/10.1016/j-.ijsolstr.2013.11.003
15. Vatulyan A. O., Nesterov S. A. On the identification problem of the thermomechanical characteristics of the finite functionally graded cylinder. Izvestiya of Saratov University. Mathematics. Mechanics. Informatics, 2021, vol. 21, iss. 1, pp. 35-47 (in Russian). https://doi.org/10.18500/1816-9791-2021-21-1-35-47
16. Tikhonov A. N., Goncharskiy A. V., Stepanov V V., Yagola A. G. Chislennye metody resheniya nekorrektnykh zadach [Numerical Methods for Solving Ill-Posed Problems]. Moscow, Nauka, 1990. 230 p. (in Russian).
Поступила в редакцию / Received 25.12.2021 Принята к публикации / Accepted 02.02.2022 Опубликована / Published 31.05.2022