Анализ напряженного состояния плоских элементов конструкций в условиях диффузии водорода

Нефедова, Ольга Анатольевна

Анализ напряженного состояния плоских элементов конструкций в условиях диффузии водорода тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Нефедова, Ольга Анатольевна АВТОР

кандидата технических наук УЧЕНАЯ СТЕПЕНЬ

Екатеринбург МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.02.04 КОД ВАК РФ

Диссертация по механике на тему «Анализ напряженного состояния плоских элементов конструкций в условиях диффузии водорода»

Автореферат диссертации на тему "Анализ напряженного состояния плоских элементов конструкций в условиях диффузии водорода"

На правах рукописи

НЕФЕДОВА ОЛЬГА АНАТОЛЬЕВНА

АНАЛИЗ НАПРЯЖЕННОГО СОСТОЯНИЯ ПЛОСКИХ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ В УСЛОВИЯХ ДИФФУЗИЙ ВОДОРОДА

Специальность 01.02.04 - Механика деформируемого твёрдого тела

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

2 8 НОЯ 2013

Екатеринбург - 2013

005541515

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте машиноведения Уральского отделения Российской академии наук

Научный руководитель:

Федотов Владимир Петрович, доктор технических наук, старший научный сотрудник

Официальные оппоненты: Емельянов Игорь Георгиевич,

доктор технических наук, профессор, Федеральное государственное бюджетное учреждение науки Институт машиноведения Уральского отделения Российской академии наук, главный научный сотрудник

Добычин Иван Александрович, доктор технических наук, доцент, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Уральский государственный университет путей сообщения», профессор

Ведущая организация: Федеральное государственное автономное

образовательное учреждение высшего профессионального образования «Уральский федеральный университет имени первого Президента России Б.Н, Ельцина»

Защита состоится 19 декабря 2013 г. в 15-00 часов на заседании диссертационного совета Д 004.023.01 при Федеральном государственном бюджетном учреждении науки Институте машиноведения Уральского отделения Российской академии наук по адресу: 620049, г. Екатеринбург, ул. Комсомольская, 34.

С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного учреждения науки Института машиноведения Уральского отделения Российской академии наук.

Автореферат разослан ноября 2013 г.

Ученый секретарь диссертационного совета доктор технических наук, профессор

Коновалов А.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследования. Цели и задачи, возникающие при анализе поведения сложных технических систем, диктуют соответствующие требования к теории и расчетному аппарату. Если для оценки прочности крупных объектов, например, моста, достаточно использовать экспериментальные данные по деформации стандартного образца, а для математических расчетов использовать схему конечных разностей или метод конечных элементов, то для соответствующего анализа конструкций в экстремальных условиях эксплуатации этого уже недостаточно. При воздействии различных физических полей могут происходить фазовые превращения, что приводит к изменению свойств, необходимых для расчетов. Путь экспериментального определения этих изменений зачастую невозможен, поскольку эксперимент требует учета только одного вида воздействия при прочих равных условиях. Если при расчете крупных объектов достаточно приписать материальной точке свойства, полученные для стандартного образца, причем только при механическом воздействии, то при наличии физических воздействий возникает проблема согласования, связанная с различными представительными объемами и характерными временами механических и физических процессов для одного и того же объекта. Среди агрессивных технологических примесей водород занимает особое место. Водород - самый подвижный из всех газов, поэтому в большинстве технологических операций возникают процессы сорбции водорода металлом, которые могут протекать при любых температурных режимах. Водородосодержащая среда, проникая в элементы конструкций, приводит к существенному ухудшению важнейших служебных характеристик металлов. В настоящее время потребление водорода в ряде областей промышленности и техники постоянно растет. Это связано с использованием водорода в качестве высокоэффективного и безопасного топлива. Поэтому разработка методов анализа напряженного состояния, роста и развития трещин в деформируемом материале в условиях диффузии водорода становится одной из основных задач для моделирования процессов производства и эксплуатации металлических изделий в химической и нефтегазовой промышленности.

тических работах А.Е. Андрейкива, A.B. Гапонцева, В.В. Кондратьева, A.A. Трещева, B.C. Харина, Р.В. Херцберга и др., а также в экспериментальных работах R.W. Davidge, H. Gao, P. Oliveira, F.L. Ri ley и др. Для моделирования процесса накопления поврежденности в металлах, находящихся под нагрузкой в агрессивных средах, используют концепции накопления поврежденности, предложенные A.B. Березиным, B.JI. Колмогоровым, В.В. Новожиловым, Ю.Н. Работновым, M. Chrzanovski, S. Taira, R.A. Oriani, E.O. Orowan и другими авторами. Однако в основе большинства расчетных моделей лежат алгоритмы, разработанные для реализации последовательных вычислений -метод сеток, метод граничных элементов, метод конечных элементов. С развитием вычислительных технологий представляется актуальным развитие методов решения сложных задач, в алгоритмы которых изначально закладывается принцип распараллеливания.

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

Задачи исследования.

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

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

Построить математическую модель связной деформационно-диффузионной задачи.

Разработать алгоритм и программное обеспечение для решения связной задачи.

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

Научная новизна результатов работы заключается в следующем.

Разработан алгоритм решения модифицированным методом граничных элементов (ММГЭ) задачи упругопластического деформирования, в структуре которого заложен принцип распараллеливания вычислений.

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

по произвольно ориентированному прямолинейному граничному элементу для произвольной точки наблюдения.

Построена математическая модель связной деформационно-диффузионной задачи.

