Научная статья на тему 'Проблемы идентификации неоднородных свойств термоупругой среды'

Проблемы идентификации неоднородных свойств термоупругой среды Текст научной статьи по специальности «Математика»

CC BY
78
28
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕОДНОРОДНЫЙ СТЕРЖЕНЬ / ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ / ТЕРМОУПРУГОСТЬ / INHOMOGENEOUS ROD / INTEGRAL EQUATION / THERMO-ELASTICITY

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

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

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

THE PROBLEMS OF IDENTIFICATION OF INHOMOGENEOUS PROPERTIES OF THERMO-ELASTIC MEDIA

The problem of identification of inhomogeneous properties of thermo-elastic media is considered. The iteration algorithm, based on step-by-step solution of Fredholm integral equation of the first kind is proposed. Results of the numerical experiments are given.

Текст научной работы на тему «Проблемы идентификации неоднородных свойств термоупругой среды»

Механика деформируемого твердого тела Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (4), с. 1657-1659

1657

УДК 539.3

ПРОБЛЕМЫ ИДЕНТИФИКАЦИИ НЕОДНОРОДНЫХ СВОЙСТВ ТЕРМОУПРУГОЙ СРЕДЫ

© 2011 г. С.А. Нестеров

Южный федеральный университет, Ростов-на-Дону

1079@list.ru

Поступила в редакцию 15.06.2011

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

Ключевые слова: неоднородный стержень, интегральные уравнения, термоупругость.

Введение

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

Постановка задачи

Рассмотрим закрепленный на торце х = 0 термоупругий стержень длиной I, в котором ко -лебания возбуждаются при помощи внезапного приложенного к торцу х = I теплового потока. Считаем, что модуль Юнга, удельная теплоемкость, плотность, коэффициент теплового расширения есть произвольные положительные функции ко -ординаты х. Тогда уравнения связанной термо-

упругости примут вид:

дсх , ч d2u

дх

_э'

дх

= Р( х)

dt2

k( х)

Э0

дх

с х = E (х)

д0

ди дх

— а( х)0

(1)

= c( х)р( х) — + Т0а( х) E (х)

д 2u

_ dt 0 ' dxdt'

а граничные и начальные условия представимы в форме

u(0, t) = 0(0, t) = 0,

d0 (2) — k(l) —(l) = qoH(t), с x (l, t) = 0,

дх

du

0( x,0) = u( x,0) = — (x,0) = 0. (3) дх

Будем считать, что известна дополнительная информация о торцевой температуре 0(l, t) = ft). Цель решения обратной задачи — реконструкция одного из коэффициентов дифференциальных операторов термоупругости при известных остальных.

Сведение прямой задачи к системе интегральных уравнений Фредгольма

При заданных произвольных значениях теплофизических и механических коэффициентов краевая задача (1)—(3) может быть исследована лишь численно. Для этого вначале вводятся безразмерные параметры

х k( ) k (zl) ) c (zl) z = _, k (z) =-, c (z) =-,

l k0 c0 p(z) = PM, E(z) = , а(z) = а(zl)

т =

p0

k0t

c0p0l

2

E0

W (z, т) =

k00

U =

а0

l

u

1658

С.А. Нестеров

, 5 =

а о т0 Е0 соРо

е = -

После этого к полученной краевой задаче в безразмерном виде применяют преобразование Лапласа и сводят ее к системе интегральных уравнений Фредгольма второго рода относительно трансформант безразмерной температуры Ж(г, р) и напряжения О(г, р):

1

Ж( г, р) = | КД г, р№ +

0

1

+1 К2( г, £ р)О& р^ + /1 (г, р), (4)

о

1

О(г,р) = |Кз(г,р) р)^+

о

1

■^(г, Ъ, р) р) ^. (5)

+

Здесь

*1( г, р) = - р (с (£, )р © + + 5а») | ^

т]п{ г,п} йп

К2( г, Ъ р) = - р5а ® | к

0 к(Л)' 1

Кз( г, р) = -е 2 р 2 а© |р(п№,

тп{ г ,п} 1

К 4

(г, р) = -е2 р2 Е(-) }р ,

/1 (г, р) =--1

тт{

1 г йг

1 ~2 1 ~ ~

|С„(г)р(г)Жл2-1(г,р)^ = — Ж -Ж^) |г=1 . ор

Аналогично были получены уравнения для нахождения поправок коэффициента теплопроводности

1

1К (г)

2

п—1

йг

1

йг = - (Жт - Жп-1) |г=1 р

и модуля Юнга

5|

1 ~ ~ Еп(г)Жп2-^г = т(Жт - Ж-1)|г=1 .

р 0к (г)

Решение обратной задачи

Обратная задача анализируется в пространстве трансформант. На первом этапе находим начальное приближение искомого коэффициента по известным его торцевым значениям. Так, например, для удельной теплоемкости начальное приближение С0(г) = с (0)(1 - г) + с(1)г . На втором этапе на основе метода линеаризации [1-3] строим итерационный процесс по уточнению реконструируемой функции, на каждом этапе которого необходимо решать интегральное уравнение Фредгольма первого рода. Так, для реконструкции удельной теплоемкости при известных других коэффициентах получено уравнение

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

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

^-1 = |(Жп-1(1, р) - Жт (1, р))2 йр.

Результаты вычислений

В настоящем исследовании натурный эксперимент заменен вычислительным, при этом все безразмерные коэффициенты, кроме искомого, полагались равными единице. Восстанавливались гладкие функции: степенные, тригонометрические, экспоненциальные. В ходе вычислительных экспериментов выяснено, что постепенное увеличение параметра связанности 5 практически не влияет на результаты реконструкции, немного увеличивая погрешность восстановления. Наибольшая погрешность восстановления удельной теплоемкости возникала в окрестности торца г = 0, а коэффициента теплопроводности - в окрестности г = 1. Монотонные функции восстанавливались гораздо лучше немонотонных.

Работа выполнена при поддержке РФФИ, грант №10-01-00194-а

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

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

2. Ватульян А.О., Нестеров С.А. Об одном подходе к восстановлению коэффициентов переноса // Изв. вузов. Северо-Кавказский регион. Естественные науки. 2009. №3. С. 39-43.

3. Ватульян А.О., Нестеров С.А. Коэффициентные обратные задачи термоупругости для неоднородных тел // Экологический вестник научных центров ЧЭС. 2009. №3. С. 24-30.

к

0

Проблемы идентификации неоднородных свойств термоупругой среды

1659

THE PROBLEMS OF IDENTIFICATION OF INHOMOGENEOUS PROPERTIES OF THERMO-ELASTIC MEDIA

S.A. Nesterov

The problem of identification of inhomogeneous properties of thermo-elastic media is considered. The iteration algorithm, based on step-by-step solution of Fredholm integral equation of the first kind is proposed. Results of the numerical experiments are given.

Keywords: inhomogeneous rod, integral equation, thermo-elasticity.

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