Механика деформируемого твердого тела Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (4), с. 1473-1475
УДК 539.3
О РЕКОНСТРУКЦИИ ОСТАТОЧНЫХ НАПРЯЖЕНИЙ В ТВЕРДЫХ ТЕЛАХ © 2011 г. В.В. Дударев, Р.Д. Недин
Южный федеральный университет, Ростов-на-Дону
dudarev_vv@mail.ru
Поступила в редакцию 15.06.2011
Представлена общая постановка задачи о реконструкции поля предварительных напряжений при анализе установившихся колебаний. Предложены методы решения прямых задач об определении компонент вектора перемещения для одномерных (интегральные уравнения Фредгольма второго рода) и двухмерных (метод конечных элементов) областей. Даны практические рекомендации по идентификации одноосного предварительного напряженного состояния. Для решения обратной задачи описан итерационный процесс на основе решения интегрального уравнения Фредгольма первого рода с использованием регуляризующей процедуры А.Н. Тихонова. Дана краткая характеристика численных экспериментов по восстановлению монотонных и немонотонных законов изменения неизвестной функции.
Ключевые слова: остаточные напряжения, обратная задача, неоднородное предварительное напряженное состояние, некорректная задача.
Введение
В настоящее время одним из важных направлений развития механики разрушения является проблема исследования остаточных напряжений в твердых телах. Остаточными или предварительными напряжениями называются такие напряжения, которые существуют в теле при отсутствии внешних нагрузок (силовых или температурных) [1]. Обычно подобные напряжения возникают в ходе различных технологических процессов: литья, сварки, скрутки, раскрутки, прокатки, закаливания, резкого изменения температурного режима, ковки и т.д. Также они могут «накапливаться» в теле при продолжительных эксплуатационных нагрузках. Область локализации таких напряжений сосредоточена в окрестности различных неоднородно стей и дефектов (трещины или полости). В силу своей природы, остаточные напряжения явным образом никак не проявляют себя. С другой стороны, при наложении критических нагрузок такие напряжения приводят к разрушениям.
Проблема идентификации и моделирования предварительно напряженного тела интересовала исследователей-механиков еще с конца XVIII века. На сегодняшний день существуют несколько методов определения остаточных напряжений: разрушающие, полуразрушающие и неразруша-ющие. Первые два класса методов хорошо разработаны и имеют достаточно высокую степень точности. Методы, относящиеся к третьему классу,
считаются наиболее востребованными, поскольку используются при диагностике объектов ответственного назначения. Одним из эффективных и хорошо зарекомендовавших себя способов нераз-рушающей идентификации является акустический метод. В его основе лежит процесс определения неизвестных данных по амплитудно-частотной характеристике (АЧХ), которая может быть снята с доступной для исследования части тела.
Общая постановка задачи
Рассмотрим упругое тело объемаом V, ограниченное поверхностью <5 = , с постоянными Ламе X, | и плотностью р. Будем считать, что в рассматриваемом объекте имеется некоторое неоднородное предварительное напряженное состояние, характеризующееся компонентами тензора напряжений с0-. Пусть колебания вызываются нагрузкой р-вш, приложенной на части 8С, а часть 8и закреплена. Линеаризованные уравнения колебаний после отделения временного множителя имеют вид [2, 3]:
Т1], 1 + Р®4 = 0 Т1] = с- + и1т с %,
с- = ик ,к81 + (и X и1 \ви = 0 Тип;К = р1,
где Е-(и) = 0.5(и,- + и- - компоненты тензора малых деформаций.
В качестве дополнительной информации в обратной задаче об определении компонент с0-
1474
В.В. Дударев, Р.Д. Недин
предлагается воспользоваться данными об измеряемых на границе компонентах поля смещений fi в зависимости от частоты колебаний ю е е [ю1; ю2]. Сформулированная обратная задача относится к классу нелинейных некорректных задач. Одним из эффективных способов решения подобных коэффициентных обратных задач является метод организации итерационных процессов.
О задачах для предварительно напряженных тел
В работе [4] была сформулирована вариационная постановка о колебаниях предварительно напряженных тел:
5(П - К + П),
где
П = 0.5^aijUi ]йу - |р1и1йя
V ^
— потенциальная энергия,
К = 0.5ю2 ¡ри2 йу
V
— кинетическая энергия,
П о = С %ии™ии }йу. На основе данной вариационной постановки получены уравнения изгибных [4], продольных и крутильных колебаний стрежня при наличии одноосного предварительного напряженного со сто -яния. В одномерном случае решение прямой задачи об определении функции смещения может быть сведено к решению интегрального уравнения Фредгольма 2-го рода. Для решения задачи об определении компонент вектора смещения в случае двумерных областей использованы конечно-элементные технологии. Одним из пакетов, позволяющих производить необходимые расчеты, является программный пакет РгееРеш++. Для описания модели предварительно напряженного тела в данной программе необходима слабая постановка, которая имеет вид
| [Я и ¿IVу + 2|огу (и)г1]- (у) + и1ту1у]о°т]- -
V
-ю2р иуг]йУ- |Р1У1йБ = 0,
•о
где у1 — пробные функции, удовлетворяющие следующему условию на границе у1 = 0.
Анализ влияния величины одноосного напряженного состояния на АЧХ тела в одномерном и двумерном случаях показал, что это влияние наиболее заметно для частот, находящихся не выше
третьей резонансной частоты (для величины
—2 0 —3
10 <Сц/Е < 10 ). На основе полученных
данных можно сказать, что процедуру реконструк-
ции следует проводить в частотном диапазоне, находящемся до третьей резонансной частоты, исключая частоты резонансов.
Решение обратной задачи об определении закона изменения одноосного предварительного напряженного состояния с°1 Ф 0 реализовано с помощью итерационных процессов. Для вычисления необходимых поправок сформулировано интегральное уравнение Фредгольма (ИУФ) 1-го рода. Решение полученного ИУФ организовано с использованием метода А.Н. Тихонова при автоматическом выборе параметра регуляризации. На первом шаге итерационного процесса начальное приближение можно выбирать в классе линейных функций вида у = ах + Ь. Коэффициенты а, Ь определяются из условия минимума функционала невязки
®2
J = |(и(ю) -[_(ю))2йю
на некотором компакте, построенном по априорной информации об ограниченности восстанавливаемой функции \ст°1 \ < С.
Эффективность предложенного метода решения обратной задачи для определения существенно неоднородного предварительно напряженного состояния в одномерных и плоских областях продемонстрирована в серии численных экспериментов [4, 5]. Для монотонных законов восстанавливаемой функции относительная погрешность составила не более 10%. В стержневых моделях немонотонные законы восстанавливались существенно хуже.
Авторы выражают благодарность своему научному руководителю проф. А.О. Ватульяну за предложенные пути решения задачи и внимание к работе.
Работа выполнена при поддержке РФФИ (грант №.10-01-00194) и в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009—2013 годы (госконтракт №П596).
Список литературы
1. Биргер И.А. Остаточные напряжения. М.: Маш-гиз, 1963. 232 с.
2. Гузь А.Н., Махорт Ф. Г., Гуща О. И. Введение в акустоупругость. Киев: Наук. думка, 1977. 152 с.
3. Ватульян А.О. Проблемы идентификации неоднородных свойств твердых тел // Вестн. Самар. унта. 2007. Вып. 54, №4. С. 93—103.
4. Ватульян А.О., Дударев В.В. О некоторых проблемах реконструкции неоднородного предварительно напряженного состояния в упругих телах // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2009. Т. 9, вып. 4, ч. 2. С. 25—32.
О реконструкции остаточных напряжений в твердых телах
1475
5. Ватульян А.О., Дударев В.В., Недин Р.Д., Саа-кян Я. Г О некоторых задачах идентификации предварительных напряжений // Современные проблемы ме-
ханики сплошной среды. Труды XIV Междунар. конф. Ростов-на-Дону, Азов, 19 - 24 июня 2010 г. Ростов н/Д: Изд-во ЮФУ, 2010. Т. 1. С. 86-90.
ON RESTORING OF RESIDUAL STRESSES IN SOLID BODIES
V.V Dudarev, R.D. Nedin
A general formulation of the prestress field reconstruction problem in the analysis of the stationary vibration is represented. The direct problem solving methods of determining the displacement vector components for one-dimensional (Fredholm integral equations of the second kind) and two-dimensional (FEM) areas are proposed. The practical recommendations of the uniaxial prestress state identification are given. The iterative process based on solving the Fredholm integral equation of the first kind using Tikhonov's regularization procedure is described for solving the inverse problem. Numerical experiments on the reconstruction of monotonic and nonmonotonic unknown function are briefly characterized.
Keywords: residual stresses, inverse problem, inhomogeneous prestressed state, ill-posed problem.