Разработан алгоритм решения на основе ММГЭ связной задачи с возможностью численной реализации на многопроцессорных вычислительных комплексах.

Теоретическая и практическая значимость работы.

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

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

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

Результаты исследования используются в курсе лекций по дисциплине «Механика пластической деформации и разрушения материалов» на кафедре «Обработка металлов давлением» Уральского федерального университета имени первого Президента России Б.Н.Ельцина.

Отдельные разделы диссертационной работы выполнялись в рамках работ по планам научно-исследовательских работ Института машиноведения УрО РАН, по госбюджетным темам «Построение математических и компьютерных моделей деформирования и разрушения материалов и конструкций и разработка методов совершенствования процессов изготовления изделий из материалов различной структуры с применением современных вычислительных методов и информационных технологий» (гос. per. № 01201354600, 2013г.); «Разработка методов оценки и управления ресурсом элементов конструкций машиностроительных систем на основе создания новых фундаментальных моделей и методов механики деформирования, прочности, разрушения, применения компьютерных методов проектирования технологии их изготовления, современного математического аппарата, алгоритмов и программ параллельного действия» (гос. per. № 01200951982, 2009-2012 гг.); Программам Президиума РАН №14 «Интеллектуальные информационные техноло-

гии, математическое моделирование, системный анализ и автоматизация», №15 «Математическое моделирование нелинейных связных физико-механических процессов на основе модифицированного метода граничных элементов»; Региональной целевой программы развития вычислительных, телекоммуникационных и информационных ресурсов УрО РАН, проекты №№ РЦП-12-П11, РЦП-13-П4.

Методология и методы исследования. В работе были использованы теоретические и вычислительные методы исследования. Решение задачи уп-ругопластического деформирования, диффузионной задачи и связной задач осуществлялось с помощью модифицированного метода граничных элементов. Компьютерная программа была написана на языках программирования c++ и Java, с использованием библиотеки параллельных вычислений MPI, и была реализована на суперкомпьютере «Уран» ИММ УрО РАН. С помощью созданной программы был решен ряд задач, результаты сравнивались с экспериментальными и расчетными данными других авторов.

Положения, выносимые на защиту.

Алгоритм решения двумерной задачи уиругопластического деформирования модифицированным методом граничных элементов.

Алгоритм решения двумерной задачи диффузии модифицированным методом граничных элементов.

Непрерывные аналитические аппроксимирующие функции для решения двумерной задачи диффузии.

Математическая модель связной деформационно-диффузионной задачи.

Алгоритм и программное обеспечение для решения двумерной связной задачи на основе ММГЭ.

Качественные оценки влияния водорода на прочность металлических деталей на примерах решения задач.

Достоверность результатов подтверждается сопоставлением результатов моделирования с экспериментами, а также с расчетами, проведенными другими авторами.

Апробация результатов работы. Основные положения и результаты диссертационной работы доложены на следующих конференциях: IV Российская научно-техническая конференция «Ресурс и диагностика материалов и конструкций», Екатеринбург, 2009 г.; X Международная конференция «Актуальные проблемы современной науки», Самара, 2009 г.; VI Российская научно-техническая конференция «Механика микронеоднородных материалов и разрушение», Екатеринбург, 2010 г.; VIII Всероссийская научная конференция с международным участием «Математическое моделирование и краевые задачи», Самара, 2011 г.; V Российская научно-техническая конференция «Ресурс и диагностика материалов и конструкций», Екатеринбург, 2011 г.;

VII Российская научно-техническая конференция «Механика микронеоднородных материалов и разрушение», Екатеринбург, 2012 г.; IX Всероссийская научная конференция с международным участием «Математическое моделирование и краевые задачи», Самара, 2013 г.

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

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

Вторая глава посвящена разработке алгоритма, наиболее адаптированного к высокому уровню распараллеливания, для решения краевой задачи деформирования при малых упругопластических деформациях. В систему уравнений краевой задачи для некоторой области £2 входят уравнения равновесия, записанные в компонентах тензора напряжений Коши а,,,

соотношения Коши, записанные для компонент тензора деформации еу-.

определяющие соотношения в рамках деформационной теории пластичности, характеризующие свойства материала,

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

О)

(2)

(3)

(4)

(5)

и граничные условия:

на Г„: и = и ;

на Г/\о1]п) = /,=/'■ (7)

Здесь гц - компоненты вектора перемещений; с - среднее давление в точке; £ — относительное изменение объема; к - коэффициент объемного сжатия, к=(1-2у)/Е\ V - коэффициент Пуассона; Е - модуль сдвига; л'/;. - компоненты девиатора напряжений; — компоненты девиатора деформации; /а - интенсивность касательных напряжений; /е - интенсивность деформации сдвига; (3 - модуль сдвига, С=Е/2/(1+г); т^- предел текучести материала при чистом сдвиге; сгг. = г5./л/3 - предел текучести материала при растяжении; а в - временное сопротивление материала на разрыв; напряжения ав и ат соответствуют деформациям ев и е7; Г = Ги и Т/ - гладкая граница области О; / - компоненты вектора поверхностных напряжений, п, - компоненты вектора внешней нормали к поверхности; по повторяющемуся индексу производится суммирование; звездочкой отмечены известные функции. Для разделения упругой и пластической областей было использовано условие текучести Ми-зеса.

В области пластического течения были приняты физические уравнения связи, разработанные В.П. Федотовым:

(8)

Для расчета прочности и прогнозирования разрушения материала в рамках модели (8) достаточно задать две экспериментальные характеристики — предел текучести при растяжении ат и временное сопротивление на разрыв од, которые имеют ясный физический смысл и зависят от параметров состояния (температуры, концентрации примеси и т.д.).

