1. Предварительные сведения
Еще в 1999 году [1] Д. Хейхарст подчеркнул первостепенную роль и значение суперкомпьютерного моделирования как инструмента создания фундамента для процесса проектирования и производства. В статье [2] обозначен путь от лабораторных испытаний материалов, выбора соответствующих действительному поведению материала определяющих уравнений до суперкомпьютерного моделирования поведения высокотемпературных инженерных компонентов. Д. Хейхарст использует континуальную механику поврежденности в качестве примера инструментария, который может быть применен для анализа и моделирования поведения повреждений широкого спектра инженерных компонентов, эксплуатируемых при высоких температурах. Подчеркивается важность использование механически обоснованных определяющих уравнений для достижения точных прогнозов и экстраполяций. Хейхарст обсуждает процедуры выбора доминирующих механизмов на основе лабораторных данных и, следовательно, соответствующих определяющих уравнений. Автор [1] приводит к утверждению, что препятствием для прогресса в использовании этих методов для создания основ и фундамента проектирования и новых технологий будет нехватка достоверных данных о материалах.
Основным технологическим изменением, произошедшим за последние годы, является появление недорогих компьютерных рабочих станций и суперкомпьютеров с доступом к современным, недорогим и быстрым средствам хранения данных. Вместе с этим численные методы и программное обеспечение для решения комбинированных задач с граничными и начальными значениями, часто с использованием метода конечных элементов (МКЭ), становятся доступными и надежными. Следовательно, становится возможным моделирование сложных физических процессов при проектировании и производстве при учете данных и информации о реальных материалах и реологических моделях. Параллельно с этими достижениями методы компьютерной визуализации и разработка средств видеоанимации достигли такого высокого уровня, что моделирование в реальном времени, иногда называемое виртуальной реальностью, становится доступным инструментом для использования при принятии решений в области проектирования и производства.
Одной из важных проблем современной механики твердого тела является внедрение вычислительной континуальной механики поврежденности в инструментарий МКЭ для анализа широкого спектра конструктивных компонентов элементов конструкций, находящихся в условиях ползучести на основе простых данных об одноосной ползучести. Преимущество метода заключается в том, что он обеспечивает прослеживаемость от определяющих уравнений, используемых в реологии процессов деформирования, повреждаемости и разрушения, связанных с фундаментальным поведением микроструктуры.
Таким образом, можно считать, что основы вычислительной механики поврежденности заложены еще в прошлом веке и сегодня ставший уже технологическим подход инкорпорирования сложных определяющих моделей в инструментарий МКЭ продолжает свое развитие и эксплуатируется множеством научных школ во всем мире.
Например, в [2] для оценки структурной целостности конструкции при повышенной температуре в реакторе на расплаве соли тория выполнен расчет повреждений при ползучести. Авторы отмечают, что модель, базирующаяся на теории повреждений при ползучести, и метод численного моделирования до сих пор не были предложены для ключевых материалов данного класса реакторов. В [2] исследуются характеристики повреждений при ползучести сплава UNS N10003 с использованием закона ползучести Нортона и модели повреждений Качанова Работнова в условиях ползучести. Во-первых, были приняты экспериментальные данные о ползучести сплава UNS N10003 при 650 °C, чтобы соответствовать материальным константам двух моделей. Во-вторых, было проанализировано и обсуждено поведение сплава UNS N10003 при повреждении при ползучести в одно- и многоосном напряженных состояниях. Результаты показали, что модель повреждения в условиях ползучести Качанова Работнова больше подходит для сплава UNS N10003, чем модель, основанная на классических уравнениях Бейли Нортона. Наконец, был разработан метод численного моделирования с помощью пользовательской подпрограммы UMAT и впоследствии проверен с помощью анализа МКЭ. Авторы показывают, что результаты конечно-элементного моделирования соответствовали имеющимся теоретическим оценкам. Авторы наглядно демонстрируют, что развитый ими подход представляет собой эффективный метод анализа повреждений в конструкциях в режиме ползучести при повышенной температуре в данном типе реакторов.
Принципы разработки пользовательских моделей материалов UMAT обсуждаются в [3]. Авторы подчеркивают важность основополагающих принципов, поскольку не все известные конституциональные модели материалов были смоделированы и подтверждены в качестве встроенных моделей материалов в рамках обычных МКЭ-решателей. Там, где обнаруживаются новые материалы или даже требуются модификации и усовершенствования существующих встроенных моделей материалов, пользователь МКЭ-пакетов обычно прибегает к вычислительному описанию своих версий конституциональных выражений, чтобы отразить такие изменения.
В [4] авторы представили прогнозирующее моделирование развития повреждений при ползучести материалов путем реализации модифицированных определяющих соотношений Робинсона Русселье и расширенного метода конечных элементов (XFEM) для решения проблемы разрушения при ползучести в задаче о росте трещины. Разработана новая модель, называемая модифицированной моделью Робинсона Русселье, для прогнозирования поведения повреждений при ползучести с точки зрения микромеханических повреждений вследствие роста пор в условиях ползучести. Интерфейс модели выполняет неявную схему интеграции в подпрограмме UMAT модуля Abaqus/Standard. Метод радиального возврата используется для интегрирования определяющего уравнения вязкопластичности в конечно-элементной формулировке. Численные модели в 2D и 3D элементах реализуются для определения корректности разработанных подпрограмм, а результаты сравниваются с точным решением для верификации. Кроме того, испытания на ползучесть при растяжении на образце из гладких прутков моделируются и испытываются при постоянной температуре 625 °C с различными уровнями напряжений. Результаты показывают, что максимальные значения напряжений, деформации ползучести и повреждения обнаруживаются вблизи центра образца при растяжении, где наблюдается шейкообразование. Кроме того, результаты сравниваются с известными литературными данными, чтобы проверить и оценить разработанную модель и показать разумное соответствие между обоими результатами. Затем проведенный анализ расширяется путем введения данных о развитии трещин в образце на основе метода XFEM. В результате предложена новая модель, получившая название модифицированной модели Робинсона Русселье, и результаты сравниваются с результатами, найденными в литературе, которые показали эволюцию роста пустот на пути распространения трещины. Таким образом, доказано [4], что модельное решение обладает потенциалом для прогнозирования поведения повреждений при ползучести с точки зрения роста трещин путем образования, роста и слияния пор на линии продолжения дефекта в конструкциях из пластичных материалов.
В [5] выполнены классический анализ ползучести и моделирование повреждений в непрерывном режиме ползучести для сварного соединения паропровода с исходным материалом 0,5Cr0,5Mo0,25В, металлом сварного шва 2,25CrMo и рассмотрена зона термического воздействия. Оценка срока службы при ползучести с использованием классического анализа деформации ползучести, включающего закон линейного накопления повреждений, часто в значительной степени отличается от практики. Поэтому в исследовании [5] выполнены три вида различного анализа, предполагающих классический анализ ползучести (без деградации материала из-за повреждений), непрерывное повреждение при ползучести при однозонном сварном пересечении и непрерывное повреждение при ползучести при многозонном сварном пересечении. Необходимые формулы для анализа континуальных повреждений при ползучести реализованы в пользовательской подпрограмме материалов (UMAT) и связаны с программным обеспечением ANSYS для конечных элементов. Представлено и обсуждено прогнозируемое распределение эквивалентных напряжений, эволюция повреждений и срок службы зарождения трещин на основе упомянутого выше моделирования. Показано, что прогнозируемый срок службы на основе классического анализа ползучести и с учетом многозонного моделирования значительно превышает срок службы, полученный на основе моделирования повреждений при ползучести. Однако прогнозируемые результаты анализа повреждений при ползучести по одно- и многозонным моделям сопоставимы с разницей примерно в 7,5 %.
Таким образом, к настоящему времени сформировалось устойчивое суждение о необходимости инкорпорирования конституциональных моделей материалов, включающих разнообразные меры поврежденности [15]. Тем не менее, несмотря на довольно сильную и достаточную изученность проблемы, которой занимается множество научных учреждений, многие вопросы остаются открытыми и требующими несомненного рассмотрения.
Задача моделирования напряженно-деформированного состояния (НДС) в условиях ползучести являлась фундаментальной проблемой нелинейной механики разрушения и остается актуальной в настоящее время [127]. Сегодня особый интерес вызывают исследование разрушения материалов и моделирование поврежденности тел с трещинами в условиях ползучести [617].
Как правило, разрушение в материале возникает из-за образования и развития множества трещин. Решения задач о телах с острой трещиной в линейной механике разрушения хорошо известны. Однако зачастую в реальных конструкциях трещина имеет конечный радиус закругления вследствие затупления ее вершины из-за пластических деформаций, образующихся при воздействии нагрузки, поэтому все чаще моделируют трещины, имеющие закругленные вершины [12]. Интерес вызывает исследование процессов разрушения в телах с трещинами именно с закругленными вершинами. Например, авторы работы [13] исследовали влияние радиуса закругления вершины выреза на поля НДС.
До настоящего времени предложено множество определяющих соотношений и кинетических уравнений для описания эволюции повреждений как совершенно новых [1416], так и базирующихся на модели поврежденности Качанова Работнова [23].
Нелинейная механика разрушения была развита с применением -интеграла. Впервые введенный Райсом в 1968 году -интеграл определяется следующим образом [24]:
(1)
В данном соотношении это нормаль к контуру интегрирования Г, плотность энергии деформации, перемещение.
Известно классическое решение Хатчинсона Райса Розенгрена [24; 25], связывающее -интеграл и поля НДС в вершине трещины при упругопластическом деформировании. Полученное решение описывается следующими соотношениями:
(2)
Эти асимптотические формулы часто называют ХРР-асимптотиками. Пользуясь аналогией Хоффа, распределения напряжений, деформаций и перемещений для случая установившейся ползучести представляются аналогично, но в терминах и вычисляются согласно следующим соотношениям:
null (3)
Известно, что -интеграл определяется как контурный интеграл по кривой , охватывающей вершину трещины, и может быть записан в следующей форме:
(4)
где
(5)
В формуле (3) в соответствии со структурой решения напряжения при будут стремиться к бесконечности, следовательно, отсутствует возможность исследования эффекта закругления вершин трещины. К тому же ХРР-решение не дает возможности учесть неоднородность материала, варианты геометриии и условия нагружения тел с трещиной.
В данном исследовании показано, что на определенном расстоянии независимо от формы вершины трещины выполняется ХРР-решение. Так актуальной и центральной становится задача извлечения асимптотического поведения напряжений, деформаций и сплошности вблизи окрестности вершины трещины в материале с введенными мерами поврежденности. Поскольку процесс накопления повреждений и эволюции НДС являются созависимыми, возникает вопрос возможности из численного решения извлечь асимптотическое поведение механических параметров задачи и является ли возможным отыскание асимптотического решения из имеющегося результата, полученного МКЭ-пакетом. Поэтому в данной статье ставится анализ радиального и углового поведения полей напряжений, деформаций и сплошности вблизи трещины в среде с учетом явления накопления повреждений.
Целью настоящего исследования является оценка влияния процесса накопления повреждений на асимптотическое поведение механических характеристик у вершины трещины средствами метода конечных элементов. Для анализа асимптотики полей и в вычислительную процедуру метода конечных элементов были с помощью UMAT [28; 29] внедрены конституциональные модели материалов, учитывающие процесс накопления повреждений.
2. Моделирование полей, ассоциированных с вершиной трещины в режиме ползучести
В пакете моделирования SIMULIA Abaqus, осуществляющем расчеты методом конечных элементов, проведено моделирование плоской пластины с центральной трещиной в условиях ползучести. Пластина имеет размеры 100 100 мм, длина трещины равна 10 мм. В первой серии расчетов трещина моделируется острой. Пластина находится в условиях одноосного растяжения.
Для моделирования ползучести избран степенной закон Бейли Нортона, имеющий следующий вид:
(6)
где компоненты скорости деформации ползучести, компоненты напряжений Коши, интенсивность напряжения, константы материала; компоненты девиатора напряжений, где символ Кронекера. Это определяющее соотношение было внедрено в расчет с помощью пользовательской подпрограммы UMAT.
Материальные константы в соотношении (6) это параметры материала, в общем случае зависящие от температуры. В качестве материала пластины выбрана сталь, имеющая следующие упругие модули: модуль Юнга , коэффициент Пуассона равен 0.3. Проведен ряд расчетов с различными материальными константами и . Принятые значения констант и условия экспериментальных исследований приведены в табл. 1.
Таблица 1 Условия проведения численных расчетов
Table 1 Conditions for making settlement calculations
модели | Модуль Юнга | Коэффициент Пуассона | , | | Нагрузка , |
.1 | 210000 | 0.3 | | 2 | 20 |
.2 | 210000 | 0.3 | | 3 | 20 |
.3 | 210000 | 0.3 | | 4 | 20 |
.4 | 210000 | 0.3 | | 5 | 20 |
Обычно показатель ползучести принимает значения от 3 до 8 для металлов и сплавов. Для отдельных современных высокопрочных сплавов достигает 18 [15; 22; 23]. Изображение модели пластины со схемой приложенной нагрузки представлено на рис. 1.
Рис. 1. Изображение модели со схемой приложенной нагрузки
Fig. 1. Scheme of the model geometry and applied loads
Сетка построена с использованием сингулярных конечных элементов в окрестности вершины трещины. Концентрические окружности поделены на 72 сектора, таким образом, раствор секторального элемента равен . Тип элементов сетки был избран CPS4. Полученная конечно-элементная сетка представлена на рис. 1.
Рис. 2. Конечно-элементная сетка в модели
Fig. 2. Typical finite element mesh of the model
В течение 35 000 часов к пластине с острой центральной трещиной была приложена одноосная растягивающая нагрузка. В результате расчета получены поля напряжений и деформаций в окрестности вершины трещины, находящейся в условиях ползучести.
Распределение компоненты тензора напряжений приведено вдоль радиальной траектории, построенной под углом , показанной на рис. 2, и построено в двойных логарифмических координатах.
Рис. 3. Путь, вдоль которого строятся радиальные распределения напряжений
Fig. 3. The path along which stress radial distributions are constructed
На рис. 4 изображены итоги вычислений в системе SIMULIA Abaqus для четырех различных моделей с острой центральной трещиной со свойствами в соответствии с табл. 1.
Поля компонент тензора напряжений приведены в различные моменты времени. Общее время моделирования равно 35 000 часов. Различными оттенками отмечены радиальные распределения поля напряжений в разные моменты времени.
Рис. 4. Результаты конечно-элементного анализа: компоненты тензора напряжений 22 в условиях ползучести c течением времени вдоль радиальной траектории ( = 0). Результаты для: a — модели 1.1, b — модели 1.2, c — модели 1.3, d — модели 1.4
Fig. 4. Results of FEM: the stress tensor components 22 under creep conditions along the radial trajectory ( = 0). Results for: a — model 1.1, b — model 1.2, c — model 1.3, d — model 1.4
Таким образом, для модели с острой центральной трещиной получено численное решение, и распределения напряжений в логарифмических координатах приведены на рис. 4.
Для каждого графика на рис. 4 построена прямая черная линия, аппроксимирующая значения компоненты тензора напряжений в последний момент времени в зоне ползучести. Полученные значения коэффициента наклона этих прямых соответствуют значениям, полученным из соотношения . Линия, аппроксимирующая значения компоненты тензора напряжений в зоне упругости, имеет при этом коэффициент наклона равный , что также соответствует аналитическому решению. Это показывает, что асимптотическое решение ХРР присутствует в расчетах для острой трещины.
Многочисленные исследования показали, что расчеты напряженно-деформированного состояния в условиях ползучести необходимо проводить в образцах с трещинами, имеющими закругленные вершины [12; 13; 27]. Ввиду этого следующая серия вычислений проведена для трещин с закругленными концами. Трещина моделируется в виде разреза с закругленными вершинами, радиус закругления составляет 0.001 мм. В настоящем исследовании проведен ряд расчетов с различными материальными константами и . Принятые значения констант и условия экспериментальных исследований приведены в табл. 2.
Таблица 2 Условия проведения численных расчетов
Table 2 Conditions for making settlement calculations
модели | Модуль Юнга | Коэффициент Пуассона | , | | Нагрузка , |
.1 | 210000 | 0.3 | | 2 | 20 |
.2 | 210000 | 0.3 | | 3 | 20 |
.3 | 210000 | 0.3 | | 4 | 20 |
Сетка в окрестности вершины трещины построена с использованием сингулярных конечных элементов. Концентрические окружности поделены на три области: перед вершиной дуга окружности величиной разбита на 36 секторов, таким образом, раствор секторального элемента равен , две другие части поделены на 17 секторов. Тип элементов сетки был избран CPS4. Полученная конечно-элементная сетка представлена на рис. 4.
В численных экспериментах в течение 35 000 часов к пластине с закругленной центральной трещиной прикладывалась одноосная растягивающая нагрузка. В результате расчета были получены поля напряжений, упругих деформаций и деформаций ползучести в окрестности вершины трещины, находящейся в условиях ползучести.
На рис. 7 изображены поля компонент тензора напряжений в логарифмических координатах с течением времени вдоль траектории, показанной на рис. 5. Различными оттенками отмечены радиальные распределения поля напряжений в разные моменты времени.
Рис. 5. Конечно-элементная сетка в модели с закругленной трещиной
Fig. 5. Typical finite element mesh of the model
Рис. 6. Путь, вдоль которого строятся радиальные распределения напряжений
Fig. 6. The path along which stress radial distributions are constructed
Рис. 7. Результаты конечно-элементного анализа: компоненты тензора напряжений 22 в условиях ползучести c течением времени. Результаты приведены для: a — модели 2.1, b — модели 2.2, c — модели 2.3
Fig. 7. Results of FEM: the stress tensor components 22 under creep conditions. Results for: a — model 2.1, b — model 2.2, c — model 2.3
Для каждого графика построена черная линия, аппроксимирующая значения компоненты тензора напряжений в последний момент времени в зоне ползучести. Полученные значения коэффициента наклона этих прямых находятся в согласии со значениями, соответствующим формуле . Можно отметить, что в области, близкой к вершине закругленной трещины, асимптотика меняется, в отличие от предыдущих расчетов с острой трещиной. Это связано с тем, что такая конфигурация трещины приводит к изменению сингулярности вблизи трещины.
Линия, аппроксимирующая значения компоненты тензора напряжений в зоне упругости, имеет при этом коэффициент наклона равный , что также соответствует аналитическому решению.
По графикам на рис. 7 можно определить границы области ползучести и зоны упругости. С течением времени область, находящаяся в условиях ползучести, увеличивается, а область, находящаяся в режиме упругости, уменьшается. Таким образом, по результатам расчета модели показаны промежутки, расстояния, на которых работает асимптотика ХРР. Определив диапазон значений расстояний, на которых работает решение ХРР, переходим к исследованию влияния на асимптотику процесса накопления поврежденности.
3. Анализ влияния прогрессивного накопления микроповреждений у кончика трещины на асимптотику
Задачей последующего исследования выступает анализ влияния развитой эволюции поврежденности вблизи вершины закругленной трещины в предположении реализации плоского напряженного состояния для степенного закона установившейся ползучести и кинетического уравнения Качанова Работнова на асимптотику напряжений.
Степенной закон ползучести, следуя классической процедуре [20], был дополнен кинетическим уравнением накопления поврежденности Качанова Работнова и математическая модель представляется в следующем виде:
(7)
где параметр сплошности Качанова [18]; гидростатическое напряжение, максимальное главное напряжение, эквивалентное напряжение, а константы находятся экспериментально. Неповрежденному, целостному материалу соответствует значение параметра , а означает, что материал полностью исчерпал несущую способность. Вместо параметра сплошности часто используется параметр поврежденности Работнова [20], связанный с параметром сплошности выражением .
Для конечно-элементного расчета необходимо задать постоянные материала фигурирующие в эволюционном уравнении 7. Обычно выбирают таким образом, чтобы выполнялось соотношение , полученное эмпирическим путем [22].
В ходе исследования проведен ряд расчетов с различными материальными константами и . Принятые значения констант и условия экспериментальных исследований приведены в табл. 3.
Таблица 3 Условия проведения численных расчетов
Table 3 Conditions for making settlement calculations
модели | Модуль Юнга | Коэффициент Пуассона | , | | , | | Нагрузка , |
.1 | 210000 | 0.3 | | 3 | | 2.5 | 10 |
.2 | 210000 | 0.3 | | 4 | | 2.8 | 10 |
.3 | 210000 | 0.3 | | 5 | | 3.5 | 10 |
Постоянные материала в эволюционном уравнении и определяющем соотношении (7) выбирались на основании экспериментальных данных, описанных в работах [15; 23], где приведены значения материальных констант для большого ряда металлов и сплавов. С использованием подпрограммы, являющейся процедурой UMAT комплекса SIMULIA Abaqus, конституциональные соотношения степенного закона ползучести и эволюционное уравнение (7) были внесены в вычислительную схему метода конечных элементов пакета Abaqus/Standart. Это дало возможность найти распределения параметра сплошности и механических полей.
Проведена серия расчетов одноосного растяжения пластины с центральной трещиной с закругленными концами в условиях ползучести с учетом процессов накопления поврежденности. В результате расчета были получены поля напряжений, деформаций и сплошности.
На рис. 7 изображены поля компоненты тензора напряжений в двойных логарифмических координатах с течением времени вдоль траектории, показанной на рис. 5 и построенной по прямой под углом к оси . Различными оттенками отмечены радиальные распределения поля напряжений в различные моменты времени.
Рис. 8. Результаты конечно-элементного анализа: компоненты тензора напряжений 22 в условиях ползучести при действии процесса накопления поврежденности c течением времени. Результаты для: a — модели 3.1, b — модели 3.2, c — модели 3.3
Fig. 8. Results of FEM: the stress tensor components 22 under creep condition under the action of damage accumulation over time. Results for: a — model 3.1, b — model 3.2, c — model 3.3
Для каждого графика построена регрессионная прямая, отмеченная на графиках черной линией, аппроксимирующая значения компоненты тензора напряжений в последний рассматриваемый момент времени в зоне ползучести. Полученные значения коэффициента наклона этой прямой существенно отличаются от значений, соответствующих формуле , а значит и асимптотике решения ХРР. Здесь также можно отметить, что в области, близкой к вершине закругленной трещины, асимптотика не наблюдается, в отличие от расчетов с острой трещиной. Результаты на рис. 7 отличаются от распределений, полученных на расчетах при действии ползучести без учета накопления поврежденности, как на рис. 7. По графикам на рис. 7 можно заметить, каким образом процесс накопления поврежденности изменяет асимптотику напряжений в окрестности вершины трещины.
На рис. 8 показано распределения сплошности с течением времени: a c состояние в 922.9, 3 923, 7 923 часа, d f состояние в 222.9, 9 923, 3 5000 часов.
На рис. 9 представлены картины распределения сплошности вдоль траектории, показанной на рис. 5, с течением времени. Различными оттенками отмечены радиальные распределения сплошности в разные моменты времени.
Рис. 9. Распределение сплошности c течением времени в a–c модели 3.1, d–f в модели 3.3
Fig. 9. Distribution of continuity over time a–c for model 3.1, d–f for model 3.3
Рис. 10. Результаты конечно-элементного анализа: распределение сплошности в условиях ползучести с учетом накопления поврежденности c течением времени. Результаты для: a — модели 3.1, b — модели 3.2, c — модели 3.3
Fig. 10. Results of FEM: distribution of continuity under creep conditions, taking into account the accumulation of damage over time. Results for: a — model 3.1, b — model 3.2, c — model 3.3
Из рис. 8, 9 можно заметить, каким образом для разных значений материальных констант меняется поле сплошности.
Заключение
В ходе исследования выполнено конечно-элементное моделирование одноосного растяжения двумерной пластины с центральной трещиной, находящейся в условиях ползучести с учетом поврежденности в связанной постановке. Моделирование выполнено в комплексе SIMULIA Abaqus с применением пользовательской процедуры UMAT. Степенной закон ползучести с помощью пользовательской процедуры был дополнен кинетическим уравнением накопления поврежденности Качанова Работнова в связанной постановке. В результате моделирования получены распределения напряжений, деформаций и сплошности в условиях ползучести с учетом накопления поврежденности с течением времени. Приведены асимптотики распределения компоненты тензора напряжений для моделей как с острой трещиной, так и с вырезом с закругленными вершинами. Моделирование выполнено для различных значений материальных констант. Приведено радиальное распределение сплошности вдоль траектории, построенной под углом к горизонтальной оси . Показано, что асимптотика в упругой области и в области развитой ползучести в первых двух сериях расчетов соответствует аналитическим решениям, однако в случае учета накопления поврежденности эта асимптотика не наблюдается.