Использование физических уравнений связи (8) в краевой задачи уп-ругопластического деформирования (1) — (7) позволило построить алгоритм решения задачи в широком диапазоне физических воздействий. Алгоритм был построен для случая плоского деформированного состояния. Решение осуществлялось по шагам деформации, деформации на шагах суммировались. На каждом расчетном шаге решается линейная задача с применением модифицированного метода граничных элементов, разработанного для решения линейных задач. Выбор ММГЭ был обусловлен тем, что в его структуре изначально заложены возможности, позволяющие распараллелить счет на всех этапах решения задачи и, следовательно, увеличить скорость вычислений с привлечением технологий параллельного программирования.

В результате применения алгоритма на каждом шаге определяются перемещения, по ним вычисляются деформации, по известным деформациям рассчитываются первый инвариант тензора деформации и интенсивность де-

формации сдвига. Первый инвариант тензора напряжений и интенсивность касательных напряжений вычисляются с использованием определяющих соотношений по точным характеристикам деформированного состояния.

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

Разработанный алгоритм дает возможность моделировать поведение упругопластического материала в процессе деформации по известным краевым условиям задачи и по известным физическим характеристикам: модулю Юнга Е, коэффициенту Пуассона V, пределу текучести при растяжении <зт и временному сопротивлению на разрыв ев-

Разработанный на основе ММГЭ алгоритм был распараллелен на каждом этапе решения. На каждом шаге коэффициенты системы линейных уравнений вычисляются независимо, следовательно, процесс построения системы линейных уравнений полностью распараллеливается. Решение системы осуществляется методом Ш-разложения, допускающим распараллеливание. Расчет напряженно-деформированного состояния во внутренних точках области также можно производить независимо, т.е. параллельно. Алгоритм был реализован в виде программы. Программа позволяет решать задачу упругопластического деформирования в рамках постановки (1) - (7) для двумерных областей любой геометрии.

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

Результаты расчетов показаны на рисунках 2, 3. На рисунке 2 приведено сравнение графиков изменения /„ вдоль оси Охз, построенных с помощью предложенного алгоритма и с использованием АЫБУЗ.

На рисунке 3 показаны результаты расчетов /о по ММГЭ при различных значениях параметров ат и а в, соответствующих различным значениям концентрации водорода. Из рисунка видно, что расчетные интенсивности касательных напряжений заметно различаются. Таким образом, при решении деформационной задачи необходимо учитывать процесс диффузии, а значит решать связную деформационно-диффузионную задачу.

Х2 Р

ь ¡1 V

//////// //////// XI

* а *

Рисунок 1 - Схема задачи упругопластического деформирования в условиях плоского деформированного состояния

300

200

100 -,-,-,-,-

0 0,002 0,004 0,006 0,008 х2, м

1 — А№ Ув; 2 - ММГЭ Рисунок 2 - Изменение /о вдоль оси Ох2

1 - ат=578 МПа, св=648 МПа; 2 - с1т=598 МПа, ав=628 МПа Рисунок 3 - Изменение /о вдоль оси Ох, при различных значениях физических параметров

Третья глава посвящена разработке алгоритма решения двумерной задачи диффузии с заложенной в нем' возможностью распараллеливания счета. Выбор метода решения диффузионной задачи вновь был сделан в пользу ММГЭ в силу его адаптированности для проведения параллельных вычислений. Краевая задача диффузии включает в себя уравнение параболического типа

Э2сМ 1 дс(х,0

эх! о а ' ^

начальные и граничные условия

ф, О = <£(*); (Ю)

на Гс: с(хо,*)=с(*о,0;

наГq{x0,t)=q{xDJ). (11)

Здесь функция с(х,() - концентрация в точке х в момент времени 10 - начальный момент времени; £) - коэффициент диффузии; Г = Гс и Г - гладкая граница области £2; х0 — граничная точка; д(х,() — диффузионный поток; с',(*), с'(х0,/), д'{хг,{) - известные функции.

Решение краевой задачи диффузии (9) - (11) осуществлялось по шагам по времени к=\,2,...К в предположении постоянства значений концентрации и диффузионного потока в узлах граничных элементов на каждом временном шаге. На каждом временном шаге решалось граничное интегральное уравнение

1 / \ \ , г \ . / \ к \

+ 1с'0(хр'(х0,х,т,(0^а(х), (12)

-с(х0,т)=О£ ) ¡с{х,^(х0,х,т) ¡д(х,1>р'{х0,х,т,()+

I 4-11,., Г Г

где т — момент наблюдения; ^ (х0^х,г,1) и С ( х0гк,т,1) - функции влияния, которые в двумерном случае определяются следующим образом:

I '

(13)

пЧ ч Ж' И |

р —- пх — -' , ч ехр

ох х 4Д(г-/)л/лО(г-7) \4£>(г-г)

где г1=г1г,\ г, х, -х01.

В соответствии с методом граничных элементов граница Г исследуемой двумерной области £2 была разбита на конечное число граничных элементов - отрезков [0/^,0/], у'=1,2,...№ Л/, на одной части которых {[о^я,], 7=1,2были заданы граничные условия для концентрации, а на другой

{[aj.i,ctj], j=N+l,N+2,...N+M } — для диффузионного потока. Тогда граничное интегральное уравнение (12) принимает вид:

\c{x0,t)=DY ) S ](c{ñ\xj)F\xB,x,Tj)~q<>\xj)G4xü,x,T,t))ixdt +

¿ It_, М с,.,

К h V-А/ .

+ Z)EÍ S \{c{'\xj)FXxa,x,Tj)-q,'^xjyr{x0,x,Tj))hdt +

(1 ч у-Л'+Ц ,

+ Jc¡(x)G,(^^,r,í0)afn(x). (14)

Здесь с(уГ(х,/)=с*(х,0, qm(x,t)=q(x,t) при хе [о,.,,a,], /=1,2,... Л'; с°\х,{)=с(х^), qtj)'(x,t)-q"(x,t) при д:е [а>у,йг;],у'=Аг+1,Л'+2,...Лг+Л/; звездочкой обозначены известные величины.

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

На каждом шаге по времени уравнение (14), записанное для всех узлов граничных элементов, порождает систему линейных уравнений, решение которой определяет значения концентрации и диффузионного потока в узлах границы, не заданные граничными условиями. Коэффициенты системы являются интегралами от произведений аппроксимирующих функций, удовлетворяющих уравнению диффузии, и функций влияния (13). Это различные комбинации произведений функции ошибок erf(x0,T,0 и функции exp(x0,T,í)-Подход ММГЭ предполагает переход от интегралов по произвольному граничному элементу к интегралам по удобному "базовому" элементу. С использованием такого преобразования были получены аналитические формулы для вычисления всех необходимых при решении двумерной задачи диффузии интегралов.

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

л: лг+л/ "•'/,. ч

+ | Е | 'К & О - 9(>)* (*, +

+ \сй{хр\Ъх,т,10)с1О.{х). (15)

Интегралы в правой части уравнения (15) также вычисляются по полученным аналитическим формулам.

На каждом шаге по времени вычисление интегралов для различных граничных точек при формировании системы линейных алгебраических уравнений может выполняться независимо. Вычисление значения концентрации во внутренних точках также можно выполнять независимо друг от друга. Это позволило распараллелить алгоритм на всех этапах решения задачи. Следовательно, как и для деформационной задачи, для задачи диффузии был построен алгоритм решения на основе ММГЭ, использующий аналитическое интегрирование и распараллеливание счета. Алгоритм был реализован в виде программы.

В качестве примера применения предложенного алгоритма и программы была рассмотрена модельная задача диффузии газа в тонкой квадратной пластине. На границе пластины была задана постоянная концентрация газа с=20л< ~3, в начальный момент времени концентрация во всей пластине принималась равной нулю. На рисунке 4 приведено сравнение результатов решения для концентрации, полученных разными методами: с помощью предложенного алгоритма (кривая 1), по методу рядов Фурье (кривая 2) и с помощью метода конечных элементов (МКЭ) (кривая 3). Все расчеты выполнялись для точки, расположенной в центре пластины.

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

1 - ММГЭ; 2 - метод рядов Фурье; 3 - МКЭ Рисунок 4 - Сравнения результатов решения для концентрации, полученных разными методами

Для рассмотренного примера было проведено сравнение времени счета по предложенному алгоритму при различных количествах задействованных процессоров и с помощью МКЭ. В таблице 1 представлены результаты при различных количествах шагов по времени. Граница пластины была разбита на 400 элементов, в методе конечных элементов область разбивалась на 200 элементов.

Таблица 1 - Время расчета модельной задачи с помощью метода конечных элементов и с помощью модифицированного метода граничных элементов при различных количествах шагов по времени

Количество шагов по времени Время расчета МКЭ, с. Время расчета ММГЭ, с.

4 процессора 80 процессоров 400 процессоров

20 155 220 23 3

80 720 1065 95 10

100 840 1264 129 15

200 5620 6260 564 60

В таблице 2 представлено время расчета задачи с помощью ММГЭ при различных количествах задействованных процессоров. Приведены результаты при различных количествах граничных элементов.

Таблица 2 - Время расчета модельной задачи ММГЭ при различных количествах граничных элементов

Количество граничных элементов Время расчета ММГЭ, с.

4 процессора 80 процессов 400 процессоров

80 60 3 3

400 1500 75 15

800 4500 225 45

Все результаты были получены на многопроцессорном вычислительном комплексе «Уран» ИММ УрО РАН. Как видно из представленных таблиц, при увеличении числа задействованных процессоров время расчета уменьшается на 1-2 порядка.

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

В работе рассматривается модель водородной хрупкости, обусловленная высоким давлением молекулярного водорода в полостях. Реальные металлические изделия никогда не обладают однородной структурой. Они содержат широкий набор полостей — пор и трещин, распределенных случайным образом. Атомарный водород, диффундируя в металле, проникает в полости, где ассоциируется в не диффундирующий молекулярный водород. В результате возникает избыточное внутреннее давление водорода внутри полостей, которое приводит к водородному охрупчиванию.

Система уравнений связной диффузионно-деформационной задачи включает в себя уравнение диффузии в поле высоких механических напряжений:

= (16) а? ят

с заданными начальными и граничными условиями:

ф,/0) = 0; (17)

на г,и, ■ с(х,/)=с"(х,0;

на Г,„: с(х,/)=0; (18)

уравнения механики деформируемого твердого тела:

= (19)

£„.=(и^+и„)/2; (20)

аи = 2Се0+^-8,етт (21)

и граничные условия:

на г„,:

на Г-:/(х) = Р;(х). (22)

Здесь Д - оператор Лапласа; Ун - парциальный молярный объем водорода в металле; Я - газовая постоянная; Т - абсолютная температура; V — градиент функции; V о- V с - скалярное произведение векторов Г = Ги и Гт, ■ Г„ -внутренняя граница области £2 (граница полости); - внешняя граница области П. Р\ - внутреннее давление в дефектах; Р2' - действующее на внешней границе области внешнее механическое напряжение. Звездочкой, как и прежде, отмечены заданные функции. Связность задачи (16) - (22) состоит в наличии градиента напряжений в правой части уравнения диффузии (16), а также обусловленного диффузией давления в первом из граничных условий (22) механической задачи.

Решение задачи (16) - (22) осуществляется по шагам по времени, на каждом шаге решение включает в себя четыре этапа. На первом этапе выполняется решение деформационной задачи в рамках линейной теории упругости, описываемой уравнениями (19) - (22). Решение осуществляется с помощью ММГЭ, в соответствии с алгоритмом уже использовавшимся в главе 2 на каждом шаге пошагового алгоритма. На первом шаге по времени полагается, что Р\ =0, на каждом последующем временном шаге Р\ определяется из решения диффузионной задачи на предыдущем шаге. В результате решения задачи (19) - (22) внутри области П определяются значения перемещений, напряжений, а также пространственных производных от напряжений, необходимые для вычисления V а.

На втором этапе рассматривается краевая задача диффузии для двумерной области, подверженной внешнему механическому нагружению. Неоднородное уравнение диффузии для этого случая имеет вид (16), краевые условия имеют вид (17), (18). Решение данной задачи также осуществляется с помощью ММГЭ по алгоритму, построенному в главе 3. Граничное интегральное уравнение, соответствующее задаче (16) - (18), имеет вид:

1 >. г

+ } | у(х, < )3' (ха,х, г,Г )сЮ.(х .

Функция у/{х,1) = Уи-Ус создает в уравнении диффузии неоднород-ЯТ

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

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

На третьем этапе вычисляется давление молекулярного водорода в каждой из полостей в предположении, что молекулярный водород удовлетворяет уравнению состояния идеального газа

где т - масса водорода в полости; М - молярная масса газа; Ум -молярный объем газа.

Для определения массы водорода в полости нужно вычислить количество водорода, поступающего в полость за время шага Лг=/г4-ь равное его полному потоку через границу полости

т = 50\ц(х,1)Ж.

'1-1

Здесь 5 - площадь поверхности полости.

На четвертом этапе проверяется условие разрушения границы каждой полости. Для этого производится сравнение внутреннего давления водорода в полости с критическим давлением на предмет выполнения неравенства:

Ркр>3<гТ.

Здесь Ркр — критическое давление, при котором наступает разрушение материала. Если критическое давление не достигнуто, расчеты продолжаются на следующем временном шаге с учетом перерасчета граничных условий деформационной задачи.

Поскольку каждый этап решения задачи допускает распараллеливание, как отмечалось в главах 2 и 3, алгоритм связной задачи также существенно распараллелен. Построенный алгоритм был реализован в виде программы. В качестве примера применения алгоритма и программы была рассмотрена модельная задача растяжения в водородосодержащей среде прямоугольной области с круговой полостью (порой) в условиях плоского напряженного состояния. Схема задачи показана на рисунке 5.

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

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

Х2 [И ' р, т

Рг

Рисунок 5 - Схема задачи растяжения пластины в водородосодержащей среде

1, (мгс)-' 4x10® 3x10® 2x10® 1x10'

1 9 17 25 33 I, с

--нет внешнего напряжения;

----приложено растягивающее напряжение

Рисунок б - Изменение потока на поверхности поры во времени

Р,, МПа------------------------------------------------

0.8 0.6 0.4 0.2

" 5 13 21 29 37 <> с

--без учета внутреннего давления в механической задаче;

-----с учетом внутреннего давления в механической задаче

Рисунок 7 - Изменение давления в поре во времени

Для рассмотренного примера был проведен анализ времени решения связной задачи ММГЭ при различных количествах задействованных процессоров (таблица 3). Граница пластины была разбита на 580 элементов, рассматриваемый интервал времени был разбит на 120 шагов.

Таблица 3 - Время расчета связной задачи с помощью ММГЭ

58 процессоров, мин. 145 процессоров, мин. 290 процессоров, мин.

260 40 25

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

ЗАКЛЮЧЕНИЕ

1. Разработан пошаговый алгоритм решения двумерной задачи упру-гопластического деформирования на основе ММГЭ, адаптированный к распараллеливанию вычислений.

2. Разработан алгоритм решения двумерной диффузионной задачи с помощью ММГЭ, в основу которого заложена возможность распараллеливания счета.

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

4. Проведено сравнение решений модельной задачи о распределении водорода в квадратной пластине, полученных при помощи ММГЭ, МКЭ и по методу рядов Фурье.

5. Построен алгоритм решения двумерной связной деформационно-диффузионной задачи на основе ММГЭ.

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

7. На примере решения задачи деформирования в водородосодержа-щей среде выполнена оценка взаимного влияния процесса диффузии и напряженного состояния в теле с дефектом.

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

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в научных изданиях, входящих в список ВАК РФ:

1. Федотов, В.П. Применение модифицированного метода граничных элементов для решения параболических задач / В.П. Федотов, О.А.Нефедова // Вестник СамГТУ. Сер. Физ.-мат. науки. - 2011. -№4 (25). - С. 93-101.

2. Федотов, В.П. Параллельные алгоритмы для анализа прочности на-водороженных конструкций / В.П. Федотов, Л.Ф. Спевак, О.А.Нефедова // Программные продукты и системы. -2012. -№3. - С. 235-239.

3.Федотов, В.П. Моделирование процессов упругопластического деформирования модифицированным методом граничных элементов / В.П. Федотов, Л.Ф. Спевак, О.А.Нефедова // Программные продукты и системы. -2013.-№4.-С. 253-257.

Основные публикации в научных изданиях, не входящих в список ВАК

РФ:

1. Федотов, В.П. Нахождение аналитического решения параболических задач с использованием фундаментального решения для аппроксимации неизвестных граничных условий / В.П.Федотов, О.А.Нефедова // Ресурс и диагностика материалов и конструкций: тезисы докладов IV Российской научно-технической конференции. - Екатеринбург, 2009. - С. 181.

2. Нефедова, O.A. Применение модифицированного метода граничных элементов для решения задач теплопроводности / О.А.Нефедова,

B.П.Федотов // Актуальные проблемы современной науки: тр. 5-го Международ. форума молодых ученых (10-й Международ, конф.). Естественные науки. Ч. 1-3: Математика. Мат. моделирование. Механика. - Самара, 2009. -

C. 143-147.

3.Федотов, В.П. Влияние диффузии на прочность элементов конструкций с дефектами / В.П.Федотов, О.А.Нефедова // Механика микронеоднородных материалов и разрушение: тезисы докладов VI Российской научно-технической конференции. - Екатеринбург, 2010. - С. 117.

4. Федотов, В.П. Оценка влияния диффузии водорода на развитие дефектов в металле / В.П. Федотов, Л.Ф. Спевак, О.А.Нефедова // Математическое моделирование и краевые задачи", 4.1: "Математические модели механики, прочности и надежности элементов конструкций: труды восьмой Всероссийской научной конференции с международным участием. - Самара, 2011.-С. 244-247

5. Федотов, В.П. Решение связной диффузионно-деформационной задачи для металлического образца с порой и трещиной / В.П. Федотов, Л.Ф. Спевак, О.А.Нефедова // Ресурс и диагностика материалов и конструкций: тезисы докладов V Российской научно-технической конференции. - Екатеринбург, 2011.-С. 143.

6. Федотов, В.П. Решение задачи диффузии с нелинейным источником методом граничных элементов / В.П. Федотов, Л.Ф. Спевак, О.А.Нефедова // Механика микронеоднородных материалов и разрушение: тезисы докладов VII Российской научно-технической конференции. - Екатеринбург, 2012. - С. 150.

7. Федотов, В.П. Применение итерационного метода для моделирования диссипативных процессов в задачах деформирования / В.П. Федотов, О.А.Нефедова // Математическое моделирование и краевые задачи: труды девятой Всероссийская научная конференция с международным участием. -Самара, 2013.

Формат 60x84 1/16 Бумага писчая. Печать офсетная. Усл. печ. л. 1. Тираж 100 экз.

Центр оперативной полиграфии «Копирус» 620075, г. Екатеринбург, ул. Мамина-Сибиряка, 137

Текст научной работы диссертации и автореферата по механике, кандидата технических наук, Нефедова, Ольга Анатольевна, Екатеринбург

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ НАУКИ ИНСТИТУТ МАШИНОВЕДЕНИЯ УРАЛЬСКОГО ОТДЕЛЕНИЯ РОССИЙСКОЙ АКАДЕМИИ НАУК

На правах рукописи

04201450641 НЕФЕДОВА Ольга Анатольевна

АНАЛИЗ НАПРЯЖЕННОГО СОСТОЯНИЯ ПЛОСКИХ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ В УСЛОВИЯХ ДИФФУЗИИ ВОДОРОДА

Специальность 01.02.04 - «Механика деформируемого твердого тела»

ДИССЕРТАЦИЯ на соискание ученой степени кандидата технических наук

Научный руководитель: Федотов Владимир Петрович, доктор технических наук, старший научный сотрудник

Екатеринбург -

2013

СОДЕРЖАНИЕ

Стр.

ВВЕДЕНИЕ........................................................4

1 ОБЗОР ИССЛЕДОВАНИЙ ПО СВЯЗНЫМ

ДЕФОРМАЦИОННО-ДИФФУЗИОННЫМ ЗАДАЧАМ..................12

2 РЕШЕНИЕ ЗАДАЧИ УПРУГОПЛАСТИЧЕСКОГО ДЕФОРМИРОВАНИЯ МОДИФИЦИРОВАННЫМ

МЕТОДОМ ГРАНИЧНЫХ ЭЛЕМЕНТОВ.............................22

2 Л Краевая задача деформирования при малых упругопластических деформациях.......................................................22

2.2 Пошаговый алгоритм решения задачи упруго пластического деформирования....................................................25

2.3 Расчет напряженно-деформированного состояния в упругопластическом теле при различных значениях

концентрации примеси...............................................30

2.4 Выводы........................................................34

3 РЕШЕНИЕ ЗАДАЧИ ДИФФУЗИИ МОДИФИЦИРОВАННЫМ МЕТОДОМ ГРАНИЧНЫХ ЭЛЕМЕНТОВ..............................35

3.1 Краевая задача диффузии.........................................35

3.2 Построение аппроксимирующих функций...........................38

3.3 Модифицированный метод граничных элементов для

двумерной задачи диффузии. Аналитическое интегрирование.............42

3.4 Модельная двумерная задача диффузии водорода.....................55

3.5 Выводы........................................................61

4 РЕШЕНИЕ СВЯЗНОЙ ДЕФОРМАЦИОННО-ДИФФУЗИОННОЙ ЗАДАЧИ МОДИФИЦИРОВАННЫМ

МЕТОДОМ ГРАНИЧНЫХ ЭЛЕМЕНТОВ.............................62

4.1 Математическая модель связной деформационно-диффузионной задачи...........................................................62

4.2 Алгоритм решения связной деформационно-диффузионной задачи .... 66

4.3 Оценка взаимного влияния диффузии и напряженного состояния

в теле с дефектом..................................................68

4.4 Выводы.......................................................73

ЗАКЛЮЧЕНИЕ...................................................75

СПИСОК ЛИТЕРАТУРЫ...........................................77

ПРИЛОЖЕНИЕ А. Акт использования результатов диссертационной работы..........................................................90

ПРИЛОЖЕНИЕ Б. Программный код для решения задачи упругопластического деформирования................................92

ПРИЛОЖЕНИЕ В. Описание компьютерной программы

115

ВВЕДЕНИЕ

Актуальность темы исследования. Цели и задачи, возникающие при анализе поведения сложных технических систем, диктуют соответствующие требования к теории и расчетному аппарату. Если для оценки прочности крупных объектов, например, моста, достаточно использовать экспериментальные данные по деформации стандартного образца, а для математических расчетов использовать схему конечных разностей или метод конечных элементов, то для соответствующего анализа конструкций в экстремальных условиях эксплуатации этого уже недостаточно. При воздействии различных физических полей, например, при сильных перепадах температур, изменении химического состава за счет диффузии примесных атомов могут происходить фазовые превращения, что приводит к изменению параметров необходимых для расчетов. Путь экспериментального определения этих изменений бесперспективен, а зачастую невозможен, поскольку эксперимент требует учета только одного вида воздействия при прочих равных условиях. Такие задачи, моделирующие взаимное влияние нескольких параллельно протекающих механических и физических процессов, будем называть связными задачами. Если при расчете крупных объектов достаточно приписать материальной точке свойства, полученные для стандартного образца, причем только при механическом воздействии, то при наличии физических воздействий возникает проблема согласования, связанная с различными представительными объемами и характерными временами механических и физических процессов для одного и того же объекта. Эта проблема порождает ряд трудностей при решении связных задач. Для математически корректной постановки связной задачи необходимо устремить представительные объемы всех процессов к нулю. Для сеточных методов, включая конечные элементы, это означает существенное увеличение числа узлов, а значит и порядка разрешающей системы линейных алгебраических уравнений. По этому пути сейчас идет развитие численных методов, полагаясь на бурное

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

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

Степень разработанности темы исследования. Различают несколько форм вредного влияния водорода на процесс разрушения деталей металлических конструкций, обусловленных разными факторами: давлением молекулярного водорода в полостях металла; ослаблением водородом межатомных связей в кристаллической решетке металла; взаимодействием атомов водорода с дислокациями; адсорбционными явлениями; хемосорбцией водорода на внутренних поверхностях микротрещин и др. Различные модели, описывающие влияние водорода на прочность металлов, приведены в теоретических работах А.Е. Андрейкива [1, 2], Г.В. Карпенко [3, 4], В.В. Кондратьева [5], A.A. Трещева [6], B.C. Харина [7], Р.В. Херцберга [8] и др., а также в экспериментальных работах О.Д. Смияна [9], Н. Gao [10], P. Oliveira [11], В.И. Швачко [12, 13] и др. Для моделирования процесса накопления поврежденности в металлах, находящихся под нагрузкой в агрессивных средах, используют концепции накопления поврежденности, предложенные A.B. Березиным [14, 15], JI.M.

Качановым [16, 17], B.JI. Колмогоровым [18 - 21], В.В. Новожиловым [22, 23], Ю.Н. Работновым [24, 25], M. Chrzanovski [26], S. Taira [27], R.A. Oriani [28, 29] и другими авторами. Однако в основе большинства расчетных моделей лежат алгоритмы, разработанные для реализации последовательных вычислений - метод сеток, метод граничных элементов, метод конечных элементов. С развитием вычислительных технологий представляется актуальным развитие методов решения сложных задач, в алгоритмы которых изначально закладывается принцип распараллеливания.

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

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

Для достижения поставленной цели необходимо было решить следующие задачи:

1) выполнить обзор исследований совместного влияния силовых воздействий и коррозионных сред на прочность деталей металлических конструкций, а также методов и существующих подходов к решению связных задач в механике деформируемого твердого тела;

2) разработать алгоритм решения задачи упругопластического деформирования, позволяющий выполнять расчеты с высокой скоростью и точностью в широком диапазоне параметров состояния материала;

3) разработать алгоритм решения диффузионной задачи, позволяющий выполнять расчеты с высокой скоростью и точностью;

4) построить математическую модель связной деформационно-диффузионной задачи;

5) разработать алгоритм и программное обеспечение для решения связной задачи;

6) получить возможные качественные и количественные оценки влияния водорода на прочность деформируемых металлических изделий.

Научная новизна результатов работы заключается в следующем:

1) разработан алгоритм решения модифицированным методом граничных элементов (ММГЭ) задачи упругопластического деформирования, в структуре которого заложен принцип распараллеливания вычислений;

2) разработан алгоритм решения двумерной задачи диффузии на основе ММГЭ, в структуре которого заложен принцип распараллеливания вычислений; для аппроксимации параметров задачи на граничном элементе построены удовлетворяющие уравнению диффузии непрерывные аналитические функции; получены аналитические формулы для вычисления интегралов от произведения аппроксимирующих функций и функций влияния по произвольно ориентированному прямолинейному граничному элементу для произвольной точки наблюдения;

3) построена математическая модель связной деформационно-диффузионной задачи;

4) разработан алгоритм решения на основе ММГЭ связной задачи с возможностью численной реализации на многопроцессорных вычислительных комплексах.

Теоретическая и практическая значимость работы.

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

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

Отдельные разделы диссертационной работы выполнялись в рамках работ по планам научно-исследовательских работ Института машиноведения УрО РАН по госбюджетным темам «Построение математических и компьютерных моделей деформирования и разрушения материалов и конструкций и разработка методов совершенствования процессов изготовления изделий из материалов различной структуры с применением современных вычислительных методов и информационных технологий» (гос. per. № 01201354600, 2013г.); «Разработка методов оценки и управления ресурсом элементов конструкций машиностроительных систем на основе создания новых фундаментальных моделей и методов механики деформирования, прочности, разрушения, применения компьютерных методов проектирования технологии их изготовления, современного математического аппарата, алгоритмов и программ параллельного действия» (гос. per. № 01200951982, 2009-2012 гг.); Программам Президиума РАН № 14 «Интеллектуальные информационные технологии, математическое

моделирование, системный анализ и автоматизация» (2009-2011гг.), № 15 «Информационные, управляющие и интеллектуальные технологии и системы» (2012-2013 гг.); Региональной целевой программы развития вычислительных, телекоммуникационных и информационных ресурсов УрО РАН, проекты №№ РЦП-12-П11, РЦП-13-П4.

Методология и методы исследований. В работе были использованы теоретические и вычислительные методы исследования. Решение задачи упругопластического деформирования, диффузионной задачи и связной задач осуществлялось с помощью модифицированного метода граничных элементов. Компьютерная программа была написана на языках программирования c++ и Java, с использованием библиотеки параллельных вычислений MPI, и была реализована на суперкомпьютере "Уран" ИММ УрО РАН. С помощью созданной программы было решено несколько задач, результаты сравнивались с экспериментальными и расчетными данными других авторов.

Положения, выносимые на защиту:

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

2) алгоритм решения двумерной задачи диффузии модифицированным методом граничных элементов;

3) непрерывные аналитические аппроксимирующие функции для решения двумерной задачи диффузии;

4) математическая модель связной деформационно-диффузионной задачи;

5) алгоритм и программное обеспечение для решения двумерной связной задачи на основе ММГЭ;

6) качественные оценки влияния водорода на прочность металлических деталей на примере решения задач.

Достоверность и обоснованность полученных результатов обеспечивается корректной постановкой задачи, использованием современных численных методов решения и подтверждается сопоставлением результатов

моделирования с экспериментами, а также с расчетами, проведенными другими авторами.

Личный вклад автора. В диссертации обобщены результаты исследований, полученные лично автором и в соавторстве. Основная роль в построении вычислительных алгоритмов и компьютерной программы, а также в тестировании программы принадлежит автору работы. Постановка задачи исследования, обсуждение и интерпретация полученных результатов проводились совместно с научным руководителем д.т.н., Федотовым В.П. Экспериментальные данные получены из опытов на растяжение серво-гидравлической испытательной машиной Ш8Т1Ю]Ч. Эксперимент проводил Д.И. Вичужанин. Расчет задачи упругопластического деформирования в пакете А^УБ выполняла М.В. Мясникова.

Апробация результатов работы. Основное содержание работы отражено в 10 публикациях, в том числе, в 3 статьях в журналах, рекомендованных ВАК России. Основные положения и результаты диссертационной работы были доложены на следующих конференциях: IV Российская научно-техническая конференция «Ресурс и диагностика материалов и конструкций», Екатеринбург, 2009 г.; X Международная конференция «Актуальные проблемы современной науки», Самара, 2009; VI Российская научно-техническая конференция «Механика микронеоднородных материалов и разрушение», Екатеринбург, 2010; VIII-aя Всероссийская научная конференция с международным участием «Математическое моделирование и краевые задачи», Самара, 2011; V Российская научно-техническая конференция «Ресурс и диагностика материалов и конструкций», Екатеринбург, 2011; VII Российская научно-техническая конференция «Механика микронеоднородных материалов и разрушение», Екатеринбург, 2012; IX Всероссийская научная конференция с международным участием «Математическое моделирование и краевые задачи», Самара, 2013.

Объем и структура диссертации. Диссертация состоит из введения, четырех глав, заключения, списка цитируемой литературы, включающего 137

наименований, и трех приложений. Содержание диссертации изложено на 121 странице, включая 21 рисунок и 6 таблиц.

1 ОБЗОР ИССЛЕДОВАНИЙ ПО СВЯЗНЫМ ДЕФОРМАЦИОННО-ДИФФУЗИОННЫМ ЗАДАЧАМ

Конструкторам и инженерам на протяжении многих веков приходится решать вопросы механики разрушения и проблемы прочности материалов, связанные с зарождением и развитием трещин. Успешные попытки решения некоторых задач были предприняты еще в XII веке при строительстве уникальных для того времени инженерных сооружений - дворцов и храмов [30]. Однако первые научные исследования в этой области появились только в эпоху Возрождения. Ими занимались великие ученые того времени Леонардо да Винчи (1452 - 1519) и Галилео Галилей (1564 - 1642). Ряд открытий в механике разрушения, выполненных в XVII веке, принадлежит Готфриду Вильгельму Лейбницу (1646 - 1716). Дальнейшими основоположниками этой науки принято считать: Р. Гука (1635 - 1703), Ш. Кулона (1736 - 1806), Л. Навье (1785 - 1836), Ж. Понселе (1788 - 1867), Г. Ламэ (1795 - 1870), Б. де